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,21 +28,21 @@ namespace util {
   // for structure functions
 
   struct inv_quant { // add more when needed
-    double nu, Q2, x, t;
+    double nu, Q2, x;
   };
 
   // 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);  
-    double t = Delta.Dot(Delta);
-    inv_quant quantities = {nu, Q2, Q2/2./P.mass()/nu, t};
+
+    double    nu         = q.Dot(P) / P.mass();
+    double    Q2         = -q.Dot(q);
+    inv_quant quantities = {nu, Q2, Q2 / 2. / P.mass() / nu};
     return quantities;
-  }*/
+  }
+
   
   // for rec
   /*inline inv_quant calc_inv_quant_rec(const std::vector<ROOT::Math::PxPyPzMVector>& parts, const double pdg_mass){
@@ -103,7 +103,7 @@ namespace util {
   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;          }
+  //inline double get_t(inv_quant quantities)  { return quantities.t;          }
 
 
   //=========================================================================================================