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