Skip to content
Snippets Groups Projects
Commit 3dbb114a authored by Ziyue Zhang's avatar Ziyue Zhang
Browse files

WIP: calc_inv_quant_both

parent 173893b5
No related branches found
No related tags found
1 merge request!39Replacing nu with y
...@@ -186,7 +186,7 @@ namespace util { ...@@ -186,7 +186,7 @@ namespace util {
return quantities; 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 //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 q(parts[0] - parts[2]);
...@@ -199,7 +199,7 @@ namespace util { ...@@ -199,7 +199,7 @@ namespace util {
double t = Delta.Dot(Delta); double t = Delta.Dot(Delta);
inv_quant quantities = {nu, Q2, Q2/2./P.mass()/nu, t}; inv_quant quantities = {nu, Q2, Q2/2./P.mass()/nu, t};
return quantities; return quantities;
}*/ }
inline double get_nu(inv_quant quantities) { return quantities.nu / 1000.; } inline double get_nu(inv_quant quantities) { return quantities.nu / 1000.; }
......
...@@ -98,9 +98,9 @@ int vm_invar(const std::string& config_name) ...@@ -98,9 +98,9 @@ int vm_invar(const std::string& config_name)
.Define("nu_sim", util::get_nu, {"invariant_quantities_sim"}) .Define("nu_sim", util::get_nu, {"invariant_quantities_sim"})
.Define("Q2_sim", util::get_Q2, {"invariant_quantities_sim"}) .Define("Q2_sim", util::get_Q2, {"invariant_quantities_sim"})
.Define("x_sim", util::get_x, {"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_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 // Define output histograms
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment