def addJetVariation(name, src, metName, metSrc):
addMetVariation(metName, metSrc)
for d in ["Up", "Down"]:
postfix = src % d
selectedJets = process.selectedJets.clone()
selectedJets.jetSelection.src = postfix
param.setJetPUIdSrc(selectedJets.jetSelection,
"selectedJets" + postfix)
setattr(process, "selectedJets" + postfix, selectedJets)
process.commonSequence += selectedJets
btagging = process.btagging.clone()
btagging.jetSrc = selectedJets.jetSelection.src.value()
setattr(process, "btagging" + postfix, btagging)
process.commonSequence += btagging
jets = ntuple.jets.jets.clone()
jets.detailsEnabled = False
jets.functions = cms.PSet()
jets.src.setModuleLabel("selectedJets" + postfix)
jets.floats.btagScaleFactor.setModuleLabel("btagging" + postfix)
del jets.floats.btagScaleFactorUncertainty
jets.bools.btagged.setModuleLabel("btagging" + postfix)
setattr(ntuple.jets, name % d, jets)
def addJetVariation(name, src, metName, metSrc):
addMetVariation(metName, metSrc)
for d in ["Up", "Down"]:
postfix = src%d
selectedJets = process.selectedJets.clone()
selectedJets.jetSelection.src = postfix
param.setJetPUIdSrc(selectedJets.jetSelection, "selectedJets"+postfix)
setattr(process, "selectedJets"+postfix, selectedJets)
process.commonSequence += selectedJets
btagging = process.btagging.clone()
btagging.jetSrc = selectedJets.jetSelection.src.value()
setattr(process, "btagging"+postfix, btagging)
process.commonSequence += btagging
jets = ntuple.jets.jets.clone()
jets.detailsEnabled = False
jets.functions = cms.PSet()
jets.src.setModuleLabel("selectedJets"+postfix)
jets.floats.btagScaleFactor.setModuleLabel("btagging"+postfix)
del jets.floats.btagScaleFactorUncertainty
jets.bools.btagged.setModuleLabel("btagging"+postfix)
setattr(ntuple.jets, name%d, jets)
process.preselectionSequence.insert(0, process.commonSequence)
addAnalysis(process,
"tauNtuple",
ntuple,
preSequence=process.preselectionSequence,
additionalCounters=preselectionCounters,
signalAnalysisCounters=True)
process.tauNtuple.eventCounter.printMainCounter = True
addSignalAnalysis = True
if addSignalAnalysis:
# Run signal analysis module on the same go with the embedding preselection without tau+MET trigger
import HiggsAnalysis.HeavyChHiggsToTauNu.signalAnalysis as signalAnalysis
module = signalAnalysis.createEDFilter(param)
param.setJetPUIdSrc(module.jetSelection, "")
module.Tree.fill = cms.untracked.bool(False)
module.histogramAmbientLevel = "Vital"
# Counters
if len(preselectionCounters) > 0:
module.eventCounter.counters = cms.untracked.VInputTag(
[cms.InputTag(c) for c in preselectionCounters])
def addModule(mod, postfix="", sequence=None):
for era, suffix, src in puEraSuffixWeights:
m = mod.clone()
m.pileupWeightReader.weightSrc = src
if era == dataEras[0] and suffix == "":
m.eventCounter.printMainCounter = cms.untracked.bool(True)
if suffix != "":
import FWCore.ParameterSet.Config as cms
import HiggsAnalysis.HeavyChHiggsToTauNu.HChSignalAnalysisParameters_cff as param
hPlusJetPtrSelectorFilter = cms.EDFilter("HPlusJetPtrSelectorFilter",
jetSelection = param.jetSelection.clone(),
tauSrc = cms.untracked.InputTag("patTausHpsPFTauTauTriggerMatched"),
removeTau = cms.bool(True),
allowEmptyTau = cms.bool(False),
histogramAmbientLevel = cms.untracked.string("Debug"),
filter = cms.bool(True),
throw = cms.bool(True),
producePt20 = cms.bool(False),
eventCounter = param.eventCounter.clone()
)
# Set Jet PU ID source
param.setJetPUIdSrc(hPlusJetPtrSelectorFilter.jetSelection,"hPlusJetPtrSelectorFilter")
import FWCore.ParameterSet.Config as cms
import HiggsAnalysis.HeavyChHiggsToTauNu.HChSignalAnalysisParameters_cff as param
hPlusJetPtrSelectorFilter = cms.EDFilter(
"HPlusJetPtrSelectorFilter",
jetSelection=param.jetSelection.clone(),
tauSrc=cms.untracked.InputTag("patTausHpsPFTauTauTriggerMatched"),
removeTau=cms.bool(True),
allowEmptyTau=cms.bool(False),
histogramAmbientLevel=cms.untracked.string("Debug"),
filter=cms.bool(True),
throw=cms.bool(True),
producePt20=cms.bool(False),
eventCounter=param.eventCounter.clone())
# Set Jet PU ID source
param.setJetPUIdSrc(hPlusJetPtrSelectorFilter.jetSelection,
"hPlusJetPtrSelectorFilter")
