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"}}
 
Gaudi::Property< std::vector< float > > m_passiveITkInnerPixelBarrelLayerRadii {this, "PassiveITkInnerPixelBarrelLayerRadii", {}}
 the specifications for building additional passive cylinders in the barrel region: for each cylinder you want to specify radius, half length in z and thickness More...
 
Gaudi::Property< std::vector< float > > m_passiveITkInnerPixelBarrelLayerHalflengthZ {this, "PassiveITkInnerPixelBarrelLayerHalflengthZ", {}}
 
Gaudi::Property< std::vector< float > > m_passiveITkInnerPixelBarrelLayerThickness {this, "PassiveITkInnerPixelBarrelLayerThickness", {}}
 
Gaudi::Property< std::vector< float > > m_passiveITkOuterPixelBarrelLayerRadii {this, "PassiveITkOuterPixelBarrelLayerRadii", {}}
 
Gaudi::Property< std::vector< float > > m_passiveITkOuterPixelBarrelLayerHalflengthZ {this, "PassiveITkOuterPixelBarrelLayerHalflengthZ", {}}
 
Gaudi::Property< std::vector< float > > m_passiveITkOuterPixelBarrelLayerThickness {this, "PassiveITkOuterPixelBarrelLayerThickness", {}}
 
Gaudi::Property< std::vector< float > > m_passiveITkStripBarrelLayerRadii {this, "PassiveITkStripBarrelLayerRadii", {}}
 
Gaudi::Property< std::vector< float > > m_passiveITkStripBarrelLayerHalflengthZ {this, "PassiveITkStripBarrelLayerHalflengthZ", {}}
 
Gaudi::Property< std::vector< float > > m_passiveITkStripBarrelLayerThickness {this, "PassiveITkStripBarrelLayerThickness", {}}
 
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 61 of file ActsTrackingGeometrySvc.cxx.

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

Member Function Documentation

◆ getNominalContext()

const ActsGeometryContext & ActsTrackingGeometrySvc::getNominalContext ( ) const
override

Definition at line 1054 of file ActsTrackingGeometrySvc.cxx.

1054 { return m_nominalContext; }

◆ initialize()

StatusCode ActsTrackingGeometrySvc::initialize ( )
override

Definition at line 69 of file ActsTrackingGeometrySvc.cxx.

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

◆ makeBeamPipeConfig()

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

Definition at line 1057 of file ActsTrackingGeometrySvc.cxx.

1058  {
1059 
1060  // adapted from InnerDetector/InDetDetDescr/InDetTrackingGeometry/src/BeamPipeBuilder.cxx
1061 
1062  PVConstLink beamPipeTopVolume = p_beamPipeMgr->getTreeTop(0);
1063 
1064  if (p_beamPipeMgr->getNumTreeTops() == 1){
1065  beamPipeTopVolume = p_beamPipeMgr->getTreeTop(0)->getChildVol(0)->getChildVol(0);
1066  }
1067 
1068  Acts::Transform3 beamPipeTransform;
1069  beamPipeTransform.setIdentity();
1070 
1071  beamPipeTransform = Acts::Translation3(beamPipeTopVolume->getX().translation());
1072 
1073  double beamPipeRadius = 20;
1074 
1075  const GeoLogVol* beamPipeLogVolume = beamPipeTopVolume->getLogVol();
1076  const GeoTube* beamPipeTube = nullptr;
1077 
1078 
1079  if (beamPipeLogVolume == nullptr) {
1080  ATH_MSG_ERROR("Beam pip volume has no log volume");
1081  throw std::runtime_error("Beam pip volume has no log volume");
1082  }
1083  // get the geoShape and translate
1084  beamPipeTube = dynamic_cast<const GeoTube*>(beamPipeLogVolume->getShape());
1085  if (beamPipeTube == nullptr){
1086  ATH_MSG_ERROR("BeamPipeLogVolume was not of type GeoTube");
1087  throw std::runtime_error{"BeamPipeLogVolume was not of type GeoTube"};
1088  }
1089 
1090  for(unsigned int i=0;i<beamPipeTopVolume->getNChildVols();i++) {
1091 
1092  if(beamPipeTopVolume->getNameOfChildVol(i) == "SectionC03"){
1093 
1094  PVConstLink childTopVolume = beamPipeTopVolume->getChildVol(i);
1095  const GeoLogVol* childLogVolume = childTopVolume->getLogVol();
1096  const GeoTube* childTube = nullptr;
1097 
1098  if (childLogVolume){
1099  childTube = dynamic_cast<const GeoTube*>(childLogVolume->getShape());
1100  if (childTube){
1101  beamPipeRadius = 0.5 * (childTube->getRMax()+childTube->getRMin());
1102  }
1103  }
1104 
1105  break; // Exit loop after SectionC03 is found
1106  }
1107 
1108  } // Loop over child volumes
1109 
1110  ATH_MSG_VERBOSE("BeamPipe constructed from Database: translation (yes) - radius "
1111  << ( beamPipeTube ? "(yes)" : "(no)") << " - r = " << beamPipeRadius );
1112 
1113  ATH_MSG_VERBOSE("BeamPipe shift estimated as : " << Amg::toString(beamPipeTransform.translation()));
1114 
1115  Acts::CylinderVolumeBuilder::Config cfg;
1116 
1117  Acts::PassiveLayerBuilder::Config bplConfig;
1118  bplConfig.layerIdentification = "BeamPipe";
1119  bplConfig.centralLayerRadii = {beamPipeRadius * 1_mm};
1120  bplConfig.centralLayerHalflengthZ = {3000_mm};
1121  bplConfig.centralLayerThickness = {1_mm};
1122  auto beamPipeBuilder = std::make_shared<const Acts::PassiveLayerBuilder>(
1123  bplConfig, makeActsAthenaLogger(this, std::string("BPLayBldr"), std::string("ActsTGSvc")));
1124 
1125  // create the volume for the beam pipe
1126  cfg.trackingVolumeHelper = cvh;
1127  cfg.volumeName = "BeamPipe";
1128  cfg.layerBuilder = beamPipeBuilder;
1129  cfg.layerEnvelopeR = {1_mm, 1_mm};
1130  cfg.buildToRadiusZero = true;
1131 
1132  return cfg;
1133 }

◆ makeHGTDLayerBuilder()

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

Definition at line 710 of file ActsTrackingGeometrySvc.cxx.

711  {
712 
713  std::string managerName = manager->getName();
714  auto matcher = [](const Acts::GeometryContext & /*gctx*/,
715  Acts::AxisDirection /*aDir*/, const Acts::Surface * /*aS*/,
716  const Acts::Surface *
717  /*bS*/) -> bool { return false; };
718 
719  Acts::SurfaceArrayCreator::Config sacCfg;
720  sacCfg.surfaceMatcher = matcher;
721  sacCfg.doPhiBinningOptimization = false;
722 
723  auto surfaceArrayCreator = std::make_shared<Acts::SurfaceArrayCreator>(
724  sacCfg,
725  makeActsAthenaLogger(this, managerName + "SrfArrCrtr", std::string("ActsTGSvc")));
726  Acts::LayerCreator::Config lcCfg;
727  lcCfg.surfaceArrayCreator = surfaceArrayCreator;
728  auto layerCreator = std::make_shared<Acts::LayerCreator>(
729  lcCfg, makeActsAthenaLogger(this, managerName + "LayCrtr", std::string("ActsTGSvc")));
730 
731  ActsHGTDLayerBuilder::Config cfg;
732  cfg.mng = static_cast<const HGTD_DetectorManager *>(manager);
733  cfg.elementStore = m_elementStore;
734  cfg.layerCreator = layerCreator;
735  cfg.idHelper = m_HGTD_idHelper;
736  return std::make_shared<const ActsHGTDLayerBuilder>(
737  cfg, makeActsAthenaLogger(this, managerName + "GMSLayBldr", std::string("ActsTGSvc")));
738 }

◆ makeLayerBuilderConfig()

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

Definition at line 740 of file ActsTrackingGeometrySvc.cxx.

741  {
742  using enum Acts::AxisDirection;
743 
744  std::string managerName = manager->getName();
745 
746  std::shared_ptr<const Acts::ILayerBuilder> gmLayerBuilder;
747  auto matcher = [](const Acts::GeometryContext & /*gctx*/,
748  Acts::AxisDirection aDir, const Acts::Surface *aS,
749  const Acts::Surface *bS) -> bool {
750  auto a = dynamic_cast<const ActsDetectorElement *>(
751  aS->associatedDetectorElement());
752  auto b = dynamic_cast<const ActsDetectorElement *>(
753  bS->associatedDetectorElement());
754  if ((not a) or (not b)) {
755  throw std::runtime_error(
756  "Cast of surface associated element to ActsDetectorElement failed "
757  "in ActsTrackingGeometrySvc::makeVolumeBuilder");
758  }
759 
760  IdentityHelper idA = a->identityHelper();
761  IdentityHelper idB = b->identityHelper();
762 
763  // check if same bec
764  // can't be same if not
765  if (idA.bec() != idB.bec())
766  return false;
767 
768  if (aDir == AxisPhi) {
769  // std::cout << idA.phi_module() << " <-> " << idB.phi_module() <<
770  // std::endl;
771  return idA.phi_module() == idB.phi_module();
772  }
773 
774  if (aDir == AxisZ) {
775  return (idA.eta_module() == idB.eta_module()) &&
776  (idA.layer_disk() == idB.layer_disk()) && (idA.bec() == idB.bec());
777  }
778 
779  if (aDir == AxisR) {
780  return (idA.eta_module() == idB.eta_module()) &&
781  (idA.layer_disk() == idB.layer_disk()) && (idB.bec() == idA.bec());
782  }
783 
784  return false;
785  };
786 
787  Acts::SurfaceArrayCreator::Config sacCfg;
788  sacCfg.surfaceMatcher = matcher;
789 
790  auto surfaceArrayCreator = std::make_shared<Acts::SurfaceArrayCreator>(
791  sacCfg,
792  makeActsAthenaLogger(this, managerName + "SrfArrCrtr", std::string("ActsTGSvc")));
793  Acts::LayerCreator::Config lcCfg;
794  lcCfg.surfaceArrayCreator = surfaceArrayCreator;
795  auto layerCreator = std::make_shared<Acts::LayerCreator>(
796  lcCfg, makeActsAthenaLogger(this, managerName + "LayCrtr", std::string("ActsTGSvc")));
797 
798  ActsLayerBuilder::Config cfg;
799  cfg.surfaceMatcher = matcher;
800 
801  // set bins from configuration
802  if (m_barrelMaterialBins.size() != 2) {
803  throw std::invalid_argument("Number of barrel material bin counts != 2");
804  }
805  std::vector<size_t> brlBins(m_barrelMaterialBins);
806  cfg.barrelMaterialBins = {brlBins.at(0), brlBins.at(1)};
807 
808  if (m_endcapMaterialBins.size() != 2) {
809  throw std::invalid_argument("Number of endcap material bin counts != 2");
810  }
811  std::vector<size_t> ecBins(m_endcapMaterialBins);
812  cfg.endcapMaterialBins = {ecBins.at(0), ecBins.at(1)};
813 
814  cfg.mng = static_cast<const InDetDD::SiDetectorManager *>(manager);
815  // use class member element store
816  cfg.elementStore = m_elementStore;
817  cfg.layerCreator = layerCreator;
818 
819  // gmLayerBuilder = std::make_shared<const ActsLayerBuilder>(
820  // cfg, makeActsAthenaLogger(this, managerName + "GMLayBldr",
821  // "ActsTGSvc"));
822 
823  // return gmLayerBuilder;
824  return cfg;
825 }

◆ 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 828 of file ActsTrackingGeometrySvc.cxx.

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

◆ makeStrawLayerBuilder()

std::shared_ptr< const Acts::ILayerBuilder > ActsTrackingGeometrySvc::makeStrawLayerBuilder ( const InDetDD::InDetDetectorManager *  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::AxisDirection /*aDir*/, 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  ActsStrawLayerBuilder::Config cfg;
701  cfg.mng = static_cast<const InDetDD::TRT_DetectorManager *>(manager);
702  cfg.elementStore = m_elementStore;
703  cfg.layerCreator = layerCreator;
704  cfg.idHelper = m_TRT_idHelper;
705  return std::make_shared<const ActsStrawLayerBuilder>(
706  cfg, makeActsAthenaLogger(this, managerName + "GMSLayBldr", std::string("ActsTGSvc")));
707 }

◆ populateAlignmentStore()

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

Definition at line 1043 of file ActsTrackingGeometrySvc.cxx.

1043  {
1044  ATH_MSG_DEBUG("Populate the alignment store with all detector elements");
1045  unsigned int nElements = 0;
1046  m_trackingGeometry->visitSurfaces([&store, &nElements](const Acts::Surface *srf) {
1047  const Acts::DetectorElementBase *detElem = srf->associatedDetectorElement();
1048  const IDetectorElement *gmde = dynamic_cast<const IDetectorElement *>(detElem);
1049  nElements += gmde->storeAlignedTransforms(store);
1050  });
1051  ATH_MSG_DEBUG("Populated with " << nElements << " elements");
1052  return nElements;
1053 }

◆ runConsistencyChecks()

bool ActsTrackingGeometrySvc::runConsistencyChecks ( ) const
private

Definition at line 445 of file ActsTrackingGeometrySvc.cxx.

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

◆ trackingGeometry()

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

Definition at line 672 of file ActsTrackingGeometrySvc.cxx.

672  {
673 
674  ATH_MSG_VERBOSE("Retrieving tracking geometry");
675  return m_trackingGeometry;
676 }

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 137 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 131 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 134 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_passiveITkInnerPixelBarrelLayerHalflengthZ

Gaudi::Property<std::vector<float> > ActsTrackingGeometrySvc::m_passiveITkInnerPixelBarrelLayerHalflengthZ {this, "PassiveITkInnerPixelBarrelLayerHalflengthZ", {}}
private

Definition at line 117 of file ActsTrackingGeometrySvc.h.

◆ m_passiveITkInnerPixelBarrelLayerRadii

Gaudi::Property<std::vector<float> > ActsTrackingGeometrySvc::m_passiveITkInnerPixelBarrelLayerRadii {this, "PassiveITkInnerPixelBarrelLayerRadii", {}}
private

the specifications for building additional passive cylinders in the barrel region: for each cylinder you want to specify radius, half length in z and thickness

Definition at line 116 of file ActsTrackingGeometrySvc.h.

◆ m_passiveITkInnerPixelBarrelLayerThickness

Gaudi::Property<std::vector<float> > ActsTrackingGeometrySvc::m_passiveITkInnerPixelBarrelLayerThickness {this, "PassiveITkInnerPixelBarrelLayerThickness", {}}
private

Definition at line 118 of file ActsTrackingGeometrySvc.h.

◆ m_passiveITkOuterPixelBarrelLayerHalflengthZ

Gaudi::Property<std::vector<float> > ActsTrackingGeometrySvc::m_passiveITkOuterPixelBarrelLayerHalflengthZ {this, "PassiveITkOuterPixelBarrelLayerHalflengthZ", {}}
private

Definition at line 121 of file ActsTrackingGeometrySvc.h.

◆ m_passiveITkOuterPixelBarrelLayerRadii

Gaudi::Property<std::vector<float> > ActsTrackingGeometrySvc::m_passiveITkOuterPixelBarrelLayerRadii {this, "PassiveITkOuterPixelBarrelLayerRadii", {}}
private

Definition at line 120 of file ActsTrackingGeometrySvc.h.

◆ m_passiveITkOuterPixelBarrelLayerThickness

Gaudi::Property<std::vector<float> > ActsTrackingGeometrySvc::m_passiveITkOuterPixelBarrelLayerThickness {this, "PassiveITkOuterPixelBarrelLayerThickness", {}}
private

Definition at line 122 of file ActsTrackingGeometrySvc.h.

◆ m_passiveITkStripBarrelLayerHalflengthZ

Gaudi::Property<std::vector<float> > ActsTrackingGeometrySvc::m_passiveITkStripBarrelLayerHalflengthZ {this, "PassiveITkStripBarrelLayerHalflengthZ", {}}
private

Definition at line 125 of file ActsTrackingGeometrySvc.h.

◆ m_passiveITkStripBarrelLayerRadii

Gaudi::Property<std::vector<float> > ActsTrackingGeometrySvc::m_passiveITkStripBarrelLayerRadii {this, "PassiveITkStripBarrelLayerRadii", {}}
private

Definition at line 124 of file ActsTrackingGeometrySvc.h.

◆ m_passiveITkStripBarrelLayerThickness

Gaudi::Property<std::vector<float> > ActsTrackingGeometrySvc::m_passiveITkStripBarrelLayerThickness {this, "PassiveITkStripBarrelLayerThickness", {}}
private

Definition at line 126 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 128 of file ActsTrackingGeometrySvc.h.

◆ m_subDetNoAlign

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

Definition at line 141 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 140 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:1057
IdentityHelper::phi_module
int phi_module() const
Definition: IdentityHelper.cxx:61
AllowedVariables::e
e
Definition: AsgElectronSelectorTool.cxx:37
SGTest::store
TestStore store
Definition: TestStore.cxx:23
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:128
get_generator_info.result
result
Definition: get_generator_info.py:21
ActsLayerBuilder::Mode::ITkPixelOuter
@ ITkPixelOuter
ATH_MSG_INFO
#define ATH_MSG_INFO(x)
Definition: AthMsgStreamMacros.h:31
ActsTrackingGeometrySvc::m_passiveITkStripBarrelLayerHalflengthZ
Gaudi::Property< std::vector< float > > m_passiveITkStripBarrelLayerHalflengthZ
Definition: ActsTrackingGeometrySvc.h:125
ActsTrackingGeometrySvc::m_consistencyCheckOutput
StringProperty m_consistencyCheckOutput
Definition: ActsTrackingGeometrySvc.h:131
ActsLayerBuilder::Mode::ITkPixelInner
@ ITkPixelInner
ActsTrk::DetectorType
DetectorType
Simple enum to Identify the Type of the ACTS sub detector.
Definition: GeometryDefs.h:17
Amg::Vector2D
Eigen::Matrix< double, 2, 1 > Vector2D
Definition: GeoPrimitives.h:48
max
constexpr double max()
Definition: ap_fixedTest.cxx:33
min
constexpr double min()
Definition: ap_fixedTest.cxx:26
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:445
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:134
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:51
python.CaloAddPedShiftConfig.type
type
Definition: CaloAddPedShiftConfig.py:42
ActsTrackingGeometrySvc::p_beamPipeMgr
const BeamPipeDetectorManager * p_beamPipeMgr
Definition: ActsTrackingGeometrySvc.h:91
ActsTrackingGeometrySvc::m_caloVolumeBuilder
ToolHandle< IActsTrackingVolumeBuilder > m_caloVolumeBuilder
Definition: ActsTrackingGeometrySvc.h:137
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:1054
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:45
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
python.BunchSpacingUtils.lb
lb
Definition: BunchSpacingUtils.py:88
lumiFormat.i
int i
Definition: lumiFormat.py:85
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::m_passiveITkStripBarrelLayerRadii
Gaudi::Property< std::vector< float > > m_passiveITkStripBarrelLayerRadii
Definition: ActsTrackingGeometrySvc.h:124
ActsTrackingGeometrySvc::m_passiveITkOuterPixelBarrelLayerRadii
Gaudi::Property< std::vector< float > > m_passiveITkOuterPixelBarrelLayerRadii
Definition: ActsTrackingGeometrySvc.h:120
ActsTrackingGeometrySvc::makeLayerBuilderConfig
ActsLayerBuilder::Config makeLayerBuilderConfig(const InDetDD::InDetDetectorManager *manager)
Definition: ActsTrackingGeometrySvc.cxx:740
ATH_CHECK
#define ATH_CHECK
Definition: AthCheckMacros.h:40
ActsDetectorElement
Definition: ActsDetectorElement.h:42
Handler::svc
AthROOTErrorHandlerSvc * svc
Definition: AthROOTErrorHandlerSvc.cxx:10
ActsTrk::IDetectorElement
Base class interface for the actual readout elements.
Definition: IDetectorElement.h:43
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
ActsTrackingGeometrySvc::m_passiveITkInnerPixelBarrelLayerThickness
Gaudi::Property< std::vector< float > > m_passiveITkInnerPixelBarrelLayerThickness
Definition: ActsTrackingGeometrySvc.h:118
IdentityHelper::eta_module
int eta_module() const
Definition: IdentityHelper.cxx:73
name
std::string name
Definition: Control/AthContainers/Root/debug.cxx:228
plotBeamSpotMon.b
b
Definition: plotBeamSpotMon.py:77
ActsTrackingGeometrySvc::m_passiveITkStripBarrelLayerThickness
Gaudi::Property< std::vector< float > > m_passiveITkStripBarrelLayerThickness
Definition: ActsTrackingGeometrySvc.h:126
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:828
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
ActsTrackingGeometrySvc::m_passiveITkOuterPixelBarrelLayerHalflengthZ
Gaudi::Property< std::vector< float > > m_passiveITkOuterPixelBarrelLayerHalflengthZ
Definition: ActsTrackingGeometrySvc.h:121
IdentityHelper::bec
int bec() const
Definition: IdentityHelper.cxx:37
python.SystemOfUnits.mm
int mm
Definition: SystemOfUnits.py:83
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
a
TList * a
Definition: liststreamerinfos.cxx:10
ActsTrackingGeometrySvc::m_subDetNoAlign
std::set< ActsTrk::DetectorType > m_subDetNoAlign
Definition: ActsTrackingGeometrySvc.h:141
ActsTrackingGeometrySvc::m_passiveITkOuterPixelBarrelLayerThickness
Gaudi::Property< std::vector< float > > m_passiveITkOuterPixelBarrelLayerThickness
Definition: ActsTrackingGeometrySvc.h:122
ATH_MSG_WARNING
#define ATH_MSG_WARNING(x)
Definition: AthMsgStreamMacros.h:32
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:679
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
DataPrepToActsConfig.pixelContainer
pixelContainer
Definition: DataPrepToActsConfig.py:9
ActsTrackingGeometrySvc::makeHGTDLayerBuilder
std::shared_ptr< const Acts::ILayerBuilder > makeHGTDLayerBuilder(const HGTD_DetectorManager *manager)
Definition: ActsTrackingGeometrySvc.cxx:710
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
Cut::all
@ all
Definition: SUSYToolsAlg.cxx:67
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
ActsTrackingGeometrySvc::m_passiveITkInnerPixelBarrelLayerRadii
Gaudi::Property< std::vector< float > > m_passiveITkInnerPixelBarrelLayerRadii
the specifications for building additional passive cylinders in the barrel region: for each cylinder ...
Definition: ActsTrackingGeometrySvc.h:116
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:140
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
ActsTrackingGeometrySvc::m_passiveITkInnerPixelBarrelLayerHalflengthZ
Gaudi::Property< std::vector< float > > m_passiveITkInnerPixelBarrelLayerHalflengthZ
Definition: ActsTrackingGeometrySvc.h:117
Generated on Tue Mar 25 2025 21:26:00 for ATLAS Offline Software by   doxygen 1.8.18