Skip to content
Snippets Groups Projects
Commit 48c4ceb1 authored by Jihee Kim's avatar Jihee Kim
Browse files

debugging

parent c4c225cd
No related branches found
No related tags found
1 merge request!149Resolve "Add true decaying particle info onto event display"
This commit is part of merge request !149. Comments created here will be created in the context of that merge request.
Hide whitespace changes
Inline Side-by-side
benchmarks/imaging_ecal/scripts/draw_cluster.py
+ 1
− 1
View file @ 48c4ceb1
...@@ -149,11 +149,11 @@ if __name__ == '__main__': ...@@ -149,11 +149,11 @@ if __name__ == '__main__':
# Read all mc particles # Read all mc particles
dfallmcp = get_all_mcp(args.file, args.iev, 'mcparticles2') dfallmcp = get_all_mcp(args.file, args.iev, 'mcparticles2')
pdgbase = ROOT.TDatabasePDG()
# Select decaying particles # Select decaying particles
dftemp = dfallmcp[dfallmcp['g4Parent'] == 1.0] dftemp = dfallmcp[dfallmcp['g4Parent'] == 1.0]
if len(dftemp) > 0: if len(dftemp) > 0:
dfdecaymcp = dftemp.copy() dfdecaymcp = dftemp.copy()
pdgbase = ROOT.TDatabasePDG()
for iptl in [0, len(dfdecaymcp) - 1]: for iptl in [0, len(dfdecaymcp) - 1]:
infoptl = pdgbase.GetParticle(int(dfdecaymcp['pid'].iloc[iptl])) infoptl = pdgbase.GetParticle(int(dfdecaymcp['pid'].iloc[iptl]))
print("{} Decaying particle = {}, pdgcode = {}, charge = {}, mass = {}"\ print("{} Decaying particle = {}, pdgcode = {}, charge = {}, mass = {}"\
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment