Add debug statements for specific events

Author: JustusRudolph

Examples/Algorithms/Fatras/src/FatrasSimulation.cpp

@@ -213,6 +213,7 @@ ActsExamples::ProcessCode ActsExamples::FatrasSimulation::execute(
     const AlgorithmContext &ctx) const {
   // read input containers
   const auto &inputParticles = m_inputParticles(ctx);
+  const unsigned nInputParticles = inputParticles.size();
 
   ACTS_DEBUG(inputParticles.size() << " input particles");
 
@@ -285,6 +286,14 @@ ActsExamples::ProcessCode ActsExamples::FatrasSimulation::execute(
 
     particlesSimulated.insert(particleSimulated);
   }
+  const unsigned nOutputParticles = particlesSimulated.size();
+  std::set<unsigned> interesting_events = {16, 374, 429, 813, 927};
+  if ( (std::find(interesting_events.begin(), interesting_events.end(),
+                  ctx.eventNumber) != interesting_events.end() ) ) {
+    ACTS_INFO("In FATRAS: event " << ctx.eventNumber << " simulated "
+              << nOutputParticles << " particles from "
+              << nInputParticles << " input particles.");
+  }
 
   // store ordered output containers
   m_outputParticles(ctx, std::move(particlesSimulated));

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

@@ -212,13 +212,20 @@ ProcessCode ParticleSelector::execute(const AlgorithmContext& ctx) const {
   }
   outputParticles.shrink_to_fit();
 
-  ACTS_DEBUG("event " << ctx.eventNumber << " selected "
-                      << outputParticles.size() << " from "
-                      << inputParticles.size() << " particles");
-  ACTS_DEBUG("filtered out because of charge: " << nInvalidCharge);
-  ACTS_DEBUG("filtered out because of hit count: " << nInvalidHitCount);
-  ACTS_DEBUG("filtered out because of measurement count: "
-             << nInvalidMeasurementCount);
+  std::set<unsigned> interesting_events = {16, 374, 429, 813, 927};
+  if ( (std::find(interesting_events.begin(), interesting_events.end(),
+                  ctx.eventNumber) != interesting_events.end() )
+        && m_cfg.rhoMax != std::numeric_limits<double>::infinity() ) {
+    ACTS_INFO("event " << ctx.eventNumber << " selected "
+                        << outputParticles.size() << " from "
+                        << inputParticles.size() << " particles");
+    ACTS_INFO("filtered out because of charge: " << nInvalidCharge);
+    ACTS_INFO("filtered out because of hit count: " << nInvalidHitCount);
+    ACTS_INFO("filtered out because of measurement count: "
+              << nInvalidMeasurementCount);
+    ACTS_INFO("filtered out because of measurement region count: "
+              << nInvalidMeasurementRegionCount);
+  }
 
   m_outputParticles(ctx, std::move(outputParticles));
 

Examples/Algorithms/TruthTracking/ActsExamples/TruthTracking/TrackTruthMatcher.cpp

@@ -67,7 +67,13 @@ ActsExamples::ProcessCode TrackTruthMatcher::execute(
 
   // For each particle within a track, how many hits did it contribute
   std::vector<ParticleHitCount> particleHitCounts;
-
+  std::set<unsigned> interesting_events = {16, 374, 429, 813, 927};
+  if ( (std::find(interesting_events.begin(), interesting_events.end(),
+                  ctx.eventNumber) != interesting_events.end() ) ) {
+    ACTS_INFO("event " << ctx.eventNumber << " has "
+                        << particles.size() << " particles and "
+                        << tracks.size() << " tracks");
+  }
   for (const auto& track : tracks) {
     // Get the majority truth particle to this track
     identifyContributingParticles(hitParticlesMap, track, particleHitCounts);
@@ -136,11 +142,12 @@ ActsExamples::ProcessCode TrackTruthMatcher::execute(
         ++particleTrackMatch.duplicates;
       }
     } else {  // not matched, i.e. fake track
-      ACTS_DEBUG("Track %u NOT MATCHED to particle %lu with %lu hits in event %lu.\
-                  Out of %lu track hits, %u were right.\n",
-                  track.tipIndex(), majorityParticleId.value(), ctx.eventNumber,
-                  particleTruthHitCount.at(majorityParticleId),
-                  nMajorityHits, track.nMeasurements());
+      ACTS_DEBUG("Track " << track.tipIndex() << " in event " << ctx.eventNumber
+                 << " NOT MATCHED to particle " << majorityParticleId.particle()
+                 << " with " << particleTruthHitCount.at(majorityParticleId) 
+                 << " hits in event " << ctx.eventNumber << ". Out of " 
+                 << track.nMeasurements() << " track hits, " 
+                 << nMajorityHits << " were right.");
       trackParticleMatching[track.index()] = {TrackMatchClassification::Fake,
                                               std::nullopt, particleHitCounts};
 

Examples/Algorithms/TruthTracking/ActsExamples/TruthTracking/TruthTrackFinder.cpp

@@ -70,7 +70,7 @@ ProcessCode TruthTrackFinder::execute(const AlgorithmContext& ctx) const {
   ACTS_VERBOSE("Create prototracks for " << particles.size() << " particles");
   for (const auto& [i, particle] : Acts::enumerate(particles)) {
     str_to_print += std::format("\tParticle {} has {} hits.",
-                                particle.particleId().value(),
+                                particle.particleId().particle(),
                                 particle.numberOfHits());
     // find the corresponding hits for this particle
     const auto& measurements =