WIP: rec_ordered(rec) and sim_ordered(sim) debug

benchmarks/dvmp/analysis/dvmp.h

@@ -54,6 +54,7 @@ namespace util {
       double e = sqrt(px*px + py*py + pz*pz + mass*mass);
       momenta[i].SetPxPyPzE(px, py, pz, e);
     }
+    cout<<Form("sim_ordered, P0 = (%f, %f, %f, %f), P1 = (%f, %f, %f, %f)", momenta[1].px(), momenta[1].py(), momenta[1].pz(), momenta[1].mass(), momenta[3].px(), momenta[3].py(), momenta[3].pz(), momenta[3].mass())<<endl;
     return momenta;
   }
   
@@ -77,6 +78,8 @@ namespace util {
       momenta[3].SetPxPyPzE(parts[i].p.x,  parts[i].p.y,  parts[i].p.z,  parts[i].energy);
     }
     
+    cout<<Form("rec_ordered, P0 = (%f, %f, %f, %f), P1 = (%f, %f, %f, %f)", momenta[1].px(), momenta[1].py(), momenta[1].pz(), momenta[1].mass(), momenta[3].px(), momenta[3].py(), momenta[3].pz(), momenta[3].mass())<<endl;
+    
     
     //search for di-lepton pair for the decay in recoil
     int daughter_pid = -1;    //unsigned
@@ -105,8 +108,6 @@ namespace util {
         
       }
     }
-    //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);
       momenta[6].SetPxPyPzE(parts[second].p.x, parts[second].p.y, parts[second].p.z, parts[second].energy);