Skip to content
Snippets Groups Projects
Commit 9a1b18eb authored by Christopher Dilks's avatar Christopher Dilks
Browse files

refactor: cleanup PDG access

parent 5a6be092
No related branches found
No related tags found
2 merge requests!309Irt algo,!293feat: dRICH benchmarks
...@@ -43,13 +43,13 @@ namespace benchmarks { ...@@ -43,13 +43,13 @@ namespace benchmarks {
if(!simpart.isAvailable()) { m_log->warn("simulated particle not available for this association"); continue; } if(!simpart.isAvailable()) { m_log->warn("simulated particle not available for this association"); continue; }
// get PDG values // get PDG values
auto simpart_pdg = simpart.getPDG(); // from MC truth auto simpart_pdg = simpart.getPDG(); // from MC truth
auto recpart_pdg = recpart.getPDG(); // PDG member of ReconstructedParticle auto recpart_pdg = recpart.getPDG(); // PDG member of ReconstructedParticle
auto pid_pdg = recpart_pdg; // ReconstructedParticle::particleIDUsed::PDG, if available (otherwise use `recpart_pdg`) auto recpart_pid_used = recpart.getParticleIDUsed();
auto recpart_pid = recpart.getParticleIDUsed(); auto pid_pdg = recpart_pid_used.isAvailable() ? recpart_pid_used.getPDG() : 0;
if(recpart_pid.isAvailable()) pid_pdg = recpart_pid.getPDG(); if(!recpart_pid_used.isAvailable())
else m_log->warn("reconstructed particle does not have particleIDUsed relation; using PDG member instead"); m_log->error("reconstructed particle does not have particleIDUsed relation");
m_log->trace("PDGs: sim | rec | pid = {:^6} | {:^6} | {:^6}", simpart_pdg, recpart_pdg, pid_pdg); m_log->trace("PDGs: sim.PDG | rec.PDG | PDG from PID = {:^6} | {:^6} | {:^6}", simpart_pdg, recpart_pdg, pid_pdg);
// check goodness of PID: if near zero, assume PID was not used for this particle and skip // check goodness of PID: if near zero, assume PID was not used for this particle and skip
if(recpart.getGoodnessOfPID()<0.01) { if(recpart.getGoodnessOfPID()<0.01) {
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment