Skimming/python/PDWG_HZZSkim_cff.py

0001 import FWCore.ParameterSet.Config as cms
0002
0003 # cuts
0004 MUON_CUT=("pt > 7 && abs(eta)<2.5 && (isGlobalMuon || isTrackerMuon)")
0005 ELECTRON_CUT=("pt > 10 && abs(eta)<2.5")
0006 DIMUON_CUT=("mass > 40 && daughter(0).pt>20 && daughter(1).pt()>7")
0007 DIELECTRON_CUT=("mass > 40 && daughter(0).pt>20 && daughter(1).pt()>10")
0008 EMU_CUT=("mass > 40 && ((daughter(0).pt>7 && daughter(1).pt()>20) || (daughter(0).pt>20 && daughter(1).pt()>10))")
0009
0010 # single lepton selectors
0011 goodHzzMuons = cms.EDFilter("MuonRefSelector",
0012                             src = cms.InputTag("muons"),
0013                             cut = cms.string(MUON_CUT)
0014                             )
0015 goodHzzElectrons = cms.EDFilter("GsfElectronRefSelector",
0016                                 src = cms.InputTag("gedGsfElectrons"),
0017                                 cut = cms.string(ELECTRON_CUT)
0018                                 )
0019
0020 # dilepton selectors
0021 diHzzMuons = cms.EDProducer("CandViewShallowCloneCombiner",
0022                             decay       = cms.string("goodHzzMuons goodHzzMuons"),
0023                             checkCharge = cms.bool(False),
0024                             cut         = cms.string(DIMUON_CUT)
0025                             )
0026 diHzzElectrons = cms.EDProducer("CandViewShallowCloneCombiner",
0027                                 decay       = cms.string("goodHzzElectrons goodHzzElectrons"),
0028                                 checkCharge = cms.bool(False),
0029                                 cut         = cms.string(DIELECTRON_CUT)
0030                                 )
0031 crossHzzLeptons  = cms.EDProducer("CandViewShallowCloneCombiner",
0032                                   decay       = cms.string("goodHzzMuons goodHzzElectrons"),
0033                                   checkCharge = cms.bool(False),
0034                                   cut         = cms.string(EMU_CUT)
0035                                   )
0036
0037 # dilepton counters
0038 diHzzMuonsFilter = cms.EDFilter("CandViewCountFilter",
0039                                 src = cms.InputTag("diHzzMuons"),
0040                                 minNumber = cms.uint32(1)
0041                                 )
0042 diHzzElectronsFilter = cms.EDFilter("CandViewCountFilter",
0043                                     src = cms.InputTag("diHzzElectrons"),
0044                                     minNumber = cms.uint32(1)
0045                                     )
0046 crossHzzLeptonsFilter = cms.EDFilter("CandViewCountFilter",
0047                                      src = cms.InputTag("crossHzzLeptons"),
0048                                      minNumber = cms.uint32(1)
0049                                      )
0050
0051 #sequences
0052 zzdiMuonSequence = cms.Sequence( goodHzzMuons * diHzzMuons * diHzzMuonsFilter )
0053 zzdiElectronSequence = cms.Sequence( goodHzzElectrons * diHzzElectrons * diHzzElectronsFilter )
0054 zzeleMuSequence = cms.Sequence( goodHzzMuons * goodHzzElectrons * crossHzzLeptons * crossHzzLeptonsFilter )
0055