TWiki> ATLAS Web>MonteCarloProduction (revision 4)EditAttach

MonteCarloProduction

This page shows how to produce Monte Carlo using MC10 settings with pathena.

There are many twiki pages explaining how to get started and set up generators.

https://twiki.cern.ch/twiki/bin/view/AtlasProtected/AtlasProductionGroup#Private_Simulation

https://twiki.cern.ch/twiki/bin/view/AtlasProtected/McProductionCommonParameters

https://twiki.cern.ch/twiki/bin/view/AtlasProtected/ExoticsMCRequestsHowTo#Setting_up_a_request_for_event_g

esp. for non-Pythia: https://twiki.cern.ch/twiki/bin/view/AtlasProtected/PreparingLesHouchesEven

If you haven't already emailed people to make sure that they approve the production and that you have the most up-to-date standard settings, this is one of the first things to do.

I followed the dashboard of a previous request http://www-f9.ijs.si/atlpy/atlprod/prodrequest/5832/

Generation and Geant

Pythia event generation and Geant simulation were done in athena release 15.6.12.9.

source setup.sh -tag=15.6.12.9,AtlasProduction,32,opt,setup

The transform scripts are:

Evgen_trf.py [options] <ecmenergy> <runnumber> <firstevent> [maxevents] <randomseed> <jobconfig> <outputevgenfile> [histogramfile] [ntuplefile] [inputgeneratorfile] [evgenjobopts]

and:

csc_atlasG4_trf.py [options] <inputevgenfile> <outputhitsfile> [maxevents] <skipevents> <randomseed> <geometryversion> [physicslist] [jobconfig] [dbrelease] [conditionstag] [dbcontent] [ignoreconfigerror] [amitag]

The command below generates 5000 events and runs Geant on the first 25. I understand that it's possible to use --nEventsPerJob and %SKIPEVENTS to run multiple subjobs on the same evgen file. https://twiki.cern.ch/twiki/bin/viewauth/Atlas/PandaAthena#example_8_How_to_run_production

Evgen_trf.py complained when I tried to generate less (or more) than 5000 events, but it's so much faster than Geant that it doesn't matter. Geant takes 4-8 mins per event

I made the random seed equal to the job number (0 in the command below)

pathena --trf \
'Evgen_trf.py \
7000 \
105587 \
0 \
5000 \
0 \
./MC10.105587.PythiaWprime800_WZ_qqee.py \
TMP.evgen.root; \
\
csc_atlasG4_trf.py \
TMP.evgen.root \
%OUT.hits.root \
25 \
0 \
0 \
ATLAS-GEO-16-00-00 \
QGSP_BERT \
VertexFromCondDB.py,CalHits.py,ParticleID.py \
%DB:ddo.000001.Atlas.Ideal.DBRelease.v120201:DBRelease-12.2.1.tar.gz \
OFLCOND-SDR-BS7T-02 \
NONE \
False \
NONE' \
--long \
--excludedSite \
RALPP,GRIF-LAL,SLAC \
--dbRelease \
ddo.000001.Atlas.Ideal.DBRelease.v120201:DBRelease-12.2.1.tar.gz \
--outDS \
user.apenson.MC10.105587.PythiaWprime800_WZ_qqee.evgen.G4.25.0

Digitization and Reconstruction

Digitization and reconstruction were done in release 16.0.2.3.

source setup.sh -tag=16.0.2.3,AtlasProduction

This command take the hits file generated by the command above and runs two more transforms on it:

pathena --trf \
'Digi_trf.py \
inputHitsFile=%IN \
outputRDOFile=TMP.RDO.root \
maxEvents=-1 \
skipEvents=0 \
DBRelease=%DB:ddo.000001.Atlas.Ideal.DBRelease.v120902:DBRelease-12.9.2.tar.gz \
conditionsTag=OFLCOND-SDR-BS7T-04-02 \
geometryVersion=ATLAS-GEO-16-00-00 \
samplingFractionDbTag=QGSP_BERT \
digiSeedOffset1=1 \
digiSeedOffset2=0; \
\
Reco_trf.py \
inputRDOFile=TMP.RDO.root \
outputESDFile=%OUT.ESD.pool.root \
outputAODFile=%OUT.AOD.pool.root \
DBRelease=%DB:ddo.000001.Atlas.Ideal.DBRelease.v120902:DBRelease-12.9.2.tar.gz \
autoConfiguration=everything \
conditionsTag=OFLCOND-SDR-BS7T-04-02 \
geometryVersion=ATLAS-GEO-16-00-00 \
preExec="rec.Commissioning.set_Value_and_Lock(True);jobproperties.Beam.energy.set_Value_and_Lock(3500*Units.GeV);muonRecFlags.writeSDOs=True" \
preInclude=RecJobTransforms/SetJetConstants-02-000.py \
postInclude_r2e=RecJobTransforms/CalibrationHitsInESDConfig.py \
triggerConfig=MCRECO:DB:TRIGGERDBMC:248,108,194' \
--long \
--dbRelease \
ddo.000001.Atlas.Ideal.DBRelease.v120902:DBRelease-12.9.2.tar.gz \
--inDS \
user.apenson.MC10.105587.PythiaWprime800_WZ_qqee.evgen.G4.25.0/ \
--outDS \
user.apenson.MC10.105587.PythiaWprime800_WZ_qqee.0

The two digiSeedOffsets must be changed for each subjob to change the pattern of noise in the calorimeter.

https://twiki.cern.ch/twiki/bin/viewauth/Atlas/RecoTrf#Special_syntax_for_preExec_postE was helpful as well as:

hn-atlas-dist-analysis-help@cern.ch and hn-atlas-job-transformations@cern.ch

Edit | Attach | Watch | Print version | History: r5 < r4 < r3 < r2 < r1 | Backlinks | Raw View | Raw edit | More topic actions...
Topic revision: r4 - 2010-12-17 - AlexPenson
 
This site is powered by the TWiki collaboration platform Powered by PerlCopyright © 2008-2021 by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
Ideas, requests, problems regarding TWiki? Send feedback