Replacing nu with y

benchmarks/dvmp/analysis/dvmp.h

@@ -40,7 +40,7 @@ namespace util {
   }
   //========================================import Geant4 and set the wanted particles in the intended order========================
   //0:e0  1:p0    2:e1    3:p1    4:recoil system (without p1)    5:l1 from 4 6:l2 from 4
-  inline auto momenta_ordered_sim(const std::vector<dd4pod::Geant4ParticleData>& parts, std::string_view mother, std::string_view daughter){//mother and daughter are not used yet; will be useful when generater is different and/or when the mcparticles doesn't follow the same order in all events
+  /*inline auto momenta_ordered_sim(const std::vector<dd4pod::Geant4ParticleData>& parts, std::string_view mother, std::string_view daughter){//mother and daughter are not used yet; will be useful when generater is different and/or when the mcparticles doesn't follow the same order in all events
     std::vector<ROOT::Math::PxPyPzMVector> momenta{7};
     int order_map[7] = {0, 3, 2, 6, 5, 7, 8};             
     for(int i = 0 ; i < 7 ; i++){
@@ -52,10 +52,10 @@ namespace util {
       momenta[i].SetPxPyPzE(px, py, pz, e);
     }
     return momenta;
-  }
+  }*/
   
   //import Reconstructed particles and set the wanted particles in the intended order========================
-  inline auto momenta_ordered_rec(const std::vector<eic::ReconstructedParticleData>& parts, std::string_view mother, std::string_view daughter){//mother and daughter are not used yet; will be useful when generater is different and/or when the mcparticles doesn't follow the same order in all events
+  /*inline auto momenta_ordered_rec(const std::vector<eic::ReconstructedParticleData>& parts, std::string_view mother, std::string_view daughter){//mother and daughter are not used yet; will be useful when generater is different and/or when the mcparticles doesn't follow the same order in all events
     std::vector<ROOT::Math::PxPyPzMVector> momenta{7};
     //0:e0  1:p0    2:e1    3:p1    4:recoil system (without p1)    5:l1 from recoil decay  6:l2 from recoil decay
     for(int i = 0 ; i < 7 ; i++) momenta[i].SetPxPyPzE(0., 0., 0., 0.);   //initialize as all 0
@@ -118,7 +118,7 @@ namespace util {
       //}
     }
     return momenta;
-  }
+  }*/
   
  //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){
@@ -186,7 +186,7 @@ namespace util {
     return quantities;
   }
   
-  inline inv_quant calc_inv_quant_both(const std::vector<ROOT::Math::PxPyPzMVector>& parts)
+  /*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]);
@@ -199,7 +199,7 @@ namespace util {
     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.; }