diff --git a/benchmarks/dvmp/analysis/dvmp.h b/benchmarks/dvmp/analysis/dvmp.h index edc6338d78caf876de8628e62178af1d77b4adb0..d64f68f5022c60ed73c948be75002f15e3b1da65 100644 --- a/benchmarks/dvmp/analysis/dvmp.h +++ b/benchmarks/dvmp/analysis/dvmp.h @@ -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.; } diff --git a/benchmarks/dvmp/analysis/vm_invar.cxx b/benchmarks/dvmp/analysis/vm_invar.cxx index f46718666ab50ebe6153dc853aaccbe25515b434..1c6d6df933eed667dcabdc84cc3753113961f604 100644 --- a/benchmarks/dvmp/analysis/vm_invar.cxx +++ b/benchmarks/dvmp/analysis/vm_invar.cxx @@ -98,9 +98,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("invariant_quantities_sim", util::calc_inv_quant_both, {"p_sim_ordered"}); //================================================================ // Define output histograms