Add pseudorapidity cut for truth track seeding to reduce failed tracking attempts

JugTrack/src/components/TrackParamTruthInit.cpp

@@ -52,6 +52,8 @@ namespace Jug::Reco {
     Gaudi::Property<double> m_maxVertexY{this, "maxVertexY", 80. * Gaudi::Units::mm};
     Gaudi::Property<double> m_maxVertexZ{this, "maxVertexZ", 200. * Gaudi::Units::mm};
     Gaudi::Property<double> m_minMomentum{this, "minMomentum", 100. * Gaudi::Units::MeV};
+    Gaudi::Property<double> m_maxEtaForward{this, "maxEtaForward", 4.0};
+    Gaudi::Property<double> m_maxEtaBackward{this, "maxEtaBackward", 4.1};
 
     SmartIF<IParticleSvc> m_pidSvc;
 
@@ -114,6 +116,14 @@ namespace Jug::Reco {
           continue;
         }
 
+        // require minimum pseudorapidity
+        if (part.ps().eta() > m_maxEtaForward || part.ps().eta() < -std::abs(m_maxEtaBackward)) {
+          if (msgLevel(MSG::DEBUG)) {
+            debug() << "ignoring particle with Eta = " << part.ps().eta() << endmsg;
+          }
+          continue;
+        }
+
         // get the particle charge
         // note that we cannot trust the mcparticles charge, as DD4hep
         // sets this value to zero! let's lookup by PDGID instead