| 1 |
import FWCore.ParameterSet.Config as cms
|
| 2 |
|
| 3 |
# HLT jet trigger
|
| 4 |
import HLTrigger.HLTfilters.hltHighLevel_cfi
|
| 5 |
hltJetHI = HLTrigger.HLTfilters.hltHighLevel_cfi.hltHighLevel.clone()
|
| 6 |
hltJetHI.HLTPaths = ["HLT_HIJet50U"]
|
| 7 |
hltJetHI.throw = False
|
| 8 |
hltJetHI.andOr = True
|
| 9 |
|
| 10 |
# selection of valid vertex
|
| 11 |
#primaryVertexFilterForJets = cms.EDFilter("VertexSelector",
|
| 12 |
# src = cms.InputTag("hiSelectedVertex"),
|
| 13 |
# cut = cms.string("!isFake && abs(z) <= 25 && position.Rho <= 2"),
|
| 14 |
# filter = cms.bool(True), # otherwise it won't filter the events
|
| 15 |
# )
|
| 16 |
|
| 17 |
from HeavyIonsAnalysis.Configuration.collisionEventSelection_cff import *
|
| 18 |
|
| 19 |
|
| 20 |
# jet energy correction (L2+L3) ??
|
| 21 |
from JetMETCorrections.Configuration.JetCorrectionServicesAllAlgos_cff import *
|
| 22 |
icPu5CaloJetsL2L3 = cms.EDProducer('CaloJetCorrectionProducer',
|
| 23 |
src = cms.InputTag('iterativeConePu5CaloJets'),
|
| 24 |
correctors = cms.vstring('ic5CaloL2L3')
|
| 25 |
)
|
| 26 |
|
| 27 |
# leading jet E_T filter
|
| 28 |
jetEtFilter = cms.EDFilter("EtMinCaloJetCountFilter",
|
| 29 |
src = cms.InputTag("icPu5CaloJetsL2L3"),
|
| 30 |
etMin = cms.double(60.0),
|
| 31 |
minNumber = cms.uint32(1)
|
| 32 |
)
|
| 33 |
|
| 34 |
# dijet E_T filter
|
| 35 |
#dijetEtFilter = cms.EDFilter("EtMinCaloJetCountFilter",
|
| 36 |
# src = cms.InputTag("icPu5CaloJetsL2L3"),
|
| 37 |
# etMin = cms.double(50.0),
|
| 38 |
# minNumber = cms.uint32(2)
|
| 39 |
# )
|
| 40 |
|
| 41 |
#from RecoHI.HiEgammaAlgos.hiEcalSpikeFilter_cfi import *
|
| 42 |
from CmsHi.PhotonAnalysis.hiEcalRecHitSpikeFilter_cfi import *
|
| 43 |
hiEcalRecHitSpikeFilter.minEt = 20.0
|
| 44 |
|
| 45 |
#HCAL cleaning
|
| 46 |
#from JetMETAnalysis.HcalReflagging.hbherechitreflaggerJETMET_cfi import *
|
| 47 |
# Broken.. commented out by Yen-Jie
|
| 48 |
# Need to update
|
| 49 |
|
| 50 |
#hbheReflagNewTimeEnv = hbherechitreflaggerJETMET.clone()
|
| 51 |
#hbheReflagNewTimeEnv.timingshapedcutsParameters.hbheTimingFlagBit=cms.untracked.int32(8)
|
| 52 |
|
| 53 |
# HCAL Timing
|
| 54 |
hcalTimingFilter = cms.EDFilter("HcalTimingFilter",
|
| 55 |
hbheHits = cms.untracked.InputTag("hbheReflagNewTimeEnv")
|
| 56 |
)
|
| 57 |
|
| 58 |
|
| 59 |
# hcal noise filter
|
| 60 |
from CommonTools.RecoAlgos.HBHENoiseFilter_cfi import *
|
| 61 |
HBHENoiseFilter.minRatio = cms.double(-99999.0)
|
| 62 |
HBHENoiseFilter.maxRatio = cms.double(99999.0)
|
| 63 |
HBHENoiseFilter.minZeros = cms.int32(100)
|
| 64 |
|
| 65 |
|
| 66 |
from Appeltel.RpPbAnalysis.PAPileUpVertexFilter_cff import *
|
| 67 |
|
| 68 |
|
| 69 |
# jet skim sequence
|
| 70 |
jetSkimSequence = cms.Sequence(hltJetHI
|
| 71 |
* collisionEventSelection
|
| 72 |
* icPu5CaloJetsL2L3
|
| 73 |
* jetEtFilter
|
| 74 |
#* dijetEtFilter
|
| 75 |
* hiEcalRecHitSpikeFilter
|
| 76 |
#* hbheReflagNewTimeEnv
|
| 77 |
* hcalTimingFilter
|
| 78 |
* HBHENoiseFilter
|
| 79 |
)
|
| 80 |
|
