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
  2. Click HERE to see the updated resolution plots with Error Bars, also eff when disable the inner-most and outer-most barrel layers (Oct 4th, 2010)

Grid submission-Athena and FTKsim

    Submitt trainning-muons to grid:
  1. Open a fresh terminal and go to /share/ftkdata1/tangjian/settings then run source setup_grid.sh
  2. go to /share/ftkdata1/tangjian/myScripts and look at a example script gridSubmission_10MSingleMuons_Oct10th2010.sh
  3. Notice the settings of --site AGLT2(or MWT2) and --outDS user.jtang.***
    Submitt FTKsim and FTKana to grid with remote input files:
  1. cd /share/ftkdata1/tangjian/settings and source setup_grid.sh
  2. go to /home/tangjian/FastTrackerAnalysis/summer2010/TrigFTKSim/scripts and use gridsub_ftk_test.sh to submit jiansplit.sh
  3. Look at jiansplit.sh, you need to submit the job in two steps:
    1. keep the command runSPLIT_reg and leave value of ftkDSblank, comment out command mergeSPLIT_reg in the end
    2. When the above job finishing running, change the script: comment out runSPLIT_reg and Un-comment mergeSPLIT_reg, use the output file name from step 1 to value ftkDS, don't forget the slash / in the end, then Re-submit the job again.
  4. If you want the grid job to keep root in addision to dgplot, you should go to file TrigFTKSim/scripts/FTKCore.py find the line with self.extra_args and delete the option --noROOT

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