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

WIP: rec_ordered(rec) and sim_ordered(sim) debug

parent 5eddef48
No related branches found
No related tags found
1 merge request!39Replacing nu with y
......@@ -81,7 +81,7 @@ namespace util {
if(parts[i].pid!=2212) continue;
//px, py, pz, e, mass are not consistent for smeared dummy rc. this is a temp solution, replacing non-smeared energy by the re-calculated energy
double energy_tmp = sqrt(parts[i].p.x*parts[i].p.x + parts[i].p.y*parts[i].p.y + parts[i].p.z*parts[i].p.z + parts[i].mass*parts[i].mass);
momenta[3].SetPxPyPzE(parts[i].p.x, parts[i].p.y, parts[i].p.z, parts[i].energy);
momenta[3].SetPxPyPzE(parts[i].p.x, parts[i].p.y, parts[i].p.z, energy_tmp);
//cout<<Form("rec_ordered, p1(px, py, pz, mass) = (%f, %f, %f, %f)", momenta[3].px(), momenta[3].py(), momenta[3].pz(), momenta[3].mass())<<endl;
}
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment