debugging

benchmarks/imaging_ecal/scripts/draw_cluster.py

@@ -149,11 +149,11 @@ if __name__ == '__main__':
 
     # Read all mc particles
     dfallmcp = get_all_mcp(args.file, args.iev, 'mcparticles2')
+    pdgbase = ROOT.TDatabasePDG()
     # Select decaying particles
     dftemp = dfallmcp[dfallmcp['g4Parent'] == 1.0]
     if len(dftemp) > 0:
         dfdecaymcp = dftemp.copy()
-        pdgbase = ROOT.TDatabasePDG()
         for iptl in [0, len(dfdecaymcp) - 1]:
             infoptl = pdgbase.GetParticle(int(dfdecaymcp['pid'].iloc[iptl]))
             print("{} Decaying particle = {}, pdgcode = {}, charge = {}, mass = {}"\