Electron Identification Based on Simple Cuts (for 2010 data)

Introduction

This twiki contains information on electron selection intended to be used with 2010 data. The selection is based on cuts on a small number of variables. Different thresholds are used for electrons found in the ECAL barrel and the ECAL endcap. Electron selection variables may be categorized in 3 groups:

• e-ID variables (shower shape, track cluster matching etc)
• isolation variables
• conversion rejection variables

The sets of cuts given here are obtained by tuning all cuts together, but the sets of cuts on each of the 3 groups of variables may be used alone quite effectively. The fake rate, and fake sources, vary with ET, and for any set of cuts, rejection power and efficiency vary with ET. However, sets of cuts optimized for ET>25 GeV are near optimal for the interval 100>ET>20 GeV, and can usefully be employed down to 15 GeV. Ultimately the most performant selection should be obtained using multi-variate techniques, likelihood fits etc. Prior to that cut-based selections can provide a useful tool to understand the data and make comparison with MC. The advantages of "Simple Cuts" are:

• Cut inversion (used in many data driven signal extraction and background subtraction methodologies) is simple
• Smallest statistics are needed for full understanding and efficiency measurement
• It is simple to cleanly separate the e-ID, isolation and conversion rejection pieces

The selection has been tuned in order to get a set of cuts with maximum background rejection for a given efficiency.

Implementing the Simple Cut Based ID in CMSSW

How to Calculate the Electron ID variables in CMSSW

Most of the variables that are currently used or planned to be used can be accessed from Reco::GsfElectron or Pat::Electron.

Isolations Calculation

Isolation variables can be accessed directly from the GsfElectron using the methods:

trackIso = electron->dr03TkSumPt()
ecalIso = electron->dr03EcalRecHitSumEt()
hcalIso = electron->dr03HcalTowerSumEt()

Notice that the cone sizes in use are DR=0.3 for all three isolations. Relative isolations are the absolute isolations, described previously, divided by the electron ET:

trackIsoRel = electron->dr03TkSumPt()/electron->p4().Pt()
ecalIsoRel = electron->dr03EcalRecHitSumEt()/electron->p4().Pt()
hcalIsoRel = electron->dr03HcalTowerSumEt()/electron->p4().Pt()

Combined Isolation as is defined with the following way:

Combined Isolation for Barrel = ( electron->dr03TkSumPt() + max(0., electron->dr03EcalRecHitSumEt() - 1.) + electron->dr03HcalTowerSumEt() ) / electron->p4().Pt()
Combined Isolation for Endcaps = ( electron->dr03TkSumPt() + electron->dr03EcalRecHitSumEt() + electron->dr03HcalTowerSumEt() ) / electron->p4().Pt()

where the -1 is the pedestal subtraction and appears only in the barrel.

ID and Conversion Rejection

Shower Shape, track-cluster matching and conversion rejection variables are calculated from the following methods of the GsfElectron:

Shower shape:

electron->hadronicOverEm()
electron->sigmaIetaIeta()

Track-cluster matching:

electron->deltaPhiSuperClusterTrackAtVtx()
electron->deltaEtaSuperClusterTrackAtVtx()

Conversion rejection:

electron->gsfTrack()->trackerExpectedHitsInner().numberOfHits() (number of missing hits)
electron->convDist() (Minimum distance between conversion tracks)
electron->convDcot() (Dcot(theta) between conversion tracks at conversion vertex)

Cut values and performance

The original tuning of electron selections based on simple cuts using the followig MC samples:

• Signal: Wenu and Zee
• Backgrounds: Wtaunu, TTbar, Wmunu, Ztautau, QCD Dijets

The cut values that appear in the following table for the electron ID variables have been tuned separately and then matched to appropriate electron isolation cuts. The user should take the conversion rejection cuts (same for both barrel and endcaps), the electron ID cuts, and then choose the isolation of their preference: either combined or relative. The efficiency number that is quoted here is approximate and refers to electrons in W events passing this selection after demanding a 25 GeV cut in the ET of the electron supercluster.

Efficiencies	95%	90%	85%	80%	70%	60%
Conversion Rejection
missing hits ≤	1	1	1	0	0	0
Dist	N/A	0.02	0.02	0.02	0.02	0.02
Δcotθ	N/A	0.02	0.02	0.02	0.02	0.02
BARREL
Combined Isolation	0.15	0.10	0.09	0.07	0.04	0.03
Relative Isolation
trackRel03	0.15	0.12	0.09	0.09	0.05	0.04
ecalRel03	2.00	0.09	0.08	0.07	0.06	0.04
hcalRel03	0.12	0.10	0.10	0.10	0.03	0.03
Electron ID
σiηiη	0.01	0.01	0.01	0.01	0.01	0.01
Δφ	0.8	0.8	0.06	0.06	0.03	0.025
Δη	0.007	0.007	0.006	0.004	0.004	0.004
HoE	0.15	0.12	0.04	0.04	0.025	0.025
ENDCAPS
Combined Isolation	0.1	0.07	0.06	0.06	0.03	0.02
Relative Isolation
trackRel03	0.08	0.05	0.05	0.04	0.025	0.025
ecalRel03	0.06	0.06	0.05	0.05	0.025	0.02
hcalRel03	0.05	0.03	0.025	0.025	0.02	0.02
Electron ID
σiηiη	0.03	0.03	0.03	0.03	0.03	0.03
Δφ	0.7	0.7	0.04	0.03	0.02	0.02
Δη	0.01	0.009	0.007	0.007	0.005	0.005
HoE	0.07	0.05	0.025	0.025	0.025	0.025

Performance of the Working Points

The comparison of the performance of the working points presented in the table of the previous session are shown here

 

 

-- MatteoSani - 20 Oct 2014