Re-enable DIS

benchmarks/dis/analysis/dis_electrons.cxx

@@ -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");