Replacing nu with y

benchmarks/dvmp/analysis/dvmp.h

@@ -120,16 +120,6 @@ namespace util {
     return momenta;
   }
   
-  
-  
-  inline double test_jpsi_pt(const std::vector<ROOT::Math::PxPyPzMVector>& momenta){
-    double px = momenta[4].px();
-    double py = momenta[4].py();
-    double pt = sqrt(px*px + py*py);
-    return pt;
-  }
-  
-  
  //for Dummy rc
   inline inv_quant calc_inv_quant_rec(const std::vector<ROOT::Math::PxPyPzMVector>& parts, const double pdg_mass, const double daughter_mass){
     int first = -1;
@@ -196,6 +186,22 @@ namespace util {
     return quantities;
   }
   
+  inline inv_quant calc_inv_quant_both(const std::vector<ROOT::Math::PxPyPzMVector>& parts)
+  {
+    //0:e0  1:p0    2:e1    3:p1    4:recoil system (without p1)    5:l1 from 4 6:l2 from 4
+    ROOT::Math::PxPyPzMVector q(parts[0] - parts[2]);
+    ROOT::Math::PxPyPzMVector P(parts[1]);
+    ROOT::Math::PxPyPzMVector Delta(parts[3] - parts[1]);//jpsi->l l' gamma, ignore gamma
+
+    //ROOT::Math::PxPyPzMVector Delta(parts[6] - parts[3]);//exact jpsi
+    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};
+    return quantities;
+  }
+  
+  
   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; }