# configure HBHEFilter according to recommendation on thi HN post
# https://hypernews.cern.ch/HyperNews/CMS/get/JetMET/1196.html
## uncomment to switch-off the isolation based filter
##HBHENoiseFilter.minIsolatedNoiseSumE = cms.double(999999.)
##HBHENoiseFilter.minNumIsolatedNoiseChannels = cms.int32(999999)
##HBHENoiseFilter.minIsolatedNoiseSumEt = cms.double(999999.)

# include this one in path to filter the event
from CommonTools.RecoAlgos.HBHENoiseFilter_cfi import *

# incelude this one to store the information without filtering the event
from CommonTools.RecoAlgos.HBHENoiseFilterResultProducer_cfi import HBHENoiseFilterResultProducer
HBHENoiseFilterRA2 = HBHENoiseFilterResultProducer.clone()

#process.patJetCorrFactorsAK5PFcorr.rho = cms.InputTag("kt6PFJetsCHS", "rho")
#for icorr in [process.patJetCorrFactorsAK5PFcorr,
#              process.patJetCorrFactorsAK7PFcorr
#              ] :
#    icorr.rho = cms.InputTag("kt6PFJets", "rho")
#    icorr.primaryVertices = cms.InputTag('goodOfflinePrimaryVertices')



#####################   Calorimeter Noise Filters ##############
#----------- hcal laser event filter 2012 -------------- #
from EventFilter.HcalRawToDigi.hcallasereventfilter2012_cfi import hcallasereventfilter2012
process.hcalLaserEventFilter2012 = hcallasereventfilter2012.clone()
#----------- HBHE noise filter ---------------------#
from CommonTools.RecoAlgos.HBHENoiseFilterResultProducer_cfi import HBHENoiseFilterResultProducer
process.HBHENoiseFilter = HBHENoiseFilterResultProducer.clone()
#----------- hcal laser event filter -------------- #
from RecoMET.METFilters.hcalLaserEventFilter_cfi import hcalLaserEventFilter
process.hcalLaserEventFilter = hcalLaserEventFilter.clone()
process.hcalLaserEventFilter.vetoByRunEventNumber=cms.untracked.bool(False)
process.hcalLaserEventFilter.vetoByHBHEOccupancy=cms.untracked.bool(True)
process.hcalLaserEventFilter.taggingMode = cms.bool(True)
#---------- ecal dead cell filter ------------#
from RecoMET.METFilters.EcalDeadCellTriggerPrimitiveFilter_cfi import EcalDeadCellTriggerPrimitiveFilter
process.EcalDeadCellTriggerPrimitiveFilter = EcalDeadCellTriggerPrimitiveFilter.clone()
process.EcalDeadCellTriggerPrimitiveFilter.tpDigiCollection = cms.InputTag("ecalTPSkimNA")
process.EcalDeadCellTriggerPrimitiveFilter.taggingMode = cms.bool(True)
#-------- The EE bad SuperCrystal filter -------#
from RecoMET.METFilters.eeBadScFilter_cfi import eeBadScFilter
process.eeBadScFilter = eeBadScFilter.clone()
process.eeBadScFilter.taggingMode = cms.bool(True)
예제 #3
0
#process.patJetCorrFactorsAK5PFcorr.primaryVertices = cms.InputTag('goodOfflinePrimaryVertices')
#process.patJetCorrFactorsAK5PFcorr.rho = cms.InputTag("kt6PFJetsCHS", "rho")
#for icorr in [process.patJetCorrFactorsAK5PFcorr,
#              process.patJetCorrFactorsAK7PFcorr
#              ] :
#    icorr.rho = cms.InputTag("kt6PFJets", "rho")
#    icorr.primaryVertices = cms.InputTag('goodOfflinePrimaryVertices')

#####################   Calorimeter Noise Filters ##############
#----------- hcal laser event filter 2012 -------------- #
from EventFilter.HcalRawToDigi.hcallasereventfilter2012_cfi import hcallasereventfilter2012
process.hcalLaserEventFilter2012 = hcallasereventfilter2012.clone()
#----------- HBHE noise filter ---------------------#
from CommonTools.RecoAlgos.HBHENoiseFilterResultProducer_cfi import HBHENoiseFilterResultProducer
process.HBHENoiseFilter = HBHENoiseFilterResultProducer.clone()
#----------- hcal laser event filter -------------- #
from RecoMET.METFilters.hcalLaserEventFilter_cfi import hcalLaserEventFilter
process.hcalLaserEventFilter = hcalLaserEventFilter.clone()
process.hcalLaserEventFilter.vetoByRunEventNumber = cms.untracked.bool(False)
process.hcalLaserEventFilter.vetoByHBHEOccupancy = cms.untracked.bool(True)
process.hcalLaserEventFilter.taggingMode = cms.bool(True)
#---------- ecal dead cell filter ------------#
from RecoMET.METFilters.EcalDeadCellTriggerPrimitiveFilter_cfi import EcalDeadCellTriggerPrimitiveFilter
process.EcalDeadCellTriggerPrimitiveFilter = EcalDeadCellTriggerPrimitiveFilter.clone(
)
process.EcalDeadCellTriggerPrimitiveFilter.tpDigiCollection = cms.InputTag(
    "ecalTPSkimNA")
process.EcalDeadCellTriggerPrimitiveFilter.taggingMode = cms.bool(True)
#-------- The EE bad SuperCrystal filter -------#
from RecoMET.METFilters.eeBadScFilter_cfi import eeBadScFilter
#
# This is a replacement for the original HBHENoiseFilter configuration.
# See https://twiki.cern.ch/twiki/bin/viewauth/CMS/HCALNoiseFilterRecipe.
# Note that this replacement relies on having the HcalNoiseSummary in the
# event record but not necessarily HBHENoiseFilterResult.
#
import FWCore.ParameterSet.Config as cms

# Module which will remake HBHENoiseFilterResult
from CommonTools.RecoAlgos.HBHENoiseFilterResultProducer_cfi import HBHENoiseFilterResultProducer
MakeHBHENoiseFilterResult = HBHENoiseFilterResultProducer.clone()

# Filter on the standard HCAL noise decision
HBHENoiseFilter = cms.EDFilter(
    'BooleanFlagFilter',
    inputLabel = cms.InputTag('MakeHBHENoiseFilterResult','HBHENoiseFilterResult'),
    reverseDecision = cms.bool(False)
)

# Customize MakeHBHENoiseFilterResult in the same manner
# as HBHENoiseFilterResultProducer
from Configuration.StandardSequences.Eras import eras
eras.run2_common.toModify(MakeHBHENoiseFilterResult, IgnoreTS4TS5ifJetInLowBVRegion=False)
eras.run2_25ns_specific.toModify(MakeHBHENoiseFilterResult, defaultDecision="HBHENoiseFilterResultRun2Loose")