Replacing nu with y

benchmarks/dvmp/analysis/vm_invar.cxx

@@ -85,9 +85,10 @@ int vm_invar(const std::string& config_name)
   // Define analysis flow
   auto d_im = d.Define("p_rec", util::momenta_RC, {"DummyReconstructedParticles"})
                   .Define("p_sim_ordered", momenta_ordered_sim, {"mcparticles2"})
+                  .Define("invariant_quantities_sim", util::calc_inv_quant_both, {"p_sim_ordered"})
                   //.Define("p_rec_ordered", momenta_ordered_rec, {"DummyReconstructedParticles"})
                   .Define("N", "p_rec.size()")
-                  .Define("p_sim", util::momenta_from_simulation, {"mcparticles2"})
+                  //.Define("p_sim", util::momenta_from_simulation, {"mcparticles2"})
                   //================================================================
                   .Define("invariant_quantities_rec", calc_inv_quant_rec, {"p_rec"})
                   //.Define("invariant_quantities_sim", util::calc_inv_quant_sim, {"p_sim"})
@@ -98,9 +99,9 @@ int vm_invar(const std::string& config_name)
                   .Define("nu_sim", util::get_nu, {"invariant_quantities_sim"})
                   .Define("Q2_sim", util::get_Q2, {"invariant_quantities_sim"})
                   .Define("x_sim", util::get_x, {"invariant_quantities_sim"})
-                  .Define("t_sim", util::get_t, {"invariant_quantities_sim"})
+                  .Define("t_sim", util::get_t, {"invariant_quantities_sim"});
                   //.Define("invariant_quantities_rec", util::calc_inv_quant_both, {"p_rec_ordered"})
-                  .Define("invariant_quantities_sim", util::calc_inv_quant_both, {"p_sim_ordered"});
+                  
   //================================================================
 
   // Define output histograms