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  DataVector< ActsTrk::ActsSeed< typename SeedingTool::external_type, 3ul > >& 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  std::vector< seed_type > seeds;
250 
251  // Space Point Grid Options
252  Acts::CylindricalSpacePointGridOptions gridOpts;
253  gridOpts.bFieldInZ = bField[2];
254  gridOpts = gridOpts.toInternalUnits();
255 
256  // Seed Finder Options
257  Acts::SeedFinderOptions finderOpts;
258  finderOpts.beamPos = Acts::Vector2(beamSpotPos[Amg::x],
259  beamSpotPos[Amg::y]);
260  finderOpts.bFieldInZ = bField[2];
261  finderOpts = finderOpts.toInternalUnits().calculateDerivedQuantities(m_finderCfg);
262 
263 
264 
265 
266  Acts::CylindricalSpacePointGrid< value_type > grid =
267  Acts::CylindricalSpacePointGridCreator::createGrid< value_type >(m_gridCfg, gridOpts, logger());
268 
269  Acts::CylindricalSpacePointGridCreator::fillGrid(m_finderCfg, finderOpts, grid,
270  spBegin, spEnd, logger());
271 
272  // Compute radius Range
273  // we rely on the fact the grid is storing the proxies
274  // with a sorting in the radius
275  float minRange = std::numeric_limits<float>::max();
276  float maxRange = std::numeric_limits<float>::lowest();
277  for (const auto& coll : grid) {
278  if (coll.empty()) {
279  continue;
280  }
281  const auto* firstEl = coll.front();
282  const auto* lastEl = coll.back();
283  minRange = std::min(firstEl->radius(), minRange);
284  maxRange = std::max(lastEl->radius(), maxRange);
285  }
286 
287  Acts::CylindricalBinnedGroup< value_type > spacePointsGrouping(std::move(grid), *m_bottomBinFinder,
288  *m_topBinFinder, m_navigation);
289 
290  // variable middle SP radial region of interest
291  const Acts::Range1D<float> rMiddleSPRange(std::floor(minRange/2)*2 + m_finderCfg.deltaRMiddleMinSPRange,
292  std::floor(maxRange/2)*2 - m_finderCfg.deltaRMiddleMaxSPRange);
293 
294  //TODO POSSIBLE OPTIMISATION come back here: see MR !52399 ( i.e. use static thread_local)
295  typename decltype(m_finder)::SeedingState state;
296 
297  // Already reserve the state vectors here.
298 
299  state.topSpVec.reserve(m_stateVectorReserveSize);
300  state.curvatures.reserve(m_stateVectorReserveSize);
301  state.impactParameters.reserve(m_stateVectorReserveSize);
302  state.linCircleTop.reserve(m_stateVectorReserveSize);
303  state.compatBottomSP.reserve(m_stateVectorReserveSize);
304  state.compatTopSP.reserve(m_stateVectorReserveSize);
305 
306  state.spacePointMutableData.resize(std::distance(spBegin, spEnd));
307 
308  for (const auto& [bottom, middle, top] : spacePointsGrouping) {
309  m_finder.createSeedsForGroup(finderOpts, state, spacePointsGrouping.grid(),
310  seeds, bottom, middle, top, rMiddleSPRange);
311  }
312 
313  if (m_seedQualitySelection) {
314  // Selection function - temporary implementation
315  // need change from ACTS for final implementation
316  // To be used only on PPP
317  auto selectionFunction = [&state] (const seed_type& seed) -> bool
318  {
319  float seed_quality = seed.seedQuality();
320  float bottom_quality = state.spacePointMutableData.quality(seed.sp()[0]->index());
321  float middle_quality = state.spacePointMutableData.quality(seed.sp()[1]->index());
322  float top_quality = state.spacePointMutableData.quality(seed.sp()[2]->index());
323 
324  if (bottom_quality > seed_quality and
325  middle_quality > seed_quality and
326  top_quality > seed_quality) {
327  return false;
328  }
329 
330  return true;
331  };
332 
333  // Select the seeds
334  std::size_t acceptedSeeds = 0;
335  for (std::size_t i(0); i<seeds.size(); ++i) {
336  const auto& seed = seeds[i];
337  if (not selectionFunction(seed)) {
338  continue;
339  }
340  // move passing seeds at the beginning of the vector/collection
341  // no need to swap them both. The seed currently at acceptedSeeds
342  // didn't make it (if i != acceptedSeeds)
343  if (acceptedSeeds != i)
344  seeds[acceptedSeeds] = std::move(seeds[i]);
345  ++acceptedSeeds;
346  }
347  // remove seeds that didn't make it
348  // they are all at the end of the collection
349  seeds.erase(seeds.begin() + acceptedSeeds, seeds.end());
350  }
351 
352 
353  // Store seeds
354  seedContainer.reserve(seeds.size());
355  for(const auto& seed: seeds) {
356  assert(seed.sp().size() == 3ul);
357  const auto bottom = seed.sp().at(0);
358  const auto middle = seed.sp().at(1);
359  const auto top = seed.sp().at(2);
360 
361  std::unique_ptr< ActsTrk::Seed > toAdd =
362  std::make_unique< ActsTrk::Seed >(bottom->externalSpacePoint(),
363  middle->externalSpacePoint(),
364  top->externalSpacePoint());
365  toAdd->setVertexZ(seed.z());
366  toAdd->setQuality(seed.seedQuality());
367  seedContainer.push_back(std::move(toAdd));
368  }
369 
370  return StatusCode::SUCCESS;
371  }
372 
373  StatusCode
374  SeedingTool::prepareConfiguration() {
375  assert(m_logger != nullptr);
376 
377  // Prepare the Acts::SeedFinderConfig object
378  // This is done only once, during initialization using the
379  // parameters set in the JO
380 
381  // Configuration for Acts::SeedFinder
382  m_finderCfg.minPt = m_minPt;
383  m_finderCfg.cotThetaMax = m_cotThetaMax;
384  m_finderCfg.impactMax = m_impactMax;
385  m_finderCfg.zMin = m_zMin;
386  m_finderCfg.zMax = m_zMax;
387  m_finderCfg.zBinEdges = m_zBinEdges;
388  m_finderCfg.rMax = m_rMax;
389  m_finderCfg.binSizeR = m_binSizeR;
390  m_finderCfg.deltaRMin = m_deltaRMin;
391  m_finderCfg.deltaRMax = m_deltaRMax;
392  m_finderCfg.deltaRMinTopSP = m_deltaRMinTopSP;
393  m_finderCfg.deltaRMaxTopSP = m_deltaRMaxTopSP;
394  m_finderCfg.deltaRMinBottomSP = m_deltaRMinBottomSP;
395  m_finderCfg.deltaRMaxBottomSP = m_deltaRMaxBottomSP;
396  m_finderCfg.deltaZMax = m_deltaZMax;
397  m_finderCfg.collisionRegionMin = m_collisionRegionMin;
398  m_finderCfg.collisionRegionMax = m_collisionRegionMax;
399  m_finderCfg.sigmaScattering = m_sigmaScattering;
400  m_finderCfg.maxPtScattering = m_maxPtScattering;
401  m_finderCfg.radLengthPerSeed = m_radLengthPerSeed;
402  m_finderCfg.maxSeedsPerSpM = m_maxSeedsPerSpM;
403  m_finderCfg.interactionPointCut = m_interactionPointCut;
404  m_finderCfg.zBinsCustomLooping = m_zBinsCustomLooping;
405  m_finderCfg.useVariableMiddleSPRange = m_useVariableMiddleSPRange;
406  m_finderCfg.deltaRMiddleMinSPRange = m_deltaRMiddleMinSPRange;
407  m_finderCfg.deltaRMiddleMaxSPRange = m_deltaRMiddleMaxSPRange;
408  m_finderCfg.seedConfirmation = m_seedConfirmation;
409  m_finderCfg.centralSeedConfirmationRange.zMinSeedConf = m_seedConfCentralZMin;
410  m_finderCfg.centralSeedConfirmationRange.zMaxSeedConf = m_seedConfCentralZMax;
411  m_finderCfg.centralSeedConfirmationRange.rMaxSeedConf = m_seedConfCentralRMax;
412  m_finderCfg.centralSeedConfirmationRange.nTopForLargeR = m_seedConfCentralNTopLargeR;
413  m_finderCfg.centralSeedConfirmationRange.nTopForSmallR = m_seedConfCentralNTopSmallR;
414  m_finderCfg.centralSeedConfirmationRange.seedConfMinBottomRadius = m_seedConfCentralMinBottomRadius;
415  m_finderCfg.centralSeedConfirmationRange.seedConfMaxZOrigin = m_seedConfCentralMaxZOrigin;
416  m_finderCfg.centralSeedConfirmationRange.minImpactSeedConf = m_seedConfCentralMinImpact;
417  m_finderCfg.forwardSeedConfirmationRange.zMinSeedConf = m_seedConfForwardZMin;
418  m_finderCfg.forwardSeedConfirmationRange.zMaxSeedConf = m_seedConfForwardZMax;
419  m_finderCfg.forwardSeedConfirmationRange.rMaxSeedConf = m_seedConfForwardRMax;
420  m_finderCfg.forwardSeedConfirmationRange.nTopForLargeR = m_seedConfForwardNTopLargeR;
421  m_finderCfg.forwardSeedConfirmationRange.nTopForSmallR = m_seedConfForwardNTopSmallR;
422  m_finderCfg.forwardSeedConfirmationRange.seedConfMinBottomRadius = m_seedConfForwardMinBottomRadius;
423  m_finderCfg.forwardSeedConfirmationRange.seedConfMaxZOrigin = m_seedConfForwardMaxZOrigin;
424  m_finderCfg.forwardSeedConfirmationRange.minImpactSeedConf = m_seedConfForwardMinImpact;
425  m_finderCfg.useDetailedDoubleMeasurementInfo = m_useDetailedDoubleMeasurementInfo;
426  m_finderCfg.toleranceParam = m_toleranceParam;
427  m_finderCfg.phiMin = m_phiMin;
428  m_finderCfg.phiMax = m_phiMax;
429  m_finderCfg.rMin = m_rMin;
430  m_finderCfg.zAlign = m_zAlign;
431  m_finderCfg.rAlign = m_rAlign;
432  m_finderCfg.sigmaError = m_sigmaError;
433 
434  // Fast tracking
435  // manually convert the two types
436  for (const auto& vec : m_rRangeMiddleSP) {
437  std::vector<float> convertedVec;
438  convertedVec.reserve(vec.size());
439  for (const auto& val : vec) {
440  convertedVec.push_back(static_cast<float>(val));
441  }
442 
443  m_finderCfg.rRangeMiddleSP.push_back(std::move(convertedVec));
444  }
445 
446 
447  // define cuts used for fast tracking configuration
448  if (m_useExperimentCuts) {
449 
450  // This function will be applied to select space points during grid filling
451  m_finderCfg.spacePointSelector
452  .connect<&ActsTrk::SeedingTool::spacePointSelectionFunction>(this);
453 
454  m_finderCfg.experimentCuts
455  .connect<&ActsTrk::SeedingTool::doubletSelectionFunction>(this);
456 
457  }
458 
459  // Configuration for Acts::SeedFilter (used by FinderCfg)
460  Acts::SeedFilterConfig filterCfg;
461  filterCfg.deltaRMin = m_deltaRMin;
462  filterCfg.maxSeedsPerSpM = m_maxSeedsPerSpM;
463  filterCfg.useDeltaRorTopRadius = m_useDeltaRorTopRadius;
464  filterCfg.seedConfirmation = m_seedConfirmationInFilter;
465  filterCfg.maxSeedsPerSpMConf = m_maxSeedsPerSpMConf;
466  filterCfg.maxQualitySeedsPerSpMConf = m_maxQualitySeedsPerSpMConf;
467  filterCfg.centralSeedConfirmationRange = m_finderCfg.centralSeedConfirmationRange;
468  filterCfg.forwardSeedConfirmationRange = m_finderCfg.forwardSeedConfirmationRange;
469  filterCfg.impactWeightFactor = m_impactWeightFactor;
470  filterCfg.zOriginWeightFactor = m_zOriginWeightFactor;
471  filterCfg.compatSeedWeight = m_compatSeedWeight;
472  filterCfg.compatSeedLimit = m_compatSeedLimit;
473  filterCfg.seedWeightIncrement = m_seedWeightIncrement;
474  filterCfg.numSeedIncrement = m_numSeedIncrement;
475  filterCfg.deltaInvHelixDiameter = m_deltaInvHelixDiameter;
476  m_finderCfg.seedFilter = std::make_unique<Acts::SeedFilter< value_type > >(filterCfg.toInternalUnits(), logger().cloneWithSuffix("Filter"));
477 
478  m_finderCfg = m_finderCfg.toInternalUnits().calculateDerivedQuantities();
479 
480  // Grid Configuration
481  m_gridCfg.minPt = m_minPt;
482  m_gridCfg.cotThetaMax = m_cotThetaMax;
483  m_gridCfg.impactMax = m_impactMax;
484  m_gridCfg.zMin = m_zMin;
485  m_gridCfg.zMax = m_zMax;
486  m_gridCfg.phiMin = m_gridPhiMin;
487  m_gridCfg.phiMax = m_gridPhiMax;
488  m_gridCfg.zBinEdges = m_zBinEdges;
489  m_gridCfg.rBinEdges = m_rBinEdges;
490  m_gridCfg.deltaRMax = m_deltaRMax;
491  m_gridCfg.rMax = m_gridRMax;
492  m_gridCfg.phiBinDeflectionCoverage = m_phiBinDeflectionCoverage;
493  m_gridCfg.maxPhiBins = m_maxPhiBins;
494  m_gridCfg = m_gridCfg.toInternalUnits();
495 
496  // Seed Finder
497  m_finder = decltype(m_finder){m_finderCfg, logger().cloneWithSuffix("Finder")};
498 
499  return StatusCode::SUCCESS;
500  }
501 
502 } // namespace ActsTrk
ActsTrk::SeedingTool::m_deltaRMaxTopSP
Gaudi::Property< float > m_deltaRMaxTopSP
Definition: SeedingTool.h:152
ActsTrk::SeedingTool::m_gridPhiMax
Gaudi::Property< float > m_gridPhiMax
Definition: SeedingTool.h:137
ActsTrk::SeedingTool::m_compatSeedLimit
Gaudi::Property< std::size_t > m_compatSeedLimit
Definition: SeedingTool.h:255
ActsTrk::SeedingTool::m_seedConfCentralMaxZOrigin
Gaudi::Property< float > m_seedConfCentralMaxZOrigin
Definition: SeedingTool.h:207
inline_hints.h
ActsTrk::SeedingTool::m_finderCfg
Acts::SeedFinderConfig< value_type > m_finderCfg
Definition: SeedingTool.h:108
ActsTrk::SeedingTool::m_seedConfCentralMinBottomRadius
Gaudi::Property< float > m_seedConfCentralMinBottomRadius
Definition: SeedingTool.h:204
ActsTrk::SeedingTool::m_compatSeedWeight
Gaudi::Property< float > m_compatSeedWeight
Definition: SeedingTool.h:253
ActsTrk::SeedingTool::m_sigmaScattering
Gaudi::Property< float > m_sigmaScattering
Definition: SeedingTool.h:164
ActsTrk::SeedingTool::m_maxSeedsPerSpM
Gaudi::Property< int > m_maxSeedsPerSpM
Definition: SeedingTool.h:170
ActsTrk::SeedingTool::m_bottomBinFinder
std::unique_ptr< Acts::GridBinFinder< 3ul > > m_bottomBinFinder
Definition: SeedingTool.h:293
ActsTrk::SeedingTool::m_seedConfForwardMaxZOrigin
Gaudi::Property< float > m_seedConfForwardMaxZOrigin
Definition: SeedingTool.h:231
ActsTrk::SeedingTool::m_maxQualitySeedsPerSpMConf
Gaudi::Property< std::size_t > m_maxQualitySeedsPerSpMConf
Definition: SeedingTool.h:265
ActsTrk::SeedingTool::m_seedConfirmation
Gaudi::Property< bool > m_seedConfirmation
Definition: SeedingTool.h:187
max
constexpr double max()
Definition: ap_fixedTest.cxx:33
ActsTrk::SeedingTool::m_phiBinDeflectionCoverage
Gaudi::Property< int > m_phiBinDeflectionCoverage
Definition: SeedingTool.h:139
ActsTrk::SeedingTool::m_phiMin
Gaudi::Property< float > m_phiMin
Definition: SeedingTool.h:242
ActsTrk::SeedingTool::m_deltaRMinTopSP
Gaudi::Property< float > m_deltaRMinTopSP
Definition: SeedingTool.h:150
ActsTrk::SeedingTool::m_deltaRMaxBottomSP
Gaudi::Property< float > m_deltaRMaxBottomSP
Definition: SeedingTool.h:156
min
constexpr double min()
Definition: ap_fixedTest.cxx:26
ActsTrk::SeedingTool::m_useVariableMiddleSPRange
Gaudi::Property< bool > m_useVariableMiddleSPRange
Definition: SeedingTool.h:178
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:87
ActsTrk::SeedingTool::m_rMin
Gaudi::Property< float > m_rMin
Definition: SeedingTool.h:244
ActsTrk::SeedingTool::m_rAlign
Gaudi::Property< float > m_rAlign
Definition: SeedingTool.h:246
ActsTrk::SeedingTool::m_seedConfCentralNTopLargeR
Gaudi::Property< size_t > m_seedConfCentralNTopLargeR
Definition: SeedingTool.h:198
athena.value
value
Definition: athena.py:124
ActsTrk::SeedingTool::m_deltaRMax
Gaudi::Property< float > m_deltaRMax
Definition: SeedingTool.h:124
ActsTrk::SeedingTool::m_deltaRMinBottomSP
Gaudi::Property< float > m_deltaRMinBottomSP
Definition: SeedingTool.h:154
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:107
ActsTrk::SeedingTool::m_useExperimentCuts
Gaudi::Property< bool > m_useExperimentCuts
Definition: SeedingTool.h:288
python.CaloAddPedShiftConfig.type
type
Definition: CaloAddPedShiftConfig.py:42
ActsTrk::SeedingTool::m_rMax
Gaudi::Property< float > m_rMax
Definition: SeedingTool.h:144
ActsTrk::SeedingTool::doubletSelectionFunction
bool doubletSelectionFunction(const value_type &middle, const value_type &other, float cotTheta, bool isBottomCandidate) const
Definition: SeedingTool.h:363
ActsTrk::SeedingTool::m_zAlign
Gaudi::Property< float > m_zAlign
Definition: SeedingTool.h:245
ActsTrk::SeedingTool::m_seedConfCentralZMax
Gaudi::Property< float > m_seedConfCentralZMax
Definition: SeedingTool.h:192
ActsTrk::SeedingTool::m_zBinNeighborsTop
Gaudi::Property< std::vector< std::pair< int, int > > > m_zBinNeighborsTop
Definition: SeedingTool.h:274
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:126
ActsTrk::SeedingTool::m_impactWeightFactor
Gaudi::Property< float > m_impactWeightFactor
Definition: SeedingTool.h:250
ActsTrk::SeedingTool::m_seedQualitySelection
Gaudi::Property< bool > m_seedQualitySelection
Definition: SeedingTool.h:112
ActsTrk::SeedingTool::m_zBinEdges
Gaudi::Property< std::vector< float > > m_zBinEdges
Definition: SeedingTool.h:128
ActsTrk::SeedingTool::m_zMax
Gaudi::Property< float > m_zMax
Definition: SeedingTool.h:122
ATH_MSG_ERROR
#define ATH_MSG_ERROR(x)
Definition: AthMsgStreamMacros.h:33
ActsTrk::SeedingTool::m_seedConfForwardZMax
Gaudi::Property< float > m_seedConfForwardZMax
Definition: SeedingTool.h:216
lumiFormat.i
int i
Definition: lumiFormat.py:85
ActsTrk::SeedingTool::m_minPt
Gaudi::Property< float > m_minPt
Definition: SeedingTool.h:116
ActsTrk::SeedingTool::m_deltaZMax
Gaudi::Property< float > m_deltaZMax
Definition: SeedingTool.h:158
ActsTrk::SeedingTool::m_maxSeedsPerSpMConf
Gaudi::Property< std::size_t > m_maxSeedsPerSpMConf
Definition: SeedingTool.h:263
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:307
ActsTrk::SeedingTool::m_binSizeR
Gaudi::Property< float > m_binSizeR
Definition: SeedingTool.h:146
ActsTrk::SeedingTool::m_topBinFinder
std::unique_ptr< Acts::GridBinFinder< 3ul > > m_topBinFinder
Definition: SeedingTool.h:294
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:219
detail::ul
unsigned long ul
Definition: PrimitiveHelpers.h:46
ActsTrk::SeedingTool::m_seedConfCentralMinImpact
Gaudi::Property< float > m_seedConfCentralMinImpact
Definition: SeedingTool.h:210
ActsTrk::SeedingTool::m_seedConfForwardMinImpact
Gaudi::Property< float > m_seedConfForwardMinImpact
Definition: SeedingTool.h:234
ActsTrk::SeedingTool::m_rBinNeighborsTop
Gaudi::Property< std::vector< std::pair< int, int > > > m_rBinNeighborsTop
Definition: SeedingTool.h:281
DataVector
Derived DataVector<T>.
Definition: DataVector.h:794
ActsTrk::SeedingTool::m_gridRMax
Gaudi::Property< float > m_gridRMax
Definition: SeedingTool.h:133
ActsTrk::SeedingTool::m_radLengthPerSeed
Gaudi::Property< float > m_radLengthPerSeed
Definition: SeedingTool.h:168
DeMoUpdate.toAdd
bool toAdd
Definition: DeMoUpdate.py:1304
ActsTrk::SeedingTool::m_deltaInvHelixDiameter
Gaudi::Property< float > m_deltaInvHelixDiameter
Definition: SeedingTool.h:269
ActsTrk::SeedingTool::logger
const Acts::Logger & logger() const
Private access to the logger.
Definition: SeedingTool.h:299
ActsTrk::SeedingTool::m_seedConfCentralRMax
Gaudi::Property< float > m_seedConfCentralRMax
Definition: SeedingTool.h:195
ActsTrk::SeedingTool::m_collisionRegionMax
Gaudi::Property< float > m_collisionRegionMax
Definition: SeedingTool.h:162
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:283
ActsTrk::SeedingTool::m_rBinsCustomLooping
Gaudi::Property< std::vector< std::size_t > > m_rBinsCustomLooping
Definition: SeedingTool.h:176
name
std::string name
Definition: Control/AthContainers/Root/debug.cxx:240
ActsTrk::SeedingTool::m_deltaRMiddleMaxSPRange
Gaudi::Property< float > m_deltaRMiddleMaxSPRange
Definition: SeedingTool.h:185
ActsTrk::SeedingTool::m_deltaRMiddleMinSPRange
Gaudi::Property< float > m_deltaRMiddleMinSPRange
Definition: SeedingTool.h:183
ActsTrk::SeedingTool::m_navigation
std::array< std::vector< std::size_t >, 3ul > m_navigation
Definition: SeedingTool.h:296
ActsTrk::SeedingTool::m_seedConfCentralZMin
Gaudi::Property< float > m_seedConfCentralZMin
Definition: SeedingTool.h:189
ActsTrk::SeedingTool::m_zBinsCustomLooping
Gaudi::Property< std::vector< size_t > > m_zBinsCustomLooping
Definition: SeedingTool.h:174
ActsTrk::SeedingTool::prepareConfiguration
StatusCode prepareConfiguration()
Definition: SeedingTool.cxx:374
ActsTrk::SeedingTool::m_phiMax
Gaudi::Property< float > m_phiMax
Definition: SeedingTool.h:243
ActsTrk::SeedingTool::m_collisionRegionMin
Gaudi::Property< float > m_collisionRegionMin
Definition: SeedingTool.h:160
ActsTrk::SeedingTool::m_seedConfCentralNTopSmallR
Gaudi::Property< size_t > m_seedConfCentralNTopSmallR
Definition: SeedingTool.h:201
ActsTrk::ActsSeed
Definition: Seed.h:18
ActsTrk::SeedingTool::m_cotThetaMax
Gaudi::Property< float > m_cotThetaMax
Definition: SeedingTool.h:118
ActsTrk::SeedingTool::m_maxPhiBins
Gaudi::Property< int > m_maxPhiBins
Definition: SeedingTool.h:141
ActsTrk::SeedingTool::initialize
virtual StatusCode initialize() override
Definition: SeedingTool.cxx:30
ActsTrk::SeedingTool::m_toleranceParam
Gaudi::Property< float > m_toleranceParam
Definition: SeedingTool.h:240
ActsTrk::SeedingTool::m_deltaRMin
Gaudi::Property< float > m_deltaRMin
Definition: SeedingTool.h:148
ActsTrk::SeedingTool::m_stateVectorReserveSize
Gaudi::Property< int > m_stateVectorReserveSize
Definition: SeedingTool.h:290
ActsTrk::SeedingTool::m_seedWeightIncrement
Gaudi::Property< float > m_seedWeightIncrement
Definition: SeedingTool.h:257
Pythia8_RapidityOrderMPI.val
val
Definition: Pythia8_RapidityOrderMPI.py:14
ActsTrk::SeedingTool::m_numPhiNeighbors
Gaudi::Property< int > m_numPhiNeighbors
Definition: SeedingTool.h:285
ActsTrk::SeedingTool::m_seedConfForwardNTopLargeR
Gaudi::Property< size_t > m_seedConfForwardNTopLargeR
Definition: SeedingTool.h:222
ActsTrk::SeedingTool::m_seedConfForwardZMin
Gaudi::Property< float > m_seedConfForwardZMin
Definition: SeedingTool.h:213
ActsTrk::SeedingTool::m_seedConfirmationInFilter
Gaudi::Property< bool > m_seedConfirmationInFilter
Definition: SeedingTool.h:261
ActsTrk::SeedingTool::m_rBinEdges
Gaudi::Property< std::vector< float > > m_rBinEdges
Definition: SeedingTool.h:131
ActsTrk::SeedingTool::m_interactionPointCut
Gaudi::Property< bool > m_interactionPointCut
Definition: SeedingTool.h:172
ActsTrk::SeedingTool::m_rRangeMiddleSP
Gaudi::Property< std::vector< std::vector< double > > > m_rRangeMiddleSP
Definition: SeedingTool.h:180
ActsTrk::SeedingTool::m_useDeltaRorTopRadius
Gaudi::Property< bool > m_useDeltaRorTopRadius
Definition: SeedingTool.h:267
SeedingTool.h
ActsTrk
The AlignStoreProviderAlg loads the rigid alignment corrections and pipes them through the readout ge...
Definition: MSTrackingVolumeBuilder.cxx:25
top
@ top
Definition: TruthClasses.h:64
ActsTrk::SeedingTool::m_seedConfForwardNTopSmallR
Gaudi::Property< size_t > m_seedConfForwardNTopSmallR
Definition: SeedingTool.h:225
ActsTrk::SeedingTool::m_gridCfg
Acts::CylindricalSpacePointGridConfig m_gridCfg
Definition: SeedingTool.h:109
ActsTrk::SeedingTool::m_zOriginWeightFactor
Gaudi::Property< float > m_zOriginWeightFactor
Definition: SeedingTool.h:252
ActsTrk::SeedingTool::m_useDetailedDoubleMeasurementInfo
Gaudi::Property< bool > m_useDetailedDoubleMeasurementInfo
Definition: SeedingTool.h:237
ActsTrk::SeedingTool::m_zBinNeighborsBottom
Gaudi::Property< std::vector< std::pair< int, int > > > m_zBinNeighborsBottom
Definition: SeedingTool.h:278
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:135
ActsTrk::SeedingTool::m_sigmaError
Gaudi::Property< float > m_sigmaError
Definition: SeedingTool.h:247
ActsTrk::SeedingTool::m_logger
std::unique_ptr< const Acts::Logger > m_logger
logging instance
Definition: SeedingTool.h:301
ActsTrk::SeedingTool::m_pixelId
const PixelID * m_pixelId
Definition: SeedingTool.h:105
ActsTrk::SeedingTool::m_maxPtScattering
Gaudi::Property< float > m_maxPtScattering
Definition: SeedingTool.h:166
ActsTrk::SeedingTool::m_numSeedIncrement
Gaudi::Property< float > m_numSeedIncrement
Definition: SeedingTool.h:259
ActsTrk::SeedingTool::m_seedConfForwardMinBottomRadius
Gaudi::Property< float > m_seedConfForwardMinBottomRadius
Definition: SeedingTool.h:228
ActsTrk::SeedingTool::m_zMin
Gaudi::Property< float > m_zMin
Definition: SeedingTool.h:120
Generated on Mon Sep 29 2025 21:18:28 for ATLAS Offline Software by   doxygen 1.8.18