Improve SimTauCPLinkProducer debug printouts.

Author: rovere

SimCalorimetry/HGCalAssociatorProducers/plugins/SimTauCPLinkProducer.cc

@@ -69,8 +69,9 @@ void SimTauProducer::buildSimTau(SimTauCPLink& t,
   bool is_leaf = (daughters.empty());
 
   if (is_leaf) {
-    LogDebug("SimTauProducer").format(" TO BE SAVED {} ", resonance_idx);
     auto const& gen_particle_barcode = gen_particle_barcodes[gen_particle_key];
+    LogDebug("SimTauProducer")
+        .format(" TO BE SAVED {}, key {}, barcode {}", resonance_idx, gen_particle_key, gen_particle_barcode);
     auto const& found_in_caloparticles = std::find_if(caloPartVec.begin(), caloPartVec.end(), [&](const auto& p) {
       return p.g4Tracks()[0].genpartIndex() == gen_particle_barcode;
     });
@@ -87,13 +88,13 @@ void SimTauProducer::buildSimTau(SimTauCPLink& t,
   } else if (generation != 0) {
     t.resonances.push_back({gen_particle.pdgId(), resonance_idx});
     resonance_idx = t.resonances.size() - 1;
-    LogDebug("SimTauProducer").format(" RESONANCE/INTERMEDIATE {} ", resonance_idx);
+    LogDebug("SimTauProducer").format(" RESONANCE/INTERMEDIATE {}", resonance_idx);
   }
 
   ++generation;
   for (auto daughter = daughters.begin(); daughter != daughters.end(); ++daughter) {
     int gen_particle_key = (*daughter).key();
-    LogDebug("SimTauProducer").format(" gen {} {} {} ", generation, gen_particle_key, (*daughter)->pdgId());
+    LogDebug("SimTauProducer").format(" gen {} {} {}", generation, gen_particle_key, (*daughter)->pdgId());
     buildSimTau(t, generation, resonance_idx, *(*daughter), gen_particle_key, calo_particle_h, gen_particle_barcodes);
   }
 }

SimDataFormats/CaloAnalysis/interface/SimTauCPLink.h

@@ -31,31 +31,38 @@ class SimTauCPLink {
 
   enum decayModes {
     kNull = -1,
-    kOneProng0PiZero,
-    kOneProng1PiZero,
-    kOneProng2PiZero,
-    kOneProng3PiZero,
-    kOneProngNPiZero,
-    kTwoProng0PiZero,
-    kTwoProng1PiZero,
-    kTwoProng2PiZero,
-    kTwoProng3PiZero,
-    kTwoProngNPiZero,
-    kThreeProng0PiZero,
-    kThreeProng1PiZero,
-    kThreeProng2PiZero,
-    kThreeProng3PiZero,
-    kThreeProngNPiZero,
-    kRareDecayMode,
-    kElectron,
-    kMuon
+    kOneProng0PiZero,    // 0
+    kOneProng1PiZero,    // 1
+    kOneProng2PiZero,    // 2
+    kOneProng3PiZero,    // 3
+    kOneProngNPiZero,    // 4
+    kTwoProng0PiZero,    // 5
+    kTwoProng1PiZero,    // 6
+    kTwoProng2PiZero,    // 7
+    kTwoProng3PiZero,    // 8
+    kTwoProngNPiZero,    // 9
+    kThreeProng0PiZero,  // 10
+    kThreeProng1PiZero,  // 11
+    kThreeProng2PiZero,  // 12
+    kThreeProng3PiZero,  // 13
+    kThreeProngNPiZero,  // 14
+    kRareDecayMode,      // 15
+    kElectron,           // 16
+    kMuon                // 17
   };
 
   void dump(void) const {
+    LogDebug("SimTauProducer")
+        .format("Decay mode: {} ", buildDecayModes())
+        .format("Leaves: {} ", leaves.size())
+        .format("Resonances: {}", resonances.size());
     for (auto const &l : leaves) {
       LogDebug("SimTauProducer")
-          .format(
-              "L {} {} CP: {} GenP idx: {}", l.pdgId(), l.resonance_idx(), l.calo_particle_idx(), l.gen_particle_idx());
+          .format("L {} {} CP: {} GenP idx: {}",
+                  l.pdgId(),
+                  l.resonance_idx(),
+                  (int)((l.calo_particle_idx() == -1) ? -1 : calo_particle_leaves[l.calo_particle_idx()].key()),
+                  l.gen_particle_idx());
     }
     for (auto const &r : resonances) {
       LogDebug("SimTauProducer").format("R {} {}", r.first, r.second);
@@ -78,7 +85,7 @@ class SimTauCPLink {
     }
   }
 
-  int buildDecayModes() {
+  int buildDecayModes() const {
     int numElectrons = 0;
     int numMuons = 0;
     int numHadrons = 0;