diff --git a/benchmarks/imaging_ecal/options/imaging_2dcluster.py b/benchmarks/imaging_ecal/options/imaging_2dcluster.py
index 6802e07d87218a43614e2ff4e0e6d09997a05cc6..d2a038ca3b92b2a3c2ea2afb23228ed78c029ef8 100644
--- a/benchmarks/imaging_ecal/options/imaging_2dcluster.py
+++ b/benchmarks/imaging_ecal/options/imaging_2dcluster.py
@@ -6,6 +6,7 @@ from GaudiKernel.SystemOfUnits import mm, MeV, rad, ns
 from Configurables import ApplicationMgr, EICDataSvc, PodioOutput, GeoSvc
 
 from Configurables import PodioInput
+from Configurables import Jug__Base__InputCopier_dd4pod__Geant4ParticleCollection_dd4pod__Geant4ParticleCollection_ as MCCopier
 from Configurables import Jug__Digi__CalorimeterHitDigi as CalHitDigi
 from Configurables import Jug__Reco__CalorimeterHitReco as CalHitReco
 from Configurables import Jug__Reco__CalorimeterHitsEtaPhiProjector as CalHitsProj
@@ -37,6 +38,12 @@ podioevent = EICDataSvc('EventDataSvc', inputs=kwargs['input'].split(','), Outpu
 podin = PodioInput('PodioReader', collections=['mcparticles', 'EcalBarrelHits'], OutputLevel=DEBUG)
 podout = PodioOutput('out', filename=kwargs['output'])
 
+
+mccopier = MCCopier('MCCopier',
+        # OutputLevel=DEBUG,
+        inputCollection='mcparticles',
+        outputCollection='mcparticles2')
+
 # use the same daq_setting for digi/reco pair
 imcal_barrel_daq = dict(
         dynamicRangeADC=3.*MeV,
@@ -86,7 +93,7 @@ imcal_barrel_clreco = RecoCoG('imcal_barrel_clreco',
 podout.outputCommands = ['keep *']
 
 ApplicationMgr(
-    TopAlg=[podin,
+    TopAlg=[podin, mccopier,
         imcal_barrel_digi, imcal_barrel_reco,
         imcal_barrel_merger, imcal_barrel_cl, imcal_barrel_clreco,
         podout],
diff --git a/benchmarks/imaging_ecal/run_emcal_barrel.sh b/benchmarks/imaging_ecal/run_emcal_barrel.sh
index 73bae02c594b960b1c55cd9d3b9c308898586b92..0dca31e54a940b5ce3ac2d4e186e02dc7cce9b8e 100644
--- a/benchmarks/imaging_ecal/run_emcal_barrel.sh
+++ b/benchmarks/imaging_ecal/run_emcal_barrel.sh
@@ -87,7 +87,7 @@ ls -lh ${CB_EMCAL_GEN_FILE}
 # Run geant4 simulations
 npsim --runType batch \
       -v WARNING \
-      --part.minimalKineticEnergy "1*TeV" \
+      --part.minimalKineticEnergy "0.5*MeV" \
       --numberOfEvents ${CB_EMCAL_NUMEV} \
       --compactFile ${CB_EMCAL_COMPACT_PATH} \
       --inputFiles ${CB_EMCAL_GEN_FILE} \
diff --git a/benchmarks/imaging_ecal/run_imcal_pion0.sh b/benchmarks/imaging_ecal/run_imcal_pion0.sh
index 760de223b81cb4fd7a994f526f2ffd44361df6ef..9506c117087195eedccf108621dce3918eb237fe 100644
--- a/benchmarks/imaging_ecal/run_imcal_pion0.sh
+++ b/benchmarks/imaging_ecal/run_imcal_pion0.sh
@@ -83,7 +83,7 @@ ls -lh ${CB_EMCAL_GEN_FILE}
 # Run geant4 simulations
 npsim --runType batch \
       -v WARNING \
-      --part.minimalKineticEnergy "1*TeV" \
+      --part.minimalKineticEnergy "0.5*MeV" \
       --numberOfEvents ${CB_EMCAL_NUMEV} \
       --compactFile ${CB_EMCAL_COMPACT_PATH} \
       --inputFiles ${CB_EMCAL_GEN_FILE} \