From 5e9986b0c4a0e3b86dff8b14f72c5fab0178527c Mon Sep 17 00:00:00 2001
From: Ziyue Zhang <Ziyue_Zhang@localhost.localdomain>
Date: Fri, 5 Mar 2021 12:08:14 -0600
Subject: [PATCH] WIP: rec_ordered(rec) and sim_ordered(sim) debug
---
benchmarks/dvmp/analysis/dvmp.h | 4 ++--
1 file changed, 2 insertions(+), 2 deletions(-)
diff --git a/benchmarks/dvmp/analysis/dvmp.h b/benchmarks/dvmp/analysis/dvmp.h
index ef1d5455..224d5b5b 100644
--- a/benchmarks/dvmp/analysis/dvmp.h
+++ b/benchmarks/dvmp/analysis/dvmp.h
@@ -105,7 +105,7 @@ namespace util {
}
}
- cout<<Form("first = %d, second = %d, jpsi new mass = %f", first, second, best_mass)<<endl;
+ //cout<<Form("first = %d, second = %d, jpsi new mass = %f", first, second, best_mass)<<endl;
if(first == -1) cout<<"jpsi not RCed"<<endl;
if(first != -1){
momenta[5].SetPxPyPzE(parts[first].p.x, parts[first].p.y, parts[first].p.z, parts[first].energy);
@@ -116,7 +116,7 @@ namespace util {
//float dp = 10.;
for(int i = 0 ; i < parts.size(); i++){
if(i==first || i==second) continue; //skip the paired leptons
- if(parts[i].pid != -11) continue;
+ if(parts[i].pid != 11) continue;
//float ptmp = 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);
//if( (k_prime.px()) * (pair_4p.px()) + (k_prime.py()) * (pair_4p.py()) + (k_prime.pz()) * (pair_4p.pz()) > 0. || ptmp >= 10.) continue; //angle between jpsi and scattered electron < pi/2, 3-momentum mag < 10.
//if(dp > 10.- ptmp){ //if there are more than one candidate of scattered electron, choose the one with highest 3-momentum mag
--
GitLab