diff --git a/benchmarks/dvmp/analysis/dvmp.h b/benchmarks/dvmp/analysis/dvmp.h
index b3375ea609a5845b49f1672fdf9127c662f04208..eb202920966863b4f60185e397bfd90b678c2129 100644
--- a/benchmarks/dvmp/analysis/dvmp.h
+++ b/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;
   }