diff --git a/ecal/emcal_electrons.sh b/ecal/emcal_electrons.sh
index eefa87ad8696f7cf2654e61362288942c8c91132..3173c35d8e1a779ad80be9eb0aaec0d11a4df289 100644
--- a/ecal/emcal_electrons.sh
+++ b/ecal/emcal_electrons.sh
@@ -5,7 +5,7 @@ if [[ ! -n  "${JUGGLER_DETECTOR}" ]] ; then
 fi
 
 if [[ ! -n  "${JUGGLER_N_EVENTS}" ]] ; then 
-  export JUGGLER_N_EVENTS=1000
+  export JUGGLER_N_EVENTS=100
 fi
 
 export JUGGLER_FILE_NAME_TAG="emcal_uniform_electrons"
diff --git a/ecal/scripts/emcal_electrons.cxx b/ecal/scripts/emcal_electrons.cxx
index ca45a209de64a25df72c5b0c628d3e2acb936090..02ac3fe0ec441971a5a48add99858e4f45b2a90f 100644
--- a/ecal/scripts/emcal_electrons.cxx
+++ b/ecal/scripts/emcal_electrons.cxx
@@ -19,7 +19,7 @@
 
 using namespace HepMC3;
 
-void emcal_electrons(int n_events = 1e6, const char* out_fname = "./data/emcal_electron_0GeVto30GeV_100kEvt.hepmc")
+void emcal_electrons(int n_events = 1e3, const char* out_fname = "./data/emcal_electron_0GeVto30GeV_100kEvt.hepmc")
 {
   WriterAscii hepmc_output(out_fname);
   int events_parsed = 0;
@@ -48,6 +48,8 @@ void emcal_electrons(int n_events = 1e6, const char* out_fname = "./data/emcal_e
     // sphere of given radius, in this case momentum.
     r1->Sphere(px, py, pz, p);
 
+    //std::cout << std::sqrt(px*px + py*py + pz*pz) - p << " is zero? \n";
+
     // type 1 is final state
     // pdgid 11 - electron 0.510 MeV/c^2
     GenParticlePtr p3 = std::make_shared<GenParticle>(