| 1 |
yilmaz |
1.2 |
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 |
|
|
|
| 48 |
|
|
hbheReflagNewTimeEnv = hbherechitreflaggerJETMET.clone()
|
| 49 |
|
|
hbheReflagNewTimeEnv.timingshapedcutsParameters.hbheTimingFlagBit=cms.untracked.int32(8)
|
| 50 |
|
|
|
| 51 |
|
|
# HCAL Timing
|
| 52 |
|
|
hcalTimingFilter = cms.EDFilter("HcalTimingFilter",
|
| 53 |
|
|
hbheHits = cms.untracked.InputTag("hbheReflagNewTimeEnv")
|
| 54 |
|
|
)
|
| 55 |
|
|
|
| 56 |
|
|
|
| 57 |
|
|
# hcal noise filter
|
| 58 |
|
|
from CommonTools.RecoAlgos.HBHENoiseFilter_cfi import *
|
| 59 |
|
|
HBHENoiseFilter.minRatio = cms.double(-99999.0)
|
| 60 |
|
|
HBHENoiseFilter.maxRatio = cms.double(99999.0)
|
| 61 |
|
|
HBHENoiseFilter.minZeros = cms.int32(100)
|
| 62 |
|
|
|
| 63 |
|
|
|
| 64 |
|
|
# jet skim sequence
|
| 65 |
|
|
jetSkimSequence = cms.Sequence(hltJetHI
|
| 66 |
|
|
* collisionEventSelection
|
| 67 |
|
|
* icPu5CaloJetsL2L3
|
| 68 |
|
|
* jetEtFilter
|
| 69 |
|
|
#* dijetEtFilter
|
| 70 |
|
|
* hiEcalRecHitSpikeFilter
|
| 71 |
|
|
* hbheReflagNewTimeEnv
|
| 72 |
|
|
* hcalTimingFilter
|
| 73 |
|
|
* HBHENoiseFilter
|
| 74 |
|
|
)
|
| 75 |
|
|
|
