Using two vertices

benchmarks/barrel_ecal/scripts/emcal_barrel_pions_electrons.cxx

@@ -76,12 +76,20 @@ void emcal_barrel_pions_electrons(int n_events = 1e6, double e_start = 0.0, doub
     GenParticlePtr p4 = std::make_shared<GenParticle>(FourVector(px, py, pz, sqrt(p * p + (0.000511 * 0.000511))), 11, 1);
 
     GenVertexPtr v1 = std::make_shared<GenVertex>();
+    GenVertexPtr v2 = std::make_shared<GenVertex>();
     v1->add_particle_in(p1);
     v1->add_particle_in(p2);
+    v1->add_particle_out(p3);
 
+    v2->add_particle_in(p1);
+    v2->add_particle_in(p2);
+    v2->add_particle_out(p4);
+/*
     if (phi <= M_PI) {v1->add_particle_out(p3);}
     if (phi > M_PI) {v1->add_particle_out(p4);}
+    */
     evt.add_vertex(v1);
+    evt.add_vertex(v2);
 
     if (events_parsed == 0) {
       std::cout << "First event: " << std::endl;