METProducers/python/METSignificanceObjects_cfi.py

0001 import FWCore.ParameterSet.Config as cms
0002
0003 ##======================================
0004 ## Muons
0005 ##======================================
0006
0007 selectedMuons = cms.EDFilter(
0008       "MuonSelector",
0009       src = cms.InputTag('muons'),
0010       cut = cms.string(    "(isTrackerMuon) && pt > 5"
0011          "&& isPFMuon"+
0012          "&& globalTrack.isNonnull"+
0013          "&& innerTrack.hitPattern.numberOfValidPixelHits > 0"+
0014          "&& innerTrack.normalizedChi2 < 10"+
0015          "&& numberOfMatches > 0"+
0016          "&& innerTrack.hitPattern.numberOfValidTrackerHits>5"+
0017          "&& globalTrack.hitPattern.numberOfValidHits>0"+
0018          "&& (pfIsolationR03.sumChargedHadronPt+max(0.,pfIsolationR03.sumNeutralHadronEt+pfIsolationR03.sumPhotonEt - 0.5*pfIsolationR03.sumPUPt))/pt < 0.3"+
0019          "&& abs(innerTrack().dxy)<2.0"
0020          ),
0021       filter = cms.bool(False)
0022       )
0023
0024
0025 ##======================================
0026 ## Electrons
0027 ##======================================
0028
0029 selectedElectrons = cms.EDFilter(
0030       "GsfElectronSelector",
0031       src = cms.InputTag('gedGsfElectrons'),
0032       cut = cms.string(
0033          "abs(eta) < 2.5 && pt > 19.5"                              +
0034          "&& (gsfTrack.hitPattern().numberOfLostHits(\'MISSING_INNER_HITS\')<=1 )" +
0035          "&& (pfIsolationVariables.sumChargedHadronPt+max(0.,pfIsolationVariables.sumNeutralHadronEt+pfIsolationVariables.sumPhotonEt - 0.5*pfIsolationVariables.sumPUPt))/et     < 0.3"  +
0036          "&& ((abs(eta) < 1.4442  "                                 +
0037          "&& abs(deltaEtaSuperClusterTrackAtVtx)            < 0.007"+
0038          "&& abs(deltaPhiSuperClusterTrackAtVtx)            < 0.8"  +
0039          "&& sigmaIetaIeta                                  < 0.01" +
0040          "&& hcalOverEcal                                   < 0.15" +
0041          "&& abs(1./superCluster.energy - 1./p)             < 0.05)"+
0042          "|| (abs(eta)  > 1.566 "+
0043          "&& abs(deltaEtaSuperClusterTrackAtVtx)            < 0.009"+
0044          "&& abs(deltaPhiSuperClusterTrackAtVtx)            < 0.10" +
0045          "&& sigmaIetaIeta                                  < 0.03" +
0046          "&& hcalOverEcal                                   < 0.10" +
0047          "&& abs(1./superCluster.energy - 1./p)             < 0.05))"
0048          ),
0049       filter = cms.bool(False)
0050       )
0051
0052
0053 ##======================================
0054 ## Photons
0055 ##======================================
0056
0057 selectedPhotons = cms.EDFilter("PhotonSelector",
0058       src = cms.InputTag("photons"),
0059       cut = cms.string(
0060          "abs(eta) < 2.5 && pt > 19.5" +
0061          "&& sigmaIetaIeta < 0.03" +
0062          "&& hadronicOverEm < 0.05" +
0063          "&& hasPixelSeed == 0" +
0064          "&& (chargedHadronIso + neutralHadronIso + photonIso)/pt < 0.2"
0065          )
0066       )
0067
0068
0069 selectionTaskForMETSig = cms.Task(
0070       selectedMuons,
0071       selectedElectrons,
0072       selectedPhotons
0073       )
0074 selectionSequenceForMETSig = cms.Sequence(selectionTaskForMETSig)
0075