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

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

parent 28817ecc
No related branches found
No related tags found
1 merge request!39Replacing nu with y
......@@ -54,7 +54,7 @@ namespace util {
double e = sqrt(px*px + py*py + pz*pz + mass*mass);
momenta[i].SetPxPyPzE(px, py, pz, e);
}
cout<<Form("sim_ordered, P0 = (%f, %f, %f, %f), P1 = (%f, %f, %f, %f)", momenta[1].px(), momenta[1].py(), momenta[1].pz(), momenta[1].mass(), momenta[3].px(), momenta[3].py(), momenta[3].pz(), momenta[3].mass())<<endl;
//cout<<Form("sim_ordered, P0 = (%f, %f, %f, %f), P1 = (%f, %f, %f, %f)", momenta[1].px(), momenta[1].py(), momenta[1].pz(), momenta[1].mass(), momenta[3].px(), momenta[3].py(), momenta[3].pz(), momenta[3].mass())<<endl;
return momenta;
}
......@@ -76,9 +76,10 @@ namespace util {
for(int i = 0 ; i < parts.size() ; i++){
if(parts[i].pid!=2212) continue;
momenta[3].SetPxPyPzE(parts[i].p.x, parts[i].p.y, parts[i].p.z, parts[i].energy);
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;
}
cout<<Form("rec_ordered, P0 = (%f, %f, %f, %f), P1 = (%f, %f, %f, %f)", momenta[1].px(), momenta[1].py(), momenta[1].pz(), momenta[1].mass(), momenta[3].px(), momenta[3].py(), momenta[3].pz(), momenta[3].mass())<<endl;
//cout<<Form("rec_ordered, P0 = (%f, %f, %f, %f), P1 = (%f, %f, %f, %f)", momenta[1].px(), momenta[1].py(), momenta[1].pz(), momenta[1].mass(), momenta[3].px(), momenta[3].py(), momenta[3].pz(), momenta[3].mass())<<endl;
//search for di-lepton pair for the decay in recoil
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment