Double counting of delta-rays suppression - sample recipe

This TWiki page is designed to produce private samples from a ParticleGun jobOptions file up to an ntuple. The gas geometry considered here is a full xenon geometry.

Event Generation

The first step is to generate some EVNT files thanks to ParticleGun. The jobOptions included to this TWiki generate 31000 muons or electrons tracks having same eta-pT spectra as it was for electrons and muons in 2013 pPb runs. In a new folder an athena release by doing :

  • setupATLAS
  • asetup 19.2.1.5,AtlasProduction,here

Then simply run athena with the corresponding jobOptions file here :

  • athena jobOption.ParticleGun_corrhist_el.py
Or here :
  • athena jobOption.ParticleGun_corrhist_mu.py
Make sure that the corresponding eta-pT histograms attached to this page are in the same folder.

Simulation

From the EVNT file produced previously we will run geant 4. In a new folder :

If you want to simulate samples with delta-rays fixed you should also modify the G4AtlasApps package :
  • sed -i '153s/0.05/30.0/' Simulation/G4Atlas/G4AtlasApps/python/atlas_idet.py
  • sed -i '154s/0.05/30.0/' Simulation/G4Atlas/G4AtlasApps/python/atlas_idet.py
Then you need to compile :
  • setupWorkArea.py
  • cd WorkArea/cmt
  • cmt br cmt config
  • source setup.sh
  • cmt br cmt make
Finally you can either submit the simulation to the GRID thanks to the script attached here that you would have to adapt a bit or locally by doing :
  • cd ../run
  • AtlasG4 _tf.py --inputEvgenFile SingleElectrons.2013.pool.root --outputHitsFile SingleElectrons.2013.HITS.pool.root --maxEvents 30000 --randomSeed 12345 --geometryVersion ATLAS-R2-2015-02-01-00_VALIDATION --conditionsTag OFLCOND-RUN12-SDR-21 --DataRunNumber 222222
In each case you would need to be sure that the previously generated EVNT files are in the run directory.

Digitization

In a clean environment in an empty directory do the following:

  • setupATLAS
  • asetup 20.1.5.12,AtlasProduction,here
  • pkgco.py TRT_Digitization-01-01-03
  • pkgco.py TRT_PAI_Process-00-01-01
  • setupWorkArea.py
  • cd WorkArea/cmt
  • cmt br cmt config
  • source setup.sh
  • cmt br cmt make
  • cd ../run
Now you are ready to run the digitization either locally or on the GRID with the script here :
  • Reco_tf.py --ignoreErrors True --maxEvents 30000 --inputHITSFile SingleElectrons.2013.HITS.pool.root --outputRDOFile SingleElectrons.2013.RDO.pool.root --conditionsTag OFLCOND-RUN12-SDR-21 --geometryVersion ATLAS-R2-2015-02-01-00_VALIDATION --DataRunNumber 222222 --postInclude 'TRT_Digitization/postInclude.OverrideTRTparameters.py,RecJobTransforms/UseFrontier.py' --autoConfiguration 'everything'
If you want to modify the digitization parameters like the high-threshold values or the TR efficiencies, this should be done at this step by modifying the file TRT_Digitization/postInclude.OverrideTRTparameters.py before running the digitization. The default values of this digitization package correspond to the delta-rays suppressed tune.

It is highly adviced when you have to generate samples with different threshold and TR efficiencies for both muons and electrons with or without the delta-rays fix to modify the script to automatically submit all the jobs.

Reconstruction and Ntuples

In a clean environment in a new directory download the file here, then do :
  • tar -xvzf ntup.tar.gz
  • setupATLAS
  • asetup 20.1.5.12,AtlasProduction,here
  • setupWorkArea.py
  • cd WorkArea/cmt
  • cmt br cmt config
  • source setup.sh
  • cmt br cmt make
  • cd ../..
Now after modifying the input and path in either the job options or the script you can run either locally :
  • athena jobOptions_RecExCommon.py
Or on the GRID :
  • ./pathena.sh

-- AnthonyEricLionti - 2016-01-25

Topic attachments
I Attachment History Action Size Date Who Comment
Unix shell scriptsh Digi.sh r1 manage 0.9 K 2016-01-25 - 17:59 AnthonyEricLionti  
Unix shell scriptsh Simu.sh r1 manage 0.8 K 2016-01-25 - 17:36 AnthonyEricLionti  
Unknown file formatroot data_histos_el_1470pt.root r1 manage 98.8 K 2016-01-25 - 16:50 AnthonyEricLionti  
Unknown file formatroot data_histos_mu_run218716pt.root r1 manage 314.5 K 2016-01-25 - 16:50 AnthonyEricLionti  
Texttxt jobOption.ParticleGun_corrhist_el.py.txt r1 manage 1.5 K 2016-01-25 - 16:50 AnthonyEricLionti  
Texttxt jobOption.ParticleGun_corrhist_mu.py.txt r1 manage 1.5 K 2016-01-25 - 16:50 AnthonyEricLionti  
Unknown file formatgz ntup.tar.gz r1 manage 208.1 K 2016-01-26 - 01:02 AnthonyEricLionti  
Edit | Attach | Watch | Print version | History: r3 < r2 < r1 | Backlinks | Raw View | WYSIWYG | More topic actions
Topic revision: r3 - 2016-01-26 - AnthonyEricLionti
 
    • 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