jetSrc = cms.untracked.InputTag("selectedPatJets"),
genSrc = cms.untracked.InputTag("genParticles"),
tauPtCut = cms.untracked.double(41),
tauEtaCut = cms.untracked.double(2.1),
pileupWeightReader = param.pileupWeightReader.clone(
enabled = False
),
eventCounter = param.eventCounter.clone(),
histogramAmbientLevel = cms.untracked.string("Informative"),
)
HChTools.addAnalysis(process, "debugAnalyzer", analyzer,
preSequence=process.commonSequence,
additionalCounters=additionalCounters)
process.debugAnalyzer.eventCounter.printMainCounter = True
for era, weight in zip(dataEras, puWeights):
m = analyzer.clone()
m.pileupWeightReader.weightSrc = weight
m.pileupWeightReader.enabled = True
HChTools.addAnalysis(process, "debugAnalyzer"+era, m,
preSequence=process.commonSequence,
additionalCounters=additionalCounters)
f = open("configDumpEmbeddingDebugTauAnalysis.py", "w")
f.write(process.dumpPython())
f.close()
muonsEnabled = cms.bool(True),
muons = Ntuple.muons.clone(
src = "tightenedMuons",
functions = cms.PSet(),
)
genParticleSrc = cms.InputTag("genParticles"),
genParticleTauSrc = cms.InputTag("genTaus"),
genTTBarEnabled = cms.bool(True),
mets = cms.PSet(),
doubles = cms.PSet(),
eventCounter = param.eventCounter.clone(),
histogramAmbientLevel = cms.untracked.string("Informative"),
)
for era, weight in zip(dataEras, puWeights):
setattr(ntuple.doubles, "weightPileup_"+era, cms.InputTag(weight))
HChTools.addAnalysis(process, "tauNtuple", ntuple,
preSequence=process.commonSequence,
additionalCounters=additionalCounters)
process.tauNtuple.eventCounter.printMainCounter = True
f = open("configDumpEmbeddingDebugTauAnalysisNtuple.py", "w")
f.write(process.dumpPython())
f.close()
def addMuonIsolationAnalyses(process, prefix, prototype, commonSequence, additionalCounters, modify=_signalAnalysisSetMuon, signalAnalysisCounters=True):
import muonAnalysis
detRelIso = muonAnalysis.isolations["sumIsoRel"]
pfRelIso = muonAnalysis.isolations["pfSumIsoRel"]
isolations = [
("RelIso05", detRelIso+" < 0.05"),
("RelIso10", detRelIso+" < 0.10"),
("RelIso15", detRelIso+" < 0.15"),
# ("RelIso20", detRelIso+" < 0.20"),
("RelIso25", detRelIso+" < 0.25"),
("RelIso50", detRelIso+" < 0.50"),
("PfRelIso05", pfRelIso+" < 0.05"),
("PfRelIso10", pfRelIso+" < 0.10"),
("PfRelIso15", pfRelIso+" < 0.15"),
# ("PfRelIso20", pfRelIso+" < 0.20"),
("PfRelIso25", pfRelIso+" < 0.25"),
("PfRelIso50", pfRelIso+" < 0.50"),
("IsoTauLikeVLoose", muonAnalysis.isolations["tauVLooseIso"]+" == 0"),
("IsoTauLikeLoose", muonAnalysis.isolations["tauLooseIso"] +" == 0"),
("IsoTauLikeMedium", muonAnalysis.isolations["tauMediumIso"]+" == 0"),
("IsoTauLikeTight", muonAnalysis.isolations["tauTightIso"] +" == 0"),
("IsoTauLikeTightSc015", muonAnalysis.isolations["tauTightSc015Iso"] +" == 0" ),
("IsoTauLikeTightSc02", muonAnalysis.isolations["tauTightSc02Iso"] +" == 0" ),
("IsoTauLikeTightIc04", muonAnalysis.isolations["tauTightIc04Iso"] +" == 0"),
("IsoTauLikeTightSc015Ic04", muonAnalysis.isolations["tauTightSc015Ic04Iso"] +" == 0" ),
("IsoTauLikeTightSc02Ic04", muonAnalysis.isolations["tauTightSc02Ic04Iso"] +" == 0" ),
("IsoTauLikeTightSumPtRel10", muonAnalysis.isolations["tauTightIso"] +" < 0.1"),
("IsoTauLikeTightSumPtRel15", muonAnalysis.isolations["tauTightIso"] +" < 0.15"),
("IsoTauLikeTightSc0SumPtRel10", muonAnalysis.isolations["tauTightSc0SumPtIsoRel"] +" < 0.1"),
("IsoTauLikeTightSc0SumPtRel15", muonAnalysis.isolations["tauTightSc0SumPtIsoRel"] +" < 0.15"),
("IsoTauLikeTightSc0Ic04SumPtRel10", muonAnalysis.isolations["tauTightSc0Ic04SumPtIsoRel"] +" < 0.1"),
("IsoTauLikeTightSc0Ic04SumPtRel15", muonAnalysis.isolations["tauTightSc0Ic04SumPtIsoRel"] +" < 0.15"),
("IsoTauLikeTightSc0Ic04NoqSumPtRel10", muonAnalysis.isolations["tauTightSc0Ic04NoqSumPtIsoRel"] +" < 0.1"),
("IsoTauLikeTightSc0Ic04NoqSumPtRel10", muonAnalysis.isolations["tauTightSc0Ic04NoqSumPtIsoRel"] +" < 0.15"),
]
tauIsolations = [
# "VLoose",
# "Loose",
# "Medium",
# "Tight"
]
for name, cut in isolations:
(sequence, counters, muons) = addMuonSelection(process, name, cut)
cseq = cms.Sequence(commonSequence*sequence)
setattr(process, prefix+name+"CommonSequence", cseq)
module = prototype.clone()
modify(module, muons)
HChTools.addAnalysis(process, prefix+name, module, cseq, additionalCounters+counters, signalAnalysisCounters)
for name in tauIsolations:
(sequence, counters, muons) = addMuonTauIsolation(process, "IsoTau"+name, "by%sIsolation"%name)
cseq = cms.Sequence(commonSequence*sequence)
setattr(process, prefix+"IsoTau"+name+"CommonSequence", cseq)
module = prototype.clone()
modify(module, muons)
HChTools.addAnalysis(process, prefix+"IsoTau"+name, module, cseq, additionalCounters+counters, signalAnalysisCounters)
def addMuonIsolationAnalyses(process, prefix, prototype, commonSequence, additionalCounters, modify=_signalAnalysisSetMuon, signalAnalysisCounters=True):
import muonAnalysis
detRelIso = muonAnalysis.isolations["sumIsoRel"]
pfRelIso = muonAnalysis.isolations["pfSumIsoRel"]
isolations = [
("RelIso05", detRelIso+" < 0.05"),
("RelIso10", detRelIso+" < 0.10"),
("RelIso15", detRelIso+" < 0.15"),
# ("RelIso20", detRelIso+" < 0.20"),
("RelIso25", detRelIso+" < 0.25"),
("RelIso50", detRelIso+" < 0.50"),
("PfRelIso05", pfRelIso+" < 0.05"),
("PfRelIso10", pfRelIso+" < 0.10"),
("PfRelIso15", pfRelIso+" < 0.15"),
# ("PfRelIso20", pfRelIso+" < 0.20"),
("PfRelIso25", pfRelIso+" < 0.25"),
("PfRelIso50", pfRelIso+" < 0.50"),
("IsoTauLikeVLoose", muonAnalysis.isolations["tauVLooseIso"]+" == 0"),
("IsoTauLikeLoose", muonAnalysis.isolations["tauLooseIso"] +" == 0"),
("IsoTauLikeMedium", muonAnalysis.isolations["tauMediumIso"]+" == 0"),
("IsoTauLikeTight", muonAnalysis.isolations["tauTightIso"] +" == 0"),
("IsoTauLikeTightSc015", muonAnalysis.isolations["tauTightSc015Iso"] +" == 0" ),
("IsoTauLikeTightSc02", muonAnalysis.isolations["tauTightSc02Iso"] +" == 0" ),
("IsoTauLikeTightIc04", muonAnalysis.isolations["tauTightIc04Iso"] +" == 0"),
("IsoTauLikeTightSc015Ic04", muonAnalysis.isolations["tauTightSc015Ic04Iso"] +" == 0" ),
("IsoTauLikeTightSc02Ic04", muonAnalysis.isolations["tauTightSc02Ic04Iso"] +" == 0" ),
("IsoTauLikeTightSumPtRel10", muonAnalysis.isolations["tauTightIso"] +" < 0.1"),
("IsoTauLikeTightSumPtRel15", muonAnalysis.isolations["tauTightIso"] +" < 0.15"),
("IsoTauLikeTightSc0SumPtRel10", muonAnalysis.isolations["tauTightSc0SumPtIsoRel"] +" < 0.1"),
("IsoTauLikeTightSc0SumPtRel15", muonAnalysis.isolations["tauTightSc0SumPtIsoRel"] +" < 0.15"),
("IsoTauLikeTightSc0Ic04SumPtRel10", muonAnalysis.isolations["tauTightSc0Ic04SumPtIsoRel"] +" < 0.1"),
("IsoTauLikeTightSc0Ic04SumPtRel15", muonAnalysis.isolations["tauTightSc0Ic04SumPtIsoRel"] +" < 0.15"),
("IsoTauLikeTightSc0Ic04NoqSumPtRel10", muonAnalysis.isolations["tauTightSc0Ic04NoqSumPtIsoRel"] +" < 0.1"),
("IsoTauLikeTightSc0Ic04NoqSumPtRel10", muonAnalysis.isolations["tauTightSc0Ic04NoqSumPtIsoRel"] +" < 0.15"),
]
tauIsolations = [
# "VLoose",
# "Loose",
# "Medium",
# "Tight"
]
for name, cut in isolations:
(sequence, counters, muons) = addMuonSelection(process, name, cut)
cseq = cms.Sequence(commonSequence*sequence)
setattr(process, prefix+name+"CommonSequence", cseq)
module = prototype.clone()
modify(module, muons)
HChTools.addAnalysis(process, prefix+name, module, cseq, additionalCounters+counters, signalAnalysisCounters)
for name in tauIsolations:
(sequence, counters, muons) = addMuonTauIsolation(process, "IsoTau"+name, "by%sIsolation"%name)
cseq = cms.Sequence(commonSequence*sequence)
setattr(process, prefix+"IsoTau"+name+"CommonSequence", cseq)
module = prototype.clone()
modify(module, muons)
HChTools.addAnalysis(process, prefix+"IsoTau"+name, module, cseq, additionalCounters+counters, signalAnalysisCounters)
# Configuration
import HiggsAnalysis.HeavyChHiggsToTauNu.HChSignalAnalysisParameters_cff as param
analyzer = cms.EDAnalyzer(
"HPlusEmbeddingDebugTauAnalyzer",
jetSrc=cms.untracked.InputTag("selectedPatJets"),
genSrc=cms.untracked.InputTag("genParticles"),
tauPtCut=cms.untracked.double(41),
tauEtaCut=cms.untracked.double(2.1),
pileupWeightReader=param.pileupWeightReader.clone(enabled=False),
eventCounter=param.eventCounter.clone(),
histogramAmbientLevel=cms.untracked.string("Informative"),
)
HChTools.addAnalysis(process,
"debugAnalyzer",
analyzer,
preSequence=process.commonSequence,
additionalCounters=additionalCounters)
process.debugAnalyzer.eventCounter.printMainCounter = True
for era, weight in zip(dataEras, puWeights):
m = analyzer.clone()
m.pileupWeightReader.weightSrc = weight
m.pileupWeightReader.enabled = True
HChTools.addAnalysis(process,
"debugAnalyzer" + era,
m,
preSequence=process.commonSequence,
additionalCounters=additionalCounters)
f = open("configDumpEmbeddingDebugTauAnalysis.py", "w")
