Kka054 at 12:39, 11 March 2009

2009-03-11T12:39:58Z

← Older revision		Revision as of 12:39, 11 March 2009
Line 202:		Line 202:
	#*Geometry Version		#*Geometry Version
	#<tt>triggerConfig</tt> (str)		#<tt>triggerConfig</tt> (str)
	#*Configuration string to use for <tt>TrigT1</tt> and <tt>HLT</tt>. Set to 'NONE'		#*Configuration string to use for <tt>TrigT1</tt> and <tt>HLT</tt>. Set to 'NONE' to switch off trigger, and set to 'DEFAULT' to use the default of the used release.
	to switch off trigger, and set to 'DEFAULT' to use the default of
	the used release.
	#[ <tt>jobConfig</tt> ] (list) default='NONE'		#[ <tt>jobConfig</tt> ] (list) default='NONE'
	#*Joboptions file with user settings, in particular the configuration settings. Default packages: <tt>RecJobTransforms</tt> , <tt>PyJobTransforms</tt>		#*Joboptions file with user settings, in particular the configuration settings. Default packages: <tt>RecJobTransforms</tt> , <tt>PyJobTransforms</tt>

Kka054 at 12:37, 11 March 2009

2009-03-11T12:37:49Z

← Older revision		Revision as of 12:37, 11 March 2009
Line 1:		Line 1:
	~~%META:TOPICINFO{author="ThomasBurgess" date="1209575726" format="1.1" reprev="1.2" version="1.2"}%~~
	~~%META:TOPICPARENT{name="WebHome"}%~~
	~~%TOC%~~

	== Walkthrough: ATLAS simulation using the CSC scripts ==		== Walkthrough: ATLAS simulation using the CSC scripts ==

Kka054: New page: %META:TOPICINFO{author="ThomasBurgess" date="1209575726" format="1.1" reprev="1.2" version="1.2"}% %META:TOPICPARENT{name="WebHome"}% %TOC% == Walkthrough: ATLAS simulation using the CSC ...

2009-03-11T12:37:32Z

New page: %META:TOPICINFO{author="ThomasBurgess" date="1209575726" format="1.1" reprev="1.2" version="1.2"}% %META:TOPICPARENT{name="WebHome"}% %TOC% == Walkthrough: ATLAS simulation using the CSC ...

New page

%META:TOPICINFO{author="ThomasBurgess" date="1209575726" format="1.1" reprev="1.2" version="1.2"}%
%META:TOPICPARENT{name="WebHome"}%
%TOC%

== Walkthrough: ATLAS simulation using the CSC scripts ==

If you follow the AthenaKlientDriftWalkThrough you can run the
examples in this guide.

Various file types of are mentioned in this text:

;'''Ntuple'''
:flat ntuples (usable directly in =root=)
;'''RDO''' '''''Raw Data Object'''''
:What comes out of the simulation
;'''ESD''' '''''Event Summary Data'''''
:Reconstructed '''RDO''' file, used to make '''AOD''' and for
;'''AOD''' '''''Analysis Object Data'''''
:Summary of an ''''ESD''' file
;'''DPD''' '''''Derived Physics Data'''''
:Analysis specific processed '''AOD'''

In general the files are <tt>root</tt> files, but not directly usable in <tt>root</tt>.

=== Event generation ===

To generate events use the script <tt>csc_evgen_trf.py</tt> that generates
events for a physics process defined in one of the <tt>jobOptions</tt> files
located in
<tt>$INSTALL/AtlasOffline/13.0.40/Generators/EvgenJobOptions/share</tt>. The
script takes the following arguments

#<tt>runNumber</tt> (int)
#*each run number corresponds to one physics process
#<tt>firstEvent</tt> (int)
#*number of the first event in the output data file
#<tt>maxEvents</tt> (int)
#*Maximum number of events to process
#<tt>randomSeed</tt> (int)
#*random seed for physics generators
#<tt>jobConfig</tt> (list)
#*<tt>jobOptions</tt> fragment containing the physics and the configuration settings
#<tt>outputEvgenFile</tt> (str)
#*Output file that contains generated events
#[ <tt>histogramFile</tt> ] (str) default='NONE'
#*Output file that contains histograms.
#[ <tt>ntupleFile</tt> ] (str) default='NONE'
#*Output file that contains ntuples.
#[ <tt>inputGeneratorFile</tt> ] (str) default='NONE'
#*Input file used by the particle generator to generate events

The following call generates a run with '''5000''' events for the process
'''Z->tau,tau''' and will save the file <tt>005188.evgen.pool.root</tt>.

<pre>
csc_evgen_trf.py\
005188 1 5000 1231243\
'CSC.005188.A3_Ztautau_filter.py'\
'005188.evgen.pool.root'
</pre>

Note that the run-number must match the =jobOptions= file given.

=== <tt>Atlfast</tt> '''AOD''' generation ===

The script <tt>csc_atlfast_trf.py</tt> truns <tt>Atlfast</tt> and produces '''AOD'''s
directly from the generated events. It takes the following arguments

#<tt>inputEvgenFile</tt> (list)
#*Input file that contains generated events
#<tt>outputAODFile</tt> (str)
#*Output file that contains '''AOD'''s
#<tt>ntupleFile</tt> (str)
#*Output file that contains ntuples.
#<tt>maxEvents</tt> (int)
#*Maximum number of events to process
#<tt>skipEvents</tt> (int)
#*Number of events to skip

To run on the file with generated events from the previous section run
<pre>
csc_atlfast_trf.py\
'005188.evgen.pool.root'\
'005188.atlfast.aod.pool.root'\
'005188.atlfast.nt.root'\
-1 0
</pre>

The '''-1''' includes all events in the input file

==== Script that makes several runs of atlfast aod ====

The following <tt>bash</tt> script generates '''nruns''' of '''nevts''' events of the type specified in ''''jobops'''.

<pre>
nruns=10
nevts="10000"
run="005188"
jobops=CSC.005188.A3_Ztautau_filter.py
for (( i=0; i<nruns; ++i));
do
evgen=${nevts}.${i}.${run}.evgen.pool.root
atlfast_aod=`echo ${evgen} | sed 's,evgen,atlfast.aod,'`
atlfast_nt=`echo ${evgen} | sed 's,evgen,atlfast.nt,'`
csc_evgen_trf.py ${run} 1 ${nevts} ${RANDOM} ${jobops} ${evgen}
csc_atlfast_trf.py ${evgen} ${atlfast_aod} ${atlfast_nt} -1 0
done
</pre>

=== Making '''ESD'''s '''AOD'''s using full ATLAS simulation ===

In these examples use the geometry definition <tt>ATLAS-CSC-01-02-00</tt> and the DEFAULT trigger. Remember that the '''Geant4''' step is very slow, run in a few events only until you know what you are doing.

==== '''GEANT4''' simulation and digitization ====

The script <tt>csc_simul_trf.py</tt> runs the '''GEANT4''' and digitization on
the generated events and produces a '''HITS''' and a '''RDO''' file. It has the
following options

#<tt>inputEvgenFile</tt> (list)
#*Input file that contains generated events
#<tt>outputHitsFile</tt> (str)
#*Output file that contains hits
#<tt>outputRDOFile</tt> (str)
#*Output file that contains RDO's
#<tt>maxEvents</tt> (int)
#*Maximum number of events to process
#<tt>skipEvents</tt> (int)
#*Number of events to skip
#<tt>randomSeed</tt> (int)
#*random seed for simulation
#<tt>geometryVersion</tt> (str)
#*Geometry Version
#<tt>digiSeedOffset1</tt> (int)
#*random seed offset for digitization
#<tt>digiSeedOffset2</tt> (int)
#*random seed offset for digitization
#[ <tt>physicsList</tt> ] (str) default='QGSP_EMV'
#*physics list to be used
#[ <tt>jobConfig</tt> ] (list) default='NONE'
#*Joboptions file with user settings, in particular the configuration settings. Default packages: .,SimuJobTransforms,PyJobTransforms
#[ <tt>DBRelease</tt> ] (str) default='NONE'
#*Tarball containing the DBRelease to use
#[ <tt>triggerConfig</tt> ] (str) default='DEFAULT'
#*Configuration string to use for =TrigT1= and =HLT=. Set to 'NONE' to switch off trigger, and set to 'DEFAULT' to use the default of the used release.

Run the full detector simulation:
<pre>
csc_simul_trf.py\
'005188.evgen.pool.root'\
'005188.hits.pool.root'\
'005188.rdo.pool.root'\
-1 0 123\
'ATLAS-CSC-01-02-00'\
12 23
</pre>

==== Reconstruction - create an '''ESD''' ====

The reconstruction can be run using the following script <tt>csc_recoESD_trf.py</tt> which takes the '''RDO''' file from the simulation and produces an '''ESD''' file. The scripts has the following options

#<tt>inputRDOFile</tt> (list)
#*Input file that contains RDOs
#<tt>outputESDFile</tt> (str)
#*Output file that contains ESDs
#<tt>ntupleFile</tt> (str)
#*Output file that contains ntuples.
#<tt>maxEvents</tt> (int)
#*Maximum number of events to process
#<tt>skipEvents</tt> (int)
#*Number of events to skip
#<tt>geometryVersion</tt> (str)
#*Geometry Version
#<tt>triggerConfig</tt> (str)
#*Configuration string to use for <tt>TrigT1</tt> and <tt>HLT</tt>. Set to 'NONE' to switch off trigger, and set to 'DEFAULT' to use the default of the used release.
#[ <tt>jobConfig</tt> ] (list) default='NONE'
#*Joboptions file with user settings, in particular the configuration settings. Default packages: RecJobTransforms, PyJobTransforms
#[ <tt>DBRelease</tt> ] (str) default='NONE'
#*Tarball containing the DBRelease to use

To run on the simulated '''RDO''' file run
<pre>
csc_recoESD_trf.py\
'005188.rdo.pool.root'\
'005188.esd.pool.root'\
'005188.rdo.ntuple.pool.root'\
-1 0\
'ATLAS-CSC-01-02-00'\
'DEFAULT'
</pre>

==== Create '''AOD''' from an '''ESD''' ====

The script <tt>csc_recoAOD_trf.py</tt> makes an '''AOD''' file from an '''ESD'''
file.

#<tt>inputESDFile</tt> (list)
#*Input file that contains ESDs
#<tt>outputAODFile</tt> (str)
#*Output file that contains AODs
#<tt>maxEvents</tt> (int)
#*Maximum number of events to process
#<tt>skipEvents</tt> (int)
#*Number of events to skip
#<tt>geometryVersion</tt> (str)
#*Geometry Version
#<tt>triggerConfig</tt> (str)
#*Configuration string to use for <tt>TrigT1</tt> and <tt>HLT</tt>. Set to 'NONE'
to switch off trigger, and set to 'DEFAULT' to use the default of
the used release.
#[ <tt>jobConfig</tt> ] (list) default='NONE'
#*Joboptions file with user settings, in particular the configuration settings. Default packages: <tt>RecJobTransforms</tt> , <tt>PyJobTransforms</tt>
#[ <tt>DBRelease</tt> ] (str) default='NONE'
#*Tarball containing the DBRelease to use

Run like this
<pre>
csc_recoAOD_trf.py\
'005188.esd.pool.root'\
'005188.aod.pool.root'\
-1 0\
'ATLAS-CSC-01-02-00'\
'DEFAULT'
</pre>

[[User:Kka054|Kka054]] 13:37, 11 March 2009 (CET)

ATLASSimulationCSCWalkThrough - Revision history

Kka054 at 12:39, 11 March 2009

Kka054 at 12:37, 11 March 2009

Kka054: New page: %META:TOPICINFO{author="ThomasBurgess" date="1209575726" format="1.1" reprev="1.2" version="1.2"}% %META:TOPICPARENT{name="WebHome"}% %TOC% == Walkthrough: ATLAS simulation using the CSC ...