Add t in dvmp; add a place holder for InvQuant in REC

WIP: Add t in dvmp; add a place holder for InvQuant with limited track selection in REC

benchmarks/dvmp/analysis/dvmp.h

@@ -28,24 +28,26 @@ namespace util {
   // for structure functions
 
   struct inv_quant { // add more when needed
-    double nu, Q2, x;
+    double nu, Q2, x, t;
   };
 
   // for simu
-  inline inv_quant calc_inv_quant_simu(const std::vector<ROOT::Math::PxPyPzMVector>& parts)
-  {
+  inline inv_quant calc_inv_quant_simu(const std::vector<ROOT::Math::PxPyPzMVector>& parts){
     ROOT::Math::PxPyPzMVector q(parts[0] - parts[2]);
     ROOT::Math::PxPyPzMVector P(parts[3]);
+    ROOT::Math::PxPyPzMVector Delta(parts[6] - parts[3]);
 
     double    nu         = q.Dot(P) / P.mass();
     double    Q2         = -q.Dot(q);
-    inv_quant quantities = {nu, Q2, Q2 / 2. / P.mass() / nu};
+    double    t          = Delta.Dot(Delta);
+    inv_quant quantities = {nu, Q2, Q2 / 2. / P.mass() / nu, t};
     return quantities;
   }
 
-  inline double get_nu_simu(inv_quant quantities) { return quantities.nu / 1000.; }
-  inline double get_Q2_simu(inv_quant quantities) { return quantities.Q2; }
-  inline double get_x_simu(inv_quant quantities) { return quantities.x; }
+  inline double get_nu(inv_quant quantities) { return quantities.nu / 1000.; }
+  inline double get_Q2(inv_quant quantities) { return quantities.Q2; }
+  inline double get_x(inv_quant quantities) { return quantities.x; }
+  inline double get_t(inv_quant quantities) { return quantities.t; }
 
   // for tracking, add later