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 |
|
|
AcdEngSc |
|
|
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.