CMS 3D CMS Logo

/data/doxygen/doxygen-1.7.3/gen/CMSSW_4_2_8/src/DPGAnalysis/Skims/python/WMuSkim_cff.py

Go to the documentation of this file.
00001 import FWCore.ParameterSet.Config as cms
00002 
00003 ### HLT filter
00004 import copy
00005 from HLTrigger.HLTfilters.hltHighLevel_cfi import *
00006 WMuHLTFilter = copy.deepcopy(hltHighLevel)
00007 WMuHLTFilter.throw = cms.bool(False)
00008 WMuHLTFilter.HLTPaths = ["HLT_Mu*","HLT_IsoMu*"]
00009 
00010 #Get muons of needed quality for Ws
00011 goodMuonsForW = cms.EDFilter("MuonViewRefSelector",
00012                              src = cms.InputTag("muons"),
00013                              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'),
00014                              filter = cms.bool(True)
00015                              )
00016 
00017 # build W->MuNu candidates using PF MET
00018 wmnPFCands = cms.EDProducer("CandViewShallowCloneCombiner",
00019                             checkCharge = cms.bool(False),
00020                             cut = cms.string('sqrt((daughter(0).pt+daughter(1).pt)*(daughter(0).pt+daughter(1).pt)-pt*pt)>50'),
00021                             decay = cms.string("goodMuonsForW pfMet")
00022                             )
00023 
00024 # W filter
00025 wmnPFFilter = cms.EDFilter("CandViewCountFilter",
00026                            src = cms.InputTag("wmnPFCands"),
00027                            minNumber = cms.uint32(1)
00028                            )
00029 
00030 # build W->MuNu candidates using TC MET
00031 wmnTCCands = cms.EDProducer("CandViewShallowCloneCombiner",
00032                             checkCharge = cms.bool(False),
00033                             cut = cms.string('sqrt((daughter(0).pt+daughter(1).pt)*(daughter(0).pt+daughter(1).pt)-pt*pt)>50'),
00034                             decay = cms.string("goodMuonsForW tcMet")
00035                             )
00036 
00037 # W filter
00038 wmnTCFilter = cms.EDFilter("CandViewCountFilter",
00039                            src = cms.InputTag("wmnTCCands"),
00040                            minNumber = cms.uint32(1)
00041                            )
00042 
00043 
00044 # W Skim sequences
00045 pfMetWMuNuSeq = cms.Sequence(WMuHLTFilter *
00046                              goodMuonsForW *
00047                              wmnPFCands *
00048                              wmnPFFilter
00049                              )
00050 
00051 
00052 tcMetWMuNuSeq = cms.Sequence(WMuHLTFilter *
00053                              goodMuonsForW *
00054                              wmnTCCands *
00055                              wmnTCFilter
00056                              )
00057 
00058