WIP: rec_ordered(rec) and sim_ordered(sim) Delta ~ e0 - e1 - l1 - l2

benchmarks/dvmp/analysis/dvmp.h

@@ -54,6 +54,10 @@ namespace util {
       double e = sqrt(px*px + py*py + pz*pz + mass*mass);
       momenta[i].SetPxPyPzE(px, py, pz, e);
     }
+    for(int i = 0 ; i < 7 ; i++){
+        cout<<Form("sim, idx = %d, P(px, py, pz, mass) = (%f, %f, %f, %f)", i, momenta[i].px(), momenta[i].py(), momenta[i].pz(), momenta[i].mass())<<endl;
+    }
+    cout<<"========================================="<<endl;
     return momenta;
   }
   
@@ -124,6 +128,10 @@ namespace util {
       //  dp = 10. - ptmp;
       //}
     }
+    for(int i = 0 ; i < 7 ; i++){
+        cout<<Form("dum, idx = %d, P(px, py, pz, mass) = (%f, %f, %f, %f)", i, momenta[i].px(), momenta[i].py(), momenta[i].pz(), momenta[i].mass())<<endl;
+    }
+    cout<<"========================================="<<endl;
     return momenta;
   }