There are various ways of getting MC produced:

1) Official Production

2) 'Organized' private production such as USATLAS RAC:

3) Private production that you do yourself. Some commands and tips below.

You should make sure that physics group conveners approve the production. They will have a quota of official production and so may suggest the other options. You should find someone who has produced MC recently (possibly via the conveners) and check with them that you have the most up-to-date standard settings.

For each datasets to be produced you will need:
- A dataset number (RunNumber)
- A cross section
- A job option (probably using pythia or herwig) that Evgen can run on
- A few validation plots on a twiki (from after the Evgen step)

You will have to run the Evgen step (presumably locally) and produce a few plots per dataset to 'validate' the setup. Then put these plots on a twiki and link to it for example here Pythia and Herwig are part of Athena, so if you have a Pythia card then you should be able to run Evgen directly on that. Otherwise you download, say, Alpgen, set it up and produce events (and possibly config) text files. Then to run Evgen you need a pythia/herwig jobOption to read the events files. Either search lxr for the job option of your favourite similar dataset: or look here:

Some twiki pages explaining how to get started:

esp. for non-Pythia:

Private Production with Pathena

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

I followed the dashboard of a previous request

Generation and Geant

Pythia event generation and Geant simulation were done in athena release

source -tag=,AtlasProduction,32,opt,setup

The transform scripts are: [options] <ecmenergy> <runnumber> <firstevent> [maxevents] <randomseed> <jobconfig> <outputevgenfile> [histogramfile] [ntuplefile] [inputgeneratorfile] [evgenjobopts]

and: [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. 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. It uses ~2GB RAM per job and is hence not allowed on xenia.

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

pathena --trf \
' \
7000 \
105587 \
0 \
5000 \
0 \
./ \
TMP.evgen.root; \
\ \
TMP.evgen.root \
%OUT.hits.root \
25 \
0 \
0 \
ATLAS-GEO-16-00-00 \
%DB:ddo.000001.Atlas.Ideal.DBRelease.v120201:DBRelease-12.2.1.tar.gz \
False \
--long \
--excludedSite \
--dbRelease \
ddo.000001.Atlas.Ideal.DBRelease.v120201:DBRelease-12.2.1.tar.gz \
--outDS \

Digitization and Reconstruction

Digitization and reconstruction were done in release

source -tag=,AtlasProduction

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

pathena --trf \
' \
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; \
\ \
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);*Units.GeV);muonRecFlags.writeSDOs=True" \
preInclude=RecJobTransforms/ \
postInclude_r2e=RecJobTransforms/ \
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 \

The two digiSeedOffsets must be changed for each subjob to change the pattern of noise in the calorimeter. was helpful as well as: and

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