Add index leaf to allow navigation to sub-ranked genealogy sim info. …

…Allow matching multiple sim particles by default
Mu2e · Nov 21, 2023 · 38c6bfa · 38c6bfa
1 parent cf3b1b3
commit 38c6bfa
Show file tree

Hide file tree

Showing 3 changed files with 7 additions and 2 deletions.
diff --git a/fcl/prolog.fcl b/fcl/prolog.fcl
@@ -97,13 +97,15 @@ dioLLWeight : {
 }
 genCountLogger: { module_type: GenEventCountReader }
 
-AllOpt : { fillMC : true
+AllOpt : {
+  fillMC : true
   trkqual : "TrkQual"
   fillTrkQual : true
   trkpid : "TrkPID"
   fillTrkPID : true
   fillHits : true
   genealogyDepth : -1
+  matchDepth : -1
 }
 
 DeM : { input : "KK"

diff --git a/inc/SimInfo.hh b/inc/SimInfo.hh
@@ -20,6 +20,7 @@ namespace mu2e
     Int_t proc = -1; // process code
     Int_t gen = -1; // generator code
     Float_t time = -1.0;  // Origin time of the SimParticle
+    Int_t index; // index into the SimInfo vector
     XYZVectorF mom = XYZVectorF();  // origin momentumof the SimParticle
     XYZVectorF pos = XYZVectorF();  // origin position of the SimParticle
     MCRelationship prirel = MCRelationship::none; // relationship to the event primary particles

diff --git a/src/InfoMCStructHelper.cc b/src/InfoMCStructHelper.cc
@@ -160,7 +160,8 @@ namespace mu2e {
           sim_info.nhits = kseedmc.simParticle(imatch)._nhits;
           sim_info.nactive = kseedmc.simParticle(imatch)._nactive;
         }
-
+        // record the index this object will have
+        sim_info.index = siminfos.size();
         siminfos.push_back(sim_info);
         if (current_sim_particle.parent().isNonnull()) {
           current_sim_particle_ptr = current_sim_particle.parent();
@@ -189,6 +190,7 @@ namespace mu2e {
         auto trkprimaryptr = kseedmc.simParticle().simParticle(_spcH);
         sim_info.trkrel = MCRelationship(spp, trkprimaryptr);
         sim_info.prirel = MCRelationship(spp, spp);
+        sim_info.index = siminfos.size();
         siminfos.push_back(sim_info);
       }
     }