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)