10 const std::string&
name,
134 return StatusCode::FAILURE;
140 return StatusCode::FAILURE;
146 return StatusCode::FAILURE;
152 return StatusCode::FAILURE;
161 "Inconsistent config zBinsCustomLooping contains bins that are not "
163 return StatusCode::FAILURE;
172 "Inconsistent config rBinsCustomLooping contains bins that are not "
174 return StatusCode::FAILURE;
194 m_gridCfg.bottomBinFinder = Acts::GridBinFinder<3ul>(
250 m_filterCfg.centralSeedConfirmationRange.nTopForLargeR =
252 m_filterCfg.centralSeedConfirmationRange.nTopForSmallR =
254 m_filterCfg.centralSeedConfirmationRange.seedConfMinBottomRadius =
256 m_filterCfg.centralSeedConfirmationRange.seedConfMaxZOrigin =
258 m_filterCfg.centralSeedConfirmationRange.minImpactSeedConf =
263 m_filterCfg.forwardSeedConfirmationRange.nTopForLargeR =
265 m_filterCfg.forwardSeedConfirmationRange.nTopForSmallR =
267 m_filterCfg.forwardSeedConfirmationRange.seedConfMinBottomRadius =
269 m_filterCfg.forwardSeedConfirmationRange.seedConfMaxZOrigin =
271 m_filterCfg.forwardSeedConfirmationRange.minImpactSeedConf =
278 Acts::Experimental::TripletSeeder(
logger().cloneWithSuffix(
"Finder"));
284 return StatusCode::SUCCESS;
289 float zabs = std::abs(sp->
z());
290 float absCotTheta = zabs /
r;
295 if (zabs > 200 &&
r < 40)
303 static constexpr
float cotThetaEta120 = 1.5095;
304 if (absCotTheta < cotThetaEta120)
308 static constexpr
float cotThetaEta340 = 14.9654;
314 static constexpr
float cotThetaEta220 = 4.4571;
315 if (absCotTheta > cotThetaEta220 &&
r > 260.)
319 static constexpr
float cotThetaEta260 = 6.6947;
320 if (absCotTheta > cotThetaEta260 &&
r > 200.)
324 static constexpr
float cotThetaEta320 = 12.2459;
325 if (absCotTheta > cotThetaEta320 &&
r > 140.)
329 static constexpr
float cotThetaEta400 = 27.2899;
330 if (absCotTheta > cotThetaEta400)
337 const Acts::Experimental::ConstSpacePointProxy2& middle,
338 const Acts::Experimental::ConstSpacePointProxy2&
other,
float cotTheta,
339 bool isBottomCandidate)
const {
345 if (std::abs(middle.zr()[0]) > 1500 and
346 middle.zr()[1] > 100 and middle.zr()[1] < 150) {
355 static constexpr
float cotThetaEta120 = 1.5095;
356 static constexpr
float cotThetaEta360 = 18.2855;
358 float absCotTheta = std::abs(
cotTheta);
360 absCotTheta > cotThetaEta120 && absCotTheta < cotThetaEta360) {
368 const Acts::Experimental::ConstSpacePointProxy2& spM,
369 const Acts::Range1D<float>& rMiddleSpRange)
const {
371 return {rMiddleSpRange.min(), rMiddleSpRange.max()};
374 throw std::runtime_error(
375 "m_rRangeMiddleSP is empty, please check the configuration.");
383 zBin == 0 ? zBin : --zBin;
388 const EventContext& ctx,
389 const std::vector<const xAOD::SpacePointContainer*>& spacePointCollections,
390 const Eigen::Vector3f& beamSpotPos,
float bFieldInZ,
395 gridCfg.bFieldInZ = bFieldInZ;
397 Acts::Experimental::CylindricalSpacePointGrid2 grid(
398 gridCfg,
logger().cloneWithSuffix(
"Grid"));
400 std::size_t totalSpacePoints = 0;
405 std::vector<const xAOD::SpacePoint*> selectedXAODSpacePoints;
406 std::vector<float> selectedSpacePointsR;
407 selectedXAODSpacePoints.reserve(totalSpacePoints);
408 selectedSpacePointsR.reserve(totalSpacePoints);
412 float x =
static_cast<float>(sp->x() - beamSpotPos[0]);
413 float y =
static_cast<float>(sp->y() - beamSpotPos[1]);
414 float z =
static_cast<float>(sp->z());
415 float r = std::hypot(
x,
y);
416 float phi = std::atan2(
y,
x);
422 grid.insert(selectedXAODSpacePoints.size(),
phi,
z,
r);
423 selectedXAODSpacePoints.push_back(sp);
424 selectedSpacePointsR.push_back(
r);
428 for (std::size_t
i = 0;
i < grid.numberOfBins(); ++
i) {
430 grid.at(
i), [&](
const Acts::Experimental::SpacePointIndex2&
a,
431 const Acts::Experimental::SpacePointIndex2&
b) {
432 return selectedSpacePointsR[a] < selectedSpacePointsR[b];
436 Acts::Experimental::SpacePointContainer2 selectedSpacePoints;
437 selectedSpacePoints.createColumns(
438 Acts::Experimental::SpacePointColumns::SourceLinks |
439 Acts::Experimental::SpacePointColumns::XY |
440 Acts::Experimental::SpacePointColumns::ZR |
441 Acts::Experimental::SpacePointColumns::VarianceZ |
442 Acts::Experimental::SpacePointColumns::VarianceR);
444 selectedSpacePoints.createColumns(
445 Acts::Experimental::SpacePointColumns::Strip);
447 selectedSpacePoints.reserve(grid.numberOfSpacePoints());
448 std::vector<Acts::Experimental::SpacePointIndex2> copyFromIndices;
449 copyFromIndices.reserve(grid.numberOfSpacePoints());
450 std::vector<Acts::Experimental::SpacePointIndexRange2> gridSpacePointRanges;
451 gridSpacePointRanges.reserve(grid.numberOfBins());
452 for (std::size_t
i = 0;
i < grid.numberOfBins(); ++
i) {
454 for (
const Acts::Experimental::SpacePointIndex2 spIndex : grid.at(
i)) {
457 auto newSp = selectedSpacePoints.createSpacePoint();
458 newSp.assignSourceLinks(
459 std::array<Acts::SourceLink, 1>{Acts::SourceLink(sp)});
460 newSp.xy() = std::array<float, 2>{
static_cast<float>(sp->
x() - beamSpotPos[0]),
461 static_cast<float>(sp->
y() - beamSpotPos[1])};
462 newSp.zr() = std::array<float, 2>{
static_cast<float>(sp->
z()),
463 selectedSpacePointsR[spIndex]};
464 newSp.varianceZ() =
static_cast<float>(sp->
varianceZ());
465 newSp.varianceR() =
static_cast<float>(sp->
varianceR());
467 Eigen::Vector3f topStripVector =
469 Eigen::Vector3f bottomStripVector =
474 newSp.topStripVector() = std::array<float, 3>{
475 topStripVector.x(), topStripVector.y(), topStripVector.z()};
476 newSp.bottomStripVector() =
477 std::array<float, 3>{bottomStripVector.x(), bottomStripVector.y(),
478 bottomStripVector.z()};
479 newSp.stripCenterDistance() = std::array<float, 3>{
480 stripCenterDistance.x(), stripCenterDistance.y(),
481 stripCenterDistance.z()};
482 newSp.topStripCenter() = std::array<float, 3>{
483 topStripCenter.x(), topStripCenter.y(), topStripCenter.z()};
486 copyFromIndices.push_back(spIndex);
489 gridSpacePointRanges.emplace_back(
begin,
end);
493 selectedXAODSpacePoints = {};
494 selectedSpacePointsR = {};
496 ACTS_VERBOSE(
"Number of space points after selection "
497 << selectedSpacePoints.size() <<
" out of " << totalSpacePoints);
501 const Acts::Range1D<float> rRange = [&]() -> Acts::Range1D<float> {
503 float maxRange = std::numeric_limits<float>::lowest();
504 for (
const Acts::Experimental::SpacePointIndexRange2&
range :
505 gridSpacePointRanges) {
509 auto first = selectedSpacePoints[
range.first];
510 auto last = selectedSpacePoints[
range.second - 1];
512 maxRange =
std::max(last.zr()[1], maxRange);
514 return {minRange, maxRange};
517 auto bottomDoubletFinder = Acts::Experimental::DoubletSeedFinder::create(
518 Acts::Experimental::DoubletSeedFinder::DerivedConfig(
520 auto topDoubletFinder = Acts::Experimental::DoubletSeedFinder::create(
521 Acts::Experimental::DoubletSeedFinder::DerivedConfig(
523 auto tripletFinder = Acts::Experimental::TripletSeedFinder::create(
528 const Acts::Range1D<float> rMiddleSpRange(
533 Acts::Experimental::BroadTripletSeedFilter::Cache filterCache;
534 Acts::Experimental::TripletSeeder::Cache cache;
536 Acts::Experimental::BroadTripletSeedFilter
filter(
539 std::vector<Acts::Experimental::SpacePointContainer2::ConstRange>
541 std::optional<Acts::Experimental::SpacePointContainer2::ConstRange>
543 std::vector<Acts::Experimental::SpacePointContainer2::ConstRange> topSpRanges;
545 Acts::Experimental::SeedContainer2 tmpSeedContainer;
546 tmpSeedContainer.reserve(seedContainer.
capacity());
548 for (
const auto [bottom, middle,
top] : grid.binnedGroup()) {
549 ACTS_VERBOSE(
"Process middle bin " << middle);
550 if (middle >= gridSpacePointRanges.size()) {
551 ATH_MSG_ERROR(
"Grid Binned Group returned an unreasonable middle bin");
552 return StatusCode::FAILURE;
555 bottomSpRanges.clear();
559 bottom, std::back_inserter(bottomSpRanges),
561 -> Acts::Experimental::SpacePointContainer2::ConstRange {
562 return selectedSpacePoints.range(gridSpacePointRanges[
b]).asConst();
565 selectedSpacePoints.range(gridSpacePointRanges[middle]).asConst();
567 top, std::back_inserter(topSpRanges),
569 -> Acts::Experimental::SpacePointContainer2::ConstRange {
570 return selectedSpacePoints.range(gridSpacePointRanges[
t]).asConst();
575 auto firstMiddleSp = middleSpRange->front();
576 auto radiusRangeForMiddle =
579 ACTS_VERBOSE(
"Validity range (radius) for the middle space point is ["
580 << radiusRangeForMiddle.first <<
", "
581 << radiusRangeForMiddle.second <<
"]");
584 cache, *bottomDoubletFinder, *topDoubletFinder, *tripletFinder,
filter,
585 selectedSpacePoints, bottomSpRanges, *middleSpRange, topSpRanges,
586 radiusRangeForMiddle, tmpSeedContainer);
592 auto selectionFunction =
594 const Acts::Experimental::MutableSeedProxy2& seed) ->
bool {
595 float seedQuality = seed.quality();
596 float bottomQuality =
597 filterState.bestSeedQualityMap.at(seed.spacePointIndices()[0]);
598 float middleQuality =
599 filterState.bestSeedQualityMap.at(seed.spacePointIndices()[1]);
601 filterState.bestSeedQualityMap.at(seed.spacePointIndices()[2]);
603 return bottomQuality <= seedQuality || middleQuality <= seedQuality ||
604 topQuality <= seedQuality;
607 seedContainer.
reserve(tmpSeedContainer.size());
610 for (Acts::Experimental::MutableSeedProxy2 seed : tmpSeedContainer) {
616 selectedSpacePoints.at(seed.spacePointIndices()[0])
620 selectedSpacePoints.at(seed.spacePointIndices()[1])
624 selectedSpacePoints.at(seed.spacePointIndices()[2])
628 auto outputSeed = std::make_unique<ActsTrk::Seed>(*bottom, *middle, *
top);
629 outputSeed->setVertexZ(seed.vertexZ());
630 outputSeed->setQuality(seed.quality());
631 seedContainer.
push_back(std::move(outputSeed));
634 return StatusCode::SUCCESS;