From f94e23fe67924d4ae550f2f6cf9361dbd631b571 Mon Sep 17 00:00:00 2001
From: "jihee.kim" <jihee.kim@anl.gov>
Date: Fri, 13 Aug 2021 17:05:01 -0500
Subject: [PATCH] reduced min kinetic E to record more mcparticles and
 mcparticles2 coll

---
 benchmarks/imaging_ecal/options/imaging_2dcluster.py | 9 ++++++++-
 benchmarks/imaging_ecal/run_emcal_barrel.sh          | 2 +-
 benchmarks/imaging_ecal/run_imcal_pion0.sh           | 2 +-
 3 files changed, 10 insertions(+), 3 deletions(-)

diff --git a/benchmarks/imaging_ecal/options/imaging_2dcluster.py b/benchmarks/imaging_ecal/options/imaging_2dcluster.py
index 6802e07d..d2a038ca 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 73bae02c..0dca31e5 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 760de223..9506c117 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} \
-- 
GitLab