Skip to content

[Simulation] EMEC simulation AF3 in G4

Rui Zhang requested to merge zhangr/athena:rui_emec_study_24.0.62 into 24.0

Implement AF3 in EMEC (e/y |eta| in 1.5-3.2, customisable) and G4 elsewhere. Study to accelerate full simulation.

It supports to run Geant4 only, AF3 only, and Geant4+AF3 with a threshold of kinematic energy of e.g. < 8GeV (for AF3)

G4 only:

AtlasG4_tf.py --inputEVNTFile $EVNTfile --outputHITSFile /tmp/HITS_$tmp --runNumber 601229 --maxEvents=$nevent --skipEvents 0 --postInclude 'default:PyJobTransforms.UseFrontier' --physicsList="FTFP_BERT_ATL" --conditionsTag "default:OFLCOND-MC23-SDR-RUN3-01" --geometryVersion="default:ATLAS-R3S-2021-03-02-00" --imf False --CA --preInclude 'sim:Campaigns.MC23cSimulationMultipleIoV'

AF3 only (including simplified Geometry):

PATH_TO_SIMPLIFIED_GEO='FastCaloSim/MC23/GeoFiles/v02/simplified_geo.gdml'

AtlasG4_tf.py --inputEVNTFile $EVNTfile --outputHITSFile /tmp/HITS_$tmp --runNumber 601229 --maxEvents=$nevent --skipEvents 0 --postInclude 'default:PyJobTransforms.UseFrontier' --physicsList="FTFP_BERT_ATL" --conditionsTag "default:OFLCOND-MC23-SDR-RUN3-01" --geometryVersion="default:ATLAS-R3S-2021-03-02-00" --imf False --CA --preInclude 'sim:Campaigns.MC23cSimulationMultipleIoV,SimuJobTransforms.FastCaloSim' --preExec "AtlasG4Tf:flags.Sim.SimplifiedGeoPath='${PATH_TO_SIMPLIFIED_GEO}';"

G4+AF3:

PATH_TO_SIMPLIFIED_GEO='FastCaloSim/MC23/GeoFiles/v02/simplified_geo.gdml'

AtlasG4_tf.py --inputEVNTFile $EVNTfile --outputHITSFile /tmp/HITS_$tmp --runNumber 601229 --maxEvents=$nevent --skipEvents 0 --postInclude 'default:PyJobTransforms.UseFrontier' --physicsList="FTFP_BERT_ATL" --conditionsTag "default:OFLCOND-MC23-SDR-RUN3-01" --geometryVersion="default:ATLAS-R3S-2021-03-02-00" --imf False --CA --preInclude 'sim:Campaigns.MC23cSimulationMultipleIoV,SimuJobTransforms.FastCaloSim' --preExec "AtlasG4Tf:flags.Sim.SimplifiedGeoPath='${PATH_TO_SIMPLIFIED_GEO}';flags.Sim.FastCalo.doEMECFCS=True"

flags.Sim.FastCalo.doEMECFCS=True will set ``Ekin < 8192, 1.5 < AbsEta < 3.2` such that e/y particles in |eta| in 1.5-3.2 and Ekin < 8GeV are simulated with AF3.

For R&D, these cuts can be altered via --postExec e.g. (note doEMECFCS=False):

PATH_TO_SIMPLIFIED_GEO='FastCaloSim/MC23/GeoFiles/v02/simplified_geo.gdml'

AtlasG4_tf.py --inputEVNTFile $EVNTfile --outputHITSFile /tmp/HITS_$tmp --runNumber 601229 --maxEvents=$nevent --skipEvents 0 --postInclude 'default:PyJobTransforms.UseFrontier' --physicsList="FTFP_BERT_ATL" --conditionsTag "default:OFLCOND-MC23-SDR-RUN3-01" --geometryVersion="default:ATLAS-R3S-2021-03-02-00" --imf False --CA --preInclude 'sim:Campaigns.MC23cSimulationMultipleIoV,SimuJobTransforms.FastCaloSim' --preExec "AtlasG4Tf:flags.Sim.SimplifiedGeoPath='${PATH_TO_SIMPLIFIED_GEO}';flags.Sim.FastCalo.doEMECFCS=False" --postExec 'sim:cfg.getPublicTool("FastSimulationMasterTool").FastSimulations["FastCaloSim"].EkinMax=9999;cfg.getPublicTool("FastSimulationMasterTool").FastSimulations["FastCaloSim"].EkinMin=99;'

For completeness, following are Reco TF command:

Reco_tf.py  --inputHITSFile $HITSfile --outputAODFile /tmp/AOD_$tmp --autoConfiguration 'everything' --maxEvents $nevent --skipEvents=0 --conditionsTag 'default:OFLCOND-MC23-SDR-RUN3-06' --geometryVersion 'default:ATLAS-R3S-2021-03-02-00' --CA "all:True" --preInclude "all:Campaigns.MC23c" --steering 'doRAWtoALL' --postInclude "default:PyJobTransforms.UseFrontier" --AMITag r15616 --AMIConfig r15616 --preExec 'all:from SimulationConfig.SimEnums import LArParameterization;flags.Sim.LArParameterization=LArParameterization.FastCaloSim'

Tagging simulation conveners @schaarsc @mbandier and FCS convener @jbeirer .

Edited by Rui Zhang

Merge request reports