Replicating Pass 6 Cuts
I've gone back to Eric's old page with the hope of picking up where he left off. Of course, this leads to the always enjoyable task of replicating the work that he has already done. The working directory for the previous analysis is:
/afs/slac/u/ek/echarles/vol2/glast_workdirs/GR_v15r50p3
From there, I have been able to more or less reconstruct the analysis that Eric was performing. Building my own version of GR v15r50p3, I used the runme.csh and runme_ag.csh scripts (with minor modifications) in order to generate a set of background and allGamma MC. Key points about these MC simulations:
- The bkg sims are on orbit using a mix of particle types...
- Other ???
Looking in gamma_ana.txt, I found the exact cuts that Eric was using. Namely, the prefilters that were being applied before the cut efficiencies were evaluated. I've listed them here for completeness.
Name |
Purpose |
Gamma Efficiency |
Bkg Efficiency |
Cut |
|---|---|---|---|---|
IrrFilter |
Remove irreducible bkg events |
|
|
((McZDir < -.2 && McId < 20) | (McZDir < .1 & McId > 20)) && McAcdZEnter > 100 && McCharge !=0 && McTHPosHitOthers < 6 |
TkrCalCore |
Take only "good" events |
|
|
TkrNumTracks>0 && CalEnergyRaw > 5 && CalCsIRLn>4 && CTBCORE>0.1 |
AcdEngFilter |
Only events with moderate signal in the ACD |
|
|
(AcdTotalEnergy+AcdRibbonEnergy<0.05) && McZDir > 0 |
Name |
Purpose |
Gamma Efficiency |
Bkg Efficiency |
Cut |
|---|---|---|---|---|
Ribbon Cut |
|
97.3% |
90.6% |
|
Corner Cut |
|
97.6% |
96.7% |
|
Basic Tile Cut |
|
90.1% |
11.2% |
|
Total Tile Energy Cut |
|
63.4% |
4.1% |
|
Veto Tile Cut |
|
90.9% |
13.4% |
|
Tile Edge Energy Cut |
|
99.2% |
88.5% |
|