python/Identification/heepElectronID_HEEPV60_80XAOD_cff.py

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