Skip to content
Snippets Groups Projects
Commit ed17fb46 authored by Ziyue Zhang's avatar Ziyue Zhang
Browse files

WIP: Test on PID instead of pmass

parent edf6f6c0
No related branches found
No related tags found
1 merge request!39Replacing nu with y
......@@ -88,9 +88,9 @@ namespace util {
if( parts[i].pid!=daughter_pid) continue;
for (int j = i + 1 ; j < parts.size() ; ++j) {
if( parts[j].pid!= -daughter_pid) continue;
ROOT::Math::PxPyPzMVector lpt_1();
ROOT::Math::PxPyPzMVector lpt_1;
lpt_1.SetPxPyPzE(parts[i].p.x, parts[i].p.y, parts[i].p.z, parts[i].energy);
ROOT::Math::PxPyPzMVector lpt_2();
ROOT::Math::PxPyPzMVector lpt_2;
lpt_2.SetPxPyPzE(parts[j].p.x, parts[j].p.y, parts[j].p.z, parts[j].energy);
const double new_mass{(lpt_1 + lpt_2).mass()};
if (fabs(new_mass - get_pdg_mass(mother)) < fabs(best_mass - get_pdg_mass(mother))) {
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment