Using new ACD merit variables in background rejection

 Current background rejection

Prefilters:

Cut

Purpose

Gamma Effic

Bkg Effic

Expression

BasicTileCut

Reject events with track pointing at struck tile:

95.1%

8.1%

Tkr1SSDVeto == 0 && AcdTkr1ActiveDist > -16 && AcdTkr1ActDistTileEnergy > .4

RibbonCut

Reject events with track pointing at struck ribbon

94.9%

7.7%

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

TotalTileEnergyCut

Reject events with excess ACD total energy

89.3%

1.2%

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

CornerCut

Reject event in the corner gap of the ACD

88.3%

1.1%

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

TileEdgeCut:

Reject events at tile edges with decreased signal

88.1%

1.1%

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

ClassificationTree:

Variable

Importance

AcdTotalTileEventEnergyRatio

2034.17

Tkr1SSDVeto

1201.47

CTBCORE

812.28

AcdTkr1ActDistNorm

575.06

AcdTkr1ActiveDistENorm

476.57

AcdTotalEnergyNorm

269.89

Tkr1ACDSideZ

256.61

AcdActDistNorm

118.93

Tkr1LATEdge

107.11

AcdTileEnergyNorm

106.51

AcdCornerDocaENorm

81.52

AcdMaxTileEnergy

74.32

AcdTileEventEnergyRatio

73.48

AcdTkr1TileEnergyNorm

25.67

AcdTileCount

24


Class

Gamma Eff

Bkg Eff

Expression

transient

83.5%

0.45%

Prefilter && CPFGamProb > 0.2

source

82.5%

0.40%

transient && CPFGamProbCPFGamProb > max(.05, .55- .3* (CTBBestLogEnergy - 1.))

diffuse

81.2%

0.36%

source && AcdDiffuseVeto


 Using New Variables

 AcdTkr1VetoSigmaHit and VetoHit1

AcdTkr1VetoSigmaHit is an estimator of the "best" hit vetoing the track, expressed in "sigma".   This is the sum in quadrature of contributions from the miss distance between the track and the tile or ribbon and the signal size being less that 1 MIP.  For MIP-like signals along the track projection both terms are 0, so this variable is 0.  If the track missed the tile  this variable is postive, likewise if the signal is less than 1 MIP.

 
Same plot, but with the X-axis remapped with log10( AcdTkr1VetoSigmaHit + 0.01)

Note the small peak in the signal distribution around 7-9 is from low energy deposits in the tile pointed to by the track.
For completeness here is a plot of  AcdTkr1VetoSigmaHit v. AcdTkr1ActiveDist


 
VetoHit1 is the same variable, including the tracker SSDVeto. 
VetoHit1 = sqrt(Acd2Tkr1VetoSigmaHit*Acd2Tkr1VetoSigmaHit+1.5*Tkr1SSDVeto*Tkr1SSDVeto)
Here is a plot of VetoHit1 

Same plot, but the X-axis remapped with log10(VetoHit1+0.01)


 
Again, note the small peak in the signal distribution around 7-9 is from low energy deposits in the tile pointed to by the track. 
 

  AcdTkr1VetoSigmaGap and VetoGap1

 AcdTkr1VetoSigmaGap is an estimator of the track passing through a gap in the ACD, expressed in "sigma".   This is the sum in quadrature of two contributions: the geometrical       calculations of the track going in the gap and the chance of seeing no signal from the gap.  For high energy tracks that pass into the corner gaps of the ACD it is 0.  For tracks going into ribbons it is 1.5 (which represents the fact that we do expect some signal from the ribbons, albiet a small one).

 Here is a plot of the low end of AcdTkr1VetoSigmaGap for all events:

 For comparison, here is a plot of AcdTkr1VetoSigmaGap v. AcdTkr1RibbonDist:

In analogy with VetoHit1 above we can define VetoGap1 to take into account the Trk1SSDVeto:
   VetoGap1 = sqrt(Acd2Tkr1VetoSigmaGap*Acd2Tkr1VetoSigmaGap+1.5*Tkr1SSDVeto*Tkr1SSDVeto)
 
Here is a plot of the low end of VetoGap1 for all events. 

If we look at only those events that pass weren't vetoed because of a track tile association or b/c of excess energy in the ACD we see this:
 

For Comparison, here is a plot of VetoGap1 v. AcdTkr1RibbonDist for the same set of events.


 

AcdTkrVetoSigmaHit and VetoHit

We can also select the best association by looping over all tracks instead of just using the best track.
Here is a plot of that variable, AcdTkrVetoSigmaHit: 


 
In analogy to VetoHit1 we can make VetoHit from  AcdTkrVetoSigmaHit. 
Here is a plot of that:
 

 
Then if we only look at events that weren't rejected because of VetoHit1 or excess energy:

 
So, we could consider cutting very loose on this and rejecting some event with a high probablity match.
 

Reproducing Existing Performance with simple cuts 

Cut

Reason

Gamma Effic

Bkg Effic

VetoHit1 < 5

Reject events with best rack pointing at struck tile/ribbon

94.8%

7.2%

TotalTileEnergyCut

Reject events with excess ACD total energy

89.3%

1.2%

VetoGap1 < 2

Reject event in gaps of the ACD

89.1%

1.2%

VetoHit < 2

Reject events if other track points at hit tile

89.0%

1.1%


Some other useful new Variables, (for Classification Tree)

 Another type of variable we consider is the scaled energy in a cone around the track.  We have three variables that reflect 15, 30 and 45 degreee cones.

Variable

AllEvents

Prefiltered Events

VetoHit1 and
VetoHit1 v. AcdEngSc

AcdEngSc
(aka  AcdTotalTileEventEnergyRatio)

AcdEng45Sc


AcdEng30Sc

AcdEng15Sc

With these variables (and some familiar ones) we can make a Classification Tree.

Variable

Importance

CTBCORE

0.1760

AcdEngSc45

0.1614

AcdEngSc15

0.1541

AcdTotalEnergy

0.1113

AcdEngSc (aka AcdTotalTileEventEnergyRatio)

0.08972

VetoHit1

0.07846

AcdEngSc30

0.06908

Tkr1SSDVeto

0.05000

VetoGap

0.04099

AcdRibbonEnergy

0.02754

VetoGap1

0.02169

VetoHit

0.01968

This plot show the performance of the CT (in red) relative to the pass 6 cuts (3 black circles).  The prefilter performances are also shown.
For comparison the CT performance without any prefilters is also show (in blue). 
 
Just to give some numbers, here are numbers for 3 different values of cutting on the CT output. 

Class

Gamma Effic

Bkg Effic

Expression

transient

86.4%

0.43%

prefilter && MVA_BDT > 0.38

source

84.1%

0.33%

prefilter && MVA_BDT > 0.43

diffuse

81.2%

0.27%

prefilter && MVA_BDT  > 0.48

 Basically, these numbers correspond to either 15% reduction in residual background at given signal efficiency or 3-4% increase in signal efficiency at given background.


  • No labels