diff --git a/benchmarks/dvmp/analysis/dvmp.h b/benchmarks/dvmp/analysis/dvmp.h index ef1d54551974b09736ebbed98a37b09d403f8373..224d5b5bccaf51401d8d79185096c0335ff09d86 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