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)
