WIP: Add y

benchmarks/dvmp/analysis/dvmp.h

@@ -21,7 +21,7 @@ namespace util {
   // for structure functions
 
   struct inv_quant { // add more when needed
-    double nu, Q2, x, t;
+    double nu, Q2, x, y, t;
   };
 
   // for simu
@@ -221,8 +221,8 @@ namespace util {
     //}else{
         t = Delta.Dot(Delta);
     //}
-    
-    inv_quant quantities = {nu, Q2, Q2/2./P.mass()/nu, t};
+    double y = q.Dot(P) / parts[0].Dot(P);
+    inv_quant quantities = {nu, Q2, Q2/2./P.mass()/nu, y, t};
     return quantities;
   }
   
@@ -230,6 +230,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_y(inv_quant quantities) { return quantities.y; }
   inline double get_t(inv_quant quantities) { return quantities.t; }
 
   // for tracking, add later