Emcal fix

emcal_electrons.cxx

@@ -40,7 +40,7 @@ void emcal_electrons(int n_events = 1e6, const char* out_fname = "./data/emcal_e
         FourVector(0.0, 0.0, 0.0, 0.938), 2212, 4);
 
     // Define momentum
-    Double_t p = 0.0 + events_parsed * 3.e-4;
+    Double_t p = r1->Uniform(0.0, 30.0);
     Double_t px;
     Double_t py;
     Double_t pz;
@@ -53,7 +53,7 @@ void emcal_electrons(int n_events = 1e6, const char* out_fname = "./data/emcal_e
     GenParticlePtr p3 = std::make_shared<GenParticle>(
         FourVector(
             px, py, pz,
-            sqrt((px * px) + (py * py) + (pz * pz) + (0.000511 * 0.000511))),
+            sqrt(p*p + (0.000511 * 0.000511))),
         11, 1);
 
     GenVertexPtr v1 = std::make_shared<GenVertex>();