![]() |
|
|||
File indexing completed on 2024-04-06 12:25:07
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 V5.1 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-HEEPV51" 0018 WP_HEEP51_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=2.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_HEEP51_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=12.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_HEEPV51 = configureHEEPElectronID_V51 ( WP_HEEP51_EB, WP_HEEP51_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_HEEPV51.idName,"b0e84f87acbc411de145ab2b8187ef1c") 0101 0102 heepElectronID_HEEPV51.isPOGApproved = cms.untracked.bool(False)
[ Source navigation ] | [ Diff markup ] | [ Identifier search ] | [ general search ] |
This page was automatically generated by the 2.2.1 LXR engine. The LXR team |
![]() ![]() |