#for histo, axis in zip(histoname, axistitle):
# for n,h in enumerate(histoname) :
foldername = s + '/' + i
if not os.path.exists(outputdir + foldername):
os.makedirs(outputdir + foldername)
if i == 'tNoCuts':
plotter.plot_mc_vs_data(foldername,
'CUT_FLOW',
rebin=1,
xaxis='CUT_FLOW',
leftside=False,
show_ratio=False,
ratio_range=1.5,
sorted=True)
plotter.save(foldername + '/CUT_FLOW')
foldername = s + '/' + i
for h in histoname:
#plotter.canvas.SetLogy(True)
#plotter.plot_data(foldername, h, rebin=rebins[n], xaxis= axistitle[n], leftside=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.plot_mc_vs_data(foldername,
h[0],
rebin=h[2],
xaxis=h[1],
leftside=False,
show_ratio=True,
ratio_range=1.5,
sorted=True)
# Override ordering
plotter.mc_samples = [
#'WZJetsTo3LNu*',
#'ZZJetsTo4L*',
'QCD',
'WW*',
'WplusJets_madgraph',
'TTplusJets_madgraph',
'Zjets_M50',
]
plotter.plot_mc_vs_data('em', 'emMass', rebin=5, leftside=False,
xaxis='m_{e#mu} (GeV)')
plotter.add_cms_blurb(sqrts)
plotter.save('mass')
plotter.plot_mc_vs_data('em', 'mPt')
plotter.save('mPt')
plotter.plot_mc_vs_data('em', 'ePt', rebin=10)
plotter.save('ePt')
plotter.plot_mc_vs_data('em', 'mAbsEta')
plotter.save('mAbsEta')
plotter.plot_mc_vs_data('em', 'eAbsEta', rebin=5)
plotter.save('eAbsEta')
plotter.plot_mc_vs_data('em', 'nvtx')
plotter.save('nvtx')
plotter.plot_mc_vs_data('em', 'bjetCSVVeto')
'WW*',
'TT*',
'WplusJets*',
"data_DoubleMu*",
]
files = []
lumifiles = []
for x in samples:
files.extend(glob.glob('results/%s/ControlZMM/%s.root' % (jobid, x)))
lumifiles.extend(glob.glob('inputs/%s/%s.lumicalc.sum' % (jobid, x)))
plotter = Plotter(files, lumifiles, output_dir)
plotter.plot_mc_vs_data('zmm', 'm1m2Mass', rebin=4)
plotter.save('mass')
plotter.plot_mc_vs_data('zmm', 'm1Pt')
plotter.save('m1Pt')
plotter.plot_mc_vs_data('zmm', 'm2Pt')
plotter.save('m2Pt')
plotter.plot_mc_vs_data('zmm', 'm1AbsEta')
plotter.save('m1AbsEta')
plotter.plot_mc_vs_data('zmm', 'm2AbsEta')
plotter.save('m2AbsEta')
plotter.plot_mc_vs_data('zmm', 'nvtx')
plotter.save('nvtx')
files = []
lumifiles = []
for x in samples:
files.extend(glob.glob('results/%s/ControlZMM/%s.root' % (jobid, x)))
lumifiles.extend(glob.glob('inputs/%s/%s.lumicalc.sum' % (jobid, x)))
plotter = Plotter(files, lumifiles, output_dir)
plotter.mc_samples = ['Zjets_M50']
sqrts = 7 if '7TeV' in jobid else 8
plotter.plot_mc_vs_data('zmm', 'm1m2Mass', rebin=2, xaxis='m_{#mu#mu} (GeV)')
plotter.add_cms_blurb(sqrts)
plotter.save('mass')
plotter.plot_mc_vs_data('zmm', 'm1m2Mass', rebin=6, xaxis='m_{#mu#mu} (GeV)')
plotter.add_cms_blurb(sqrts)
plotter.save('mass_rebin')
plotter.plot_mc_vs_data('zmm', 'm1Pt')
plotter.save('m1Pt')
plotter.plot_mc_vs_data('zmm', 'm1Pt', 5)
plotter.save('m1Pt_rebin')
plotter.plot_mc_vs_data('zmm', 'm2Pt')
plotter.save('m2Pt')
plotter.plot_mc_vs_data('zmm', 'm2Pt', 5)
plotter.save('m2Pt_rebin')
plotter.plot_mc_vs_data('zmm', 'm1AbsEta')
# Override ordering
plotter.mc_samples = [
'TTplusJets_madgraph',
'WplusJets_madgraph',
'Zjets_M50',
'WZJetsTo3LNu*',
'WW*',
'ZZJetsTo4L*',
]
sqrts = 7 if '7TeV' in jobid else 8
plotter.plot_mc_vs_data('em', 'emMass', rebin=10, leftside=False,
xaxis='m_{e#mu} (GeV)')
plotter.add_cms_blurb(sqrts)
plotter.save('mass')
plotter.plot_mc_vs_data('em', 'mPt')
plotter.save('mPt')
plotter.plot_mc_vs_data('em', 'ePt')
plotter.save('ePt')
plotter.plot_mc_vs_data('em', 'mAbsEta')
plotter.save('mAbsEta')
plotter.plot_mc_vs_data('em', 'eAbsEta')
plotter.save('eAbsEta')
plotter.plot_mc_vs_data('em', 'nvtx')
plotter.save('nvtx')
plotter.plot_mc_vs_data('em', 'bjetCSVVeto')
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
##
##files=[]
##lumifiles=[]
##channel = 'em'
##for x in filelist:
'TTplusJets_madgraph-no-fakes',
#'TTplusJets_madgraph',
#'WZJetsTo3LNu*',
#'ZZJetsTo4L*',
'Zjets_M50-no-fakes',
#'Zjets_M50',
]
for suffix, samples in [('', os_ss_samples), ('-fakes', fakes_samples)]:
plotter.mc_samples = samples
#from pdb import set_trace; set_trace()
plotter.plot_mc_vs_data('em', 'emMass', rebin=5, leftside=False,
xaxis='m_{e#mu} (GeV)')#, show_ratio=True)
plotter.add_cms_blurb(sqrts)
plotter.save('mass' + suffix)
plotter.plot_mc_vs_data('em', 'emMass', rebin=10, leftside=False,
xaxis='m_{e#mu} (GeV)')#, show_ratio=True)
plotter.add_cms_blurb(sqrts)
plotter.save('mass_rebin' + suffix)
plotter.plot_mc_vs_data('em', 'mPt', rebin=10, show_ratio=True)
plotter.save('mPt' + suffix)
plotter.plot_mc_vs_data('em', 'ePt', rebin=10, show_ratio=True)
plotter.save('ePt' + suffix)
plotter.plot_mc_vs_data('em', 'mAbsEta', rebin=5, show_ratio=True)
plotter.save('mAbsEta' + suffix)
plotter.plot_mc_vs_data('em', 'eAbsEta', rebin=5, show_ratio=True)
plotter.save('eAbsEta' + suffix)
#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'],
h[0],
rebin=h[2],
xaxis=h[1],
leftside=False,
show_ratio=False,
ratio_range=3,
rescale=10)
if not os.path.exists(outputdir + foldername):
os.makedirs(outputdir + foldername)
