weight=total_weight) cfg_example.cut = cuts[cut_name] cfg_example.vars = variables channel = "l#nujj" if cut_name.find("_mu") >= 0: channel = "#mu#nujj" elif cut_name.find("_e") >= 0: channel = "e#nujj" if cut_name.find("HP") >= 0: channel += " HP" elif cut_name.find("LP") >= 0: channel += " LP" plots = createHistograms(cfg_example, verbose=False) for variable in variables: plot = plots[variable.name] plot.Group('Diboson', ['WWTo1L1Nu2Q', 'WWTo1L1Nu2Q', 'WZTo1L1Nu2Q']) plot.Group('Top', [ 'TT_pow_W', 'TT_pow_nonW', 'TToLeptons_tch_powheg', 'TBarToLeptons_tch_powheg', 'TToLeptons_sch', 'TBar_tWch', 'T_tWch' ]) plot.Group('VJets', [ 'WJetsToLNu_HT100to200', 'WJetsToLNu_HT200to400', 'WJetsToLNu_HT400to600', 'WJetsToLNu_HT600to800', 'WJetsToLNu_HT800to1200', 'WJetsToLNu_HT1200to2500', 'WJetsToLNu_HT2500toInf', 'DYJetsToLL_M50_HT100to200', 'DYJetsToLL_M50_HT200to400', 'DYJetsToLL_M50_HT400to600', 'DYJetsToLL_M50_HT600to800', 'DYJetsToLL_M50_HT800to1200', 'DYJetsToLL_M50_HT1200to2500', 'DYJetsToLL_M50_HT2500toInf'
# VariableCfg(name='mvis', binning={'nbinsx':35, 'xmin':0, 'xmax':350}, unit='GeV', xtitle='m_{vis}') # ] for cut_name in cuts: cfg_example = HistogramCfg(name='example', var=None, cfgs=samples, cut=inc_cut, lumi=int_lumi, weight=total_weight) cfg_example.cut = cuts[cut_name] cfg_example.vars = variables plots = createHistograms(cfg_example, verbose=True) for variable in variables: plot = plots[variable.name] plot.Group('QCD', [ 'QCD_HT2000toInf', 'QCD_HT1500to2000', 'QCD_HT1000to1500', 'QCD_HT700to1000', 'QCD_HT500to700', 'QCD_HT300to500', 'QCD_HT200to300' ]) plot.Group('TT', ['TTJets']) # plot.Group('ZLL', ['ZL', 'ZJ'], style=plot.Hist('ZL').style) # plot.Group('Electroweak', ['W', 'VV']) HistDrawer.draw(plot, plot_dir='plots/' + cut_name, channel='VV') # plot.WriteDataCard(filename='datacard_mm.root', dir='mm_' + cut_name)