Versions Compared

Old Version 22

changes.mady.by.user Alex Drlica-Wagner

Saved on

compared with

New Version 23

changes.mady.by.user Alex Drlica-Wagner

Saved on

Key

  • This line was added.
  • This line was removed.
  • Formatting was changed.

...

 

nSig

nBkg

Before Prefilters

~9x10^5

~4x10^6

After Prefilters

55015

48614

Pass 6 Performance

First defining some ACD variables:

Code Block

    "AcdCornerDocaENorm" : "AcdCornerDoca*(min(1000, max(30,CTBBestEnergy)))^.5/10.",
    "AcdTileEventEnergyRatio" : "AcdTkr1ActDistTileEnergy/max(1., CTBBestEnergy)*100",
    "AcdTkr1ActiveDistENorm" : "AcdTkr1ActiveDist * sqrt(min(3000., max(10, CTBBestEnergy)))/10.",
    "AcdTotalTileEventEnergyRatio" : "AcdTotalEnergy/max(1., CTBBestEnergy) * 100",

Then, After applying the prefilters (which cut out a large portion of un-reconstructable events), I finally apply the ACD background rejection cuts, omitting the CT analysis to the files that I generated with Eric's scripts. I've also listed, in parentheses, the background efficiency removing the Earth10 photon contamination. The efficiency for each cut is listed stacked with the preceding cuts.

Name

Purpose

Gamma Efficiency

Bkg Efficiency

Cut

BasicTileCut

Reject events with track pointing at struck tile.

95.2%

7.85 (7.37)%

Code Block
(Tkr1SSDVeto == 0 && AcdTkr1ActiveDist > -16 && 
AcdTkr1ActDistTileEnergy > .4)

RibbonCut

Reject events with track pointing at struck ribbon.

94.9%

7.43 (6.95)%

Code Block
(AcdRibbonActDist > -(2 +350/sqrt(max(20,CTBBestEnergy))) && 
Tkr1SSDVeto < 3 && AcdRibbonEnergy > .05 )

TotalTileEnergyCut

Reject events with excess ACD total energy.

89.3%

1.17 (1.04)%

Code Block
(AcdTotalTileEventEnergyRatio > .8 || (AcdTkr1ActiveDistENorm > -300 && 
AcdTotalTileEventEnergyRatio > max(.005, .1 - .0001*AcdTkr1ActiveDistENorm)))

CornerCut

Reject events in the corner gap of the ACD.

88.4%

1.15 (1.03)%

Code Block
((Tkr1LATEdge/1.5) ^ 2 + (AcdCornerDocaENorm - 10)^ 2 < 3800 || 
(Tkr1LATEdge < 80 && abs(AcdCornerDocaENorm-2) < 4)) && Tkr1SSDVeto < 3

TileEdgeCut

Reject events at tile edges with decreased signal.

88.3%

1.15 (1.03)%

Code Block
(abs(AcdTkr1ActiveDistENorm) < 15 && AcdTotalTileEventEnergyRatio > .005)

...

Moving on to Pass7, the ACD analysis has be modified a fair amount. I rebuilt GR v17r35p10 and generated roughly the same number of signal and background events (again using Eric allGamma.txt and background.txt scripts). Since Pass7 was just a re-design of the worksheet, the prefilters should have the same efficiency (Irreducible 13%, AcdEngCut 47%) with a slight change in CalTkrCore since it depends on CTBCore (2.2% for bkg and 6.4% for sig) leaving a total of 1.3% of the bkg and 6.4% of the gamma-ray events.

Again defining some ACD variables:

Code Block

    "AcdCornerDocaENorm" : "AcdCornerDoca*(min(1000, max(30, CTBBestEnergy)))^.5/10." ,
    "AcdTkr1RibbonActDistENorm":"AcdTkr1RibbonActDist * sqrt(min(3000., max(10, CTBBestEnergy)))/10.",
    "AcdTkr1ActiveDistENorm":"AcdTkr1ActiveDist * sqrt(min(3000., max(10, CTBBestEnergy)))/10." ,
    "Tkr1ACDTopX":"Tkr1X0 + Tkr1XDir*(755-Tkr1Z0)/Tkr1ZDir" ,
    "Tkr1ACDTopY":"Tkr1Y0 + Tkr1YDir*(755-Tkr1Z0)/Tkr1ZDir" ,
    "AcdTileEventEnergyRatio":"100*AcdTkr1ActDistTileEnergy/max(10., CTBBestEnergy)" ,
    "AcdTotalTileEventEnergyRatio":"100.*AcdTotalEnergy/max(10., CTBBestEnergy)" ,
    "Tkr1ACDSideZ":"min((Tkr1Z0 + abs(Tkr1ZDir*(840 - abs(Tkr1Y0))/Tkr1YDir)),(Tkr1Z0 + abs(Tkr1ZDir*(840 - abs(Tkr1X0))/Tkr1XDir)))" ,
    "AcdTkrVActiveDistENorm":"AcdActiveDist3D  * sqrt(min(3000., max(10, CTBBestEnergy)))/10." ,
    # Why this 1* is necessary, we may never know...
    "AcdTkr1RibbonActDistMaxTileEnergy":"1*max(AcdTkr1RibbonActEnergyPmtA, AcdTkr1RibbonActEnergyPmtB)",

Name

Purpose

Gamma Efficiency

Bkg Efficiency

Cut

BasicTileCut

Reject events with track pointing at struck tile.

94.2%

6.18%

Code Block
Tkr1SSDVeto< 5 && AcdTkr1ActDistTileEnergy > .7 && AcdTkr1ActiveDistENorm> -350

RibbonCut

Reject events with track pointing at struck ribbon.

93%

5.85%

Code Block
(AcdTkr1RibbonActDistENorm > -40  && 
Tkr1SSDVeto < 3 &&  AcdTkr1RibbonActDistMaxTileEnergy  > .04) || 
(AcdTkr1RibbonDist > -1/(CTBBestEnergy/100)  && Tkr1SSDVeto < 2)

TotalTileEnergyCut

Reject events with excess ACD total energy.

88.1%

2.08%

Code Block
AcdTotalTileEventEnergyRatio > .8 || AcdTkr1ActiveDistENorm > -200 && 
AcdTotalTileEventEnergyRatio > max(.005, .1 -  .0001*AcdTkr1ActiveDistENorm) * 
max(1., CTBBestLogEnergy/2.5)

CornerCut

Reject events in the corner gap of the ACD.

86.9%

2.06%

Code Block
((Tkr1LATEdge/1.5)^2 + (AcdCornerDocaENorm - 10)^2 < 6400  && 
Tkr1SSDVeto < 3)
|| (Tkr1LATEdge < 300 && abs(AcdCornerDocaENorm-2) < 4)

TileEdgeCut

Reject events at tile edges with decreased signal.

86.9%

2.06%

Code Block
Tkr1SSDVeto == 0 & abs(AcdTkr1ActiveDistENorm) < 10 & 
AcdTkr1ActDistTileEnergy > .025

VetoTileCut

Reject events with ...

86.6%

1.71%

Code Block
(AcdTkrVActiveDistENorm > -100 && AcdActDistTileEnergy /
sqrt(max(1., CTBBestLogEnergy-3.5)) > .9 +.15* TkrVSSDVeto) ||
 (abs(AcdTkrVActiveDistENorm) < 15 && AcdActDistTileEnergy > .25 && TkrVSSDVeto < 2)

...