Single Particle Study (new)

Event Generation

This topic will probably be covered later once someone needs to produce new single particle EVNT files. For now there's a technique which makes use of the discontinued ParticleGenerator package described here. There are also EVNT files with 50k muon and electron tracks having same eta-pT spectra as it was for electrons and muons in 2013 pPb runs:

/afs/cern.ch/user/m/maevskiy/public/singleParticleStudyNew/EvtGen/el.pPb.EVNT.pool.root
/afs/cern.ch/user/m/maevskiy/public/singleParticleStudyNew/EvtGen/mu.pPb.EVNT.pool.root

Simulation

In a clean environment in an empty directory do the following:
asetup 17.3.12.3,AtlasProduction,here,slc5
cmt co -r TRT_GeoModel-00-02-57 InnerDetector/InDetDetDescr/TRT_GeoModel
cmt co -r G4AtlasApps-00-07-43-27 Simulation/G4Atlas/G4AtlasApps
setupWorkArea.py
cd WorkArea/cmt
cmt br cmt config
source setup.sh
cmt br cmt make
cd ../../

To run it simply do:

AtlasG4_trf.py \
inputEvgenFile=/afs/cern.ch/user/m/maevskiy/public/singleParticleStudyNew/EvtGen/mu.pPb.EVNT.pool.root \
outputHitsFile=outputFile.HITS.pool.root \
maxEvents=5 skipEvents=0 randomSeed=1237 \
geometryVersion=ATLAS-R2-2015-01-05-00_VALIDATION \
conditionsTag=OFLCOND-RUN12-SDR-12 \
postInclude="RecJobTransforms/UseFrontier.py" \
postExec='
from IOVDbSvc.CondDB import conddb
conddb.addFolderSplitOnline("TRT","/TRT/Onl/Cond/StatusHT","/TRT/Cond/StatusHT")
' \
DataRunNumber=1000000

This will run within the "baseline" Argon scenario. In order to change it to the "worstcase" scenario set OFLCOND-RUN12-SDR-13 instead of OFLCOND-RUN12-SDR-12.

Alternative way (yet to check):

ln -s /afs/cern.ch/atlas/project/database/DBREL/packaging/DBRelease/current/geomDB geomDB
ln -s /afs/cern.ch/atlas/project/database/DBREL/packaging/DBRelease/current/triggerDB triggerDB
ln -s /afs/cern.ch/user/a/atlcond/coolrep/sqlite200 sqlite200

AtlasG4_trf.py \
inputEvgenFile=/afs/cern.ch/user/m/maevskiy/public/singleParticleStudyNew/EvtGen/mu.pPb.EVNT.pool.root \
outputHitsFile=outputFile.HITS.pool.root \
maxEvents=5 skipEvents=0 randomSeed=1237 \
geometryVersion=ATLAS-GEO-20-00-01 \
conditionsTag=OFLCOND-MC12-SDR-06 \
postExec='
from IOVDbSvc.CondDB import conddb
conddb.blockFolder("/TRT/Cond/StatusHT")
conddb.addFolderWithTag("TRT_OFL","/TRT/Cond/StatusHT","TrtStrawStatusHT-MC-run2-scenario1_00-00",force=True,forceMC=True)
' \
preExec='
from GeoModelSvc.GeoModelSvcConf import GeoModelSvc
GeoModelSvc.TRT_VersionOverride="TRT-GEO-02"'

# to switch from baseline argon scenario to worst case scenario change
# "TrtStrawStatusHT-MC-run2-scenario1_00-00" to "TrtStrawStatusHT-MC-run2-scenario2_00-00"

Digitization

In a clean environment in an empty directory do the following:
asetup 17.3.12.3,AtlasProduction,here,slc5
cmt co -r TRT_Digitization-00-11-08-04 InnerDetector/InDetDigitization/TRT_Digitization
setupWorkArea.py
cp /afs/cern.ch/user/m/maevskiy/public/singleParticleStudyNew/Digi/TRTProcessingOfStraw.cxx InnerDetector/InDetDigitization/TRT_Digitization/src/
cd WorkArea/cmt/
cmt br cmt config
source setup.sh
cmt br cmt make
cd ../../

Whenever you perform any change with the source files you should repeat these steps:

cd WorkArea/cmt/
cmt br cmt make
cd ../../
And make sure that compilation was successful.

To run digi do:

cd WorkArea/run/
ln -s /afs/cern.ch/atlas/project/database/DBREL/packaging/DBRelease/current/geomDB geomDB
ln -s /afs/cern.ch/atlas/project/database/DBREL/packaging/DBRelease/current/triggerDB triggerDB
ln -s /afs/cern.ch/user/a/atlcond/coolrep/sqlite200 sqlite200

Digi_trf.py \
maxEvents=5 \
outputRDOFile=OUT.RDO.pool.root \
digiRndmSvc='AtDSFMTGenSvc' \
inputHitsFile=INPUT.HITS.FILENAME.root \
conditionsTag='OFLCOND-RUN12-SDR-12' \
autoConfiguration='everything' \
postExec='
from IOVDbSvc.CondDB import conddb
conddb.addFolderSplitOnline("TRT","/TRT/Onl/Cond/StatusHT","/TRT/Cond/StatusHT")
from TRT_Digitization.TRT_DigitizationConf import TRTDigitizationTool
TRTDigitizationTool.UseArgonStraws = True
TRTDigitizationTool.UseConditionsHTStatus = True' \
geometryVersion='ATLAS-R2-2015-01-05-00_VALIDATION'

cd ../../

Again, change OFLCOND-RUN12-SDR-12 to OFLCOND-RUN12-SDR-13 for the worst case scenario.

Alternative way (not checked yet):

cd WorkArea/run/
ln -s /afs/cern.ch/atlas/project/database/DBREL/packaging/DBRelease/current/geomDB geomDB
ln -s /afs/cern.ch/atlas/project/database/DBREL/packaging/DBRelease/current/triggerDB triggerDB
ln -s /afs/cern.ch/user/a/atlcond/coolrep/sqlite200 sqlite200

Digi_trf.py \
maxEvents=5 \
outputRDOFile=OUT.RDO.pool.root \
digiRndmSvc='AtDSFMTGenSvc' \
inputHitsFile=INPUT.HITS.FILENAME.root \
preExec='
from GeoModelSvc.GeoModelSvcConf import GeoModelSvc
GeoModelSvc.TRT_VersionOverride="TRT-GEO-02"
' \
conditionsTag=OFLCOND-MC12-SDR-06 \
autoConfiguration='everything' \
postExec='
from IOVDbSvc.CondDB import conddb
conddb.blockFolder("/TRT/Cond/StatusHT")
conddb.addFolderWithTag("TRT_OFL","/TRT/Cond/StatusHT","TrtStrawStatusHT-MC-run2-scenario1_00-00",force=True,forceMC=True)
from TRT_Digitization.TRT_DigitizationConf import TRTDigitizationTool
TRTDigitizationTool.UseArgonStraws = True
TRTDigitizationTool.UseConditionsHTStatus = True' \
geometryVersion=ATLAS-GEO-20-00-01

cd ../../

# to switch from baseline argon scenario to worst case scenario change
# "TrtStrawStatusHT-MC-run2-scenario1_00-00" to "TrtStrawStatusHT-MC-run2-scenario2_00-00"

-- ArtemMaevskiy - 2014-11-07

Edit | Attach | Watch | Print version | History: r4 < r3 < r2 < r1 | Backlinks | Raw View | WYSIWYG | More topic actions
Topic revision: r4 - 2014-11-22 - 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