diff --git a/benchmarks/dvmp/analysis/dvmp.h b/benchmarks/dvmp/analysis/dvmp.h index 4d95048c1fdc8ba0d8c57b3d5d48a041ed413324..20afc7aa45b1952707ecd3b36b21bb4364e003f9 100644 --- a/benchmarks/dvmp/analysis/dvmp.h +++ b/benchmarks/dvmp/analysis/dvmp.h @@ -58,7 +58,7 @@ namespace util { } //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}; //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 @@ -79,6 +79,7 @@ namespace util { //search for di-lepton pair for the decay in recoil int daughter_pid = -1; //unsigned + cout<<"_daughter_ = "<<daughter<<endl; if(daughter == "electron"){ daughter_pid = 11; }else if(daughter == "muon"){