Replacing nu with y

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