Configuration/python/HI_DiJetSkim_cff.py

0001 import FWCore.ParameterSet.Config as cms
0002
0003 # HLT jet trigger
0004 import HLTrigger.HLTfilters.hltHighLevel_cfi
0005 hltJetHI = HLTrigger.HLTfilters.hltHighLevel_cfi.hltHighLevel.clone()
0006 hltJetHI.HLTPaths = ["HLT_HIJet50U"]
0007 hltJetHI.throw = False
0008 hltJetHI.andOr = True
0009
0010 # selection of valid vertex
0011 primaryVertexFilterForJets = cms.EDFilter("VertexSelector",
0012     src = cms.InputTag("hiSelectedVertex"),
0013     cut = cms.string("!isFake && abs(z) <= 25 && position.Rho <= 2"),
0014     filter = cms.bool(True),   # otherwise it won't filter the events
0015     )
0016
0017 # jet energy correction (L2+L3)
0018 from JetMETCorrections.Configuration.JetCorrectorsAllAlgos_cff import ic5CaloL2RelativeCorrector, ic5CaloL3AbsoluteCorrector, ic5CaloL2L3Corrector, ic5CaloL2L3CorrectorTask
0019 icPu5CaloJetsL2L3 = cms.EDProducer('CorrectedCaloJetProducer',
0020     src = cms.InputTag('iterativeConePu5CaloJets'),
0021     correctors = cms.VInputTag('ic5CaloL2L3Corrector')
0022     )
0023
0024 # leading jet E_T filter
0025 jetEtFilter = cms.EDFilter("EtMinCaloJetCountFilter",
0026     src = cms.InputTag("icPu5CaloJetsL2L3"),
0027     etMin = cms.double(110.0),
0028     minNumber = cms.uint32(1)
0029     )
0030
0031 # Dijet requirement
0032 leadingCaloJet = cms.EDFilter( "LargestEtCaloJetSelector",
0033     src = cms.InputTag( "icPu5CaloJetsL2L3" ),
0034     filter = cms.bool( False ),
0035     maxNumber = cms.uint32( 1 )
0036     )
0037
0038 goodLeadingJet = cms.EDFilter("CaloJetSelector",
0039     src = cms.InputTag("leadingCaloJet"),
0040     cut = cms.string("et > 130")
0041     )
0042
0043 goodSecondJet = cms.EDFilter("CaloJetSelector",
0044     src = cms.InputTag("icPu5CaloJetsL2L3"),
0045     cut = cms.string("et > 50")
0046     )
0047
0048 backToBackDijets = cms.EDProducer("CandViewShallowCloneCombiner",
0049     checkCharge = cms.bool(False),
0050     cut = cms.string('abs(deltaPhi(daughter(0).phi,daughter(1).phi)) > 2.5'),
0051     decay = cms.string("goodLeadingJet goodSecondJet")
0052     )
0053
0054 dijetFilter = cms.EDFilter("CandViewCountFilter",
0055     src = cms.InputTag("backToBackDijets"),
0056     minNumber = cms.uint32(1)
0057     )
0058
0059 # dijet skim sequence
0060 diJetSkimSequence = cms.Sequence(hltJetHI
0061                                  * primaryVertexFilterForJets
0062                                  * icPu5CaloJetsL2L3
0063                                  * jetEtFilter
0064                                  * leadingCaloJet
0065                                  * goodLeadingJet
0066                                  * goodSecondJet
0067                                  * backToBackDijets
0068                                  * dijetFilter,
0069                                  ic5CaloL2L3CorrectorTask
0070                                  )