Skip to content
Snippets Groups Projects

Replacing nu with y

Merged Ziyue Zhang requested to merge ziyue_work_branch into master
1 file
+ 1
1
Compare changes
  • Side-by-side
  • Inline
@@ -81,7 +81,7 @@ namespace util {
@@ -81,7 +81,7 @@ namespace util {
if(parts[i].pid!=2212) continue;
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
//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);
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;
//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;
}
}
Loading