cfgTools.addCategories(process.vbfTagDumper,
cats,
variables = all_variables,
histograms = []
)
print cats
#process.vbfTagDumper.nameTemplate ="$PROCESS_$SQRTS_$LABEL_$SUBCAT_$CLASSNAME"
process.vbfTagDumper.nameTemplate = "$PROCESS_$SQRTS_$CLASSNAME_$SUBCAT_$LABEL"
customize.setDefault("maxEvents" , -1 ) # max-number of events
customize.setDefault("targetLumi", 1.00e+3 ) # define integrated lumi
customize(process)
from HLTrigger.HLTfilters.hltHighLevel_cfi import hltHighLevel
process.hltHighLevel = hltHighLevel.clone(HLTPaths = cms.vstring("HLT_Diphoton30_18_R9Id_OR_IsoCaloId_AND_HE_R9Id_Mass95_v1") )
process.options = cms.untracked.PSet( wantSummary = cms.untracked.bool(True) )
process.hltRequirement = cms.Sequence()
if customize.processId == "Data" and requireTriggerOnData:
process.hltRequirement += process.hltHighLevel
if doSystematics:
process.p1 = cms.Path(
process.hltRequirement*
process.flashggDiPhotonSystematics*
(process.flashggUnpackedJets*process.jetSystematicsSequence)*
(process.flashggTagSequence+process.systematicsTagSequences)*
0), ## remaining EB is low R9
("EEHighR9", "min(leadingPhoton.r9,subLeadingPhoton.r9)>0.94",
0), ## then EE high R9
("EELowR9", "1", 0), ## evereything elese is EE low R9
],
## variables to be dumped in trees/datasets. Same variables for all categories
## if different variables wanted for different categories, can add categorie one by one with cfgTools.addCategory
variables=[
"CMS_hgg_mass[320,100,180]:=mass",
"leadPt :=leadingPhoton.pt",
"subleadPt :=subLeadingPhoton.pt",
"minR9 :=min(leadingPhoton.r9,subLeadingPhoton.r9)",
"maxEta :=max(abs(leadingPhoton.superCluster.eta),abs(leadingPhoton.superCluster.eta))",
"leadIDMVA :=leadingView.phoIdMvaWrtChosenVtx",
"subleadIDMVA :=subLeadingView.phoIdMvaWrtChosenVtx",
],
## histograms to be plotted.
## the variables need to be defined first
histograms=[
"CMS_hgg_mass>>mass(320,100,180)",
"subleadPt:leadPt>>ptSubVsLead(180,20,200:180,20,200)",
"minR9>>minR9(110,0,1.1)",
"maxEta>>maxEta[0.,0.1,0.2,0.3,0.4,0.6,0.8,1.0,1.2,1.4442,1.566,1.7,1.8,2.,2.2,2.3,2.5]"
])
process.p1 = cms.Path(process.kinPreselDiPhotons * process.diphotonDumper)
from flashgg.MetaData.JobConfig import customize
customize.setDefault("maxEvents", 100)
customize(process)