MCatNLO Step 3

On this page we will use the <>.event files that were produced by MCatNLO (see step 2) as input to Herwig. Using the Athena program we will use Herwig to have the top quarks decay and finally we will produce a combined ntuple (CBNT) and Pool file as was done for Rome.

1) Create a few directories

cd /project/atlas/users/ivov/MyMCatNLO/
mkdir MyMCatNLO_BASICS
mkdir MyMCatNLO_OUTPUT
mkdir MyMCatNLO_OUTPUT/InputAndLogFiles

2) Get necessary files required for Athena running

Since we are using MCatNLO we require a joboptions file and some other Athena-specific python files. In addition a file called inparmMcAtNlo.dat is expected that provides some input variables.

cd MyMCatNLO_BASICS
get_files PDGTABLE.MeV
get_files PartPropSvc.py

There are 2 'basic' input-files required. Each of them will be changed by a script that is talked about in step 4.

cp /afs/cern.ch/users/i/ivo/public/MCatNLO_scripts/joboptions_MyMCatNLO_Step3.py.BASIC .

The lines that will be changed by the script are related to the number of events, the random number seeds for the generators and the output files.

theApp.EvtMax = XNEVENTSX
AtRndmGenSvc.Seeds = ["JIMMY  XRNDJIMMY0X XRNDJIMMY1X", "JIMMY_INIT  XRNDJIMMY2X XRNDJIMMY3X",
                      "TAUOLA XRNDTAUOLA0X XRNDTAUOLA1X", "TAUOLA_INIT XRNDTAUOLA2X XRNDTAUOLA3X", 
                      "PHOTOS XRNDPHOTOS0X XRNDPHOTOS1X", "PHOTOS_INIT XRNDPHOTOS2X XRNDPHOTOS3X"]
NTupleSvc.Output     = [ "FILE1 DATAFILE='./XOutputDirCBNTX/MCatNLO.CBNT.JobXJOBNUMBERX.root' OPT='NEW'" ]
Stream1.OutputFile = "./XOutputDirPOOLX/MCatNLO.POOL.JobXJOBNUMBERX.root"

When using MCatNLO we require an additional file is present with some McatNLO settings. Also copy this file to the basics directory.

cp /afs/cern.ch/users/i/ivo/public/MCatNLO_scripts/joboptions_MyMCatNLO_Step3.py.BASIC .

 './XINPUTDIRX/XINPUTNAMEX'     ! event file
  XNEVENTSX                     ! number of events

3) Get script and taylor it to your needs

cd /project/atlas/users/ivov/MyMCatNLO/
cp /afs/cern.ch/users/i/ivo/public/MCatNLO_scripts/ShipOff_MyMCatNLO_Step3.py .

In the script we need to define a few variables of which the input and output directories are most important.

input_dir           = "/project/atlas/users/ivov/MyMCatNLO/MyMCatNLO_MATRIXELEMENT/"
output_dir          = "/project/atlas/users/ivov/MyMCatNLO/MyMCatNLO_OUTPUT/"
steering_files_dir  = "/project/atlas/users/ivov/MyMCatNLO/MyMCatNLO_BASICS/"

Nevents             =   12500
Nevents_joboptions  =       1
Njobs               =       1
f_LogFile           =       0

Important:

The variable called Nevents is linked to the number of events that are in the <>.events files. Make sure therefore that this number is the same as in the header of the <>.events file.

4)Run the script

When running the script a directory called Job1 filled with the files from the 'basic' directory and changed so that the variables are unique. The output Ntuple is stored in the output directory.