process.load("FWCore.MessageService.MessageLogger_cfi")
process.maxEvents = cms.untracked.PSet(input = cms.untracked.int32(-1))
process.MessageLogger.cerr.FwkReport.reportEvery = cms.untracked.int32( 1 )
process.source = cms.Source ("PoolSource",
fileNames = cms.untracked.vstring(options.inputFiles))
if options.maxEvents > 0:
process.source.eventsToProcess = cms.untracked.VEventRange('1:1-1:'+str(options.maxEvents))
process.TFileService = cms.Service("TFileService",
fileName = cms.string(options.outputFile))
from flashgg.Taggers.tagsDumpers_cfi import createTagDumper
import flashgg.Taggers.dumperConfigTools as cfgTools
process.vhEtTagDumper = createTagDumper("VHEtTag")
process.vhEtTagDumper.dumpTrees = True
process.vhEtTagDumper.dumpHistos = True
dipho_variables=["dipho_sumpt := diPhoton.sumPt",
"dipho_cosphi := abs(cos(diPhoton.leadingPhoton.phi - diPhoton.subLeadingPhoton.phi))",
"mass := diPhoton.mass",
"leadPt := diPhoton.leadingPhoton.pt",
"leadEt := diPhoton.leadingPhoton.et",
"leadEta := diPhoton.leadingPhoton.eta",
"leadPhi := diPhoton.leadingPhoton.phi",
"lead_sieie := diPhoton.leadingPhoton.sigmaIetaIeta",
"lead_hoe := diPhoton.leadingPhoton.hadronicOverEm",
"lead_sigmaEoE := diPhoton.leadingPhoton.sigEOverE",
"lead_ptoM := diPhoton.leadingPhoton.pt/diPhoton.mass",
"leadR9 := diPhoton.leadingPhoton.r9",
if systlabel == "":
continue
newseq = cloneProcessingSnippet(process,process.flashggTagSequence,systlabel)
if systlabel in jetsystlabels:
for i in range(len(jetSystematicsInputTags)):
massSearchReplaceAnyInputTag(newseq,jetSystematicsInputTags[i],cms.InputTag(jetSystematicsInputTags[i].moduleLabel,systlabel))
for name in newseq.moduleNames():
module = getattr(process,name)
if hasattr(module,"SystLabel"):
module.SystLabel = systlabel
process.systematicsTagSequences += newseq
process.flashggVBFTagMerger.src.append(cms.InputTag("flashggVBFTag" + systlabel))
# set the VBF dumper
process.vbfTagDumper = createTagDumper("VBFTag")
process.vbfTagDumper.dumpTrees = True
process.vbfTagDumper.dumpHistos = True
process.vbfTagDumper.dumpWorkspace = False
if doSystematics:
process.vbfTagDumper.src = "flashggVBFTagMerger"
else:
# use the trigger-diphoton-preselection
massSearchReplaceAnyInputTag(process.flashggTagSequence,cms.InputTag("flashggDiPhotons"),cms.InputTag("flashggPreselectedDiPhotons"))
# get the variable list
import flashgg.Taggers.VBFTagVariables as var
new_variables = [
"prompt_pho_1 := diPhoton.leadingPhoton.genMatchType()",
# "file:myMicroAODOutputFile.root"
# "root://eoscms.cern.ch//eos/cms/store/group/phys_higgs/cmshgg/sethzenz/flashgg/RunIISpring15-ReMiniAOD-1_1_0-25ns/1_1_0/VBFHToGG_M-125_13TeV_powheg_pythia8/RunIISpring15-ReMiniAOD-1_1_0-25ns-1_1_0-v0-RunIISpring15MiniAODv2-74X_mcRun2_asymptotic_v2-v1/160105_224017/0000/myMicroAODOutputFile_1.root"
# "/store/group/phys_higgs/cmshgg/szenz/flashgg/RunIISpring15-ReReco74X-Rerun-1_1_0-25ns/1_2_0/DoubleEG/RunIISpring15-ReReco74X-Rerun-1_1_0-25ns-1_2_0-v0-Run2015D-04Dec2015-v2/160117_214114/0000/myMicroAODOutputFile_10.root"
# "root://eoscms.cern.ch//eos/cms/store/group/phys_higgs/cmshgg/sethzenz/flashgg/RunIISpring15-ReMiniAOD-1_1_0-25ns/1_1_0/ttHJetToGG_M125_13TeV_amcatnloFXFX_madspin_pythia8/RunIISpring15-ReMiniAOD-1_1_0-25ns-1_1_0-v0-RunIISpring15MiniAODv2-74X_mcRun2_asymptotic_v2-v1/160105_224456/0000/myMicroAODOutputFile_2.root"
#"root://eoscms.cern.ch//eos/cms//store/group/phys_higgs/cmshgg/sethzenz/flashgg/RunIISpring15-ReMiniAOD-1_1_0-25ns/1_1_0/VHToGG_M120_13TeV_amcatnloFXFX_madspin_pythia8/RunIISpring15-ReMiniAOD-1_1_0-25ns-1_1_0-v0-RunIISpring15MiniAODv2-74X_mcRun2_asymptotic_v2-v1/160105_224138/0000/myMicroAODOutputFile_1.root"
"root://eoscms.cern.ch//eos/cms//store/group/phys_higgs/cmshgg/ferriff/flashgg/RunIIFall15DR76-1_3_0-25ns_ext1/1_3_1/DYJetsToLL_M-50_TuneCUETP8M1_13TeV-madgraphMLM-pythia8/RunIIFall15DR76-1_3_0-25ns_ext1-1_3_1-v0-RunIIFall15MiniAODv2-PU25nsData2015v1_76X_mcRun2_asymptotic_v12-v1/160127_112132/0000/myMicroAODOutputFile_156.root"
# "/store/group/phys_higgs/cmshgg/ferriff/flashgg/RunIIFall15DR76-1_3_0-25ns_ext1/1_3_1/ttHJetToGG_M120_13TeV_amcatnloFXFX_madspin_pythia8/RunIIFall15DR76-1_3_0-25ns_ext1-1_3_1-v0-RunIIFall15MiniAODv2-PU25nsData2015v1_76X_mcRun2_asymptotic_v12-v1/160127_024939/0000/myMicroAODOutputFile_1.root"
))
process.TFileService = cms.Service("TFileService",
fileName=cms.string("test.root"))
import flashgg.Taggers.dumperConfigTools as cfgTools
from flashgg.Taggers.tagsDumpers_cfi import createTagDumper
process.ZPlusJetTagDumper = createTagDumper("ZPlusJetTag")
process.ZPlusJetTagDumper.dumpTrees = True
process.ZPlusJetTagDumper.dumpHistos = True
process.ZPlusJetTagDumper.dumpWorkspace = False
if doJetSystTrees:
process.ZPlusJetTagDumper.src = cms.InputTag("flashggSystTagMerger")
# DY
process.flashggPreselectedDiPhotons.variables = cms.vstring(
'pfPhoIso03', 'trkSumPtHollowConeDR03', 'full5x5_sigmaIetaIeta',
'full5x5_r9', '1-passElectronVeto')
if doJetSystTrees:
for syst in jetsystlabels:
getattr(process, "flashggPreselectedDiPhotons%s" %
syst).variables = process.flashggPreselectedDiPhotons.variables
if customize.runOnZ != '':
process.flashggPreselectedDiPhotons.variables[-1] = "-(passElectronVeto - 1)"
if customize.runOnZ == 'double':
process.hltHighLevel.HLTPaths = ["HLT_Diphoton30_18_R9Id_OR_IsoCaloId_AND_HE_R9Id_DoublePixelSeedMatch_Mass70_v*"]
if customize.runOnZ == 'single' and customize.processType == 'data':
process.hltHighLevel.HLTPaths = ["HLT_Ele27_eta2p1_WPLoose_Gsf_v*"]
if customize.runOnZ == 'single' and customize.processType != 'data':
process.hltHighLevel.HLTPaths = ["HLT_Ele27_eta2p1_WP75_Gsf_v*"]
import flashgg.Taggers.dumperConfigTools as cfgTools
from flashgg.Taggers.tagsDumpers_cfi import createTagDumper
# ## FIXME switch to preselected diphotons
#process.flashggDiPhotonMVANew.DiPhotonTag = "flashggPreselectedDiPhotons"
process.tagDumper = createTagDumper("UntaggedTag")
process.tagDumper.src = "flashggUntagged"
#process.tagDumper.processIndex=cms.int32(1)
process.tagDumper.splitLumiWeight=cms.untracked.bool(True)
#process.tagDumper.throwOnUnclassified= False
process.tagDumper.dumpTrees = True
process.tagDumper.dumpWorkspace = False
process.tagDumper.quietRooFit = True
process.flashggDiPhotonMVA.diphotonMVAweightfile = customize.diphoxml
minmass=100.
maxmass=180.
if customize.runOnZ:
minmass=70.
-1] = "-(passElectronVeto - 1)"
if customize.runOnZ == 'double':
process.hltHighLevel.HLTPaths = [
"HLT_Diphoton30_18_R9Id_OR_IsoCaloId_AND_HE_R9Id_DoublePixelSeedMatch_Mass70_v*"
]
if customize.runOnZ == 'single' and customize.processType == 'data':
process.hltHighLevel.HLTPaths = ["HLT_Ele27_eta2p1_WPLoose_Gsf_v*"]
if customize.runOnZ == 'single' and customize.processType != 'data':
process.hltHighLevel.HLTPaths = ["HLT_Ele27_eta2p1_WP75_Gsf_v*"]
import flashgg.Taggers.dumperConfigTools as cfgTools
from flashgg.Taggers.tagsDumpers_cfi import createTagDumper
# ## FIXME switch to preselected diphotons
#process.flashggDiPhotonMVANew.DiPhotonTag = "flashggPreselectedDiPhotons"
process.tagDumper = createTagDumper("UntaggedTag")
process.tagDumper.src = "flashggUntagged"
#process.tagDumper.processIndex=cms.int32(1)
process.tagDumper.splitLumiWeight = cms.untracked.bool(True)
#process.tagDumper.throwOnUnclassified= False
process.tagDumper.dumpTrees = True
process.tagDumper.dumpWorkspace = False
process.tagDumper.quietRooFit = True
process.flashggDiPhotonMVA.diphotonMVAweightfile = customize.diphoxml
minmass = 100.
maxmass = 180.
if customize.runOnZ:
minmass = 70.
process.source = cms.Source(
"PoolSource",
fileNames=cms.untracked.vstring(
"root://eoscms.cern.ch//eos/cms/store/group/phys_higgs/cmshgg/sethzenz/flashgg/RunIISpring16DR80X-2_2_0-25ns_DYWithPDF/2_2_0/DYJetsToLL_M-50_TuneCUETP8M1_13TeV-amcatnloFXFX-pythia8/RunIISpring16DR80X-2_2_0-25ns_DYWithPDF-2_2_0-v0-RunIISpring16MiniAODv2-PUSpring16_80X_mcRun2_asymptotic_2016_miniAODv2_v0-v1/160718_161741/0000/myMicroAODOutputFile_1.root"
# "root://eoscms.cern.ch//eos/cms/store/group/phys_higgs/cmshgg/sethzenz/flashgg/RunIISpring16DR80X-2_1_2-25ns_ICHEP16/2_1_2/VBFHToGG_M-125_13TeV_powheg_pythia8/RunIISpring16DR80X-2_1_2-25ns_ICHEP16-2_1_2-v0-RunIISpring16MiniAODv1-PUSpring16RAWAODSIM_80X_mcRun2_asymptotic_2016_v3-v1/160628_143400/0000/myMicroAODOutputFile_5.root"
#"/store/group/phys_higgs/cmshgg/ferriff/flashgg/RunIISpring16DR80X-2_0_0-25ns/2_0_0/VBFHToGG_M-120_13TeV_powheg_pythia8/RunIISpring16DR80X-2_0_0-25ns-2_0_0-v0-RunIISpring16MiniAODv1-PUSpring16RAWAODSIM_80X_mcRun2_asymptotic_2016_v3-v1/160524_093617/0000/myMicroAODOutputFile_1.root"
#"root://eoscms.cern.ch//eos/cms/store/group/phys_higgs/cmshgg/ferriff/flashgg/RunIISpring16DR80X-2_0_0-25ns/2_0_0/DYToEE_NNPDF30_13TeV-powheg-pythia8/RunIISpring16DR80X-2_0_0-25ns-2_0_0-v0-RunIISpring16MiniAODv1-PUSpring16_80X_mcRun2_asymptotic_2016_v3-v1/160524_084452/0000/myMicroAODOutputFile_110.root"
))
process.TFileService = cms.Service("TFileService",
fileName=cms.string("test.root"))
import flashgg.Taggers.dumperConfigTools as cfgTools
from flashgg.Taggers.tagsDumpers_cfi import createTagDumper
process.vbfTagDumper = createTagDumper("VBFTag")
process.vbfTagDumper.dumpTrees = True
process.vbfTagDumper.dumpHistos = True
process.vbfTagDumper.dumpWorkspace = False
if dumpJetSysTrees:
process.vbfTagDumper.src = "flashggSystTagMerger"
# use custum TMVA weights
# process.flashggVBFMVA.vbfMVAweightfile = cms.FileInPath("flashgg/Taggers/data/TMVAClassification_dijetMVA_76x_17_02_15_BDTG.weights.xml")
# process.flashggVBFMVA.MVAMethod = cms.string("BDTG")
# process.flashggDiPhotonMVA.diphotonMVAweightfile = cms.FileInPath("flashgg/Taggers/data/TMVAClassification_BDT_QCDeroded_v100_rereco.weights.xml")
#process.flashggVBFMVA.vbfMVAweightfile = cms.FileInPath("flashgg/Taggers/data/TMVA_classification_dijet-mva-80x-ICHEP-12-07-2016-v03_BDTG.weights.xml")
#process.flashggVBFMVA.vbfDiPhoDiJetMVAweightfile = cms.FileInPath("flashgg/Taggers/data/TMVA_classification_combined-mva-80x-ICHEP-v02_BDTG.weights.xml")
# QCD Recovery
if switchToPuppi:
process.flashggUnpackedJets.JetsTag = cms.InputTag("flashggFinalPuppiJets")
from flashgg.Taggers.flashggTagOutputCommands_cff import tagDefaultOutputCommand
#process.out = cms.OutputModule("PoolOutputModule", fileName = cms.untracked.string('/tmp/nancy/vhEtTagOutputFile_Data.root'),
# outputCommands = tagDefaultOutputCommand
# )
process.TFileService = cms.Service("TFileService",
fileName=cms.string("test.root"))
from flashgg.Taggers.tagsDumpers_cfi import createTagDumper
import flashgg.Taggers.dumperConfigTools as cfgTools
process.vhEtTagDumper = createTagDumper("VHEtTag")
process.vhEtTagDumper.dumpTrees = True
process.vhEtTagDumper.dumpHistos = True
dipho_variables = [
"dipho_sumpt := diPhoton.sumPt",
"dipho_pt := diPhoton.pt",
"dipho_cosphi := cos(diPhoton.leadingPhoton.phi - diPhoton.subLeadingPhoton.phi)",
"mass := diPhoton.mass",
"phi := diPhoton.momentum.phi",
"leadPt := diPhoton.leadingPhoton.pt",
"leadEt := diPhoton.leadingPhoton.et",
"leadEta := diPhoton.leadingPhoton.eta",
"leadPhi := diPhoton.leadingPhoton.phi",
"lead_sieie := diPhoton.leadingPhoton.full5x5_sigmaIetaIeta",
"lead_oldhoe := diPhoton.leadingPhoton.hadronicOverEm",
from flashgg.MetaData.samples_utils import SamplesManager
process.source = cms.Source ("PoolSource",
fileNames = cms.untracked.vstring(
#"/store/group/phys_higgs/cmshgg/sethzenz/flashgg/RunIISummer16-2_4_1-25ns_Moriond17/2_4_1/DYJetsToLL_M-50_TuneCUETP8M1_13TeV-amcatnloFXFX-pythia8/RunIISummer16-2_4_1-25ns_Moriond17-2_4_1-v0-RunIISummer16MiniAODv2-PUMoriond17_HCALDebug_80X_mcRun2_asymptotic_2016_TrancheIV_v6-v1/170114_082421/0000/myMicroAODOutputFile_9.root"
"/store/group/phys_higgs/cmshgg/ferriff/flashgg/RunIISpring16DR80X-2_3_0-25ns_Moriond17_MiniAODv2/2_3_0/DoubleEG/RunIISpring16DR80X-2_3_0-25ns_Moriond17_MiniAODv2-2_3_0-v0-Run2016D-23Sep2016-v1/161117_082833/0000/myMicroAODOutputFile_99.root"
)
)
process.TFileService = cms.Service("TFileService",
fileName = cms.string("test.root"))
import flashgg.Taggers.dumperConfigTools as cfgTools
from flashgg.Taggers.tagsDumpers_cfi import createTagDumper
process.ZPlusJetTagDumper = createTagDumper("ZPlusJetTag")
process.ZPlusJetTagDumper.dumpTrees = True
process.ZPlusJetTagDumper.dumpHistos = True
process.ZPlusJetTagDumper.dumpWorkspace = False
if customize.doJetSystTrees:
process.ZPlusJetTagDumper.src = cms.InputTag("flashggSystTagMerger")
# DY
process.flashggPreselectedDiPhotons.variables = cms.vstring('pfPhoIso03',
'trkSumPtHollowConeDR03',
'full5x5_sigmaIetaIeta',
'full5x5_r9',
'1-passElectronVeto')
if customize.doJetSystTrees:
for syst in jetsystlabels_:
# ee bad supercluster filter on data
process.load('RecoMET.METFilters.eeBadScFilter_cfi')
process.eeBadScFilter.EERecHitSource = cms.InputTag("reducedEgamma","reducedEERecHits") # Saved MicroAOD Collection (data only)
process.dataRequirements = cms.Sequence()
if customize.processId == "Data":
process.dataRequirements += process.hltHighLevel
process.dataRequirements += process.eeBadScFilter
#process.load("flashgg.MicroAOD.flashggMets_cfi")
process.load("flashgg/Taggers/flashggTagTester_cfi")
from flashgg.Taggers.tagsDumpers_cfi import createTagDumper
import flashgg.Taggers.dumperConfigTools as cfgTools
process.TTHDiLeptonTagDumper = createTagDumper("TTHDiLeptonTag")
#process.load("flashgg.Taggers.tthDumper_cfi")
#process.flashggMuMuGamma.PhotonTag=cms.InputTag('flashggUpdatedIdMVAPhotons')
process.TTHDiLeptonTagDumper.dumpTrees = True
process.TTHDiLeptonTagDumper.dumpHistos = False
process.TTHDiLeptonTagDumper.dumpWorkspace = False
process.TTHDiLeptonTagDumper.nameTemplate = cms.untracked.string("$PROCESS_$SQRTS_$CLASSNAME_$SUBCAT_$LABEL")
#from PhysicsTools.PatAlgos.tools.helpers import massSearchReplaceAnyInputTag
#massSearchReplaceAnyInputTag(process.flashggTagSequence,cms.InputTag("flashggDiPhotons"),cms.InputTag("flashggPreselectedDiPhotons"))
import flashgg.Taggers.ttHTagVariables as var
if customize.processId == "Data":
variables_ = var.generic_variables + var.dipho_variables
else:
variables_ = var.generic_variables + var.dipho_variables
