Skip to content
Snippets Groups Projects
Commit a08c64f1 authored by Vakhtang Tsulaia's avatar Vakhtang Tsulaia
Browse files

Merge branch 'GeoModelTestRemSagTest' into 'main'

MuonGeoModelTest - Remove scan test along the sagged surface

See merge request !75789
parents 6cc5390b 1bcb84ad
No related branches found
No related tags found
9 merge requests!77731Draft: Updates to ZDC reconstruction,!77728Draft: updates to ZDC reconstruction,!77522Draft: sTGC Pad Trigger Emulator,!76725ZdcNtuple: Fix cppcheck warning.,!76611L1CaloFEXByteStream: Fix out-of-bounds array accesses.,!76475Punchthrough AF3 implementation in FastG4,!76474Punchthrough AF3 implementation in FastG4,!75789MuonGeoModelTest - Remove scan test along the sagged surface,!75127Draft: Seeder types rebase2 grid fix
......@@ -183,7 +183,6 @@ StatusCode GeoModelMdtTest::dumpToTree(const EventContext& ctx, const MdtReadout
m_tubeTransform.push_back(layTransf);
const Amg::Vector3D tubePos = layTransf.translation();
const Amg::Vector3D roPos = readoutEle->ROPos(tube_id);
m_roPos.push_back(roPos);
m_activeTubeLength.push_back(readoutEle->getActiveTubeLength(lay,tube));
......@@ -200,22 +199,6 @@ StatusCode GeoModelMdtTest::dumpToTree(const EventContext& ctx, const MdtReadout
m_cablingMROD.push_back(translation.mrod);
m_cablingTdcId.push_back(translation.tdcId);
m_cablingTdcCh.push_back(translation.channelId);
}
if (!m_dumpSurfaces) continue;
const Amg::Vector3D globalDir {(tubePos - roPos).unit()};
const Trk::SaggedLineSurface& surf{readoutEle->surface(tube_id)};
for (double l = readoutEle->tubeLength(tube_id) /2; l > 0; l = l -100. ) {
Amg::Vector2D lPos{Amg::Vector2D::Zero()};
surf.globalToLocal(roPos + l * globalDir,Amg::Vector3D::Zero(),lPos);
std::unique_ptr<Trk::StraightLineSurface> sagged{surf.correctedSurface(lPos)};
m_layDistTubeLay.push_back(lay);
m_layDistTubeNum.push_back(tube);
m_layDistPosAlongWire.push_back(l);
const Amg::Transform3D layTransf{sagged->transform()};
m_layDist.push_back(layTransf);
}
}
}
......
......@@ -47,7 +47,6 @@ class GeoModelMdtTest : public AthHistogramAlgorithm {
Gaudi::Property<std::vector<std::string>> m_selectStat{
this, "TestStations", {"BIL1A3"}};
Gaudi::Property<bool> m_dumpSurfaces{this, "dumpSurfaces", false, "Adds the bounds and surfaces of each tube to the dump"};
/// Write a TTree for validation purposes
......@@ -91,13 +90,6 @@ class GeoModelMdtTest : public AthHistogramAlgorithm {
/// Position of the readout
MuonVal::ThreeVectorBranch m_roPos{m_tree, "readOutPos"};
/// Distorted transformations
MuonVal::VectorBranch<unsigned short>& m_layDistTubeLay{m_tree.newVector<unsigned short>("DistTubeLayer")};
MuonVal::VectorBranch<unsigned short>& m_layDistTubeNum{m_tree.newVector<unsigned short>("DistTubeNumber")};
MuonVal::VectorBranch<double>& m_layDistPosAlongWire{m_tree.newVector<double>("DistPosAlongTube")};
MuonVal::CoordSystemsBranch m_layDist{m_tree, "DistLayer"};
/// Alignment parameters
MuonVal::ScalarBranch<float>& m_ALineTransS{m_tree.newScalar<float>("ALineTransS", 0.)};
MuonVal::ScalarBranch<float>& m_ALineTransT{m_tree.newScalar<float>("ALineTransT", 0.)};
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment