python/Identification/heepElectronID_HEEPV60_cff.py

0001 from PhysicsTools.SelectorUtils.centralIDRegistry import central_id_registry
0002
0003 # Common functions and classes for ID definition are imported here:
0004 from RecoEgamma.ElectronIdentification.Identification.heepElectronID_tools import *
0005
0006 #
0007 # The HEEP ID cuts V6.0 below are optimized IDs for PHYS14 Scenario PU 20, bx 25ns
0008 # The cut values are taken from the twiki:
0009 #       https://twiki.cern.ch/twiki/bin/view/CMS/HEEPElectronIdentificationRun2#Selection_Cuts_HEEP_V5_1
0010 #       (where they may not stay, if a newer version of cuts becomes available for these
0011 #        conditions)
0012 # See also the presentation explaining these working points (this will not change):
0013 #   [ ... none available for this particular version ... ]
0014
0015
0016 # The cut values for the  Barrel and Endcap
0017 idName = "heepElectronID-HEEPV60"
0018 WP_HEEP60_EB = HEEP_WorkingPoint_V1(
0019     idName=idName,
0020     dEtaInSeedCut=0.004,
0021     dPhiInCut=0.06,
0022     full5x5SigmaIEtaIEtaCut=9999,
0023     # Two constants for the GsfEleFull5x5E2x5OverE5x5Cut
0024     minE1x5OverE5x5Cut=0.83,
0025     minE2x5OverE5x5Cut=0.94,
0026     # Three constants for the GsfEleHadronicOverEMLinearCut
0027     #     cut = constTerm if value < slopeStart
0028     #     cut = slopeTerm * (value - slopeStart) + constTerm if value >= slopeStart
0029     hOverESlopeTerm=0.05,
0030     hOverESlopeStart=0.00,
0031     hOverEConstTerm=1.00,
0032     # Three constants for the GsfEleTrkPtIsoCut:
0033     #     cut = constTerm if value < slopeStart
0034     #     cut = slopeTerm * (value - slopeStart) + constTerm if value >= slopeStart
0035     trkIsoSlopeTerm=0.00,
0036     trkIsoSlopeStart=0.00,
0037     trkIsoConstTerm=5.00,
0038     # Three constants for the GsfEleEmHadD1IsoRhoCut:
0039     #     cut = constTerm if value < slopeStart
0040     #     cut = slopeTerm * (value - slopeStart) + constTerm if value >= slopeStart
0041     # Also for the same cut, the effective area for the rho correction of the isolation
0042     ehIsoSlopeTerm=0.03,
0043     ehIsoSlopeStart=0.00,
0044     ehIsoConstTerm=2.00,
0045     effAreaForEHIso=0.28,
0046     # other cuts
0047     dxyCut=0.02,
0048     maxMissingHitsCut=1,
0049     ecalDrivenCut=1
0050     )
0051
0052 WP_HEEP60_EE = HEEP_WorkingPoint_V1(
0053      idName=idName,
0054     dEtaInSeedCut=0.006,
0055     dPhiInCut=0.06,
0056     full5x5SigmaIEtaIEtaCut=0.03,
0057     # Two constants for the GsfEleFull5x5E2x5OverE5x5Cut
0058     minE1x5OverE5x5Cut=-1.0,
0059     minE2x5OverE5x5Cut=-1.0,
0060     # Three constants for the GsfEleHadronicOverEMLinearCut
0061     #     cut = constTerm if value < slopeStart
0062     #     cut = slopeTerm * (value - slopeStart) + constTerm if value >= slopeStart
0063     hOverESlopeTerm=0.05,
0064     hOverESlopeStart=0.00,
0065     hOverEConstTerm=5,
0066     # Three constants for the GsfEleTrkPtIsoCut:
0067     #     cut = constTerm if value < slopeStart
0068     #     cut = slopeTerm * (value - slopeStart) + constTerm if value >= slopeStart
0069     trkIsoSlopeTerm=0.00,
0070     trkIsoSlopeStart=0.00,
0071     trkIsoConstTerm=5.00,
0072     # Three constants for the GsfEleEmHadD1IsoRhoCut:
0073     #     cut = constTerm if value < slopeStart
0074     #     cut = slopeTerm * (value - slopeStart) + constTerm if value >= slopeStart
0075     # Also for the same cut, the effective area for the rho correction of the isolation
0076     ehIsoSlopeTerm=0.03,
0077     ehIsoSlopeStart=50.0,
0078     ehIsoConstTerm=2.50,
0079     effAreaForEHIso=0.28,
0080     # other cuts
0081     dxyCut=0.05,
0082     maxMissingHitsCut=1,
0083     ecalDrivenCut=1
0084
0085     )
0086
0087 #
0088 # Finally, set up VID configuration for all cuts
0089 #
0090 heepElectronID_HEEPV60  = configureHEEPElectronID_V60 ( WP_HEEP60_EB, WP_HEEP60_EE )
0091
0092 #
0093 # The MD5 sum numbers below reflect the exact set of cut variables
0094 # and values above. If anything changes, one has to
0095 # 1) comment out the lines below about the registry,
0096 # 2) run "calculateMD5 <this file name> <one of the VID config names just above>
0097 # 3) update the MD5 sum strings below and uncomment the lines again.
0098 #
0099
0100 central_id_registry.register(heepElectronID_HEEPV60.idName,"df10ac7e3a9c22f63fa7936573beaafb")
0101
0102 heepElectronID_HEEPV60.isPOGApproved = cms.untracked.bool(True)