compact top

Core/include/Acts/Seeding/SeedFinder.ipp

@@ -123,6 +123,9 @@ void SeedFinder<external_spacepoint_t, platform_t>::createSeedsForGroup(
       seedFilterState.nTopSeedConf = rM > seedConfRange.rMaxSeedConf
                                          ? seedConfRange.nTopForLargeR
                                          : seedConfRange.nTopForSmallR;
+      // set max bottom radius for seed confirmation
+      seedFilterState.rMaxSeedConf = seedConfRange.rMaxSeedConf;
+      // continue if number of top SPs is smaller than minimum
       if (state.compatTopSP.size() < seedFilterState.nTopSeedConf) {
         continue;
       }
@@ -319,6 +322,18 @@ void SeedFinder<external_spacepoint_t, platform_t>::filterCandidates(
       rotationTermsUVtoXY[1] = spM.y() * sinTheta / spM.radius();
     }
 
+    // minimum number of compatible top SPs to trigger the filter for a certain
+    // middle bottom pair if seedConfirmation is false we always ask for at
+    // least one compatible top to trigger the filter
+    size_t minCompatibleTopSPs = 2;
+    if (!m_config.seedConfirmation or
+        state.compatBottomSP[b]->radius() > seedFilterState.rMaxSeedConf) {
+      minCompatibleTopSPs = 1;
+    }
+    if (m_config.seedConfirmation and seedFilterState.numQualitySeeds) {
+      minCompatibleTopSPs++;
+    }
+
     for (size_t index_t = t0; index_t < numTopSP; index_t++) {
       const std::size_t& t = sorted_tops[index_t];
 
@@ -538,7 +553,8 @@ void SeedFinder<external_spacepoint_t, platform_t>::filterCandidates(
       state.ptVec.push_back(pT);
     }  // loop on tops
 
-    if (state.topSpVec.empty()) {
+    // continue if number of top SPs is smaller than minimum required for filter
+    if (state.topSpVec.size() < minCompatibleTopSPs) {
       continue;
     }
 

Core/include/Acts/Seeding/SeedFinderConfig.hpp

@@ -71,6 +71,11 @@ struct SeedFinderConfig {
   // difference in tan(theta)
   bool arithmeticAverageCotTheta = false;
 
+  // radius of bottom component of seed that is used to define the number of
+  // compatible top required
+  float rMaxSeedConf =
+      std::numeric_limits<float>::max();  // Acts::UnitConstants::mm
+
   // non equidistant binning in z
   std::vector<float> zBinEdges;