Skip to content
Snippets Groups Projects

Replacing nu with y

Merged Ziyue Zhang requested to merge ziyue_work_branch into master
1 file
+ 2
1
Compare changes
  • Side-by-side
  • Inline
@@ -58,7 +58,7 @@ namespace util {
@@ -58,7 +58,7 @@ namespace util {
}
}
//import Reconstructed particles and set the wanted particles in the intended order========================
//import Reconstructed particles and set the wanted particles in the intended order========================
inline auto momenta_ordered_rec(const std::vector<eic::ReconstructedParticleData>& parts, std::string_view mother, std::string_view daughter){//mother and daughter are not used yet; will be useful when generater is different and/or when the mcparticles doesn't follow the same order in all events
inline auto momenta_ordered_rec(const std::vector<eic::ReconstructedParticleData>& parts, std::string_view mother, std::string_view daughter){
std::vector<ROOT::Math::PxPyPzMVector> momenta{7};
std::vector<ROOT::Math::PxPyPzMVector> momenta{7};
//0:e0 1:p0 2:e1 3:p1 4:recoil system (without p1) 5:l1 from recoil decay 6:l2 from recoil decay
//0:e0 1:p0 2:e1 3:p1 4:recoil system (without p1) 5:l1 from recoil decay 6:l2 from recoil decay
for(int i = 0 ; i < 7 ; i++) momenta[i].SetPxPyPzE(0., 0., 0., 0.); //initialize as all 0
for(int i = 0 ; i < 7 ; i++) momenta[i].SetPxPyPzE(0., 0., 0., 0.); //initialize as all 0
@@ -79,6 +79,7 @@ namespace util {
@@ -79,6 +79,7 @@ namespace util {
//search for di-lepton pair for the decay in recoil
//search for di-lepton pair for the decay in recoil
int daughter_pid = -1; //unsigned
int daughter_pid = -1; //unsigned
 
cout<<"_daughter_ = "<<daughter<<endl;
if(daughter == "electron"){
if(daughter == "electron"){
daughter_pid = 11;
daughter_pid = 11;
}else if(daughter == "muon"){
}else if(daughter == "muon"){
Loading