Skip to content
Snippets Groups Projects
Commit 821efdb6 authored by Jihee Kim's avatar Jihee Kim
Browse files

reduced min kinetic E to record more mcparticles and mcparticles2 coll

parent 5d98d8e0
Branches
Tags
1 merge request!149Resolve "Add true decaying particle info onto event display"
...@@ -6,6 +6,7 @@ from GaudiKernel.SystemOfUnits import mm, MeV, rad, ns ...@@ -6,6 +6,7 @@ from GaudiKernel.SystemOfUnits import mm, MeV, rad, ns
from Configurables import ApplicationMgr, EICDataSvc, PodioOutput, GeoSvc from Configurables import ApplicationMgr, EICDataSvc, PodioOutput, GeoSvc
from Configurables import PodioInput 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__Digi__CalorimeterHitDigi as CalHitDigi
from Configurables import Jug__Reco__CalorimeterHitReco as CalHitReco from Configurables import Jug__Reco__CalorimeterHitReco as CalHitReco
from Configurables import Jug__Reco__CalorimeterHitsEtaPhiProjector as CalHitsProj from Configurables import Jug__Reco__CalorimeterHitsEtaPhiProjector as CalHitsProj
...@@ -37,6 +38,12 @@ podioevent = EICDataSvc('EventDataSvc', inputs=kwargs['input'].split(','), Outpu ...@@ -37,6 +38,12 @@ podioevent = EICDataSvc('EventDataSvc', inputs=kwargs['input'].split(','), Outpu
podin = PodioInput('PodioReader', collections=['mcparticles', 'EcalBarrelHits'], OutputLevel=DEBUG) podin = PodioInput('PodioReader', collections=['mcparticles', 'EcalBarrelHits'], OutputLevel=DEBUG)
podout = PodioOutput('out', filename=kwargs['output']) podout = PodioOutput('out', filename=kwargs['output'])
mccopier = MCCopier('MCCopier',
# OutputLevel=DEBUG,
inputCollection='mcparticles',
outputCollection='mcparticles2')
# use the same daq_setting for digi/reco pair # use the same daq_setting for digi/reco pair
imcal_barrel_daq = dict( imcal_barrel_daq = dict(
dynamicRangeADC=3.*MeV, dynamicRangeADC=3.*MeV,
...@@ -86,7 +93,7 @@ imcal_barrel_clreco = RecoCoG('imcal_barrel_clreco', ...@@ -86,7 +93,7 @@ imcal_barrel_clreco = RecoCoG('imcal_barrel_clreco',
podout.outputCommands = ['keep *'] podout.outputCommands = ['keep *']
ApplicationMgr( ApplicationMgr(
TopAlg=[podin, TopAlg=[podin, mccopier,
imcal_barrel_digi, imcal_barrel_reco, imcal_barrel_digi, imcal_barrel_reco,
imcal_barrel_merger, imcal_barrel_cl, imcal_barrel_clreco, imcal_barrel_merger, imcal_barrel_cl, imcal_barrel_clreco,
podout], podout],
......
...@@ -87,7 +87,7 @@ ls -lh ${CB_EMCAL_GEN_FILE} ...@@ -87,7 +87,7 @@ ls -lh ${CB_EMCAL_GEN_FILE}
# Run geant4 simulations # Run geant4 simulations
npsim --runType batch \ npsim --runType batch \
-v WARNING \ -v WARNING \
--part.minimalKineticEnergy "1*TeV" \ --part.minimalKineticEnergy "0.5*MeV" \
--numberOfEvents ${CB_EMCAL_NUMEV} \ --numberOfEvents ${CB_EMCAL_NUMEV} \
--compactFile ${CB_EMCAL_COMPACT_PATH} \ --compactFile ${CB_EMCAL_COMPACT_PATH} \
--inputFiles ${CB_EMCAL_GEN_FILE} \ --inputFiles ${CB_EMCAL_GEN_FILE} \
......
...@@ -83,7 +83,7 @@ ls -lh ${CB_EMCAL_GEN_FILE} ...@@ -83,7 +83,7 @@ ls -lh ${CB_EMCAL_GEN_FILE}
# Run geant4 simulations # Run geant4 simulations
npsim --runType batch \ npsim --runType batch \
-v WARNING \ -v WARNING \
--part.minimalKineticEnergy "1*TeV" \ --part.minimalKineticEnergy "0.5*MeV" \
--numberOfEvents ${CB_EMCAL_NUMEV} \ --numberOfEvents ${CB_EMCAL_NUMEV} \
--compactFile ${CB_EMCAL_COMPACT_PATH} \ --compactFile ${CB_EMCAL_COMPACT_PATH} \
--inputFiles ${CB_EMCAL_GEN_FILE} \ --inputFiles ${CB_EMCAL_GEN_FILE} \
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment