Resolve "pi0 resolution for ECal barrel"

Closes #19 (closed)

benchmarks/barrel_ecal/scripts/emcal_barrel_pions.cxx

@@ -35,7 +35,7 @@ void emcal_barrel_pions(int n_events = 1e6, double e_start = 0.0, double e_end =
     // FourVector(px,py,pz,e,pdgid,status)
     // type 4 is beam
     // pdgid 11 - electron
-    // pdgid 111 - pi0
+    // pdgid 111 - pi0 
     // pdgid 2212 - proton
     GenParticlePtr p1 = std::make_shared<GenParticle>(FourVector(0.0, 0.0, 10.0, 10.0), 11, 4);
     GenParticlePtr p2 = std::make_shared<GenParticle>(FourVector(0.0, 0.0, 0.0, 0.938), 2212, 4);
@@ -55,13 +55,16 @@ void emcal_barrel_pions(int n_events = 1e6, double e_start = 0.0, double e_end =
 
     // type 1 is final state
     // pdgid 211 - pion+ 139.570 MeV/c^2
-    GenParticlePtr p3 = std::make_shared<GenParticle>(FourVector(px, py, pz, sqrt(p * p + (0.139570 * 0.139570))), 211, 1);
+    // pdgid 111 - pion0 134.977 MeV/c^2
+    //GenParticlePtr p3 = std::make_shared<GenParticle>(FourVector(px, py, pz, sqrt(p * p + (0.139570 * 0.139570))), 211, 1);
+    GenParticlePtr p4 = std::make_shared<GenParticle>(FourVector(px, py, pz, sqrt(p * p + (0.134977 * 0.134977))), 111, 1);
 
     GenVertexPtr v1 = std::make_shared<GenVertex>();
     v1->add_particle_in(p1);
     v1->add_particle_in(p2);
 
-    v1->add_particle_out(p3);
+    //v1->add_particle_out(p3);
+    v1->add_particle_out(p4);
     evt.add_vertex(v1);
 
     if (events_parsed == 0) {