Difference between revisions of "MCatNLO Step 3"

From Atlas Wiki
Jump to navigation Jump to search
 
(4 intermediate revisions by the same user not shown)
Line 14: Line 14:
 
mkdir ATHENA_BASICS
 
mkdir ATHENA_BASICS
 
mkdir POOLCBNT
 
mkdir POOLCBNT
mkdir POOLCBNT/InputAndLogFiles
+
mkdir POOLCBNT/InputAndLogfiles
 
</pre>
 
</pre>
 
</font>
 
</font>
 
  
 
==2) Get the basic input files ==
 
==2) Get the basic input files ==
Line 40: Line 39:
 
</font>
 
</font>
  
This is a good check to see if you have set-up athena correctly.
+
This is a good check to see if you have set-up athena correctly. If this fails please make sure you have the Athena envoronment set-up correctly.
  
 
<b>(b) Get the basic joboptions file</b>
 
<b>(b) Get the basic joboptions file</b>
Line 100: Line 99:
 
<b>(b) Tailor the script to your environment</b>
 
<b>(b) Tailor the script to your environment</b>
  
In the script you need to define the variables related to the input and output directories. Look for the tag "Main Program".
+
In the script you need to define the variables related to the input and output directories. These are defined in the first "routine" in the script.
  
 
<font color=blue size="-1">  
 
<font color=blue size="-1">  
Line 110: Line 109:
 
</font>
 
</font>
  
There are however also parameters related to the number of events, the number of jobs etc etc. Look again for the tag "Main Program".
+
At the same place in the script also parameters related to the number of events, the number of jobs etc have to be set.
  
 
<font color=blue size="-1">  
 
<font color=blue size="-1">  
Line 134: Line 133:
 
</font>  
 
</font>  
  
Once you get it running for a single job, set the parameter f_Logfile to 1 and produce many files where the logfile is automatically written to the output directory.
+
Once you get it running for a single job, set the parameter <b>f_Logfile=1</b> and <b>Njobs=2</b> to produce two files  
 +
where the logfile is automatically written to the output  
 +
directory.
  
 
<u>You are finished</u>
 
<u>You are finished</u>
  
To POOL files you have produced can then be transformed to AOD after full reconstruction or Atlfast.
 
  
<b> What did we produce </b>
+
<b> What did we produce ??</b>
  
 
<font color=blue size="-1">  
 
<font color=blue size="-1">  
 
<pre>
 
<pre>
 
[ivov@njord MyMCatNLO]$ ls -lrt ./POOLCBNT/*
 
[ivov@njord MyMCatNLO]$ ls -lrt ./POOLCBNT/*
-rw-r--r--    1 ivov    atlas     224799 May 10 00:21 ./POOLCBNT/MCatNLO.POOL.Job1.root
+
-rw-r--r--    1 ivov    atlas     2511240 May 10 11:38 ./POOLCBNT/MCatNLO.POOL.Job1.root
-rw-r--r--    1 ivov    atlas     103891 May 10 00:21 ./POOLCBNT/MCatNLO.CBNT.Job1.root
+
-rw-r--r--    1 ivov    atlas     1044531 May 10 11:38 ./POOLCBNT/MCatNLO.CBNT.Job1.root
-rw-r--r--    1 ivov    atlas     218281 May 10 00:22 ./POOLCBNT/MCatNLO.POOL.Job2.root
+
-rw-r--r--    1 ivov    atlas     2429914 May 10 11:38 ./POOLCBNT/MCatNLO.POOL.Job2.root
-rw-r--r--    1 ivov    atlas     102282 May 10 00:22 ./POOLCBNT/MCatNLO.CBNT.Job2.root
+
-rw-r--r--    1 ivov    atlas     1014065 May 10 11:38 ./POOLCBNT/MCatNLO.CBNT.Job2.root
  
 
./POOLCBNT/InputAndLogfiles:
 
./POOLCBNT/InputAndLogfiles:
total 568
+
total 540
-rw-r--r--    1 ivov    atlas      276179 May 10 00:21 Logfile.Job1.log
+
-rw-r--r--    1 ivov    atlas      250943 May 10 11:38 Logfile.Job1.log
-rw-r--r--    1 ivov    atlas        5113 May 10 00:21 joboptions_Job1.py
+
-rw-r--r--    1 ivov    atlas        5114 May 10 11:38 joboptions_Job1.py
-rw-r--r--    1 ivov    atlas        492 May 10 00:21 inparmMcAtNlo.dat.Job1
+
-rw-r--r--    1 ivov    atlas        492 May 10 11:38 inparmMcAtNlo.dat.Job1
-rw-r--r--    1 ivov    atlas      268323 May 10 00:22 Logfile.Job2.log
+
-rw-r--r--    1 ivov    atlas      263851 May 10 11:38 Logfile.Job2.log
-rw-r--r--    1 ivov    atlas        5113 May 10 00:22 joboptions_Job2.py
+
-rw-r--r--    1 ivov    atlas        5114 May 10 11:38 joboptions_Job2.py
-rw-r--r--    1 ivov    atlas        492 May 10 00:22 inparmMcAtNlo.dat.Job2
+
-rw-r--r--    1 ivov    atlas        492 May 10 11:38 inparmMcAtNlo.dat.Job2
 
</pre>
 
</pre>
 
</font>
 
</font>
 +
 +
To POOL files you have produced can then be transformed to AOD after full reconstruction or Atlfast.

Latest revision as of 09:39, 10 May 2005

On this page we will use the <>.event files that were produced by MCatNLO and use Herwig (using the Athena framework) to decay the tops and produce the parton shower. We will produce a combined ntuple CBNT and POOL file.

Main idea:

The main idea is to call athena with a joboptions file where you specify the location of the event file. In addition it is required to have an additional input file. We will again use a script that will take the BASIC forms for these input files and prodice unique input files for each job and take care of the bookkeeping.

1) Set-up Athena and output directories

We will use Athena to do the decay. Make sure that you have correctly set up the athena environment. To see how to do this please check xxx.

cd /MyMCatNLO
mkdir ATHENA_BASICS
mkdir POOLCBNT
mkdir POOLCBNT/InputAndLogfiles

2) Get the basic input files

Since we have used MCatNLO to produce the and momenta, we have to signal to Herwig that it must these from your external <>.event-file. In Athena, this is done inside a joboptions file. Again we start at again to your MCatNLO working directory area. Create two new directories where we'll store the basic input files (basic joboptions file and other athena related files) and the output files (the combined ntuple for the generator information and POOL files).

cd /MyMCatNLO/ATHENA_BASICS

(a) Get some Athena python-files and check set-up

Get these athena python files into your BASICS area using get_files.

get_files PDGTABLE.MeV
get_files PartPropSvc.py

This is a good check to see if you have set-up athena correctly. If this fails please make sure you have the Athena envoronment set-up correctly.

(b) Get the basic joboptions file

Get the basic joboptions file: joboptions_MCatNLO_Rome_POOLCBNT.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"


(c) Get the inparmMcAtNlo.dat (MCatNLO parameter) file

When starting the athena job, a file called inparmMcAtNlo.dat is expected in your run directory. It contains the location of the <>.event file and some additional input variables related to the generated MCatNLO sample.

Get the additional MCatNLO parameter file: inparmMcAtNlo.dat.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.

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

3) Get the script and taylor it to your needs

The last step is to get the script and start the jobs. For each job it will:

o create a directory called Job<job#>
o create two unique input files
o copy all otehr nexessary files here and start the job
o copy output files (CBNT,POOL and Logfiles) to an output directory
o remove directory Job<job#>

(a) Get the script

Go to your working area:

cd /MyMCatNLO/

and get the python script: ShipOff_MCatNLO_Rome_POOLCBNT.py


(b) Tailor the script to your environment

In the script you need to define the variables related to the input and output directories. These are defined in the first "routine" in the script.

input_dir           = "/data/atlas/users/ivov/MyMCatNLO/MATRIXELEMENT/"
output_dir          = "/data/atlas/users/ivov/MyMCatNLO/POOLCBNT/"
steering_files_dir  = "/data/atlas/users/ivov/MyMCatNLO/ATHENA_BASICS/"

At the same place in the script also parameters related to the number of events, the number of jobs etc have to be set.

Nevents_eventfile   =     500
Nevents_joboptions  =      20
Njobs               =       2
f_LogFile           =       0

Important: The variable called Nevents_eventfile 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 (should be made automatic). The variable Nevents_joboptions controls the number of events that will be processed.

4) Run the script

Once the number of events and the number of jobs are defined start the run:

./ShipOff_MCatNLO_Rome_POOLCBNT.py

Once you get it running for a single job, set the parameter f_Logfile=1 and Njobs=2 to produce two files where the logfile is automatically written to the output directory.

You are finished


What did we produce ??

[ivov@njord MyMCatNLO]$ ls -lrt ./POOLCBNT/*
-rw-r--r--    1 ivov     atlas     2511240 May 10 11:38 ./POOLCBNT/MCatNLO.POOL.Job1.root
-rw-r--r--    1 ivov     atlas     1044531 May 10 11:38 ./POOLCBNT/MCatNLO.CBNT.Job1.root
-rw-r--r--    1 ivov     atlas     2429914 May 10 11:38 ./POOLCBNT/MCatNLO.POOL.Job2.root
-rw-r--r--    1 ivov     atlas     1014065 May 10 11:38 ./POOLCBNT/MCatNLO.CBNT.Job2.root

./POOLCBNT/InputAndLogfiles:
total 540
-rw-r--r--    1 ivov     atlas      250943 May 10 11:38 Logfile.Job1.log
-rw-r--r--    1 ivov     atlas        5114 May 10 11:38 joboptions_Job1.py
-rw-r--r--    1 ivov     atlas         492 May 10 11:38 inparmMcAtNlo.dat.Job1
-rw-r--r--    1 ivov     atlas      263851 May 10 11:38 Logfile.Job2.log
-rw-r--r--    1 ivov     atlas        5114 May 10 11:38 joboptions_Job2.py
-rw-r--r--    1 ivov     atlas         492 May 10 11:38 inparmMcAtNlo.dat.Job2

To POOL files you have produced can then be transformed to AOD after full reconstruction or Atlfast.