From 94aeeb080a2f2670f629830f5c85e93dfa16c94f Mon Sep 17 00:00:00 2001
From: Ziyue Zhang <Ziyue_Zhang@localhost.localdomain>
Date: Fri, 5 Mar 2021 11:37:57 -0600
Subject: [PATCH] WIP: rec_ordered(rec) and sim_ordered(sim) debug
---
benchmarks/dvmp/analysis/dvmp.h | 3 ++-
1 file changed, 2 insertions(+), 1 deletion(-)
diff --git a/benchmarks/dvmp/analysis/dvmp.h b/benchmarks/dvmp/analysis/dvmp.h
index 4d95048c..20afc7aa 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"){
--
GitLab