Skip to content
Snippets Groups Projects
Commit 358b5dc9 authored by Sylvester Joosten's avatar Sylvester Joosten
Browse files

don't call unmatched charged clusters 'photons'

parent 80c7bd0d
No related branches found
No related tags found
1 merge request!359don't call unmatched charged clusters 'photons'
...@@ -119,14 +119,14 @@ public: ...@@ -119,14 +119,14 @@ public:
<< endmsg; << endmsg;
} }
// get mass/PID from mcparticles, photon in case the matched particle is charged. // get mass/PID from mcparticles, 0 (unidentified) in case the matched particle is charged.
const auto& mc = mcparticles[mcID.value]; const auto& mc = mcparticles[mcID.value];
const double mass = (!mc.charge()) ? mc.mass() : 0; const double mass = (!mc.charge()) ? mc.mass() : 0;
const int32_t pid = (!mc.charge()) ? mc.pdgID() : 22; const int32_t pid = (!mc.charge()) ? mc.pdgID() : 0;
if (msgLevel(MSG::DEBUG)) { if (msgLevel(MSG::DEBUG)) {
if (mc.charge()) { if (mc.charge()) {
debug() << " --> associated mcparticle is not a neutral (PID: " << mc.pdgID() debug() << " --> associated mcparticle is not a neutral (PID: " << mc.pdgID()
<< "), setting the reconstructed particle ID to photon" << endmsg; << "), setting the reconstructed particle ID to 0 (unidentified)" << endmsg;
} }
debug() << " --> found matching associated mcparticle with PID: " << pid << ", energy: " << mc.energy() debug() << " --> found matching associated mcparticle with PID: " << pid << ", energy: " << mc.energy()
<< endmsg; << endmsg;
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment