Back to home page

Project CMSSW displayed by LXR

 
 

    


File indexing completed on 2021-02-14 12:49:25

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.JetCorrectionServicesAllAlgos_cff import *
0019 icPu5CaloJetsL2L3 = cms.EDProducer('CaloJetCorrectionProducer',
0020     src = cms.InputTag('iterativeConePu5CaloJets'),
0021     correctors = cms.vstring('ic5CaloL2L3')
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                                  )