Skip to content
Snippets Groups Projects
Commit cb3f30a9 authored by Wouter Deconinck's avatar Wouter Deconinck
Browse files

p_diff to mom_diff

parent 2f5aa2f1
No related branches found
No related tags found
1 merge request!81Re-enable DIS
This commit is part of merge request !81. Comments created here will be created in the context of that merge request.
......@@ -55,7 +55,7 @@ inline auto convertMtoE(const std::vector<ROOT::Math::PxPyPzMVector>& mom)
}
// Momentum
auto momentum = [](std::vector<ROOT::Math::PxPyPzEVector> const& in) {
auto momentum(std::vector<ROOT::Math::PxPyPzEVector> const& in) {
std::vector<double> result;
for (size_t i = 0; i < in.size(); ++i) {
result.push_back(in[i].P());
......@@ -236,20 +236,21 @@ int dis_electrons(const std::string& config_name)
.Define("elec_Q2_sim", "Q2_sim[0]")
.Define("Q2_diff", "(elec_Q2_rec_rand - elec_Q2_sim)/elec_Q2_sim")
.Define("p_diff", sub, {"p_rec","p_sim"})
.Define("mom_rec", momentum, {"p_rec"})
.Define("mom_sim", momentum, {"p_sim"})
.Define("p_diff", sub, {"p_rec_sort","p_sim_sort"})
.Define("mom_rec", momentum, {"p_rec_sort"})
.Define("mom_sim", momentum, {"p_sim_sort"})
.Define("mom_diff", momentum, {"p_diff"})
/*
.Define("electrons_sim", elec_PID_sim, {"sorted_sim"})
.Define("Q2_sim_elec_pid", Q2_from_sim, {"electrons_sim"})
.Define("elec_Q2_sim_pid", "Q2_sim_elec_pid[0]");
.Define("elec_Q2_sim_pid", "Q2_sim_elec_pid[0]")
*/
;
//Momentum
auto h_mom_sim = d0.Histo1D({"h_mom_sim", "; GeV; counts", 100, 0, 50}, "mom_sim");
auto h_mom_rec = d0.Histo1D({"h_mom_rec", "; GeV; counts", 100, 0, 50}, "mom_rec");
auto h_mom_res = d0.Histo1D({"h_mom_res", "; ; counts", 100, -10, 10}, "p_diff");
auto h_mom_res = d0.Histo1D({"h_mom_res", "; ; counts", 100, -10, 10}, "mom_diff");
//Q2
auto h_Q2_sim = d0.Histo1D({"h_Q2_sim", "; GeV; counts", 100, -5, 25}, "elec_Q2_sim");
auto h_Q2_rec = d0.Histo1D({"h_Q2_rec", "; GeV; counts", 100, -5, 25}, "elec_Q2_rec_rand");
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment