Back to home page

Project CMSSW displayed by LXR

 
 

    


File indexing completed on 2021-02-14 13:34:27

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_V70
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-HEEPV70"
0019 WP_HEEP70_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_HEEP70_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_HEEPV70  = configureHEEPElectronID_V70 (idName, WP_HEEP70_EB, WP_HEEP70_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_HEEPV70.idName,"49b6b60e9f16727f241eb34b9d345a8f")
0102 heepElectronID_HEEPV70.isPOGApproved = cms.untracked.bool(True)
0103