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

WIP: calc_inv_quant_both (utilize)

parent 8eb153f1
No related branches found
No related tags found
1 merge request!39Replacing nu with y
......@@ -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
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment