CMS 3D CMS Logo

 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
heepElectronID_HEEPV51_cff.py
Go to the documentation of this file.
1 from PhysicsTools.SelectorUtils.centralIDRegistry import central_id_registry
2 
3 # Common functions and classes for ID definition are imported here:
5 
6 #
7 # The HEEP ID cuts V5.1 below are optimized IDs for PHYS14 Scenario PU 20, bx 25ns
8 # The cut values are taken from the twiki:
9 # https://twiki.cern.ch/twiki/bin/view/CMS/HEEPElectronIdentificationRun2#Selection_Cuts_HEEP_V5_1
10 # (where they may not stay, if a newer version of cuts becomes available for these
11 # conditions)
12 # See also the presentation explaining these working points (this will not change):
13 # [ ... none available for this particular version ... ]
14 
15 
16 # The cut values for the Barrel and Endcap
17 idName = "heepElectronID-HEEPV51"
18 WP_HEEP51_EB = HEEP_WorkingPoint_V1(
19  idName=idName,
20  dEtaInSeedCut=0.004,
21  dPhiInCut=0.06,
22  full5x5SigmaIEtaIEtaCut=9999,
23  # Two constants for the GsfEleFull5x5E2x5OverE5x5Cut
24  minE1x5OverE5x5Cut=0.83,
25  minE2x5OverE5x5Cut=0.94,
26  # Three constants for the GsfEleHadronicOverEMLinearCut
27  # cut = constTerm if value < slopeStart
28  # cut = slopeTerm * (value - slopeStart) + constTerm if value >= slopeStart
29  hOverESlopeTerm=0.05,
30  hOverESlopeStart=0.00,
31  hOverEConstTerm=2.00,
32  # Three constants for the GsfEleTrkPtIsoCut:
33  # cut = constTerm if value < slopeStart
34  # cut = slopeTerm * (value - slopeStart) + constTerm if value >= slopeStart
35  trkIsoSlopeTerm=0.00,
36  trkIsoSlopeStart=0.00,
37  trkIsoConstTerm=5.00,
38  # Three constants for the GsfEleEmHadD1IsoRhoCut:
39  # cut = constTerm if value < slopeStart
40  # cut = slopeTerm * (value - slopeStart) + constTerm if value >= slopeStart
41  # Also for the same cut, the effective area for the rho correction of the isolation
42  ehIsoSlopeTerm=0.03,
43  ehIsoSlopeStart=0.00,
44  ehIsoConstTerm=2.00,
45  effAreaForEHIso=0.28,
46  # other cuts
47  dxyCut=0.02,
48  maxMissingHitsCut=1,
49  ecalDrivenCut=1
50  )
51 
52 WP_HEEP51_EE = HEEP_WorkingPoint_V1(
53  idName=idName,
54  dEtaInSeedCut=0.006,
55  dPhiInCut=0.06,
56  full5x5SigmaIEtaIEtaCut=0.03,
57  # Two constants for the GsfEleFull5x5E2x5OverE5x5Cut
58  minE1x5OverE5x5Cut=-1.0,
59  minE2x5OverE5x5Cut=-1.0,
60  # Three constants for the GsfEleHadronicOverEMLinearCut
61  # cut = constTerm if value < slopeStart
62  # cut = slopeTerm * (value - slopeStart) + constTerm if value >= slopeStart
63  hOverESlopeTerm=0.05,
64  hOverESlopeStart=0.00,
65  hOverEConstTerm=12.5,
66  # Three constants for the GsfEleTrkPtIsoCut:
67  # cut = constTerm if value < slopeStart
68  # cut = slopeTerm * (value - slopeStart) + constTerm if value >= slopeStart
69  trkIsoSlopeTerm=0.00,
70  trkIsoSlopeStart=0.00,
71  trkIsoConstTerm=5.00,
72  # Three constants for the GsfEleEmHadD1IsoRhoCut:
73  # cut = constTerm if value < slopeStart
74  # cut = slopeTerm * (value - slopeStart) + constTerm if value >= slopeStart
75  # Also for the same cut, the effective area for the rho correction of the isolation
76  ehIsoSlopeTerm=0.03,
77  ehIsoSlopeStart=50.0,
78  ehIsoConstTerm=2.50,
79  effAreaForEHIso=0.28,
80  # other cuts
81  dxyCut=0.05,
82  maxMissingHitsCut=1,
83  ecalDrivenCut=1
84 
85  )
86 
87 #
88 # Finally, set up VID configuration for all cuts
89 #
90 heepElectronID_HEEPV51 = configureHEEPElectronID_V51 ( WP_HEEP51_EB, WP_HEEP51_EE )
91 
92 #
93 # The MD5 sum numbers below reflect the exact set of cut variables
94 # and values above. If anything changes, one has to
95 # 1) comment out the lines below about the registry,
96 # 2) run "calculateMD5 <this file name> <one of the VID config names just above>
97 # 3) update the MD5 sum strings below and uncomment the lines again.
98 #
99 
100 central_id_registry.register(heepElectronID_HEEPV51.idName,"b0e84f87acbc411de145ab2b8187ef1c")
101 
102 heepElectronID_HEEPV51.isPOGApproved = cms.untracked.bool(False)