cfg_example.cut = cuts[cut_name]
cfg_example.vars = variables
channel = "l#nujj"
if cfg_example.cut.find("mu_"):
channel = "#mu#nujj"
elif cfg_example.cut.find("e_"):
channel = "e#nujj"
if cfg_example.cut.find("HP"):
channel += " HP"
elif cfg_example.cut.find("LP"):
channel += " LP"
plots = createHistograms(cfg_example, verbose=False)
for variable in variables:
plot = plots[variable.name]
plot.Group('Diboson', ['WWTo1L1Nu2Q', 'WZTo1L1Nu2Q'])
# plot.Group('TT_W', ['TTJets_W'])
plot.Group('Top', ['TT_pow_ext3_W', 'TT_pow_ext3_nonW', 'TToLeptons_tch_powheg', 'TBarToLeptons_tch_powheg', 'TToLeptons_sch', 'TBar_tWch', 'T_tWch'])
# plot.Group('W', ['WJetsToLNu_HT100to200', 'WJetsToLNu_HT200to400', 'WJetsToLNu_HT400to600', 'WJetsToLNu_HT600to800', 'WJetsToLNu_HT800to1200', 'WJetsToLNu_HT1200to2500', 'WJetsToLNu_HT2500toInf'])
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_HT600toInf'])
plot.Group('QCD', ['QCD_HT2000toInf', 'QCD_HT1500to2000', 'QCD_HT1000to1500', 'QCD_HT700to1000', 'QCD_HT500to700']) # , 'QCD_HT300to500'
# plot.Group('TT_nonW', ['TTJets_nonW'])
# plot.Group('Single t', ['TToLeptons_tch_powheg', 'TBarToLeptons_tch_powheg', 'TToLeptons_sch', 'TBar_tWch', 'T_tWch'])
# plot.Group('ZLL', ['DYJetsToLL_M50_HT100to200', 'DYJetsToLL_M50_HT200to400', 'DYJetsToLL_M50_HT400to600', 'DYJetsToLL_M50_HT600toInf'])
plot.Group('data_obs', ['data_SingleMuon', 'data_SingleElectron', 'data_MET']) #, 'data_JetHT'
#['WpWpJJ', 'ZGTo2LG', 'ZGJets', 'WGToLNuG', 'WGJets', 'WW', 'WWDouble', 'WWTo1L1Nu2Q', 'WWToLNuQQ_ext', 'WWToLNuQQ', 'WWTo2L2Nu', 'WZ', 'WZTo3LNu_amcatnlo', 'WZTo3LNu', 'WZTo2L2Q', 'WZTo1L1Nu2Q', 'WZTo1L3Nu', 'ZZ', 'VVTo2L2Nu', 'ZZTo4L', 'ZZTo2Q2Nu', 'ZZTo2L2Q', 'ZZTo2L2Nu'])
HistDrawer.draw(plot, plot_dir='plots_lnujj/'+cut_name, channel=channel)
# plot.WriteDataCard(filename='datacard_mm.root', dir='mm_' + cut_name)
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'
])
plot.Group('QCD', [
'QCD_HT2000toInf', 'QCD_HT1500to2000', 'QCD_HT1000to1500',
'QCD_HT700to1000', 'QCD_HT500to700'
]) # , 'QCD_HT300to500'
plot.Group('data_obs',
['data_SingleMuon', 'data_SingleElectron', 'data_MET'
]) #, 'data_JetHT'
#['WpWpJJ', 'ZGTo2LG', 'ZGJets', 'WGToLNuG', 'WGJets', 'WW', 'WWDouble', 'WWTo1L1Nu2Q', 'WWToLNuQQ_ext', 'WWToLNuQQ', 'WWTo2L2Nu', 'WZ', 'WZTo3LNu_amcatnlo', 'WZTo3LNu', 'WZTo2L2Q', 'WZTo1L1Nu2Q', 'WZTo1L3Nu', 'ZZ', 'VVTo2L2Nu', 'ZZTo4L', 'ZZTo2Q2Nu', 'ZZTo2L2Q', 'ZZTo2L2Nu'])
HistDrawer.draw(plot,
plot_dir='plots_lnujj_VJetsReweight/' + cut_name,
channel=channel)
# plot.WriteDataCard(filename='datacard_mm.root', dir='mm_' + cut_name)
weight=total_weight)
cfg_example.cut = cuts[cut_name]
print cfg_example.cut
cfg_example.vars = variables
channel = ""
# channel = "jj"
# 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', ['TTHad_pow'])
plot.Group('WJets', ['WJetsToQQ_HT800toInf'])
plot.Group('ZJets', ['ZJetsToQQ_HT800toInf'])
plot.Group('QCD', [
'QCD_HT100to200', 'QCD_HT200to300', 'QCD_HT300to500',
'QCD_HT500to700', 'QCD_HT700to1000', 'QCD_HT1000to1500',
'QCD_HT1500to2000', 'QCD_HT2000toInf'
])
plot.Group('data_obs', ['data_JetHT_Run2017F'])
HistDrawer.draw(plot,
plot_dir='plots_jj_Run2017F/' + cut_name,
channel=channel)
# plot.WriteDataCard(filename='datacard_mm.root', dir='mm_' + cut_name)
# 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)
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'
])
plot.Group('QCD', [
'QCD_HT2000toInf', 'QCD_HT1500to2000', 'QCD_HT1000to1500',
'QCD_HT700to1000', 'QCD_HT500to700'
]) # , 'QCD_HT300to500'
plot.Group('data_obs',
['data_SingleMuon', 'data_SingleElectron', 'data_MET'
]) #, 'data_JetHT'
#['WpWpJJ', 'ZGTo2LG', 'ZGJets', 'WGToLNuG', 'WGJets', 'WW', 'WWDouble', 'WWTo1L1Nu2Q', 'WWToLNuQQ_ext', 'WWToLNuQQ', 'WWTo2L2Nu', 'WZ', 'WZTo3LNu_amcatnlo', 'WZTo3LNu', 'WZTo2L2Q', 'WZTo1L1Nu2Q', 'WZTo1L3Nu', 'ZZ', 'VVTo2L2Nu', 'ZZTo4L', 'ZZTo2Q2Nu', 'ZZTo2L2Q', 'ZZTo2L2Nu'])
HistDrawer.draw(plot,
plot_dir='plots_lnujj_vbf_withSignal/' + cut_name,
channel=channel)
# plot.WriteDataCard(filename='datacard_mm.root', dir='mm_' + cut_name)
cfg_example = HistogramCfg(name='example',
var=None,
cfgs=samples,
cut='',
lumi=int_lumi,
weight=total_weight)
cfg_example.cut = cuts[cut_name]
cfg_example.vars = variables
plots = createHistograms(cfg_example, verbose=False)
for variable in variables:
plot = plots[variable.name]
# plot.Group('Diboson', ['WWTo1L1Nu2Q', 'WZTo1L1Nu2Q'])
# plot.Group('TT_W', ['TTJets_W'])
# plot.Group('W', ['WJetsToLNu_HT100to200', 'WJetsToLNu_HT200to400', 'WJetsToLNu_HT400to600', 'WJetsToLNu_HT600to800', 'WJetsToLNu_HT800to1200', 'WJetsToLNu_HT1200to2500', 'WJetsToLNu_HT2500toInf'])
plot.Group('QCD', [
'QCD_HT2000toInf', 'QCD_HT1500to2000', 'QCD_HT1000to1500',
'QCD_HT700to1000', 'QCD_HT500to700', 'QCD_HT300to500'
])
# plot.Group('TT_nonW', ['TTJets_nonW'])
# plot.Group('Single t', ['TToLeptons_tch_powheg', 'TBarToLeptons_tch_powheg', 'TToLeptons_sch', 'TBar_tWch', 'T_tWch'])
# plot.Group('ZLL', ['DYJetsToLL_M50_HT100to200', 'DYJetsToLL_M50_HT200to400', 'DYJetsToLL_M50_HT400to600', 'DYJetsToLL_M50_HT600toInf'])
plot.Group('data_obs', ['data_JetHT']) #, 'data_JetHT'
#['WpWpJJ', 'ZGTo2LG', 'ZGJets', 'WGToLNuG', 'WGJets', 'WW', 'WWDouble', 'WWTo1L1Nu2Q', 'WWToLNuQQ_ext', 'WWToLNuQQ', 'WWTo2L2Nu', 'WZ', 'WZTo3LNu_amcatnlo', 'WZTo3LNu', 'WZTo2L2Q', 'WZTo1L1Nu2Q', 'WZTo1L3Nu', 'ZZ', 'VVTo2L2Nu', 'ZZTo4L', 'ZZTo2Q2Nu', 'ZZTo2L2Q', 'ZZTo2L2Nu'])
HistDrawer.draw(plot,
plot_dir='plots_VV/' + cut_name,
channel='VV #rightarrow jj')
# plot.WriteDataCard(filename='datacard_mm.root', dir='mm_' + cut_name)
lumi=int_lumi,
weight=total_weight)
cfg_example.cut = cuts[cut_name]
cfg_example.vars = variables
plots = createHistograms(cfg_example, verbose=False)
for variable in variables:
plot = plots[variable.name]
plot.Group('Diboson', ['WWTo1L1Nu2Q', 'WZTo1L1Nu2Q'])
plot.Group('TT_W', ['TTJets_W'])
plot.Group('W', [
'WJetsToLNu_HT100to200', 'WJetsToLNu_HT200to400',
'WJetsToLNu_HT400to600', 'WJetsToLNu_HT600to800',
'WJetsToLNu_HT800to1200', 'WJetsToLNu_HT1200to2500',
'WJetsToLNu_HT2500toInf'
])
# plot.Group('QCD', ['QCD_HT2000toInf', 'QCD_HT1500to2000', 'QCD_HT1000to1500', 'QCD_HT700to1000', 'QCD_HT500to700']) # , 'QCD_HT300to500'
plot.Group('TT_nonW', ['TTJets_nonW'])
# plot.Group('Single t', ['TToLeptons_tch_powheg', 'TBarToLeptons_tch_powheg', 'TToLeptons_sch', 'TBar_tWch', 'T_tWch'])
# plot.Group('ZLL', ['DYJetsToLL_M50_HT100to200', 'DYJetsToLL_M50_HT200to400', 'DYJetsToLL_M50_HT400to600', 'DYJetsToLL_M50_HT600toInf'])
plot.Group('data_obs',
['data_SingleMuon', 'data_SingleElectron', 'data_MET'
]) #, 'data_JetHT'
#['WpWpJJ', 'ZGTo2LG', 'ZGJets', 'WGToLNuG', 'WGJets', 'WW', 'WWDouble', 'WWTo1L1Nu2Q', 'WWToLNuQQ_ext', 'WWToLNuQQ', 'WWTo2L2Nu', 'WZ', 'WZTo3LNu_amcatnlo', 'WZTo3LNu', 'WZTo2L2Q', 'WZTo1L1Nu2Q', 'WZTo1L3Nu', 'ZZ', 'VVTo2L2Nu', 'ZZTo4L', 'ZZTo2Q2Nu', 'ZZTo2L2Q', 'ZZTo2L2Nu'])
HistDrawer.draw(plot,
plot_dir='plots/' + cut_name,
channel='WV #rightarrow l#nujj')
# plot.WriteDataCard(filename='datacard_mm.root', dir='mm_' + cut_name)