Skims/python/WMuSkim_cff.py

0001 import FWCore.ParameterSet.Config as cms
0002
0003 ### HLT filter
0004 import copy
0005 from HLTrigger.HLTfilters.hltHighLevel_cfi import *
0006 WMuHLTFilter = copy.deepcopy(hltHighLevel)
0007 WMuHLTFilter.throw = cms.bool(False)
0008 WMuHLTFilter.HLTPaths = ["HLT_Mu*","HLT_IsoMu*"]
0009
0010 #Get muons of needed quality for Ws
0011 goodMuonsForW = cms.EDFilter("MuonViewRefSelector",
0012                              src = cms.InputTag("muons"),
0013                              cut = cms.string('isGlobalMuon=1 && isTrackerMuon=1 && abs(eta)<2.1 && abs(globalTrack().dxy)<0.2 && pt>20. && globalTrack().normalizedChi2<10 && globalTrack().hitPattern().numberOfValidTrackerHits>10 && globalTrack().hitPattern().numberOfValidMuonHits>0 && globalTrack().hitPattern().numberOfValidPixelHits>0 && numberOfMatches>1 && (isolationR03().sumPt+isolationR03().emEt+isolationR03().hadEt)<0.15*pt'),
0014                              filter = cms.bool(True)
0015                              )
0016
0017 # build W->MuNu candidates using PF MET
0018 wmnPFCands = cms.EDProducer("CandViewShallowCloneCombiner",
0019                             checkCharge = cms.bool(False),
0020                             cut = cms.string('sqrt((daughter(0).pt+daughter(1).pt)*(daughter(0).pt+daughter(1).pt)-pt*pt)>50'),
0021                             decay = cms.string("goodMuonsForW pfMet")
0022                             )
0023
0024 # W filter
0025 wmnPFFilter = cms.EDFilter("CandViewCountFilter",
0026                            src = cms.InputTag("wmnPFCands"),
0027                            minNumber = cms.uint32(1)
0028                            )
0029
0030 # build W->MuNu candidates using TC MET
0031 wmnTCCands = cms.EDProducer("CandViewShallowCloneCombiner",
0032                             checkCharge = cms.bool(False),
0033                             cut = cms.string('sqrt((daughter(0).pt+daughter(1).pt)*(daughter(0).pt+daughter(1).pt)-pt*pt)>50'),
0034                             decay = cms.string("goodMuonsForW tcMet")
0035                             )
0036
0037 # W filter
0038 wmnTCFilter = cms.EDFilter("CandViewCountFilter",
0039                            src = cms.InputTag("wmnTCCands"),
0040                            minNumber = cms.uint32(1)
0041                            )
0042
0043
0044 # W Skim sequences
0045 pfMetWMuNuSeq = cms.Sequence(WMuHLTFilter *
0046                              goodMuonsForW *
0047                              wmnPFCands *
0048                              wmnPFFilter
0049                              )
0050
0051
0052 tcMetWMuNuSeq = cms.Sequence(WMuHLTFilter *
0053                              goodMuonsForW *
0054                              wmnTCCands *
0055                              wmnTCFilter
0056                              )
0057
0058