Single Particle Study

Event Generation

For generating the events with single particles I use ParticleGenerator. To run event generation using Athena one can use Generate_trf.py, which also requires a job options file. Here is an example of a job options file for running ParticleGenerator:

from AthenaCommon.AlgSequence import AlgSequence 
topAlg = AlgSequence("TopAlg") 

from ParticleGenerator.ParticleGeneratorConf import ParticleGenerator
topAlg += ParticleGenerator()


ParticleGenerator = topAlg.ParticleGenerator

# For VERBOSE output from ParticleGenerator.
ParticleGenerator.OutputLevel = 1

ParticleGenerator.orders = [
 "PDGcode: sequence -13 13",
 "p: constant 20000",
 "eta: flat -0.7 0.7",
 "phi: flat -3.14159 3.14159"
 ]


from EvgenJobTransforms.EvgenConfig import evgenConfig

evgenConfig.generators += [ "ParticleGenerator" ]
evgenConfig.description = "simpleJobOption"

These job options will in each event create one muon with 20 GeV momentum (not pt!), all muons being flatly distributed over -0.7<|eta|<0.7. In order to run such a job options file I login to lxplus, cd to a convenient directory and type following lines:

asetup 17.2.8.8,AtlasProduction,here,slc5
Generate_trf.py jobConfig=../MyJobOptionsFile.py firstEvent=1 maxEvents=10 randomSeed=12345 runNumber=54321 outputEVNTFile=outputFile.EVNT.pool.root ecmEnergy=8000

Simulation

To run simulation with full Xenon configuration I did the following:

asetup 17.2.5.2,AtlasProduction,here,slc5
cmt co -r TRT_GeoModel-00-02-42 InnerDetector/InDetDetDescr/TRT_GeoModel
setupWorkArea.py
cd WorkArea/cmt
cmt br cmt config
source setup.sh
cmt br cmt make
cd ../run
AtlasG4_trf.py inputEvgenFile=../eventsFile.EVNT.pool.root outputHitsFile=outputFile.HITS.pool.root maxEvents=10 skipEvents=0 randomSeed=1237 geometryVersion=ATLAS-GEO-20-00-01 conditionsTag=OFLCOND-MC12-SDR-06

To use implementation of Argon in layer0 made by Alex everything is the same (except for the cmt co line, of course). You can find the code in /afs/cern.ch/user/m/maevskiy/public/singleParticleStudy/AlexCode.

Digitization

Sample text

D3PD production

Sample text

D3PDs processing

-- ArtemMaevskiy - 17 Jan 2014

Edit | Attach | Watch | Print version | History: r5 | r4 < r3 < r2 < r1 | Backlinks | Raw View | Raw edit | More topic actions...
Topic revision: r1 - 2014-01-17 - ArtemMaevskiy
 
    • Cern Search Icon Cern Search
    • TWiki Search Icon TWiki Search
    • Google Search Icon Google Search

    Main All webs login

This site is powered by the TWiki collaboration platform Powered by PerlCopyright & 2008-2019 by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
Ideas, requests, problems regarding TWiki? Send feedback