From 3dbb114a30db8a82a30e310cc70b0af8c0a32412 Mon Sep 17 00:00:00 2001
From: Ziyue Zhang <Ziyue_Zhang@localhost.localdomain>
Date: Thu, 4 Mar 2021 18:15:20 -0600
Subject: [PATCH] WIP: calc_inv_quant_both

---
 benchmarks/dvmp/analysis/dvmp.h       | 4 ++--
 benchmarks/dvmp/analysis/vm_invar.cxx | 4 ++--
 2 files changed, 4 insertions(+), 4 deletions(-)

diff --git a/benchmarks/dvmp/analysis/dvmp.h b/benchmarks/dvmp/analysis/dvmp.h
index edc6338d..d64f68f5 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 f4671866..1c6d6df9 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
-- 
GitLab