Back to home page

Project CMSSW displayed by LXR

 
 

    


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

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)