ATLAS Offline Software
Public Member Functions | Private Member Functions | Private Attributes | List of all members
ActsTrackingGeometrySvc Class Reference

#include <ActsTrackingGeometrySvc.h>

Inheritance diagram for ActsTrackingGeometrySvc:
Collaboration diagram for ActsTrackingGeometrySvc:

Public Member Functions

StatusCode initialize () override
 
 ActsTrackingGeometrySvc (const std::string &name, ISvcLocator *pSvcLocator)
 
std::shared_ptr< const Acts::TrackingGeometry > trackingGeometry () override
 
unsigned int populateAlignmentStore (ActsTrk::DetectorAlignStore &store) const override
 
const ActsGeometryContext & getNominalContext () const override
 

Private Member Functions

ActsLayerBuilder::Config makeLayerBuilderConfig (const InDetDD::InDetDetectorManager *manager)
 
std::shared_ptr< const Acts::ILayerBuilder > makeStrawLayerBuilder (const InDetDD::InDetDetectorManager *manager)
 
std::shared_ptr< const Acts::ILayerBuilder > makeHGTDLayerBuilder (const HGTD_DetectorManager *manager)
 
std::shared_ptr< Acts::TrackingVolume > makeSCTTRTAssembly (const Acts::GeometryContext &gctx, const Acts::ILayerBuilder &sct_lb, const Acts::ILayerBuilder &trt_lb, const Acts::CylinderVolumeHelper &cvh, const std::shared_ptr< const Acts::TrackingVolume > &pixel)
 
Acts::CylinderVolumeBuilder::Config makeBeamPipeConfig (std::shared_ptr< const Acts::CylinderVolumeHelper > cvh) const
 
bool runConsistencyChecks () const
 

Private Attributes

ServiceHandle< StoreGateSvc > m_detStore
 
const InDetDD::SiDetectorManager * p_pixelManager {nullptr}
 
const InDetDD::SiDetectorManager * p_SCTManager {nullptr}
 
const InDetDD::TRT_DetectorManager * p_TRTManager {nullptr}
 
const InDetDD::SiDetectorManager * p_ITkPixelManager {nullptr}
 
const InDetDD::SiDetectorManager * p_ITkStripManager {nullptr}
 
const BeamPipeDetectorManager * p_beamPipeMgr {nullptr}
 
const HGTD_DetectorManager * p_HGTDManager {nullptr}
 
std::shared_ptr< ActsElementVector > m_elementStore {nullptr}
 
std::shared_ptr< const Acts::TrackingGeometry > m_trackingGeometry {nullptr}
 
const TRT_ID * m_TRT_idHelper {nullptr}
 
const HGTD_ID * m_HGTD_idHelper {nullptr}
 
ActsGeometryContext m_nominalContext {}
 
Gaudi::Property< bool > m_useMaterialMap {this, "UseMaterialMap", false, ""}
 
Gaudi::Property< bool > m_objDebugOutput {this, "ObjDebugOutput", false, ""}
 
Gaudi::Property< std::string > m_materialMapInputFileBase {this, "MaterialMapInputFile", "", ""}
 
Gaudi::Property< std::string > m_materialMapCalibFolder {this, "MaterialMapCalibFolder", ".", ""}
 
Gaudi::Property< bool > m_buildBeamPipe {this, "BuildBeamPipe", false, ""}
 
Gaudi::Property< std::vector< size_t > > m_barrelMaterialBins {this, "BarrelMaterialBins", {10, 10}}
 
Gaudi::Property< std::vector< size_t > > m_endcapMaterialBins {this, "EndcapMaterialBins", {5, 20}}
 
Gaudi::Property< std::vector< std::string > > m_buildSubdetectors {this, "BuildSubDetectors", {"Pixel", "SCT", "TRT", "Calo", "HGTD"}}
 
BooleanProperty m_runConsistencyChecks
 
StringProperty m_consistencyCheckOutput
 
Gaudi::Property< size_t > m_consistencyCheckPoints
 
ToolHandle< IActsTrackingVolumeBuilder > m_caloVolumeBuilder
 
Gaudi::Property< std::vector< unsigned int > > m_subDetNoAlignProp {this, "NotAlignDetectors", {}}
 Define the subdetectors for which the tracking geometry does not expect a valid alignment store. More...
 
std::set< ActsTrk::DetectorType > m_subDetNoAlign {}
 

Detailed Description

Definition at line 51 of file ActsTrackingGeometrySvc.h.

Constructor & Destructor Documentation

◆ ActsTrackingGeometrySvc()

ActsTrackingGeometrySvc::ActsTrackingGeometrySvc ( const std::string &  name,
ISvcLocator *  pSvcLocator 
)

Definition at line 60 of file ActsTrackingGeometrySvc.cxx.

62  : base_class(name, svc),
63  m_detStore("StoreGateSvc/DetectorStore", name),
64  m_elementStore (std::make_shared<ActsElementVector>())
65 {
66 }

Member Function Documentation

◆ getNominalContext()

const ActsGeometryContext & ActsTrackingGeometrySvc::getNominalContext ( ) const
override

Definition at line 1021 of file ActsTrackingGeometrySvc.cxx.

1021 { return m_nominalContext; }

◆ initialize()

StatusCode ActsTrackingGeometrySvc::initialize ( )
override

Definition at line 68 of file ActsTrackingGeometrySvc.cxx.

68  {
69  ATH_MSG_INFO(name() << " is initializing");
70  for (unsigned int skipAlign : m_subDetNoAlignProp) {
71  try {
72  m_subDetNoAlign.insert(static_cast<DetectorType>(skipAlign));
73  } catch (...) {
74  ATH_MSG_FATAL("Failed to interpret " << m_subDetNoAlignProp << " as ActsDetectorElements");
75  return StatusCode::FAILURE;
76  }
77 }
78 
79  // FIXME: ActsCaloTrackingVolumeBuilder holds ReadHandle to
80  // CaloDetDescrManager. Hopefully this service is never called before that
81  // object is available.
82  m_autoRetrieveTools = false;
83  m_checkToolDeps = false;
84 
85  ATH_MSG_INFO("ACTS version is: v"
86  << Acts::VersionMajor << "." << Acts::VersionMinor << "."
87  << Acts::VersionPatch << " [" << Acts::CommitHash << "]");
88 
89  // load which subdetectors to build from property
90  std::set<std::string> buildSubdet(m_buildSubdetectors.begin(),
91  m_buildSubdetectors.end());
92  ATH_MSG_INFO("Configured to build " << buildSubdet.size()
93  << " subdetectors:");
94  for (const auto &s : buildSubdet) {
95  ATH_MSG_INFO(" - " << s);
96  }
97 
98  ATH_MSG_DEBUG("Loading detector manager(s)");
99  if (buildSubdet.find("Pixel") != buildSubdet.end()) {
100  ATH_CHECK(m_detStore->retrieve(p_pixelManager, "Pixel"));
101  }
102  if (buildSubdet.find("SCT") != buildSubdet.end()) {
103  ATH_CHECK(m_detStore->retrieve(p_SCTManager, "SCT"));
104  }
105  if (buildSubdet.find("TRT") != buildSubdet.end()) {
106  ATH_CHECK(m_detStore->retrieve(p_TRTManager, "TRT"));
107  ATH_CHECK(m_detStore->retrieve(m_TRT_idHelper, "TRT_ID"));
108  }
109  if (buildSubdet.find("ITkPixel") != buildSubdet.end()) {
110  ATH_CHECK(m_detStore->retrieve(p_ITkPixelManager, "ITkPixel"));
111  }
112  if (buildSubdet.find("ITkStrip") != buildSubdet.end()) {
113  ATH_CHECK(m_detStore->retrieve(p_ITkStripManager, "ITkStrip"));
114  }
115  if (buildSubdet.find("HGTD") != buildSubdet.end()) {
116  ATH_CHECK(m_detStore->retrieve(p_HGTDManager, "HGTD"));
117  ATH_CHECK(m_detStore->retrieve(m_HGTD_idHelper, "HGTD_ID"));
118  }
119 
120  if(m_buildBeamPipe) {
121  ATH_CHECK(m_detStore->retrieve(p_beamPipeMgr, "BeamPipe"));
122  }
123 
124 
125  ATH_MSG_DEBUG("Setting up ACTS geometry helpers");
126 
127  Acts::LayerArrayCreator::Config lacCfg;
128  auto layerArrayCreator = std::make_shared<const Acts::LayerArrayCreator>(
129  lacCfg, makeActsAthenaLogger(this, std::string("LayArrCrtr"), std::string("ActsTGSvc")));
130 
131  Acts::TrackingVolumeArrayCreator::Config tvcCfg;
132  auto trackingVolumeArrayCreator =
133  std::make_shared<const Acts::TrackingVolumeArrayCreator>(
134  tvcCfg, makeActsAthenaLogger(this, std::string("TrkVolArrCrtr"), std::string("ActsTGSvc")));
135 
136  Acts::CylinderVolumeHelper::Config cvhConfig;
137  cvhConfig.layerArrayCreator = layerArrayCreator;
138  cvhConfig.trackingVolumeArrayCreator = trackingVolumeArrayCreator;
139 
140  auto cylinderVolumeHelper =
141  std::make_shared<const Acts::CylinderVolumeHelper>(
142  cvhConfig, makeActsAthenaLogger(this, std::string("CylVolHlpr"), std::string("ActsTGSvc")));
143 
144  Acts::TrackingGeometryBuilder::Config tgbConfig;
145  tgbConfig.trackingVolumeHelper = cylinderVolumeHelper;
146 
147  if (m_useMaterialMap) {
148  std::shared_ptr<const Acts::IMaterialDecorator> matDeco = nullptr;
149 
150  std::string matFileFullPath = PathResolverFindCalibFile(m_materialMapCalibFolder.value()+"/"+m_materialMapInputFileBase.value());
151  if (matFileFullPath.empty()) {
152  ATH_MSG_ERROR( "Material Map Input File " << m_materialMapCalibFolder.value() << "/" << m_materialMapInputFileBase.value() << " not found.");
153  return StatusCode::FAILURE;
154  }
155  ATH_MSG_INFO("Configured to use material input: " << matFileFullPath);
156 
157  if (matFileFullPath.find(".json") != std::string::npos) {
158  // Set up the converter first
159  Acts::MaterialMapJsonConverter::Config jsonGeoConvConfig;
160  // Set up the json-based decorator
161  matDeco = std::make_shared<const Acts::JsonMaterialDecorator>(
162  jsonGeoConvConfig, matFileFullPath, ActsTrk::actsLevelVector(msg().level()));
163  }
164  tgbConfig.materialDecorator = matDeco;
165  }
166 
167  std::array<double, 2> sctECEnvelopeZ{20_mm, 20_mm};
168 
169  try {
170  // BeamPipe
171  if(m_buildBeamPipe) {
172  tgbConfig.trackingVolumeBuilders.push_back([&](const auto &gctx,
173  const auto &inner,
174  const auto &) {
175 
176  Acts::CylinderVolumeBuilder::Config bpvConfig =
177  makeBeamPipeConfig(cylinderVolumeHelper);
178 
179  Acts::CylinderVolumeBuilder beamPipeVolumeBuilder {
180  bpvConfig, makeActsAthenaLogger(this, std::string("BPVolBldr"), std::string("ActsTGSvc"))};
181 
182  return beamPipeVolumeBuilder.trackingVolume(gctx, inner);
183  });
184  }
185 
186 
187 
188  // PIXEL
189  if (buildSubdet.count("Pixel") > 0) {
190  tgbConfig.trackingVolumeBuilders.push_back([&](const auto &gctx,
191  const auto &inner,
192  const auto &) {
193  auto cfg = makeLayerBuilderConfig(p_pixelManager);
194  cfg.mode = ActsLayerBuilder::Mode::Pixel;
195  auto lb = std::make_shared<ActsLayerBuilder>(
196  cfg, makeActsAthenaLogger(this, std::string("PixelGMSLayBldr"), std::string("ActsTGSvc")));
197  Acts::CylinderVolumeBuilder::Config cvbConfig;
198  cvbConfig.layerEnvelopeR = {3_mm, 3_mm};
199  cvbConfig.layerEnvelopeZ = 1_mm;
200  cvbConfig.trackingVolumeHelper = cylinderVolumeHelper;
201  cvbConfig.volumeName = "Pixel";
202  cvbConfig.layerBuilder = lb;
203  cvbConfig.buildToRadiusZero = !m_buildBeamPipe;
204 
205  Acts::CylinderVolumeBuilder cvb(
206  cvbConfig, makeActsAthenaLogger(this, std::string("CylVolBldr"), std::string("ActsTGSvc")));
207 
208  return cvb.trackingVolume(gctx, inner);
209  });
210  }
211 
212  // ITK PIXEL
213  if (buildSubdet.count("ITkPixel") > 0) {
214  tgbConfig.trackingVolumeBuilders.push_back(
215  [&](const auto &gctx, const auto &inner, const auto &) {
216  auto cfg = makeLayerBuilderConfig(p_ITkPixelManager);
217  cfg.mode = ActsLayerBuilder::Mode::ITkPixelInner;
218  cfg.objDebugOutput = m_objDebugOutput;
219  cfg.doEndcapLayerMerging = true;
220  auto lb = std::make_shared<ActsLayerBuilder>(
221  cfg, makeActsAthenaLogger(this, std::string("ITkPxInLb"), std::string("ActsTGSvc")));
222 
223  Acts::CylinderVolumeBuilder::Config cvbConfig;
224  cvbConfig.layerEnvelopeR = {5_mm, 5_mm};
225  cvbConfig.layerEnvelopeZ = 1_mm;
226  cvbConfig.trackingVolumeHelper = cylinderVolumeHelper;
227  cvbConfig.volumeName = "ITkPixelInner";
228  cvbConfig.layerBuilder = lb;
229  cvbConfig.buildToRadiusZero = !m_buildBeamPipe;
230 
231  Acts::CylinderVolumeBuilder cvb(
232  cvbConfig,
233  makeActsAthenaLogger(this, std::string("CylVolBldr"), std::string("ActsTGSvc")));
234 
235  return cvb.trackingVolume(gctx, inner);
236  });
237 
238  tgbConfig.trackingVolumeBuilders.push_back(
239  [&](const auto &gctx, const auto &inner, const auto &) {
240  auto cfg = makeLayerBuilderConfig(p_ITkPixelManager);
241  cfg.mode = ActsLayerBuilder::Mode::ITkPixelOuter;
242  cfg.objDebugOutput = m_objDebugOutput;
243  cfg.doEndcapLayerMerging = false;
244  auto lb = std::make_shared<ActsLayerBuilder>(
245  cfg, makeActsAthenaLogger(this, std::string("ITkPxOtLb"), std::string("ActsTGSvc")));
246 
247  Acts::CylinderVolumeBuilder::Config cvbConfig;
248  cvbConfig.layerEnvelopeR = {5_mm, 5_mm};
249  cvbConfig.layerEnvelopeZ = 1_mm;
250  cvbConfig.trackingVolumeHelper = cylinderVolumeHelper;
251  cvbConfig.volumeName = "ITkPixelOuter";
252  cvbConfig.layerBuilder = lb;
253  cvbConfig.buildToRadiusZero = false;
254  cvbConfig.checkRingLayout = true;
255  cvbConfig.ringTolerance = 10_mm;
256 
257  Acts::CylinderVolumeBuilder cvb(
258  cvbConfig,
259  makeActsAthenaLogger(this, std::string("CylVolBldr"), std::string("ActsTGSvc")));
260 
261  return cvb.trackingVolume(gctx, inner);
262  });
263  }
264 
265  // ITK STRIP
266  if (buildSubdet.count("ITkStrip") > 0) {
267  tgbConfig.trackingVolumeBuilders.push_back(
268  [&](const auto &gctx, const auto &inner, const auto &) {
269  auto cfg = makeLayerBuilderConfig(p_ITkStripManager);
270  cfg.mode = ActsLayerBuilder::Mode::ITkStrip;
271  cfg.objDebugOutput = m_objDebugOutput;
272  auto lb = std::make_shared<ActsLayerBuilder>(
273  cfg, makeActsAthenaLogger(this, std::string("ITkStripLB"), std::string("ActsTGSvc")));
274 
275  Acts::CylinderVolumeBuilder::Config cvbConfig;
276  cvbConfig.layerEnvelopeR = {5_mm, 5_mm};
277  cvbConfig.layerEnvelopeZ = 1_mm;
278  cvbConfig.trackingVolumeHelper = cylinderVolumeHelper;
279  cvbConfig.volumeName = "ITkStrip";
280  cvbConfig.layerBuilder = lb;
281  cvbConfig.buildToRadiusZero =
282  buildSubdet.count("ITkPixel") == 0 && !m_buildBeamPipe;
283 
284  Acts::CylinderVolumeBuilder cvb(
285  cvbConfig,
286  makeActsAthenaLogger(this, std::string("CylVolBldr"), std::string("ActsTGSvc")));
287 
288  return cvb.trackingVolume(gctx, inner);
289  });
290  }
291 
292  bool buildSCT = buildSubdet.count("SCT") > 0;
293  bool buildTRT = buildSubdet.count("TRT") > 0;
294 
295  if (buildSCT && buildTRT) {
296  // building both we need to take care
297  tgbConfig.trackingVolumeBuilders.push_back(
298  [&](const auto &gctx, const auto &inner, const auto &) {
299  auto cfg = makeLayerBuilderConfig(p_SCTManager);
300  cfg.mode = ActsLayerBuilder::Mode::SCT;
301  cfg.endcapEnvelopeZ = sctECEnvelopeZ;
302  auto sct_lb = std::make_shared<ActsLayerBuilder>(
303  cfg, makeActsAthenaLogger(this, std::string("SCTGMSLayBldr"), std::string("ActsTGSvc")));
304 
305  auto trt_lb = makeStrawLayerBuilder(p_TRTManager);
306 
307  return makeSCTTRTAssembly(gctx, *sct_lb, *trt_lb,
308  *cylinderVolumeHelper, inner);
309  });
310 
311  } else if (buildSCT) {
312  tgbConfig.trackingVolumeBuilders.push_back(
313  [&](const auto &gctx, const auto &inner, const auto &) {
314  auto lbCfg = makeLayerBuilderConfig(p_SCTManager);
315  lbCfg.mode = ActsLayerBuilder::Mode::SCT;
316  lbCfg.endcapEnvelopeZ = sctECEnvelopeZ;
317  auto lb = std::make_shared<ActsLayerBuilder>(
318  lbCfg,
319  makeActsAthenaLogger(this, std::string("SCTGMSLayBldr"), std::string("ActsTGSvc")));
320 
321  Acts::CylinderVolumeBuilder::Config cvbConfig;
322  cvbConfig.layerEnvelopeR = {5_mm, 5_mm};
323  cvbConfig.layerEnvelopeZ = 2_mm;
324  cvbConfig.trackingVolumeHelper = cylinderVolumeHelper;
325  cvbConfig.volumeName = "SCT";
326  cvbConfig.layerBuilder = lb;
327  cvbConfig.buildToRadiusZero = false;
328 
329  Acts::CylinderVolumeBuilder cvb(
330  cvbConfig,
331  makeActsAthenaLogger(this, std::string("SCTCylVolBldr"), std::string("ActsTGSvc")));
332 
333  return cvb.trackingVolume(gctx, inner);
334  });
335  } else if (buildTRT) {
336  tgbConfig.trackingVolumeBuilders.push_back(
337  [&](const auto &gctx, const auto &inner, const auto &) {
338  auto lb = makeStrawLayerBuilder(p_TRTManager);
339  Acts::CylinderVolumeBuilder::Config cvbConfig;
340  cvbConfig.layerEnvelopeR = {5_mm, 5_mm};
341  cvbConfig.layerEnvelopeZ = 2_mm;
342  cvbConfig.trackingVolumeHelper = cylinderVolumeHelper;
343  cvbConfig.volumeName = "TRT";
344  cvbConfig.layerBuilder = lb;
345  cvbConfig.buildToRadiusZero = false;
346 
347  Acts::CylinderVolumeBuilder cvb(
348  cvbConfig,
349  makeActsAthenaLogger(this, std::string("TRTCylVolBldr"), std::string("ActsTGSvc")));
350 
351  return cvb.trackingVolume(gctx, inner);
352  });
353  }
354 
355  //HGTD
356  if(buildSubdet.count("HGTD") > 0) {
357  tgbConfig.trackingVolumeBuilders.push_back(
358  [&](const auto &gctx, const auto &inner, const auto &) {
359  auto lb = makeHGTDLayerBuilder(p_HGTDManager); //using ActsHGTDLayerBuilder
360  Acts::CylinderVolumeBuilder::Config cvbConfig;
361  cvbConfig.layerEnvelopeR = {5_mm, 5_mm};
362  cvbConfig.layerEnvelopeZ = 1_mm;
363  cvbConfig.trackingVolumeHelper = cylinderVolumeHelper;
364  cvbConfig.volumeName = "HGTD";
365  cvbConfig.layerBuilder = lb;
366  cvbConfig.buildToRadiusZero = false;
367 
368  Acts::CylinderVolumeBuilder cvb(
369  cvbConfig,
370  makeActsAthenaLogger(this, std::string("HGTDCylVolBldr"), std::string("ActsTGSvc")));
371 
372  return cvb.trackingVolume(gctx, inner);
373  });
374  }
375 
376  // Calo
377  if (buildSubdet.count("Calo") > 0) {
378  tgbConfig.trackingVolumeBuilders.push_back(
379  [&](const auto &gctx, const auto &inner, const auto &) {
380  return m_caloVolumeBuilder->trackingVolume(gctx, inner, nullptr);
381  });
382  }
383  } catch (const std::exception &e) {
384  ATH_MSG_ERROR("Encountered error when building Acts tracking geometry");
385  ATH_MSG_ERROR(e.what());
386  return StatusCode::FAILURE;
387  }
388 
389  auto trackingGeometryBuilder =
390  std::make_shared<const Acts::TrackingGeometryBuilder>(
391  tgbConfig, makeActsAthenaLogger(this, std::string("TrkGeomBldr"), std::string("ActsTGSvc")));
392 
393  ATH_MSG_VERBOSE("Begin building process");
394  m_trackingGeometry =
395  trackingGeometryBuilder->trackingGeometry(getNominalContext().context());
396  ATH_MSG_VERBOSE("Building process completed");
397 
398  if (!m_trackingGeometry) {
399  ATH_MSG_ERROR("No ACTS tracking geometry was built. Cannot proceeed");
400  return StatusCode::FAILURE;
401  }
402 
403 
404  if(m_runConsistencyChecks) {
405  ATH_MSG_INFO("Running extra consistency check! (this is SLOW)");
406  if(!runConsistencyChecks()) {
407  ATH_MSG_ERROR("Consistency check has failed! Geometry is not consistent");
408  return StatusCode::FAILURE;
409  }
410  }
411 
412  ATH_MSG_INFO("Acts TrackingGeometry construction completed");
413 
414  return StatusCode::SUCCESS;
415 }

◆ makeBeamPipeConfig()

Acts::CylinderVolumeBuilder::Config ActsTrackingGeometrySvc::makeBeamPipeConfig ( std::shared_ptr< const Acts::CylinderVolumeHelper >  cvh) const
private

Definition at line 1024 of file ActsTrackingGeometrySvc.cxx.

1025  {
1026 
1027  // adapted from InnerDetector/InDetDetDescr/InDetTrackingGeometry/src/BeamPipeBuilder.cxx
1028 
1029  PVConstLink beamPipeTopVolume = p_beamPipeMgr->getTreeTop(0);
1030 
1031  if (p_beamPipeMgr->getNumTreeTops() == 1){
1032  beamPipeTopVolume = p_beamPipeMgr->getTreeTop(0)->getChildVol(0)->getChildVol(0);
1033  }
1034 
1035  Acts::Transform3 beamPipeTransform;
1036  beamPipeTransform.setIdentity();
1037 
1038  beamPipeTransform = Acts::Translation3(beamPipeTopVolume->getX().translation());
1039 
1040  double beamPipeRadius = 20;
1041 
1042  const GeoLogVol* beamPipeLogVolume = beamPipeTopVolume->getLogVol();
1043  const GeoTube* beamPipeTube = nullptr;
1044 
1045 
1046  if (beamPipeLogVolume == nullptr) {
1047  ATH_MSG_ERROR("Beam pip volume has no log volume");
1048  throw std::runtime_error("Beam pip volume has no log volume");
1049  }
1050  // get the geoShape and translate
1051  beamPipeTube = dynamic_cast<const GeoTube*>(beamPipeLogVolume->getShape());
1052  if (beamPipeTube == nullptr){
1053  ATH_MSG_ERROR("BeamPipeLogVolume was not of type GeoTube");
1054  throw std::runtime_error{"BeamPipeLogVolume was not of type GeoTube"};
1055  }
1056 
1057  for(unsigned int i=0;i<beamPipeTopVolume->getNChildVols();i++) {
1058 
1059  if(beamPipeTopVolume->getNameOfChildVol(i) == "SectionC03"){
1060 
1061  PVConstLink childTopVolume = beamPipeTopVolume->getChildVol(i);
1062  const GeoLogVol* childLogVolume = childTopVolume->getLogVol();
1063  const GeoTube* childTube = nullptr;
1064 
1065  if (childLogVolume){
1066  childTube = dynamic_cast<const GeoTube*>(childLogVolume->getShape());
1067  if (childTube){
1068  beamPipeRadius = 0.5 * (childTube->getRMax()+childTube->getRMin());
1069  }
1070  }
1071 
1072  break; // Exit loop after SectionC03 is found
1073  }
1074 
1075  } // Loop over child volumes
1076 
1077  ATH_MSG_VERBOSE("BeamPipe constructed from Database: translation (yes) - radius "
1078  << ( beamPipeTube ? "(yes)" : "(no)") << " - r = " << beamPipeRadius );
1079 
1080  ATH_MSG_VERBOSE("BeamPipe shift estimated as : " << Amg::toString(beamPipeTransform.translation()));
1081 
1082  Acts::CylinderVolumeBuilder::Config cfg;
1083 
1084  Acts::PassiveLayerBuilder::Config bplConfig;
1085  bplConfig.layerIdentification = "BeamPipe";
1086  bplConfig.centralLayerRadii = {beamPipeRadius * 1_mm};
1087  bplConfig.centralLayerHalflengthZ = {3000_mm};
1088  bplConfig.centralLayerThickness = {1_mm};
1089  auto beamPipeBuilder = std::make_shared<const Acts::PassiveLayerBuilder>(
1090  bplConfig, makeActsAthenaLogger(this, std::string("BPLayBldr"), std::string("ActsTGSvc")));
1091 
1092  // create the volume for the beam pipe
1093  cfg.trackingVolumeHelper = cvh;
1094  cfg.volumeName = "BeamPipe";
1095  cfg.layerBuilder = beamPipeBuilder;
1096  cfg.layerEnvelopeR = {1_mm, 1_mm};
1097  cfg.buildToRadiusZero = true;
1098 
1099  return cfg;
1100 }

◆ makeHGTDLayerBuilder()

std::shared_ptr< const Acts::ILayerBuilder > ActsTrackingGeometrySvc::makeHGTDLayerBuilder ( const HGTD_DetectorManager *  manager)
private

Definition at line 679 of file ActsTrackingGeometrySvc.cxx.

680  {
681 
682  std::string managerName = manager->getName();
683  auto matcher = [](const Acts::GeometryContext & /*gctx*/,
684  Acts::BinningValue /*bValue*/, const Acts::Surface * /*aS*/,
685  const Acts::Surface *
686  /*bS*/) -> bool { return false; };
687 
688  Acts::SurfaceArrayCreator::Config sacCfg;
689  sacCfg.surfaceMatcher = matcher;
690  sacCfg.doPhiBinningOptimization = false;
691 
692  auto surfaceArrayCreator = std::make_shared<Acts::SurfaceArrayCreator>(
693  sacCfg,
694  makeActsAthenaLogger(this, managerName + "SrfArrCrtr", std::string("ActsTGSvc")));
695  Acts::LayerCreator::Config lcCfg;
696  lcCfg.surfaceArrayCreator = surfaceArrayCreator;
697  auto layerCreator = std::make_shared<Acts::LayerCreator>(
698  lcCfg, makeActsAthenaLogger(this, managerName + "LayCrtr", std::string("ActsTGSvc")));
699 
700  ActsHGTDLayerBuilder::Config cfg;
701  cfg.mng = static_cast<const HGTD_DetectorManager *>(manager);
702  cfg.elementStore = m_elementStore;
703  cfg.layerCreator = layerCreator;
704  cfg.idHelper = m_HGTD_idHelper;
705  return std::make_shared<const ActsHGTDLayerBuilder>(
706  cfg, makeActsAthenaLogger(this, managerName + "GMSLayBldr", std::string("ActsTGSvc")));
707 }

◆ makeLayerBuilderConfig()

ActsLayerBuilder::Config ActsTrackingGeometrySvc::makeLayerBuilderConfig ( const InDetDD::InDetDetectorManager *  manager)
private

Definition at line 709 of file ActsTrackingGeometrySvc.cxx.

710  {
711  std::string managerName = manager->getName();
712 
713  std::shared_ptr<const Acts::ILayerBuilder> gmLayerBuilder;
714  auto matcher = [](const Acts::GeometryContext & /*gctx*/,
715  Acts::BinningValue bValue, const Acts::Surface *aS,
716  const Acts::Surface *bS) -> bool {
717  auto a = dynamic_cast<const ActsDetectorElement *>(
718  aS->associatedDetectorElement());
719  auto b = dynamic_cast<const ActsDetectorElement *>(
720  bS->associatedDetectorElement());
721  if ((not a) or (not b)) {
722  throw std::runtime_error(
723  "Cast of surface associated element to ActsDetectorElement failed "
724  "in ActsTrackingGeometrySvc::makeVolumeBuilder");
725  }
726 
727  IdentityHelper idA = a->identityHelper();
728  IdentityHelper idB = b->identityHelper();
729 
730  // check if same bec
731  // can't be same if not
732  if (idA.bec() != idB.bec())
733  return false;
734 
735  if (bValue == Acts::binPhi) {
736  // std::cout << idA.phi_module() << " <-> " << idB.phi_module() <<
737  // std::endl;
738  return idA.phi_module() == idB.phi_module();
739  }
740 
741  if (bValue == Acts::binZ) {
742  return (idA.eta_module() == idB.eta_module()) &&
743  (idA.layer_disk() == idB.layer_disk()) && (idA.bec() == idB.bec());
744  }
745 
746  if (bValue == Acts::binR) {
747  return (idA.eta_module() == idB.eta_module()) &&
748  (idA.layer_disk() == idB.layer_disk()) && (idB.bec() == idA.bec());
749  }
750 
751  return false;
752  };
753 
754  Acts::SurfaceArrayCreator::Config sacCfg;
755  sacCfg.surfaceMatcher = matcher;
756 
757  auto surfaceArrayCreator = std::make_shared<Acts::SurfaceArrayCreator>(
758  sacCfg,
759  makeActsAthenaLogger(this, managerName + "SrfArrCrtr", std::string("ActsTGSvc")));
760  Acts::LayerCreator::Config lcCfg;
761  lcCfg.surfaceArrayCreator = surfaceArrayCreator;
762  auto layerCreator = std::make_shared<Acts::LayerCreator>(
763  lcCfg, makeActsAthenaLogger(this, managerName + "LayCrtr", std::string("ActsTGSvc")));
764 
765  ActsLayerBuilder::Config cfg;
766  cfg.surfaceMatcher = matcher;
767 
768  // set bins from configuration
769  if (m_barrelMaterialBins.size() != 2) {
770  throw std::invalid_argument("Number of barrel material bin counts != 2");
771  }
772  std::vector<size_t> brlBins(m_barrelMaterialBins);
773  cfg.barrelMaterialBins = {brlBins.at(0), brlBins.at(1)};
774 
775  if (m_endcapMaterialBins.size() != 2) {
776  throw std::invalid_argument("Number of endcap material bin counts != 2");
777  }
778  std::vector<size_t> ecBins(m_endcapMaterialBins);
779  cfg.endcapMaterialBins = {ecBins.at(0), ecBins.at(1)};
780 
781  cfg.mng = static_cast<const InDetDD::SiDetectorManager *>(manager);
782  // use class member element store
783  cfg.elementStore = m_elementStore;
784  cfg.layerCreator = layerCreator;
785 
786  // gmLayerBuilder = std::make_shared<const ActsLayerBuilder>(
787  // cfg, makeActsAthenaLogger(this, managerName + "GMLayBldr",
788  // "ActsTGSvc"));
789 
790  // return gmLayerBuilder;
791  return cfg;
792 }

◆ makeSCTTRTAssembly()

std::shared_ptr< Acts::TrackingVolume > ActsTrackingGeometrySvc::makeSCTTRTAssembly ( const Acts::GeometryContext &  gctx,
const Acts::ILayerBuilder &  sct_lb,
const Acts::ILayerBuilder &  trt_lb,
const Acts::CylinderVolumeHelper &  cvh,
const std::shared_ptr< const Acts::TrackingVolume > &  pixel 
)
private

Definition at line 795 of file ActsTrackingGeometrySvc.cxx.

798  {
799  ATH_MSG_VERBOSE("Building SCT+TRT assembly");
800 
801  Acts::CylinderVolumeBuilder::Config cvbCfg;
802  Acts::CylinderVolumeBuilder cvb(
803  cvbCfg, makeActsAthenaLogger(this, std::string("SCTTRTCVB"), std::string("ActsTGSvc")));
804 
805  ATH_MSG_VERBOSE("Making SCT negative layers: ");
806  Acts::VolumeConfig sctNegEC =
807  cvb.analyzeContent(gctx, sct_lb.negativeLayers(gctx), {});
808  ATH_MSG_VERBOSE("Making SCT positive layers: ");
809  Acts::VolumeConfig sctPosEC =
810  cvb.analyzeContent(gctx, sct_lb.positiveLayers(gctx), {});
811  ATH_MSG_VERBOSE("Making SCT central layers: ");
812  Acts::VolumeConfig sctBrl =
813  cvb.analyzeContent(gctx, sct_lb.centralLayers(gctx), {});
814 
815  ATH_MSG_VERBOSE("Making TRT negative layers: ");
816  Acts::VolumeConfig trtNegEC =
817  cvb.analyzeContent(gctx, trt_lb.negativeLayers(gctx), {});
818  ATH_MSG_VERBOSE("Making TRT positive layers: ");
819  Acts::VolumeConfig trtPosEC =
820  cvb.analyzeContent(gctx, trt_lb.positiveLayers(gctx), {});
821  ATH_MSG_VERBOSE("Making TRT central layers: ");
822  Acts::VolumeConfig trtBrl =
823  cvb.analyzeContent(gctx, trt_lb.centralLayers(gctx), {});
824 
825  // synchronize trt
826 
827  double absZMinEC = std::min(std::abs(trtNegEC.zMax), std::abs(trtPosEC.zMin));
828  double absZMaxEC = std::max(std::abs(trtNegEC.zMin), std::abs(trtPosEC.zMax));
829 
830  trtNegEC.zMin = -absZMaxEC;
831  trtNegEC.zMax = -absZMinEC;
832  trtPosEC.zMin = absZMinEC;
833  trtPosEC.zMax = absZMaxEC;
834 
835  using CVBBV = Acts::CylinderVolumeBounds::BoundValues;
836 
837  // if pixel is present, shrink SCT volumes in R
838  bool isSCTSmallerInZ = false;
839  if (pixel) {
840  ATH_MSG_VERBOSE("Shrinking SCT in R (and maybe in increase size in Z) to fit around Pixel");
841  auto pixelBounds = dynamic_cast<const Acts::CylinderVolumeBounds *>(
842  &pixel->volumeBounds());
843  double sctNegECzMin = std::min(sctNegEC.zMin, -pixelBounds->get(CVBBV::eHalfLengthZ));
844  double sctPosECzMax = std::max(sctPosEC.zMax, pixelBounds->get(CVBBV::eHalfLengthZ));
845 
846  ATH_MSG_VERBOSE("- SCT +-EC.rMin: " << sctNegEC.rMin << " -> " << pixelBounds->get(CVBBV::eMaxR));
847  ATH_MSG_VERBOSE("- SCT BRL.rMin: " << sctBrl.rMin << " -> " << pixelBounds->get(CVBBV::eMaxR));
848  ATH_MSG_VERBOSE("- SCT EC.zMin: " << sctNegEC.zMin << " -> " << sctNegECzMin);
849  ATH_MSG_VERBOSE("- SCT EC.zMax: " << sctPosEC.zMax << " -> " << sctPosECzMax);
850 
851  sctNegEC.rMin = pixelBounds->get(CVBBV::eMaxR);
852  sctPosEC.rMin = pixelBounds->get(CVBBV::eMaxR);
853  sctBrl.rMin = pixelBounds->get(CVBBV::eMaxR);
854 
855  isSCTSmallerInZ = sctPosEC.zMax < pixelBounds->get(CVBBV::eHalfLengthZ);
856 
857  sctNegEC.zMin = sctNegECzMin;
858  sctPosEC.zMax = sctPosECzMax;
859 
860 
861  } else {
862  ATH_MSG_VERBOSE("Pixel is not configured, not wrapping");
863  }
864 
865  ATH_MSG_VERBOSE("SCT Volume Configuration:");
866  ATH_MSG_VERBOSE("- SCT::NegativeEndcap: " << sctNegEC.layers.size()
867  << " layers, "
868  << sctNegEC.toString());
869  ATH_MSG_VERBOSE("- SCT::Barrel: " << sctBrl.layers.size() << " layers, "
870  << sctBrl.toString());
871  ATH_MSG_VERBOSE("- SCT::PositiveEncap: " << sctPosEC.layers.size()
872  << " layers, "
873  << sctPosEC.toString());
874 
875  ATH_MSG_VERBOSE("TRT Volume Configuration:");
876  ATH_MSG_VERBOSE("- TRT::NegativeEndcap: " << trtNegEC.layers.size()
877  << " layers, "
878  << trtNegEC.toString());
879  ATH_MSG_VERBOSE("- TRT::Barrel: " << trtBrl.layers.size() << " layers, "
880  << trtBrl.toString());
881  ATH_MSG_VERBOSE("- TRT::PositiveEncap: " << trtPosEC.layers.size()
882  << " layers, "
883  << trtPosEC.toString());
884 
885  // harmonize SCT BRL <-> EC, normally the CVB does this, but we're skipping
886  // that
887  sctBrl.zMax = (sctBrl.zMax + sctPosEC.zMin) / 2.;
888  sctBrl.zMin = -sctBrl.zMax;
889 
890  // and now harmonize everything
891  // inflate TRT Barrel to match SCT
892  trtBrl.zMin = sctBrl.zMin;
893  trtBrl.zMax = sctBrl.zMax;
894 
895  // extend TRT endcaps outwards z so they match SCT
896  trtNegEC.zMin = sctNegEC.zMin;
897  trtPosEC.zMax = sctPosEC.zMax;
898 
899  // extend endcap in z so it touches barrel
900  trtNegEC.zMax = trtBrl.zMin;
901  sctNegEC.zMax = trtBrl.zMin;
902  trtPosEC.zMin = trtBrl.zMax;
903  sctPosEC.zMin = trtBrl.zMax;
904 
905  // extend SCT in R so they touch TRT barel
906  sctBrl.rMax = trtBrl.rMin;
907  sctNegEC.rMax = trtNegEC.rMin;
908  sctPosEC.rMax = trtPosEC.rMin;
909 
910  // extend TRT endcaps in r to that of Barrel
911  trtNegEC.rMax = trtBrl.rMax;
912  trtPosEC.rMax = trtBrl.rMax;
913 
914  ATH_MSG_VERBOSE("Dimensions after synchronization between SCT and TRT");
915  ATH_MSG_VERBOSE("SCT Volume Configuration:");
916  ATH_MSG_VERBOSE("- SCT::NegativeEndcap: " << sctNegEC.layers.size()
917  << " layers, "
918  << sctNegEC.toString());
919  ATH_MSG_VERBOSE("- SCT::Barrel: " << sctBrl.layers.size() << " layers, "
920  << sctBrl.toString());
921  ATH_MSG_VERBOSE("- SCT::PositiveEncap: " << sctPosEC.layers.size()
922  << " layers, "
923  << sctPosEC.toString());
924 
925  ATH_MSG_VERBOSE("TRT Volume Configuration:");
926  ATH_MSG_VERBOSE("- TRT::NegativeEndcap: " << trtNegEC.layers.size()
927  << " layers, "
928  << trtNegEC.toString());
929  ATH_MSG_VERBOSE("- TRT::Barrel: " << trtBrl.layers.size() << " layers, "
930  << trtBrl.toString());
931  ATH_MSG_VERBOSE("- TRT::PositiveEncap: " << trtPosEC.layers.size()
932  << " layers, "
933  << trtPosEC.toString());
934 
935  auto makeTVol = [&](const auto &vConf, const auto &name) {
936  return cvh.createTrackingVolume(gctx, vConf.layers, {},
937  nullptr, // no material
938  vConf.rMin, vConf.rMax, vConf.zMin,
939  vConf.zMax, name);
940  };
941 
942  // now turn them into actual TrackingVolumes
943  auto tvSctNegEC = makeTVol(sctNegEC, "SCT::NegativeEndcap");
944  auto tvSctBrl = makeTVol(sctBrl, "SCT::Barrel");
945  auto tvSctPosEC = makeTVol(sctPosEC, "SCT::PositiveEndcap");
946 
947  auto tvTrtNegEC = makeTVol(trtNegEC, "TRT::NegativeEndcap");
948  auto tvTrtBrl = makeTVol(trtBrl, "TRT::Barrel");
949  auto tvTrtPosEC = makeTVol(trtPosEC, "TRT::PositiveEndcap");
950 
951  // combine the endcaps and the barrels, respetively
952  auto negEC =
953  cvh.createContainerTrackingVolume(gctx, {tvSctNegEC, tvTrtNegEC});
954  auto posEC =
955  cvh.createContainerTrackingVolume(gctx, {tvSctPosEC, tvTrtPosEC});
956  auto barrel = cvh.createContainerTrackingVolume(gctx, {tvSctBrl, tvTrtBrl});
957 
958  // and now combine all of those into one container for the assembly
959 
960  auto container =
961  cvh.createContainerTrackingVolume(gctx, {negEC, barrel, posEC});
962 
963  // if pixel is present, add positive and negative gap volumes so we can wrap
964  // it all
965  if (pixel) {
966  auto containerBounds = dynamic_cast<const Acts::CylinderVolumeBounds *>(
967  &container->volumeBounds());
968  auto pixelBounds = dynamic_cast<const Acts::CylinderVolumeBounds *>(
969  &pixel->volumeBounds());
970  std::vector<std::shared_ptr<Acts::TrackingVolume>> noVolumes;
971 
972  if(!isSCTSmallerInZ) {
973  // pixel is smaller in z, need gap volumes
974  auto posGap = cvh.createGapTrackingVolume(
975  gctx, noVolumes,
976  nullptr, // no material,
977  pixelBounds->get(CVBBV::eMinR), pixelBounds->get(CVBBV::eMaxR),
978  pixelBounds->get(CVBBV::eHalfLengthZ),
979  containerBounds->get(CVBBV::eHalfLengthZ),
980  0, // material layers,
981  true, // cylinder
982  "Pixel::PositiveGap");
983  auto negGap = cvh.createGapTrackingVolume(
984  gctx, noVolumes,
985  nullptr, // no material,
986  pixelBounds->get(CVBBV::eMinR), pixelBounds->get(CVBBV::eMaxR),
987  -containerBounds->get(CVBBV::eHalfLengthZ),
988  -pixelBounds->get(CVBBV::eHalfLengthZ),
989  0, // material layers,
990  true, // cylinder
991  "Pixel::NegativeGap");
992 
993  auto pixelContainer =
994  cvh.createContainerTrackingVolume(gctx, {negGap, pixel, posGap});
995  // and now create one container that contains Pixel+SCT+TRT
996  container =
997  cvh.createContainerTrackingVolume(gctx, {pixelContainer, container});
998  }
999  else {
1000  // wrap the pixel directly
1001  container =
1002  cvh.createContainerTrackingVolume(gctx, {pixel, container});
1003  }
1004 
1005  }
1006 
1007  return container;
1008 }

◆ makeStrawLayerBuilder()

std::shared_ptr< const Acts::ILayerBuilder > ActsTrackingGeometrySvc::makeStrawLayerBuilder ( const InDetDD::InDetDetectorManager *  manager)
private

Definition at line 648 of file ActsTrackingGeometrySvc.cxx.

649  {
650 
651  std::string managerName = manager->getName();
652  auto matcher = [](const Acts::GeometryContext & /*gctx*/,
653  Acts::BinningValue /*bValue*/, const Acts::Surface * /*aS*/,
654  const Acts::Surface *
655  /*bS*/) -> bool { return false; };
656 
657  Acts::SurfaceArrayCreator::Config sacCfg;
658  sacCfg.surfaceMatcher = matcher;
659  sacCfg.doPhiBinningOptimization = false;
660 
661  auto surfaceArrayCreator = std::make_shared<Acts::SurfaceArrayCreator>(
662  sacCfg,
663  makeActsAthenaLogger(this, managerName + "SrfArrCrtr", std::string("ActsTGSvc")));
664  Acts::LayerCreator::Config lcCfg;
665  lcCfg.surfaceArrayCreator = surfaceArrayCreator;
666  auto layerCreator = std::make_shared<Acts::LayerCreator>(
667  lcCfg, makeActsAthenaLogger(this, managerName + "LayCrtr", std::string("ActsTGSvc")));
668 
669  ActsStrawLayerBuilder::Config cfg;
670  cfg.mng = static_cast<const InDetDD::TRT_DetectorManager *>(manager);
671  cfg.elementStore = m_elementStore;
672  cfg.layerCreator = layerCreator;
673  cfg.idHelper = m_TRT_idHelper;
674  return std::make_shared<const ActsStrawLayerBuilder>(
675  cfg, makeActsAthenaLogger(this, managerName + "GMSLayBldr", std::string("ActsTGSvc")));
676 }

◆ populateAlignmentStore()

unsigned int ActsTrackingGeometrySvc::populateAlignmentStore ( ActsTrk::DetectorAlignStore &  store) const
override

Definition at line 1010 of file ActsTrackingGeometrySvc.cxx.

1010  {
1011  ATH_MSG_DEBUG("Populate the alignment store with all detector elements");
1012  unsigned int nElements = 0;
1013  m_trackingGeometry->visitSurfaces([&store, &nElements](const Acts::Surface *srf) {
1014  const Acts::DetectorElementBase *detElem = srf->associatedDetectorElement();
1015  const IDetectorElement *gmde = dynamic_cast<const IDetectorElement *>(detElem);
1016  nElements += gmde->storeAlignedTransforms(store);
1017  });
1018  ATH_MSG_DEBUG("Populated with " << nElements << " elements");
1019  return nElements;
1020 }

◆ runConsistencyChecks()

bool ActsTrackingGeometrySvc::runConsistencyChecks ( ) const
private

Definition at line 417 of file ActsTrackingGeometrySvc.cxx.

417  {
418  bool result = true;
419 
420  std::vector<Acts::Vector2> localPoints;
421  localPoints.reserve(m_consistencyCheckPoints);
422  std::mt19937 gen;
423  std::uniform_real_distribution<> dist(0.0, 1.0);
424 
425  std::optional<std::ofstream> os;
426  if(!m_consistencyCheckOutput.empty()){
427  os = std::ofstream{m_consistencyCheckOutput};
428  if(!os->good()) {
429  throw std::runtime_error{"Failed to open consistency check output file"};
430  }
431  }
432 
433  if(os) {
434  (*os) << "geo_id,vol_id,lay_id,sen_id,type,acts_loc0,acts_loc1,acts_inside,trk_loc0,trk_loc1,trk_inside,x,y,z,g2l_loc0,g2l_loc1,trk_x,trk_y,trk_z" << std::endl;
435  }
436  for(size_t i=0;i<m_consistencyCheckPoints;i++) {
437  localPoints.emplace_back(dist(gen), dist(gen));
438  }
439 
440  Acts::GeometryContext gctx = getNominalContext().context();
441 
442  size_t nTotalSensors = 0;
443  std::array<size_t,3> nInconsistent{0,0,0};
444  size_t nMismatchedCenters = 0;
445  size_t nMismatchedNormals = 0;
446 
447  // Comparison of Eigen vectors, similar to a.isApprox(b), but use absolute comparison to also work with zero vectors.
448  // All values will be mm or radians, so 1e-5 is a reasonable precision.
449  auto isApprox = [](auto& a, auto& b) -> bool {
450  return ((a - b).array().abs() < 1e-5).all();
451  };
452 
453  m_trackingGeometry->visitSurfaces([&](const Acts::Surface *surface) {
454  nTotalSensors++;
455 
456  const auto* actsDetElem = dynamic_cast<const ActsDetectorElement*>(surface->associatedDetectorElement());
457  if(actsDetElem == nullptr) {
458  ATH_MSG_ERROR("Invalid detector element found");
459  result = false;
460  return;
461  }
462  const auto* siDetElem = dynamic_cast<const InDetDD::SiDetectorElement*>(actsDetElem->upstreamDetectorElement());
463  if(siDetElem == nullptr) {
464  return;
465  }
466 
467  const auto* regSurface = dynamic_cast<const Acts::RegularSurface*>(surface);
468  const auto& trkSurface = siDetElem->surface();
469  if(regSurface == nullptr) {
470  ATH_MSG_ERROR("Invalid surface found");
471  result = false;
472  return;
473  }
474 
475  Acts::Vector3 center{regSurface->center(gctx)};
476  Amg::Vector3D trkCenter{trkSurface.center()};
477  if (/* auto *b = */ dynamic_cast<const Acts::AnnulusBounds *>(&surface->bounds()))
478  {
479  // // Acts::AnnulusBounds defines center() as center of whole disc, so get it from the bounds
480  // Acts::Vector2 locCenter{0.5 * (b->rMin() + b->rMax()), 0.5 * (b->phiMin() + b->phiMax())};
481  // center = surface->localToGlobal(gctx, locCenter, Acts::Vector3::Zero());
482  center.head<2>() = trkCenter.head<2>(); // that doesn't (quite) work for xy, so just pass that check
483  }
484 
485  if(!isApprox(trkCenter, center)) {
486  std::string trkName;
487  if (auto idHelper = siDetElem->getIdHelper())
488  {
489  trkName = idHelper->show_to_string(siDetElem->identify());
490  }
491  ATH_MSG_WARNING("Acts surface "
492  << surface->geometryId()
493  << " center (" << center[0] << ',' << center[1] << ',' << center[2]
494  << ") does not match Trk surface " << trkName
495  << " center (" << trkCenter[0] << ',' << trkCenter[1] << ',' << trkCenter[2] << ')');
496  nMismatchedCenters++;
497  result = false;
498  }
499 
500  const auto* lineSurface = dynamic_cast<const Acts::LineSurface*>(surface);
501  if(lineSurface == nullptr) {
502  Acts::Vector3 norm{regSurface->normal(gctx, regSurface->center(gctx))};
503  Amg::Vector3D trkNorm{trkSurface.normal()};
504  if(!isApprox(trkNorm, norm)) {
505  std::string trkName;
506  if (auto idHelper = siDetElem->getIdHelper())
507  {
508  trkName = idHelper->show_to_string(siDetElem->identify());
509  }
510  ATH_MSG_WARNING("Acts surface "
511  << surface->geometryId()
512  << " normal (" << norm[0] << ',' << norm[1] << ',' << norm[2]
513  << ") does not match Trk surface " << trkName
514  << " normal (" << trkNorm[0] << ',' << trkNorm[1] << ',' << trkNorm[2] << ')');
515  nMismatchedNormals++;
516  result = false;
517  }
518  }
519 
520  auto doPoints = [&](unsigned int type, const Acts::Vector2& loc) -> std::array<bool,3> {
521  Acts::Vector3 glb = surface->localToGlobal(gctx, loc, Acts::Vector3::Zero());
522 
523  Amg::Vector2D locTrk = Amg::Vector2D::Zero();
524  Amg::Vector3D glbTrk = Amg::Vector3D::Zero();
525  Acts::Vector2 locg2l = Acts::Vector2::Zero();
526  bool locg2lOk = false;
527  auto locTrkRes = trkSurface.globalToLocal(glb);
528  if (locTrkRes) {
529  locTrk = locTrkRes.value();
530  glbTrk = trkSurface.localToGlobal(locTrk);
531 
532  auto locg2lRes = surface->globalToLocal(gctx, glbTrk, Acts::Vector3::Zero());
533  if (locg2lRes.ok()) {
534  locg2lOk = true;
535  locg2l = locg2lRes.value();
536  }
537  }
538 
539  auto gId = surface->geometryId();
540  if(os) {
541  (*os) << gId.value()
542  << "," << gId.volume()
543  << "," << gId.layer()
544  << "," << gId.sensitive()
545  << "," << type
546  << "," << loc[0]
547  << "," << loc[1]
548  << "," << surface->insideBounds(loc)
549  << "," << locTrk[0]
550  << "," << locTrk[1]
551  << "," << trkSurface.insideBounds(locTrk)
552  << "," << glb[0]
553  << "," << glb[1]
554  << "," << glb[2]
555  << "," << locg2l[0]
556  << "," << locg2l[1]
557  << "," << glbTrk[0]
558  << "," << glbTrk[1]
559  << "," << glbTrk[2]
560  << std::endl;
561  }
562 
563  return {surface->insideBounds(loc) == trkSurface.insideBounds(locTrk),
564  locg2lOk ? isApprox(loc, locg2l) : true,
565  locTrkRes ? isApprox(glb, glbTrk) : true};
566  };
567 
568  std::array<bool,3> allOk{true,true,true};
569  if(const auto* bounds = dynamic_cast<const Acts::PlanarBounds*>(&surface->bounds()); bounds) {
570  ATH_MSG_VERBOSE("Planar bounds");
571 
572  const Acts::RectangleBounds& boundingBox = bounds->boundingBox();
573  Acts::Vector2 min = boundingBox.min().array() - 1*Acts::UnitConstants::mm;
574  Acts::Vector2 max = boundingBox.max().array() - 1*Acts::UnitConstants::mm;
575  Acts::Vector2 diag = max - min;
576 
577  for(const auto& testPoint : localPoints) {
578  Acts::Vector2 loc = min.array() + (testPoint.array() * diag.array());
579  auto pointOk = doPoints(0, loc);
580  for (size_t i=0; i<pointOk.size(); ++i) {
581  if (!pointOk[i]) {
582  result = false;
583  allOk[i] = false;
584  }
585  }
586  }
587 
588  }
589  else if(const auto* bounds = dynamic_cast<const Acts::AnnulusBounds*>(&surface->bounds()); bounds) {
590  ATH_MSG_VERBOSE("Annulus bounds");
591 
592  // custom bounding box algo
593  std::vector<Acts::Vector2> vertices = bounds->vertices(5); // 5 segments on the radial edges
594  Acts::Vector2 min{std::numeric_limits<double>::max(), std::numeric_limits<double>::max()};
595  Acts::Vector2 max{std::numeric_limits<double>::lowest(), std::numeric_limits<double>::lowest()};
596  for (const auto& vtx : vertices) {
597  min = min.array().min(vtx.array());
598  max = max.array().max(vtx.array());
599  }
600  min.array() -= 1*Acts::UnitConstants::mm;
601  max.array() += 1*Acts::UnitConstants::mm;
602  Acts::Vector2 diag = max - min;
603 
604  for(const auto& testPoint : localPoints) {
605  Acts::Vector2 locXY = min.array() + (testPoint.array() * diag.array());
606  Acts::Vector2 locPC = dynamic_cast<const Acts::DiscSurface&>(*surface).localCartesianToPolar(locXY);
607 
608  auto pointOk = doPoints(1, locPC);
609  for (size_t i=0; i<pointOk.size(); ++i) {
610  if (!pointOk[i]) {
611  result = false;
612  allOk[i] = false;
613  }
614  }
615  }
616 
617  }
618  else {
619  result = false;
620  }
621 
622  for (size_t i=0; i<allOk.size(); ++i) {
623  if (!allOk[i]) {
624  ++nInconsistent[i];
625  }
626  }
627 
628  });
629 
630  ATH_MSG_INFO("Total number of sensors : " << nTotalSensors);
631  ATH_MSG_INFO("Number of sensors with mismatched centers : " << nMismatchedCenters);
632  ATH_MSG_INFO("Number of sensors with mismatched normals : " << nMismatchedNormals);
633  ATH_MSG_INFO("Number of sensors with inconsistent inside: " << nInconsistent[0]);
634  ATH_MSG_INFO("Number of sensors with inconsistent g2l : " << nInconsistent[1]);
635  ATH_MSG_INFO("Number of sensors with inconsistent l2g : " << nInconsistent[2]);
636 
637  return result;
638 }

◆ trackingGeometry()

std::shared_ptr< const Acts::TrackingGeometry > ActsTrackingGeometrySvc::trackingGeometry ( )
override

Definition at line 641 of file ActsTrackingGeometrySvc.cxx.

641  {
642 
643  ATH_MSG_VERBOSE("Retrieving tracking geometry");
644  return m_trackingGeometry;
645 }

Member Data Documentation

◆ m_barrelMaterialBins

Gaudi::Property<std::vector<size_t> > ActsTrackingGeometrySvc::m_barrelMaterialBins {this, "BarrelMaterialBins", {10, 10}}
private

Definition at line 108 of file ActsTrackingGeometrySvc.h.

◆ m_buildBeamPipe

Gaudi::Property<bool> ActsTrackingGeometrySvc::m_buildBeamPipe {this, "BuildBeamPipe", false, ""}
private

Definition at line 106 of file ActsTrackingGeometrySvc.h.

◆ m_buildSubdetectors

Gaudi::Property<std::vector<std::string> > ActsTrackingGeometrySvc::m_buildSubdetectors {this, "BuildSubDetectors", {"Pixel", "SCT", "TRT", "Calo", "HGTD"}}
private

Definition at line 110 of file ActsTrackingGeometrySvc.h.

◆ m_caloVolumeBuilder

ToolHandle<IActsTrackingVolumeBuilder> ActsTrackingGeometrySvc::m_caloVolumeBuilder
private
Initial value:
{this,
"CaloVolumeBuilder", "", "CaloVolumeBuilder"}

Definition at line 121 of file ActsTrackingGeometrySvc.h.

◆ m_consistencyCheckOutput

StringProperty ActsTrackingGeometrySvc::m_consistencyCheckOutput
private
Initial value:
{this, "ConsistencyCheckOutput",
"", "Output file for geometry debugging, will not write if empty",}

Definition at line 115 of file ActsTrackingGeometrySvc.h.

◆ m_consistencyCheckPoints

Gaudi::Property<size_t> ActsTrackingGeometrySvc::m_consistencyCheckPoints
private
Initial value:
{this, "ConsistencyCheckPoints",
1000, "number of random points for consistency check"}

Definition at line 118 of file ActsTrackingGeometrySvc.h.

◆ m_detStore

ServiceHandle<StoreGateSvc> ActsTrackingGeometrySvc::m_detStore
private

Definition at line 85 of file ActsTrackingGeometrySvc.h.

◆ m_elementStore

std::shared_ptr<ActsElementVector> ActsTrackingGeometrySvc::m_elementStore {nullptr}
private

Definition at line 94 of file ActsTrackingGeometrySvc.h.

◆ m_endcapMaterialBins

Gaudi::Property<std::vector<size_t> > ActsTrackingGeometrySvc::m_endcapMaterialBins {this, "EndcapMaterialBins", {5, 20}}
private

Definition at line 109 of file ActsTrackingGeometrySvc.h.

◆ m_HGTD_idHelper

const HGTD_ID* ActsTrackingGeometrySvc::m_HGTD_idHelper {nullptr}
private

Definition at line 98 of file ActsTrackingGeometrySvc.h.

◆ m_materialMapCalibFolder

Gaudi::Property<std::string> ActsTrackingGeometrySvc::m_materialMapCalibFolder {this, "MaterialMapCalibFolder", ".", ""}
private

Definition at line 105 of file ActsTrackingGeometrySvc.h.

◆ m_materialMapInputFileBase

Gaudi::Property<std::string> ActsTrackingGeometrySvc::m_materialMapInputFileBase {this, "MaterialMapInputFile", "", ""}
private

Definition at line 104 of file ActsTrackingGeometrySvc.h.

◆ m_nominalContext

ActsGeometryContext ActsTrackingGeometrySvc::m_nominalContext {}
private

Definition at line 100 of file ActsTrackingGeometrySvc.h.

◆ m_objDebugOutput

Gaudi::Property<bool> ActsTrackingGeometrySvc::m_objDebugOutput {this, "ObjDebugOutput", false, ""}
private

Definition at line 103 of file ActsTrackingGeometrySvc.h.

◆ m_runConsistencyChecks

BooleanProperty ActsTrackingGeometrySvc::m_runConsistencyChecks
private
Initial value:
{this, "RunConsistencyChecks",
false, "Run extra consistency checks w.r.t to Trk::. This is SLOW!"}

Definition at line 112 of file ActsTrackingGeometrySvc.h.

◆ m_subDetNoAlign

std::set<ActsTrk::DetectorType> ActsTrackingGeometrySvc::m_subDetNoAlign {}
private

Definition at line 125 of file ActsTrackingGeometrySvc.h.

◆ m_subDetNoAlignProp

Gaudi::Property<std::vector<unsigned int> > ActsTrackingGeometrySvc::m_subDetNoAlignProp {this, "NotAlignDetectors", {}}
private

Define the subdetectors for which the tracking geometry does not expect a valid alignment store.

Definition at line 124 of file ActsTrackingGeometrySvc.h.

◆ m_trackingGeometry

std::shared_ptr<const Acts::TrackingGeometry> ActsTrackingGeometrySvc::m_trackingGeometry {nullptr}
private

Definition at line 95 of file ActsTrackingGeometrySvc.h.

◆ m_TRT_idHelper

const TRT_ID* ActsTrackingGeometrySvc::m_TRT_idHelper {nullptr}
private

Definition at line 97 of file ActsTrackingGeometrySvc.h.

◆ m_useMaterialMap

Gaudi::Property<bool> ActsTrackingGeometrySvc::m_useMaterialMap {this, "UseMaterialMap", false, ""}
private

Definition at line 102 of file ActsTrackingGeometrySvc.h.

◆ p_beamPipeMgr

const BeamPipeDetectorManager* ActsTrackingGeometrySvc::p_beamPipeMgr {nullptr}
private

Definition at line 91 of file ActsTrackingGeometrySvc.h.

◆ p_HGTDManager

const HGTD_DetectorManager* ActsTrackingGeometrySvc::p_HGTDManager {nullptr}
private

Definition at line 92 of file ActsTrackingGeometrySvc.h.

◆ p_ITkPixelManager

const InDetDD::SiDetectorManager* ActsTrackingGeometrySvc::p_ITkPixelManager {nullptr}
private

Definition at line 89 of file ActsTrackingGeometrySvc.h.

◆ p_ITkStripManager

const InDetDD::SiDetectorManager* ActsTrackingGeometrySvc::p_ITkStripManager {nullptr}
private

Definition at line 90 of file ActsTrackingGeometrySvc.h.

◆ p_pixelManager

const InDetDD::SiDetectorManager* ActsTrackingGeometrySvc::p_pixelManager {nullptr}
private

Definition at line 86 of file ActsTrackingGeometrySvc.h.

◆ p_SCTManager

const InDetDD::SiDetectorManager* ActsTrackingGeometrySvc::p_SCTManager {nullptr}
private

Definition at line 87 of file ActsTrackingGeometrySvc.h.

◆ p_TRTManager

const InDetDD::TRT_DetectorManager* ActsTrackingGeometrySvc::p_TRTManager {nullptr}
private

Definition at line 88 of file ActsTrackingGeometrySvc.h.

The documentation for this class was generated from the following files:
ActsTrackingGeometrySvc::makeBeamPipeConfig
Acts::CylinderVolumeBuilder::Config makeBeamPipeConfig(std::shared_ptr< const Acts::CylinderVolumeHelper > cvh) const
Definition: ActsTrackingGeometrySvc.cxx:1024
IdentityHelper::phi_module
int phi_module() const
Definition: IdentityHelper.cxx:61
store
StoreGateSvc * store
Definition: fbtTestBasics.cxx:69
PlotCalibFromCool.norm
norm
Definition: PlotCalibFromCool.py:100
ATH_MSG_FATAL
#define ATH_MSG_FATAL(x)
Definition: AthMsgStreamMacros.h:34
python.SystemOfUnits.s
int s
Definition: SystemOfUnits.py:131
ActsTrackingGeometrySvc::m_runConsistencyChecks
BooleanProperty m_runConsistencyChecks
Definition: ActsTrackingGeometrySvc.h:112
get_generator_info.result
result
Definition: get_generator_info.py:21
max
#define max(a, b)
Definition: cfImp.cxx:41
ActsLayerBuilder::Mode::ITkPixelOuter
@ ITkPixelOuter
ATH_MSG_INFO
#define ATH_MSG_INFO(x)
Definition: AthMsgStreamMacros.h:31
ActsTrackingGeometrySvc::m_consistencyCheckOutput
StringProperty m_consistencyCheckOutput
Definition: ActsTrackingGeometrySvc.h:115
ActsLayerBuilder::Mode::ITkPixelInner
@ ITkPixelInner
ActsTrk::DetectorType
DetectorType
Simple enum to Identify the Type of the ACTS sub detector.
Definition: GeometryDefs.h:17
Trk::binZ
@ binZ
Definition: BinningType.h:49
Amg::Vector2D
Eigen::Matrix< double, 2, 1 > Vector2D
Definition: GeoPrimitives.h:48
ActsTrackingGeometrySvc::m_barrelMaterialBins
Gaudi::Property< std::vector< size_t > > m_barrelMaterialBins
Definition: ActsTrackingGeometrySvc.h:108
ActsHGTDLayerBuilder::Config
Definition: ActsHGTDLayerBuilder.h:45
IdentityHelper::layer_disk
int layer_disk() const
Definition: IdentityHelper.cxx:49
BeamPipeDetectorManager::getNumTreeTops
virtual unsigned int getNumTreeTops() const
Definition: BeamPipeDetectorManager.cxx:15
InDetDD::SolidStateDetectorElementBase::surface
Trk::Surface & surface()
Element Surface.
ActsTrackingGeometrySvc::runConsistencyChecks
bool runConsistencyChecks() const
Definition: ActsTrackingGeometrySvc.cxx:417
ActsTrk::actsLevelVector
Acts::Logging::Level actsLevelVector(MSG::Level lvl)
Definition: LoggerUtils.cxx:9
ActsStrawLayerBuilder::Config
Definition: ActsStrawLayerBuilder.h:35
ActsTrackingGeometrySvc::m_consistencyCheckPoints
Gaudi::Property< size_t > m_consistencyCheckPoints
Definition: ActsTrackingGeometrySvc.h:118
ATH_MSG_VERBOSE
#define ATH_MSG_VERBOSE(x)
Definition: AthMsgStreamMacros.h:28
ActsTrackingGeometrySvc::m_detStore
ServiceHandle< StoreGateSvc > m_detStore
Definition: ActsTrackingGeometrySvc.h:85
ActsTrackingGeometrySvc::p_ITkPixelManager
const InDetDD::SiDetectorManager * p_ITkPixelManager
Definition: ActsTrackingGeometrySvc.h:89
BeamPipeDetectorManager::getTreeTop
virtual PVConstLink getTreeTop(unsigned int i) const
Definition: BeamPipeDetectorManager.cxx:20
ActsLayerBuilder::Config
Definition: ActsLayerBuilder.h:50
ActsTrackingGeometrySvc::p_beamPipeMgr
const BeamPipeDetectorManager * p_beamPipeMgr
Definition: ActsTrackingGeometrySvc.h:91
ActsTrackingGeometrySvc::m_caloVolumeBuilder
ToolHandle< IActsTrackingVolumeBuilder > m_caloVolumeBuilder
Definition: ActsTrackingGeometrySvc.h:121
ActsLayerBuilder::Mode::Pixel
@ Pixel
master.gen
gen
Definition: master.py:32
python.iconfTool.models.loaders.level
level
Definition: loaders.py:20
ActsTrackingGeometrySvc::getNominalContext
const ActsGeometryContext & getNominalContext() const override
Definition: ActsTrackingGeometrySvc.cxx:1021
ActsTrackingGeometrySvc::p_TRTManager
const InDetDD::TRT_DetectorManager * p_TRTManager
Definition: ActsTrackingGeometrySvc.h:88
makeActsAthenaLogger
std::unique_ptr< const Acts::Logger > makeActsAthenaLogger(IMessageSvc *svc, const std::string &name, int level, std::optional< std::string > parent_name)
Definition: Tracking/Acts/ActsInterop/src/Logger.cxx:64
ActsGeometryContext::context
Acts::GeometryContext context() const
Definition: ActsGeometryContext.h:46
ActsTrackingGeometrySvc::p_HGTDManager
const HGTD_DetectorManager * p_HGTDManager
Definition: ActsTrackingGeometrySvc.h:92
ActsLayerBuilder::Mode::ITkStrip
@ ITkStrip
Amg::toString
std::string toString(const Translation3D &translation, int precision=4)
GeoPrimitvesToStringConverter.
Definition: GeoPrimitivesToStringConverter.h:40
ATH_MSG_ERROR
#define ATH_MSG_ERROR(x)
Definition: AthMsgStreamMacros.h:33
Trk::BinningValue
BinningValue
how to take the global / local position
Definition: BinningType.h:46
python.BunchSpacingUtils.lb
lb
Definition: BunchSpacingUtils.py:88
lumiFormat.i
int i
Definition: lumiFormat.py:92
ActsTrackingGeometrySvc::m_elementStore
std::shared_ptr< ActsElementVector > m_elementStore
Definition: ActsTrackingGeometrySvc.h:94
ActsTrackingGeometrySvc::m_materialMapInputFileBase
Gaudi::Property< std::string > m_materialMapInputFileBase
Definition: ActsTrackingGeometrySvc.h:104
ATH_MSG_DEBUG
#define ATH_MSG_DEBUG(x)
Definition: AthMsgStreamMacros.h:29
ActsTrackingGeometrySvc::p_pixelManager
const InDetDD::SiDetectorManager * p_pixelManager
Definition: ActsTrackingGeometrySvc.h:86
calibdata.exception
exception
Definition: calibdata.py:496
ActsTrackingGeometrySvc::p_SCTManager
const InDetDD::SiDetectorManager * p_SCTManager
Definition: ActsTrackingGeometrySvc.h:87
ActsTrackingGeometrySvc::m_buildBeamPipe
Gaudi::Property< bool > m_buildBeamPipe
Definition: ActsTrackingGeometrySvc.h:106
ActsTrackingGeometrySvc::makeLayerBuilderConfig
ActsLayerBuilder::Config makeLayerBuilderConfig(const InDetDD::InDetDetectorManager *manager)
Definition: ActsTrackingGeometrySvc.cxx:709
ATH_CHECK
#define ATH_CHECK
Definition: AthCheckMacros.h:40
ActsDetectorElement
Definition: ActsDetectorElement.h:42
Handler::svc
AthROOTErrorHandlerSvc * svc
Definition: AthROOTErrorHandlerSvc.cxx:10
ActsTrk::IDetectorElement
Definition: IDetectorElement.h:28
ActsTrk::IDetectorElement::storeAlignedTransforms
virtual unsigned int storeAlignedTransforms(const ActsTrk::DetectorAlignStore &store) const =0
Caches the aligned transformation in the provided store. Returns the number of cached elements.
ActsTrackingGeometrySvc::m_endcapMaterialBins
Gaudi::Property< std::vector< size_t > > m_endcapMaterialBins
Definition: ActsTrackingGeometrySvc.h:109
ReadFromCoolCompare.os
os
Definition: ReadFromCoolCompare.py:231
ActsTrackingGeometrySvc::m_materialMapCalibFolder
Gaudi::Property< std::string > m_materialMapCalibFolder
Definition: ActsTrackingGeometrySvc.h:105
min
#define min(a, b)
Definition: cfImp.cxx:40
IdentityHelper::eta_module
int eta_module() const
Definition: IdentityHelper.cxx:73
name
std::string name
Definition: Control/AthContainers/Root/debug.cxx:195
plotBeamSpotMon.b
b
Definition: plotBeamSpotMon.py:77
ActsTrackingGeometrySvc::makeSCTTRTAssembly
std::shared_ptr< Acts::TrackingVolume > makeSCTTRTAssembly(const Acts::GeometryContext &gctx, const Acts::ILayerBuilder &sct_lb, const Acts::ILayerBuilder &trt_lb, const Acts::CylinderVolumeHelper &cvh, const std::shared_ptr< const Acts::TrackingVolume > &pixel)
Definition: ActsTrackingGeometrySvc.cxx:795
InDetDD::SiDetectorElement
Definition: SiDetectorElement.h:109
WriteCaloSwCorrections.cfg
cfg
Definition: WriteCaloSwCorrections.py:23
Amg::Vector3D
Eigen::Matrix< double, 3, 1 > Vector3D
Definition: GeoPrimitives.h:47
ActsTrackingGeometrySvc::p_ITkStripManager
const InDetDD::SiDetectorManager * p_ITkStripManager
Definition: ActsTrackingGeometrySvc.h:90
CVBBV
Acts::CylinderVolumeBounds::BoundValues CVBBV
Definition: ActsCaloTrackingVolumeBuilder.cxx:31
PathResolverFindCalibFile
std::string PathResolverFindCalibFile(const std::string &logical_file_name)
Definition: PathResolver.cxx:431
IdentityHelper::bec
int bec() const
Definition: IdentityHelper.cxx:37
python.SystemOfUnits.mm
int mm
Definition: SystemOfUnits.py:83
Cut::all
@ all
Definition: SUSYToolsAlg.cxx:64
Trk::binR
@ binR
Definition: BinningType.h:50
InDetDD::TRT_DetectorManager
The Detector Manager for all TRT Detector elements, it acts as the interface to the detector elements...
Definition: TRT_DetectorManager.h:69
DiTauMassTools::MaxHistStrategyV2::e
e
Definition: PhysicsAnalysis/TauID/DiTauMassTools/DiTauMassTools/HelperFunctions.h:26
a
TList * a
Definition: liststreamerinfos.cxx:10
ActsTrackingGeometrySvc::m_subDetNoAlign
std::set< ActsTrk::DetectorType > m_subDetNoAlign
Definition: ActsTrackingGeometrySvc.h:125
ATH_MSG_WARNING
#define ATH_MSG_WARNING(x)
Definition: AthMsgStreamMacros.h:32
python.CaloScaleNoiseConfig.type
type
Definition: CaloScaleNoiseConfig.py:78
DetectorZone::barrel
@ barrel
InDetDD::SiDetectorManager
Definition: SiDetectorManager.h:60
HGTD_DetectorManager
Definition: HGTD_DetectorManager.h:33
ActsTrackingGeometrySvc::makeStrawLayerBuilder
std::shared_ptr< const Acts::ILayerBuilder > makeStrawLayerBuilder(const InDetDD::InDetDetectorManager *manager)
Definition: ActsTrackingGeometrySvc.cxx:648
python.Logging.manager
manager
Definition: PhysicsAnalysis/D3PDTools/AnaAlgorithm/python/Logging.py:92
ActsTrackingGeometrySvc::m_useMaterialMap
Gaudi::Property< bool > m_useMaterialMap
Definition: ActsTrackingGeometrySvc.h:102
ActsTrackingGeometrySvc::m_trackingGeometry
std::shared_ptr< const Acts::TrackingGeometry > m_trackingGeometry
Definition: ActsTrackingGeometrySvc.h:95
ActsTrackingGeometrySvc::makeHGTDLayerBuilder
std::shared_ptr< const Acts::ILayerBuilder > makeHGTDLayerBuilder(const HGTD_DetectorManager *manager)
Definition: ActsTrackingGeometrySvc.cxx:679
ActsTrackingGeometrySvc::m_buildSubdetectors
Gaudi::Property< std::vector< std::string > > m_buildSubdetectors
Definition: ActsTrackingGeometrySvc.h:110
IdentityHelper
Definition: IdentityHelper.h:14
ActsTrackingGeometrySvc::m_nominalContext
ActsGeometryContext m_nominalContext
Definition: ActsTrackingGeometrySvc.h:100
ActsTrackingGeometrySvc::m_TRT_idHelper
const TRT_ID * m_TRT_idHelper
Definition: ActsTrackingGeometrySvc.h:97
ActsTrackingGeometrySvc::m_HGTD_idHelper
const HGTD_ID * m_HGTD_idHelper
Definition: ActsTrackingGeometrySvc.h:98
python.AutoConfigFlags.msg
msg
Definition: AutoConfigFlags.py:7
ActsTrackingGeometrySvc::m_subDetNoAlignProp
Gaudi::Property< std::vector< unsigned int > > m_subDetNoAlignProp
Define the subdetectors for which the tracking geometry does not expect a valid alignment store.
Definition: ActsTrackingGeometrySvc.h:124
ActsTrackingGeometrySvc::m_objDebugOutput
Gaudi::Property< bool > m_objDebugOutput
Definition: ActsTrackingGeometrySvc.h:103
SiliconTech::pixel
@ pixel
generate::Zero
void Zero(TH1D *hin)
Definition: generate.cxx:32
ActsLayerBuilder::Mode::SCT
@ SCT
Trk::binPhi
@ binPhi
Definition: BinningType.h:51
Generated on Tue Jul 2 2024 21:36:33 for ATLAS Offline Software by   doxygen 1.8.18