How to Use Athena

Codes

check out Athena package
Digitization Codes

Run Athena Locally

    How to use the script to submit the muon trainning jobs:
  1. Refer to /home/tangjian/submissions/submit_training_set1_set2_bkg.sh or THIS
  2. go to your athena dir: /home/tangjian/FastTrackerAnalysis/summer2010 and type source ./setup.sh
  3. run commands from submit_training_set1_set2_bkg.sh separatly:
    1. athena.py -c 'RandomSeed1 = 78200 ; RandomSeed2=8000200 ' jo.TrainingMuons.ParticleGeneratorFullEta.py (this generate four-momentum for all event. Note: remove %RNDM from the .sh script to run locally)
    2. ./csc_atlasG4_trf.py ParticleGenerator.pool.root hits.pool.root -1 0 112200 ATLAS-GEO-02-01-00 QGSP_BERT NONE NONE OFLCOND-SIM-00-00-06 jobConfig=jo.IDonly.py,jo.CaloMuonOffSimulDigitConfig.py,jo.ZeroedBeamWidth.py (Note: replace the "%DB=..." option with "NONE", and 112200 is the random number and you have to remove '%RNDM=' in all the scripts )
    3. ./csc_digi_trf.py hits.pool.root rdo.pool.root -1 0 ATLAS-GEO-02-01-00 98200 84200 NONE NONE jo.IDonly.py,jo.DigiIDonly.py,jo.DigiAllNoiseOff.py,jo.DigiZeroPileup.py,jo.ZeroedBeamWidth.py NONE AtRanluxGenSvc QGSP_BERT default NONE NONE NONE False OFLCOND-SIM-00-00-06 (This is the digitization process. Note: again replace the "%DB=..." option with "NONE", also -1 option is for the default 20,000 events, you can change it to whatever positive number you want, say 100)
    4. athena.py jo.RecoWrapper.py (this generate wrapper files)
    5. bzip2 ftksim_wrap_raw.dat (tarball the wrapper file)

    Generate tracks from wrapper files--How to compile correctly:
  1. cd ~tangjian/FastTrackerAnalysis/summer2010/
  2. source ~antonk/setup/setup_all_public.sh
  3. echo $ROOTSYS and echo $SITEROOT to check the version
  4. cd TrigFTKSim/standalone
  5. source bootstrap.sh
  6. gmake clean
  7. gmake -f Makefile.grid.mk -j4 (-j4 uses 4 CPUs)
  8. [to check if it works:] $./road_finder (this test the compiled codes)

  9. cd ../../TrigFTKAna
  10. source bootstrap.sh
  11. gmake clean
  12. gmake -j4
  13. [to check if it works:] $ ana/ana_standard_overview
  14. ** you can refer to HERE and HERE
  15. To Run SingleMuon events:
    ./ftk.py -a fullchain,STOP,ana(optional) -d raw_11L_50x64x144_20KSingleMuonsJuly23rd -n100 -o MAXPATTERNS=8000000 -r /share/ftkdata1/tangjian/deadchannels -t /share/ftkdata/antonk/bases/v14 -p -m4 -q pbs(or use 'threaded_4')
  16. [see THIS error log, jobs in tier3 are deleted for some user randomly]

    What to do after change the digitization codes:
  1. go to the Digitization code dir: /home/tangjian/FastTrackerAnalysis/summer2010/InnerDetector/InDetDigitization
  2. apply the code changes you want (will be detailed later)
  3. go to the "cmt" dir of each of the subdirs of InDetDigitization and type "cmt make"
  4. then submit muon trainning jobs locally following the instructions above

    Results:
  1. Click HERE to see the result presented on Sep 20th, 2010

Disable Channels

    I: Disable All Pixels
  1. In InnerDetector/InDetDigitization/SiDigitization/src/SiRandomDisabledCellGenerator.cxx [Line 87: replace 'if(RandFlat::**)' with if(true), which results in only sct hits in the wrapper file, so this is where we can work on dead-channel enabling stuff.]
    II: Disable SCT channels (add parameter in jobOption file)
  1. In folder './SCT_Digitization/share/SCT_Digitization_jobOptions.py' add 'sct.DisableProbability = 0.5555' to set the disable prob of SCT
  2. In constructor of 'SCT_Digitization' add 'declareProperty("DisableProbability", m_disableProbability=0.0, "Testing new disable prob variable");'
  3. In 'SCT_Digitization.cxx' method 'SCT_Digitization::initDisabledCells()', add 'm_sct_RandomDisabledCellGenerator->setDisableProbability(m_disableProbability);'

Documents

How Silicon Detector Works
Tracts in Modern Physics, 2008 Hartmann

Transmitters
Slides I
Slides II
Slides III