Line Code
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55
import FWCore.ParameterSet.Config as cms

# cuts
MUON_CUT=("pt > 7 && abs(eta)<2.5 && (isGlobalMuon || isTrackerMuon)")
ELECTRON_CUT=("pt > 10 && abs(eta)<2.5")
DIMUON_CUT=("mass > 40 && daughter(0).pt>20 && daughter(1).pt()>7")
DIELECTRON_CUT=("mass > 40 && daughter(0).pt>20 && daughter(1).pt()>10")
EMU_CUT=("mass > 40 && ((daughter(0).pt>7 && daughter(1).pt()>20) || (daughter(0).pt>20 && daughter(1).pt()>10))")

# single lepton selectors
goodHzzMuons = cms.EDFilter("MuonRefSelector",
                            src = cms.InputTag("muons"),
                            cut = cms.string(MUON_CUT)
                            )
goodHzzElectrons = cms.EDFilter("GsfElectronRefSelector",
                                src = cms.InputTag("gedGsfElectrons"),
                                cut = cms.string(ELECTRON_CUT)
                                )

# dilepton selectors
diHzzMuons = cms.EDProducer("CandViewShallowCloneCombiner",
                            decay       = cms.string("goodHzzMuons goodHzzMuons"),
                            checkCharge = cms.bool(False),
                            cut         = cms.string(DIMUON_CUT)
                            )
diHzzElectrons = cms.EDProducer("CandViewShallowCloneCombiner",
                                decay       = cms.string("goodHzzElectrons goodHzzElectrons"),
                                checkCharge = cms.bool(False),
                                cut         = cms.string(DIELECTRON_CUT)
                                )
crossHzzLeptons  = cms.EDProducer("CandViewShallowCloneCombiner",
                                  decay       = cms.string("goodHzzMuons goodHzzElectrons"),
                                  checkCharge = cms.bool(False),
                                  cut         = cms.string(EMU_CUT)
                                  )

# dilepton counters
diHzzMuonsFilter = cms.EDFilter("CandViewCountFilter",
                                src = cms.InputTag("diHzzMuons"),
                                minNumber = cms.uint32(1)
                                )
diHzzElectronsFilter = cms.EDFilter("CandViewCountFilter",
                                    src = cms.InputTag("diHzzElectrons"),
                                    minNumber = cms.uint32(1)
                                    )
crossHzzLeptonsFilter = cms.EDFilter("CandViewCountFilter",
                                     src = cms.InputTag("crossHzzLeptons"),
                                     minNumber = cms.uint32(1)
                                     )

#sequences
zzdiMuonSequence = cms.Sequence( goodHzzMuons * diHzzMuons * diHzzMuonsFilter )
zzdiElectronSequence = cms.Sequence( goodHzzElectrons * diHzzElectrons * diHzzElectronsFilter )
zzeleMuSequence = cms.Sequence( goodHzzMuons * goodHzzElectrons * crossHzzLeptons * crossHzzLeptonsFilter )