feat: Add ability to remove secondaries in ParticleSelector (#2265)

Allows to filter secondary particles with `ParticleSelector`

blocked by
- #2266

Examples/Algorithms/TruthTracking/ActsExamples/TruthTracking/ParticleSelector.cpp

@@ -49,6 +49,7 @@ ActsExamples::ParticleSelector::ParticleSelector(const Config& config,
                                       << ")");
   ACTS_DEBUG("remove charged particles " << m_cfg.removeCharged);
   ACTS_DEBUG("remove neutral particles " << m_cfg.removeNeutral);
+  ACTS_DEBUG("remove secondary particles " << m_cfg.removeSecondaries);
 }
 
 ActsExamples::ProcessCode ActsExamples::ParticleSelector::execute(
@@ -65,7 +66,8 @@ ActsExamples::ProcessCode ActsExamples::ParticleSelector::execute(
     const bool validNeutral = (p.charge() == 0) and not m_cfg.removeNeutral;
     const bool validCharged = (p.charge() != 0) and not m_cfg.removeCharged;
     const bool validCharge = validNeutral or validCharged;
-    return validCharge and
+    const bool validSecondary = not m_cfg.removeSecondaries or !p.isSecondary();
+    return validCharge and validSecondary and
            within(p.transverseMomentum(), m_cfg.ptMin, m_cfg.ptMax) and
            within(std::abs(eta), m_cfg.absEtaMin, m_cfg.absEtaMax) and
            within(eta, m_cfg.etaMin, m_cfg.etaMax) and

Examples/Algorithms/TruthTracking/ActsExamples/TruthTracking/ParticleSelector.hpp

@@ -58,6 +58,9 @@ class ParticleSelector final : public IAlgorithm {
     bool removeCharged = false;
     /// Remove neutral particles.
     bool removeNeutral = false;
+
+    /// Remove secondaries.
+    bool removeSecondaries = false;
   };
 
   ParticleSelector(const Config& config, Acts::Logging::Level level);

Examples/Python/python/acts/examples/simulation.py

@@ -42,8 +42,9 @@
         "m",  # (min,max)
         "removeCharged",  # bool
         "removeNeutral",  # bool
+        "removeSecondaries",  # bool
     ],
-    defaults=[(None, None)] * 8 + [None] * 2,
+    defaults=[(None, None)] * 8 + [None] * 3,
 )
 
 
@@ -371,6 +372,7 @@ def addParticleSelection(
                 mMax=config.m[1],
                 removeCharged=config.removeCharged,
                 removeNeutral=config.removeNeutral,
+                removeSecondaries=config.removeSecondaries,
             ),
             level=customLogLevel(),
             inputParticles=inputParticles,

Examples/Python/src/TruthTracking.cpp

@@ -126,6 +126,7 @@ void addTruthTracking(Context& ctx) {
     ACTS_PYTHON_MEMBER(ptMax);
     ACTS_PYTHON_MEMBER(removeCharged);
     ACTS_PYTHON_MEMBER(removeNeutral);
+    ACTS_PYTHON_MEMBER(removeSecondaries);
     ACTS_PYTHON_STRUCT_END();
 
     pythonRangeProperty(c, "rho", &Config::rhoMin, &Config::rhoMax);