axistitle = [
'm1 p_{T} (GeV)', 'm2 p_{T} (GeV)', 'm1-m2 Inv Mass (GeV)',
'm1 Gen p_{T} (GeV)', 'm2 Gen p_{T} (GeV)'
]
foldername = channel #+"/fromHiggs"
if not os.path.exists(outputdir + '/' + foldername):
os.makedirs(outputdir + '/' + foldername)
for n, h in enumerate(histoname):
# plotter.simpleplot_mc(foldername, ['VBF_LFV_HToMuTau_M125*', 'GluGlu_LFV_HToMuTau_M125*'], h[0], rebin=h[2], xaxis=h[1], leftside=False, xrange=None , preprocess=None, sort=True, forceLumi=1000)
# plotter.simpleplot_mc(foldername, 'WplusHToTauTau_M125*', h[0], rebin=h[2], xaxis=h[1], leftside=False, xrange=None , preprocess=None, sort=True, forceLumi=1000)
plotter.plot_mc(foldername,
h,
1,
xaxis=axistitle[n],
leftside=False,
show_ratio=False,
ratio_range=1.5,
sort=False)
# plotter.plot_mc_vs_data(foldername, h, rebin=rebins[n], xaxis= axistitle[n], leftside=False, show_ratio=True, ratio_range=1.5, sort=True)
plotter.save(foldername + '/' + h)
##canvas = ROOT.TCanvas("canvas","canvas",800,800)
##legend = ROOT.TLegend(0.2,0.8, 0.4, 0.7)
##LFVStack = ROOT.THStack("stack","")
##
##mypath = 'results/%s/LFVEMuAnalyserGen/' %jobid
##filelist = [f for f in listdir(mypath) if (isfile(join(mypath, f)) and 'data' not in f)]
##
###print lfvfilelist
###print smfilelist
#plotter.mc_samples = mc_samples
for i in sign:
for j in process:
for k in ptcut:
for nj in njets:
for n, h in enumerate(histoname):
foldername = i + '/' + j + '/ept' + str(
int(k)) + '/' + str(int(nj))
#plotter.canvas.SetLogy(True)
plotter.plot_mc(foldername,
['ggHiggsToETau', 'vbfHiggToETau'],
h[0],
rebin=h[2],
xaxis=h[1],
leftside=False,
show_ratio=False,
ratio_range=1.5,
rescale=10)
#plotter.simpleplot_mc(foldername,h[0], rebin=h[2], xaxis= h[1], leftside=False)
if not os.path.exists(outputdir + foldername):
os.makedirs(outputdir + foldername)
plotter.save(foldername + '/mc_' + h)
foldername = i + '/' + j + '/ept' + str(k) + '/' + str(
nj) + '/selected'
#plotter.canvas.SetLogy(True)
plotter.plot_mc(foldername,
['ggHiggsToETau', 'vbfHiggsToETau'],
