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")