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

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

parent dd07ceac
No related branches found
No related tags found
1 merge request!39Replacing nu with y
This commit is part of merge request !39. Comments created here will be created in the context of that merge request.
...@@ -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"){
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment