ATLAS Offline Software
SeedingTool.cxx
Go to the documentation of this file.
1 /* Copyright (C) 2002-2023 CERN for the benefit of the ATLAS collaboration
2 */
3 
4 #include "src/SeedingTool.h"
5 
6 // ACTS
7 #include "Acts/Seeding/SeedFilterConfig.hpp"
8 #include "Acts/Seeding/BinnedGroup.hpp"
9 #include "Acts/Seeding/SeedFilter.hpp"
10 #include "Acts/Seeding/SeedFinder.hpp"
11 #include "Acts/Seeding/SeedFinderConfig.hpp"
12 #include "Acts/Definitions/Units.hpp"
13 #include "Acts/Seeding/SeedConfirmationRangeConfig.hpp"
14 
15 using namespace Acts::HashedStringLiteral;
16 
17 namespace ActsTrk {
18  SeedingTool::SeedingTool(const std::string& type,
19  const std::string& name,
20  const IInterface* parent)
21  : base_class(type, name, parent)
22  {}
23 
24  StatusCode SeedingTool::initialize() {
25  ATH_MSG_DEBUG("Initializing " << name() << "...");
26 
27  ATH_MSG_DEBUG("Properties Summary:");
28  ATH_MSG_DEBUG(" " << m_zBinNeighborsTop);
29  ATH_MSG_DEBUG(" " << m_zBinNeighborsBottom);
30  ATH_MSG_DEBUG(" " << m_rBinNeighborsTop);
31  ATH_MSG_DEBUG(" " << m_rBinNeighborsBottom);
32  ATH_MSG_DEBUG(" " << m_numPhiNeighbors);
33 
34  ATH_MSG_DEBUG(" * Used by SpacePointGridConfig");
35  ATH_MSG_DEBUG(" " << m_minPt);
36  ATH_MSG_DEBUG(" " << m_cotThetaMax);
37  ATH_MSG_DEBUG(" " << m_impactMax);
38  ATH_MSG_DEBUG(" " << m_zMin);
39  ATH_MSG_DEBUG(" " << m_zMax);
40  ATH_MSG_DEBUG(" " << m_gridPhiMin);
41  ATH_MSG_DEBUG(" " << m_gridPhiMax);
42  ATH_MSG_DEBUG(" " << m_zBinEdges);
43  ATH_MSG_DEBUG(" " << m_rBinEdges);
44  ATH_MSG_DEBUG(" " << m_deltaRMax);
45  ATH_MSG_DEBUG(" " << m_gridRMax);
46  ATH_MSG_DEBUG(" " << m_phiBinDeflectionCoverage);
47 
48  ATH_MSG_DEBUG(" * Used by SeedFinderConfig:");
49  ATH_MSG_DEBUG(" " << m_minPt);
50  ATH_MSG_DEBUG(" " << m_cotThetaMax);
51  ATH_MSG_DEBUG(" " << m_impactMax);
52  ATH_MSG_DEBUG(" " << m_zMin);
53  ATH_MSG_DEBUG(" " << m_zMax);
54  ATH_MSG_DEBUG(" " << m_zBinEdges);
55  ATH_MSG_DEBUG(" " << m_rMax);
56  ATH_MSG_DEBUG(" " << m_deltaRMin);
57  ATH_MSG_DEBUG(" " << m_deltaRMax);
58  ATH_MSG_DEBUG(" " << m_deltaRMinTopSP);
59  ATH_MSG_DEBUG(" " << m_deltaRMaxTopSP);
60  ATH_MSG_DEBUG(" " << m_deltaRMinBottomSP);
61  ATH_MSG_DEBUG(" " << m_deltaRMaxBottomSP);
62  ATH_MSG_DEBUG(" " << m_deltaZMax);
63  ATH_MSG_DEBUG(" " << m_collisionRegionMin);
64  ATH_MSG_DEBUG(" " << m_collisionRegionMax);
65  ATH_MSG_DEBUG(" " << m_sigmaScattering);
66  ATH_MSG_DEBUG(" " << m_maxPtScattering);
67  ATH_MSG_DEBUG(" " << m_radLengthPerSeed);
68  ATH_MSG_DEBUG(" " << m_maxSeedsPerSpM);
69  ATH_MSG_DEBUG(" " << m_interactionPointCut);
70  ATH_MSG_DEBUG(" " << m_zBinsCustomLooping);
71  ATH_MSG_DEBUG(" " << m_rBinsCustomLooping);
72  ATH_MSG_DEBUG(" " << m_useVariableMiddleSPRange);
73  if ( m_useVariableMiddleSPRange ) {
74  ATH_MSG_DEBUG(" " << m_deltaRMiddleMinSPRange);
75  ATH_MSG_DEBUG(" " << m_deltaRMiddleMaxSPRange);
76  } else if ( not m_rRangeMiddleSP.empty() )
77  ATH_MSG_DEBUG(" " << m_rRangeMiddleSP);
78  ATH_MSG_DEBUG(" " << m_seedConfirmation);
79  if ( m_seedConfirmation ) {
80  ATH_MSG_DEBUG(" " << m_seedConfCentralZMin);
81  ATH_MSG_DEBUG(" " << m_seedConfCentralZMax);
82  ATH_MSG_DEBUG(" " << m_seedConfCentralRMax);
83  ATH_MSG_DEBUG(" " << m_seedConfCentralNTopLargeR);
84  ATH_MSG_DEBUG(" " << m_seedConfCentralNTopSmallR);
85  ATH_MSG_DEBUG(" " << m_seedConfCentralMinBottomRadius);
86  ATH_MSG_DEBUG(" " << m_seedConfCentralMaxZOrigin);
87  ATH_MSG_DEBUG(" " << m_seedConfCentralMinImpact);
88  ATH_MSG_DEBUG(" " << m_seedConfForwardZMin);
89  ATH_MSG_DEBUG(" " << m_seedConfForwardZMax);
90  ATH_MSG_DEBUG(" " << m_seedConfForwardRMax);
91  ATH_MSG_DEBUG(" " << m_seedConfForwardNTopLargeR);
92  ATH_MSG_DEBUG(" " << m_seedConfForwardNTopSmallR);
93  ATH_MSG_DEBUG(" " << m_seedConfForwardMinBottomRadius);
94  ATH_MSG_DEBUG(" " << m_seedConfForwardMaxZOrigin);
95  ATH_MSG_DEBUG(" " << m_seedConfForwardMinImpact);
96  }
97  ATH_MSG_DEBUG(" " << m_useDetailedDoubleMeasurementInfo);
98  ATH_MSG_DEBUG(" " << m_toleranceParam);
99  ATH_MSG_DEBUG(" " << m_phiMin);
100  ATH_MSG_DEBUG(" " << m_phiMax);
101  ATH_MSG_DEBUG(" " << m_rMin);
102  ATH_MSG_DEBUG(" " << m_zAlign);
103  ATH_MSG_DEBUG(" " << m_rAlign);
104  ATH_MSG_DEBUG(" " << m_sigmaError);
105 
106  ATH_MSG_DEBUG(" * Used by SeedFilterConfig:");
107  ATH_MSG_DEBUG(" " << m_deltaRMin);
108  ATH_MSG_DEBUG(" " << m_maxSeedsPerSpM);
109  ATH_MSG_DEBUG(" " << m_useDeltaRorTopRadius);
110  ATH_MSG_DEBUG(" " << m_seedConfirmationInFilter);
111  if (m_seedConfirmationInFilter) {
112  ATH_MSG_DEBUG(" " << m_maxSeedsPerSpMConf);
113  ATH_MSG_DEBUG(" " << m_maxQualitySeedsPerSpMConf);
114  ATH_MSG_DEBUG(" " << m_seedConfCentralZMin);
115  ATH_MSG_DEBUG(" " << m_seedConfCentralZMax);
116  ATH_MSG_DEBUG(" " << m_seedConfCentralRMax);
117  ATH_MSG_DEBUG(" " << m_seedConfCentralNTopLargeR);
118  ATH_MSG_DEBUG(" " << m_seedConfCentralNTopSmallR);
119  ATH_MSG_DEBUG(" " << m_seedConfCentralMinBottomRadius);
120  ATH_MSG_DEBUG(" " << m_seedConfCentralMaxZOrigin);
121  ATH_MSG_DEBUG(" " << m_seedConfCentralMinImpact);
122  ATH_MSG_DEBUG(" " << m_seedConfForwardZMin);
123  ATH_MSG_DEBUG(" " << m_seedConfForwardZMax);
124  ATH_MSG_DEBUG(" " << m_seedConfForwardRMax);
125  ATH_MSG_DEBUG(" " << m_seedConfForwardNTopLargeR);
126  ATH_MSG_DEBUG(" " << m_seedConfForwardNTopSmallR);
127  ATH_MSG_DEBUG(" " << m_seedConfForwardMinBottomRadius);
128  ATH_MSG_DEBUG(" " << m_seedConfForwardMaxZOrigin);
129  ATH_MSG_DEBUG(" " << m_seedConfForwardMinImpact);
130  }
131  ATH_MSG_DEBUG(" " << m_impactWeightFactor);
132  ATH_MSG_DEBUG(" " << m_compatSeedWeight);
133  ATH_MSG_DEBUG(" " << m_compatSeedLimit);
134  ATH_MSG_DEBUG(" " << m_seedWeightIncrement);
135  ATH_MSG_DEBUG(" " << m_numSeedIncrement);
136  ATH_MSG_DEBUG(" " << m_deltaInvHelixDiameter);
137 
138  // Make the logger And Propagate to ACTS routines
139  m_logger = makeActsAthenaLogger(this, "Acts");
140 
141  if (m_zBinEdges.size() - 1 !=
142  m_zBinNeighborsTop.size() and
143  not m_zBinNeighborsTop.empty()) {
144  ATH_MSG_ERROR("Inconsistent config zBinNeighborsTop");
145  return StatusCode::FAILURE;
146  }
147 
148  if (m_zBinEdges.size() - 1 !=
149  m_zBinNeighborsBottom.size() and
150  not m_zBinNeighborsBottom.empty()) {
151  ATH_MSG_ERROR("Inconsistent config zBinNeighborsBottom");
152  return StatusCode::FAILURE;
153  }
154 
155  if (m_rBinEdges.size() - 1 !=
156  m_rBinNeighborsTop.size() and
157  not m_rBinNeighborsTop.empty()) {
158  ATH_MSG_ERROR("Inconsistent config rBinNeighborsTop");
159  return StatusCode::FAILURE;
160  }
161 
162  if (m_rBinEdges.size() - 1 !=
163  m_rBinNeighborsBottom.size() and
164  not m_rBinNeighborsBottom.empty()) {
165  ATH_MSG_ERROR("Inconsistent config rBinNeighborsBottom");
166  return StatusCode::FAILURE;
167  }
168 
169  if (m_zBinsCustomLooping.size() != 0) {
170  // zBinsCustomLooping can contain a number of elements <= to the total number
171  // of bin in zBinEdges
172  for (std::size_t i : m_zBinsCustomLooping) {
173  if (i >= m_zBinEdges.size()) {
174  ATH_MSG_ERROR("Inconsistent config zBinsCustomLooping contains bins that are not in zBinEdges");
175  return StatusCode::FAILURE;
176  }
177  }
178  }
179 
180  if (m_rBinsCustomLooping.size() != 0) {
181  for (std::size_t i : m_rBinsCustomLooping) {
182  if (i >= m_rBinEdges.size()) {
183  ATH_MSG_ERROR("Inconsistent config rBinsCustomLooping contains bins that are not in rBinEdges");
184  return StatusCode::FAILURE;
185  }
186  }
187  }
188 
189  ATH_CHECK( prepareConfiguration() );
190 
191  m_bottomBinFinder = std::make_unique< Acts::GridBinFinder< 3ul > >(m_numPhiNeighbors.value(),
192  m_zBinNeighborsBottom.value(),
193  m_rBinNeighborsBottom.value());
194  m_topBinFinder = std::make_unique< Acts::GridBinFinder< 3ul > >(m_numPhiNeighbors.value(),
195  m_zBinNeighborsTop.value(),
196  m_rBinNeighborsTop.value());
197 
198  m_navigation[0ul] = {};
199  m_navigation[1ul] = m_finderCfg.zBinsCustomLooping;
200  m_navigation[2ul] = m_rBinsCustomLooping.value();
201 
202  return StatusCode::SUCCESS;
203  }
204 
205  StatusCode
206  SeedingTool::createSeeds(const EventContext& /*ctx*/,
207  const Acts::SpacePointContainer<ActsTrk::SpacePointCollector, Acts::detail::RefHolder>& spContainer,
208  const Acts::Vector3& beamSpotPos,
209  const Acts::Vector3& bField,
210  ActsTrk::SeedContainer& seedContainer ) const
211  {
212  // Create Seeds
213  //TODO POSSIBLE OPTIMISATION come back here: see MR !52399 ( i.e. use static thread_local)
214  ATH_CHECK(createSeeds(spContainer.begin(),
215  spContainer.end(),
216  beamSpotPos,
217  bField,
218  seedContainer));
219 
220  return StatusCode::SUCCESS;
221  }
222 
223  template< typename external_iterator_t >
224  StatusCode
225  SeedingTool::createSeeds(external_iterator_t spBegin,
226  external_iterator_t spEnd,
227  const Acts::Vector3& beamSpotPos,
228  const Acts::Vector3& bField,
229  DataVector< Acts::Seed< typename SeedingTool::external_type, 3ul > >& seedContainer) const
230  {
231  static_assert(std::is_same<typename external_spacepoint< external_iterator_t >::type, const value_type&>::value,
232  "Inconsistent type");
233 
234  if (spBegin == spEnd)
235  return StatusCode::SUCCESS;
236 
237  std::vector< seed_type > seeds;
238 
239  // Space Point Grid Options
240  Acts::CylindricalSpacePointGridOptions gridOpts;
241  gridOpts.bFieldInZ = bField[2];
242  gridOpts = gridOpts.toInternalUnits();
243 
244  // Seed Finder Options
245  Acts::SeedFinderOptions finderOpts;
246  finderOpts.beamPos = Acts::Vector2(beamSpotPos[Amg::x],
247  beamSpotPos[Amg::y]);
248  finderOpts.bFieldInZ = bField[2];
249  finderOpts = finderOpts.toInternalUnits().calculateDerivedQuantities(m_finderCfg);
250 
251 
252 
253 
254  Acts::CylindricalSpacePointGrid< value_type > grid =
255  Acts::CylindricalSpacePointGridCreator::createGrid< value_type >(m_gridCfg, gridOpts, logger());
256 
257  Acts::CylindricalSpacePointGridCreator::fillGrid(m_finderCfg, finderOpts, grid,
258  spBegin, spEnd, logger());
259 
260  // Compute radius Range
261  // we rely on the fact the grid is storing the proxies
262  // with a sorting in the radius
263  float minRange = std::numeric_limits<float>::max();
264  float maxRange = std::numeric_limits<float>::lowest();
265  for (const auto& coll : grid) {
266  if (coll.empty()) {
267  continue;
268  }
269  const auto* firstEl = coll.front();
270  const auto* lastEl = coll.back();
271  minRange = std::min(firstEl->radius(), minRange);
272  maxRange = std::max(lastEl->radius(), maxRange);
273  }
274 
275  Acts::CylindricalBinnedGroup< value_type > spacePointsGrouping(std::move(grid), *m_bottomBinFinder,
276  *m_topBinFinder, m_navigation);
277 
278  // variable middle SP radial region of interest
279  const Acts::Range1D<float> rMiddleSPRange(std::floor(minRange/2)*2 + m_finderCfg.deltaRMiddleMinSPRange,
280  std::floor(maxRange/2)*2 - m_finderCfg.deltaRMiddleMaxSPRange);
281 
282  //TODO POSSIBLE OPTIMISATION come back here: see MR !52399 ( i.e. use static thread_local)
283  typename decltype(m_finder)::SeedingState state;
284  state.spacePointMutableData.resize(std::distance(spBegin, spEnd));
285 
286  for (const auto [bottom, middle, top] : spacePointsGrouping) {
287  m_finder.createSeedsForGroup(finderOpts, state, spacePointsGrouping.grid(),
288  seeds, bottom, middle, top, rMiddleSPRange);
289  }
290 
291  if (m_seedQualitySelection) {
292  // Selection function - temporary implementation
293  // need change from ACTS for final implementation
294  // To be used only on PPP
295  auto selectionFunction = [&state] (const seed_type& seed) -> bool
296  {
297  float seed_quality = seed.seedQuality();
298  float bottom_quality = state.spacePointMutableData.quality(seed.sp()[0]->index());
299  float middle_quality = state.spacePointMutableData.quality(seed.sp()[1]->index());
300  float top_quality = state.spacePointMutableData.quality(seed.sp()[2]->index());
301 
302  if (bottom_quality > seed_quality and
303  middle_quality > seed_quality and
304  top_quality > seed_quality) {
305  return false;
306  }
307 
308  return true;
309  };
310 
311  // Select the seeds
312  std::size_t acceptedSeeds = 0;
313  for (std::size_t i(0); i<seeds.size(); ++i) {
314  const auto& seed = seeds[i];
315  if (not selectionFunction(seed)) {
316  continue;
317  }
318  // move passing seeds at the beginning of the vector/collection
319  // no need to swap them both. The seed currently at acceptedSeeds
320  // didn't make it (if i != acceptedSeeds)
321  if (acceptedSeeds != i)
322  seeds[acceptedSeeds] = std::move(seeds[i]);
323  ++acceptedSeeds;
324  }
325  // remove seeds that didn't make it
326  // they are all at the end of the collection
327  seeds.erase(seeds.begin() + acceptedSeeds, seeds.end());
328  }
329 
330 
331  // Store seeds
332  seedContainer.reserve(seeds.size());
333  for(const auto& seed: seeds) {
334  const auto [bottom, middle, top] = seed.sp();
335 
336  std::unique_ptr< ActsTrk::Seed > toAdd =
337  std::make_unique< ActsTrk::Seed >(bottom->externalSpacePoint(),
338  middle->externalSpacePoint(),
339  top->externalSpacePoint());
340  toAdd->setVertexZ(seed.z());
341  toAdd->setQuality(seed.seedQuality());
342  seedContainer.push_back(std::move(toAdd));
343  }
344 
345  return StatusCode::SUCCESS;
346  }
347 
348  StatusCode
349  SeedingTool::prepareConfiguration() {
350  assert(m_logger != nullptr);
351 
352  // Prepare the Acts::SeedFinderConfig object
353  // This is done only once, during initialization using the
354  // parameters set in the JO
355 
356  // Configuration for Acts::SeedFinder
357  m_finderCfg.minPt = m_minPt;
358  m_finderCfg.cotThetaMax = m_cotThetaMax;
359  m_finderCfg.impactMax = m_impactMax;
360  m_finderCfg.zMin = m_zMin;
361  m_finderCfg.zMax = m_zMax;
362  m_finderCfg.zBinEdges = m_zBinEdges;
363  m_finderCfg.rMax = m_rMax;
364  m_finderCfg.binSizeR = m_binSizeR;
365  m_finderCfg.deltaRMin = m_deltaRMin;
366  m_finderCfg.deltaRMax = m_deltaRMax;
367  m_finderCfg.deltaRMinTopSP = m_deltaRMinTopSP;
368  m_finderCfg.deltaRMaxTopSP = m_deltaRMaxTopSP;
369  m_finderCfg.deltaRMinBottomSP = m_deltaRMinBottomSP;
370  m_finderCfg.deltaRMaxBottomSP = m_deltaRMaxBottomSP;
371  m_finderCfg.deltaZMax = m_deltaZMax;
372  m_finderCfg.collisionRegionMin = m_collisionRegionMin;
373  m_finderCfg.collisionRegionMax = m_collisionRegionMax;
374  m_finderCfg.sigmaScattering = m_sigmaScattering;
375  m_finderCfg.maxPtScattering = m_maxPtScattering;
376  m_finderCfg.radLengthPerSeed = m_radLengthPerSeed;
377  m_finderCfg.maxSeedsPerSpM = m_maxSeedsPerSpM;
378  m_finderCfg.interactionPointCut = m_interactionPointCut;
379  m_finderCfg.zBinsCustomLooping = m_zBinsCustomLooping;
380  m_finderCfg.useVariableMiddleSPRange = m_useVariableMiddleSPRange;
381  m_finderCfg.deltaRMiddleMinSPRange = m_deltaRMiddleMinSPRange;
382  m_finderCfg.deltaRMiddleMaxSPRange = m_deltaRMiddleMaxSPRange;
383  m_finderCfg.seedConfirmation = m_seedConfirmation;
384  m_finderCfg.centralSeedConfirmationRange.zMinSeedConf = m_seedConfCentralZMin;
385  m_finderCfg.centralSeedConfirmationRange.zMaxSeedConf = m_seedConfCentralZMax;
386  m_finderCfg.centralSeedConfirmationRange.rMaxSeedConf = m_seedConfCentralRMax;
387  m_finderCfg.centralSeedConfirmationRange.nTopForLargeR = m_seedConfCentralNTopLargeR;
388  m_finderCfg.centralSeedConfirmationRange.nTopForSmallR = m_seedConfCentralNTopSmallR;
389  m_finderCfg.centralSeedConfirmationRange.seedConfMinBottomRadius = m_seedConfCentralMinBottomRadius;
390  m_finderCfg.centralSeedConfirmationRange.seedConfMaxZOrigin = m_seedConfCentralMaxZOrigin;
391  m_finderCfg.centralSeedConfirmationRange.minImpactSeedConf = m_seedConfCentralMinImpact;
392  m_finderCfg.forwardSeedConfirmationRange.zMinSeedConf = m_seedConfForwardZMin;
393  m_finderCfg.forwardSeedConfirmationRange.zMaxSeedConf = m_seedConfForwardZMax;
394  m_finderCfg.forwardSeedConfirmationRange.rMaxSeedConf = m_seedConfForwardRMax;
395  m_finderCfg.forwardSeedConfirmationRange.nTopForLargeR = m_seedConfForwardNTopLargeR;
396  m_finderCfg.forwardSeedConfirmationRange.nTopForSmallR = m_seedConfForwardNTopSmallR;
397  m_finderCfg.forwardSeedConfirmationRange.seedConfMinBottomRadius = m_seedConfForwardMinBottomRadius;
398  m_finderCfg.forwardSeedConfirmationRange.seedConfMaxZOrigin = m_seedConfForwardMaxZOrigin;
399  m_finderCfg.forwardSeedConfirmationRange.minImpactSeedConf = m_seedConfForwardMinImpact;
400  m_finderCfg.useDetailedDoubleMeasurementInfo = m_useDetailedDoubleMeasurementInfo;
401  m_finderCfg.toleranceParam = m_toleranceParam;
402  m_finderCfg.phiMin = m_phiMin;
403  m_finderCfg.phiMax = m_phiMax;
404  m_finderCfg.rMin = m_rMin;
405  m_finderCfg.zAlign = m_zAlign;
406  m_finderCfg.rAlign = m_rAlign;
407  m_finderCfg.sigmaError = m_sigmaError;
408 
409  // Fast tracking
410  // manually convert the two types
411  for (const auto& vec : m_rRangeMiddleSP) {
412  std::vector<float> convertedVec;
413 
414  for (const auto& val : vec) {
415  convertedVec.push_back(static_cast<float>(val));
416  }
417 
418  m_finderCfg.rRangeMiddleSP.push_back(convertedVec);
419  }
420 
421  // define cuts used for fast tracking configuration
422  if (m_useExperimentCuts) {
423 
424  // This function will be applied to select space points during grid filling
425  m_finderCfg.spacePointSelector
426  .connect<itkFastTrackingSPselect>();
427 
428  m_finderCfg.experimentCuts
429  .connect<itkFastDoubletCut>();
430 
431  }
432 
433  // Configuration for Acts::SeedFilter (used by FinderCfg)
434  Acts::SeedFilterConfig filterCfg;
435  filterCfg.deltaRMin = m_deltaRMin;
436  filterCfg.maxSeedsPerSpM = m_maxSeedsPerSpM;
437  filterCfg.useDeltaRorTopRadius = m_useDeltaRorTopRadius;
438  filterCfg.seedConfirmation = m_seedConfirmationInFilter;
439  filterCfg.maxSeedsPerSpMConf = m_maxSeedsPerSpMConf;
440  filterCfg.maxQualitySeedsPerSpMConf = m_maxQualitySeedsPerSpMConf;
441  filterCfg.centralSeedConfirmationRange = m_finderCfg.centralSeedConfirmationRange;
442  filterCfg.forwardSeedConfirmationRange = m_finderCfg.forwardSeedConfirmationRange;
443  filterCfg.impactWeightFactor = m_impactWeightFactor;
444  filterCfg.zOriginWeightFactor = m_zOriginWeightFactor;
445  filterCfg.compatSeedWeight = m_compatSeedWeight;
446  filterCfg.compatSeedLimit = m_compatSeedLimit;
447  filterCfg.seedWeightIncrement = m_seedWeightIncrement;
448  filterCfg.numSeedIncrement = m_numSeedIncrement;
449  filterCfg.deltaInvHelixDiameter = m_deltaInvHelixDiameter;
450  m_finderCfg.seedFilter = std::make_unique<Acts::SeedFilter< value_type > >(filterCfg.toInternalUnits(), logger().cloneWithSuffix("Filter"));
451 
452  m_finderCfg = m_finderCfg.toInternalUnits().calculateDerivedQuantities();
453 
454  // Grid Configuration
455  m_gridCfg.minPt = m_minPt;
456  m_gridCfg.cotThetaMax = m_cotThetaMax;
457  m_gridCfg.impactMax = m_impactMax;
458  m_gridCfg.zMin = m_zMin;
459  m_gridCfg.zMax = m_zMax;
460  m_gridCfg.phiMin = m_gridPhiMin;
461  m_gridCfg.phiMax = m_gridPhiMax;
462  m_gridCfg.zBinEdges = m_zBinEdges;
463  m_gridCfg.rBinEdges = m_rBinEdges;
464  m_gridCfg.deltaRMax = m_deltaRMax;
465  m_gridCfg.rMax = m_gridRMax;
466  m_gridCfg.phiBinDeflectionCoverage = m_phiBinDeflectionCoverage;
467  m_gridCfg.maxPhiBins = m_maxPhiBins;
468  m_gridCfg = m_gridCfg.toInternalUnits();
469 
470  // Seed Finder
471  m_finder = {m_finderCfg, logger().cloneWithSuffix("Finder")};
472 
473  return StatusCode::SUCCESS;
474  }
475 
476 } // namespace ActsTrk
ActsTrk::SeedingTool::m_deltaRMaxTopSP
Gaudi::Property< float > m_deltaRMaxTopSP
Definition: SeedingTool.h:135
ActsTrk::SeedingTool::m_gridPhiMax
Gaudi::Property< float > m_gridPhiMax
Definition: SeedingTool.h:120
ActsTrk::SeedingTool::m_compatSeedLimit
Gaudi::Property< std::size_t > m_compatSeedLimit
Definition: SeedingTool.h:238
ActsTrk::SeedingTool::m_seedConfCentralMaxZOrigin
Gaudi::Property< float > m_seedConfCentralMaxZOrigin
Definition: SeedingTool.h:190
ActsTrk::SeedingTool::m_finderCfg
Acts::SeedFinderConfig< value_type > m_finderCfg
Definition: SeedingTool.h:91
ActsTrk::SeedingTool::m_seedConfCentralMinBottomRadius
Gaudi::Property< float > m_seedConfCentralMinBottomRadius
Definition: SeedingTool.h:187
ActsTrk::SeedingTool::m_compatSeedWeight
Gaudi::Property< float > m_compatSeedWeight
Definition: SeedingTool.h:236
ActsTrk::SeedingTool::m_sigmaScattering
Gaudi::Property< float > m_sigmaScattering
Definition: SeedingTool.h:147
ActsTrk::SeedingTool::m_maxSeedsPerSpM
Gaudi::Property< int > m_maxSeedsPerSpM
Definition: SeedingTool.h:153
ActsTrk::SeedingTool::m_bottomBinFinder
std::unique_ptr< Acts::GridBinFinder< 3ul > > m_bottomBinFinder
Definition: SeedingTool.h:274
ActsTrk::SeedingTool::m_seedConfForwardMaxZOrigin
Gaudi::Property< float > m_seedConfForwardMaxZOrigin
Definition: SeedingTool.h:214
ActsTrk::SeedingTool::m_maxQualitySeedsPerSpMConf
Gaudi::Property< std::size_t > m_maxQualitySeedsPerSpMConf
Definition: SeedingTool.h:248
ActsTrk::SeedingTool::m_seedConfirmation
Gaudi::Property< bool > m_seedConfirmation
Definition: SeedingTool.h:170
max
constexpr double max()
Definition: ap_fixedTest.cxx:33
ActsTrk::SeedingTool::m_phiBinDeflectionCoverage
Gaudi::Property< int > m_phiBinDeflectionCoverage
Definition: SeedingTool.h:122
ActsTrk::SeedingTool::m_phiMin
Gaudi::Property< float > m_phiMin
Definition: SeedingTool.h:225
ActsTrk::SeedingTool::m_deltaRMinTopSP
Gaudi::Property< float > m_deltaRMinTopSP
Definition: SeedingTool.h:133
ActsTrk::SeedingTool::m_deltaRMaxBottomSP
Gaudi::Property< float > m_deltaRMaxBottomSP
Definition: SeedingTool.h:139
min
constexpr double min()
Definition: ap_fixedTest.cxx:26
ActsTrk::SeedingTool::itkFastTrackingSPselect
static bool itkFastTrackingSPselect(const value_type &sp)
Definition: SeedingTool.h:286
ActsTrk::SeedingTool::m_useVariableMiddleSPRange
Gaudi::Property< bool > m_useVariableMiddleSPRange
Definition: SeedingTool.h:161
Amg::y
@ y
Definition: GeoPrimitives.h:35
ActsTrk::SeedingTool::external_spacepoint::type
typename std::conditional< std::is_pointer< typename spacepoint_iterator_t::value_type >::value, typename std::remove_const< typename std::remove_pointer< typename spacepoint_iterator_t::value_type >::type >::type, typename std::remove_const< typename spacepoint_iterator_t::value_type >::type >::type type
Definition: SeedingTool.h:73
ActsTrk::SeedingTool::m_rMin
Gaudi::Property< float > m_rMin
Definition: SeedingTool.h:227
ActsTrk::SeedingTool::m_rAlign
Gaudi::Property< float > m_rAlign
Definition: SeedingTool.h:229
ActsTrk::SeedingTool::m_seedConfCentralNTopLargeR
Gaudi::Property< size_t > m_seedConfCentralNTopLargeR
Definition: SeedingTool.h:181
athena.value
value
Definition: athena.py:124
ActsTrk::SeedingTool::m_deltaRMax
Gaudi::Property< float > m_deltaRMax
Definition: SeedingTool.h:107
ActsTrk::SeedingTool::m_deltaRMinBottomSP
Gaudi::Property< float > m_deltaRMinBottomSP
Definition: SeedingTool.h:137
vec
std::vector< size_t > vec
Definition: CombinationsGeneratorTest.cxx:12
ActsTrk::SeedingTool::createSeeds
virtual StatusCode createSeeds(const EventContext &ctx, const Acts::SpacePointContainer< ActsTrk::SpacePointCollector, Acts::detail::RefHolder > &spContainer, const Acts::Vector3 &beamSpotPos, const Acts::Vector3 &bField, ActsTrk::SeedContainer &seedContainer) const override
Definition: SeedingTool.cxx:206
ActsTrk::SeedingTool::m_finder
Acts::SeedFinder< value_type, Acts::CylindricalSpacePointGrid< value_type > > m_finder
Definition: SeedingTool.h:90
ActsTrk::SeedingTool::m_useExperimentCuts
Gaudi::Property< bool > m_useExperimentCuts
Definition: SeedingTool.h:271
ActsTrk::SeedingTool::m_rMax
Gaudi::Property< float > m_rMax
Definition: SeedingTool.h:127
ActsTrk::SeedingTool::m_zAlign
Gaudi::Property< float > m_zAlign
Definition: SeedingTool.h:228
ActsTrk::SeedingTool::m_seedConfCentralZMax
Gaudi::Property< float > m_seedConfCentralZMax
Definition: SeedingTool.h:175
ActsTrk::SeedingTool::m_zBinNeighborsTop
Gaudi::Property< std::vector< std::pair< int, int > > > m_zBinNeighborsTop
Definition: SeedingTool.h:257
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
ActsTrk::SeedingTool::m_impactMax
Gaudi::Property< float > m_impactMax
Definition: SeedingTool.h:109
ActsTrk::SeedingTool::m_impactWeightFactor
Gaudi::Property< float > m_impactWeightFactor
Definition: SeedingTool.h:233
ActsTrk::SeedingTool::m_seedQualitySelection
Gaudi::Property< bool > m_seedQualitySelection
Definition: SeedingTool.h:95
ActsTrk::SeedingTool::m_zBinEdges
Gaudi::Property< std::vector< float > > m_zBinEdges
Definition: SeedingTool.h:111
ActsTrk::SeedingTool::m_zMax
Gaudi::Property< float > m_zMax
Definition: SeedingTool.h:105
ATH_MSG_ERROR
#define ATH_MSG_ERROR(x)
Definition: AthMsgStreamMacros.h:33
ActsTrk::SeedingTool::m_seedConfForwardZMax
Gaudi::Property< float > m_seedConfForwardZMax
Definition: SeedingTool.h:199
lumiFormat.i
int i
Definition: lumiFormat.py:85
ActsTrk::SeedingTool::m_minPt
Gaudi::Property< float > m_minPt
Definition: SeedingTool.h:99
ActsTrk::SeedingTool::m_deltaZMax
Gaudi::Property< float > m_deltaZMax
Definition: SeedingTool.h:141
ActsTrk::SeedingTool::m_maxSeedsPerSpMConf
Gaudi::Property< std::size_t > m_maxSeedsPerSpMConf
Definition: SeedingTool.h:246
EL::StatusCode
::StatusCode StatusCode
StatusCode definition for legacy code.
Definition: PhysicsAnalysis/D3PDTools/EventLoop/EventLoop/StatusCode.h:22
ATH_MSG_DEBUG
#define ATH_MSG_DEBUG(x)
Definition: AthMsgStreamMacros.h:29
Amg::x
@ x
Definition: GeoPrimitives.h:34
ActsTrk::SeedingTool::m_binSizeR
Gaudi::Property< float > m_binSizeR
Definition: SeedingTool.h:129
ActsTrk::SeedingTool::m_topBinFinder
std::unique_ptr< Acts::GridBinFinder< 3ul > > m_topBinFinder
Definition: SeedingTool.h:275
test_pyathena.parent
parent
Definition: test_pyathena.py:15
ATH_CHECK
#define ATH_CHECK
Definition: AthCheckMacros.h:40
ActsTrk::SeedingTool::m_seedConfForwardRMax
Gaudi::Property< float > m_seedConfForwardRMax
Definition: SeedingTool.h:202
detail::ul
unsigned long ul
Definition: PrimitiveHelpers.h:46
ActsTrk::SeedingTool::m_seedConfCentralMinImpact
Gaudi::Property< float > m_seedConfCentralMinImpact
Definition: SeedingTool.h:193
ActsTrk::SeedingTool::m_seedConfForwardMinImpact
Gaudi::Property< float > m_seedConfForwardMinImpact
Definition: SeedingTool.h:217
ActsTrk::SeedingTool::m_rBinNeighborsTop
Gaudi::Property< std::vector< std::pair< int, int > > > m_rBinNeighborsTop
Definition: SeedingTool.h:264
DataVector
Derived DataVector<T>.
Definition: DataVector.h:794
ActsTrk::SeedingTool::m_gridRMax
Gaudi::Property< float > m_gridRMax
Definition: SeedingTool.h:116
ActsTrk::SeedingTool::m_radLengthPerSeed
Gaudi::Property< float > m_radLengthPerSeed
Definition: SeedingTool.h:151
DeMoUpdate.toAdd
bool toAdd
Definition: DeMoUpdate.py:1304
ActsTrk::SeedingTool::m_deltaInvHelixDiameter
Gaudi::Property< float > m_deltaInvHelixDiameter
Definition: SeedingTool.h:252
ActsTrk::SeedingTool::logger
const Acts::Logger & logger() const
Private access to the logger.
Definition: SeedingTool.h:280
ActsTrk::SeedingTool::m_seedConfCentralRMax
Gaudi::Property< float > m_seedConfCentralRMax
Definition: SeedingTool.h:178
ActsTrk::SeedingTool::m_collisionRegionMax
Gaudi::Property< float > m_collisionRegionMax
Definition: SeedingTool.h:145
ActsTrk::SeedingTool::m_rBinNeighborsBottom
Gaudi::Property< std::vector< std::pair< int, int > > > m_rBinNeighborsBottom
Definition: SeedingTool.h:266
ActsTrk::SeedingTool::m_rBinsCustomLooping
Gaudi::Property< std::vector< std::size_t > > m_rBinsCustomLooping
Definition: SeedingTool.h:159
name
std::string name
Definition: Control/AthContainers/Root/debug.cxx:228
ActsTrk::SeedingTool::m_deltaRMiddleMaxSPRange
Gaudi::Property< float > m_deltaRMiddleMaxSPRange
Definition: SeedingTool.h:168
ActsTrk::SeedingTool::m_deltaRMiddleMinSPRange
Gaudi::Property< float > m_deltaRMiddleMinSPRange
Definition: SeedingTool.h:166
ActsTrk::SeedingTool::m_navigation
std::array< std::vector< std::size_t >, 3ul > m_navigation
Definition: SeedingTool.h:277
ActsTrk::SeedingTool::m_seedConfCentralZMin
Gaudi::Property< float > m_seedConfCentralZMin
Definition: SeedingTool.h:172
ActsTrk::SeedingTool::itkFastDoubletCut
static bool itkFastDoubletCut(float bottomRadius, float cotTheta)
Definition: SeedingTool.h:304
ActsTrk::SeedingTool::m_zBinsCustomLooping
Gaudi::Property< std::vector< size_t > > m_zBinsCustomLooping
Definition: SeedingTool.h:157
ActsTrk::SeedingTool::prepareConfiguration
StatusCode prepareConfiguration()
Definition: SeedingTool.cxx:349
ActsTrk::SeedingTool::m_phiMax
Gaudi::Property< float > m_phiMax
Definition: SeedingTool.h:226
ActsTrk::SeedingTool::m_collisionRegionMin
Gaudi::Property< float > m_collisionRegionMin
Definition: SeedingTool.h:143
ActsTrk::SeedingTool::m_seedConfCentralNTopSmallR
Gaudi::Property< size_t > m_seedConfCentralNTopSmallR
Definition: SeedingTool.h:184
ActsTrk::SeedingTool::m_cotThetaMax
Gaudi::Property< float > m_cotThetaMax
Definition: SeedingTool.h:101
ActsTrk::SeedingTool::m_maxPhiBins
Gaudi::Property< int > m_maxPhiBins
Definition: SeedingTool.h:124
ActsTrk::SeedingTool::initialize
virtual StatusCode initialize() override
Definition: SeedingTool.cxx:24
ActsTrk::SeedingTool::m_toleranceParam
Gaudi::Property< float > m_toleranceParam
Definition: SeedingTool.h:223
ActsTrk::SeedingTool::seed_type
Acts::Seed< value_type, 3ul > seed_type
Definition: SeedingTool.h:43
ActsTrk::SeedingTool::m_deltaRMin
Gaudi::Property< float > m_deltaRMin
Definition: SeedingTool.h:131
ActsTrk::SeedingTool::m_seedWeightIncrement
Gaudi::Property< float > m_seedWeightIncrement
Definition: SeedingTool.h:240
Pythia8_RapidityOrderMPI.val
val
Definition: Pythia8_RapidityOrderMPI.py:14
python.CaloScaleNoiseConfig.type
type
Definition: CaloScaleNoiseConfig.py:78
ActsTrk::SeedingTool::m_numPhiNeighbors
Gaudi::Property< int > m_numPhiNeighbors
Definition: SeedingTool.h:268
ActsTrk::SeedingTool::m_seedConfForwardNTopLargeR
Gaudi::Property< size_t > m_seedConfForwardNTopLargeR
Definition: SeedingTool.h:205
ActsTrk::SeedingTool::m_seedConfForwardZMin
Gaudi::Property< float > m_seedConfForwardZMin
Definition: SeedingTool.h:196
ActsTrk::SeedingTool::m_seedConfirmationInFilter
Gaudi::Property< bool > m_seedConfirmationInFilter
Definition: SeedingTool.h:244
ActsTrk::SeedingTool::m_rBinEdges
Gaudi::Property< std::vector< float > > m_rBinEdges
Definition: SeedingTool.h:114
ActsTrk::SeedingTool::m_interactionPointCut
Gaudi::Property< bool > m_interactionPointCut
Definition: SeedingTool.h:155
ActsTrk::SeedingTool::m_rRangeMiddleSP
Gaudi::Property< std::vector< std::vector< double > > > m_rRangeMiddleSP
Definition: SeedingTool.h:163
ActsTrk::SeedingTool::m_useDeltaRorTopRadius
Gaudi::Property< bool > m_useDeltaRorTopRadius
Definition: SeedingTool.h:250
SeedingTool.h
ActsTrk
The AlignStoreProviderAlg loads the rigid alignment corrections and pipes them through the readout ge...
Definition: MuonDetectorBuilderTool.cxx:54
top
@ top
Definition: TruthClasses.h:64
ActsTrk::SeedingTool::m_seedConfForwardNTopSmallR
Gaudi::Property< size_t > m_seedConfForwardNTopSmallR
Definition: SeedingTool.h:208
ActsTrk::SeedingTool::m_gridCfg
Acts::CylindricalSpacePointGridConfig m_gridCfg
Definition: SeedingTool.h:92
ActsTrk::SeedingTool::m_zOriginWeightFactor
Gaudi::Property< float > m_zOriginWeightFactor
Definition: SeedingTool.h:235
ActsTrk::SeedingTool::m_useDetailedDoubleMeasurementInfo
Gaudi::Property< bool > m_useDetailedDoubleMeasurementInfo
Definition: SeedingTool.h:220
ActsTrk::SeedingTool::m_zBinNeighborsBottom
Gaudi::Property< std::vector< std::pair< int, int > > > m_zBinNeighborsBottom
Definition: SeedingTool.h:261
value_type
Definition: EDM_MasterSearch.h:11
Amg::distance
float distance(const Amg::Vector3D &p1, const Amg::Vector3D &p2)
calculates the distance between two point in 3D space
Definition: GeoPrimitivesHelpers.h:54
ActsTrk::SeedingTool::m_gridPhiMin
Gaudi::Property< float > m_gridPhiMin
Definition: SeedingTool.h:118
ActsTrk::SeedingTool::m_sigmaError
Gaudi::Property< float > m_sigmaError
Definition: SeedingTool.h:230
ActsTrk::SeedingTool::m_logger
std::unique_ptr< const Acts::Logger > m_logger
logging instance
Definition: SeedingTool.h:282
ActsTrk::SeedingTool::m_maxPtScattering
Gaudi::Property< float > m_maxPtScattering
Definition: SeedingTool.h:149
ActsTrk::SeedingTool::m_numSeedIncrement
Gaudi::Property< float > m_numSeedIncrement
Definition: SeedingTool.h:242
ActsTrk::SeedingTool::m_seedConfForwardMinBottomRadius
Gaudi::Property< float > m_seedConfForwardMinBottomRadius
Definition: SeedingTool.h:211
ActsTrk::SeedingTool::m_zMin
Gaudi::Property< float > m_zMin
Definition: SeedingTool.h:103
Generated on Fri Jan 10 2025 21:17:27 for ATLAS Offline Software by   doxygen 1.8.18