-
Notifications
You must be signed in to change notification settings - Fork 2
/
bkg.py
433 lines (366 loc) · 26.6 KB
/
bkg.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
433
from os import system
import ROOT
from ROOT import TFile,THStack,TCanvas,TLegend,kRed,kYellow,kWhite,TLatex,kBlue,TIter,TList
from sys import argv
#setup
print 'setup'
ROOT.gROOT.SetBatch()
ROOT.gROOT.SetBatch()
ROOT.gROOT.LoadMacro("tdrstyle.C");
ROOT.gROOT.ProcessLine( 'setTDRStyle();')
#setTDRStyle()
ROOT.gROOT.LoadMacro("CMS_lumi.C+");
path_base='/nfs/dust/cms/user/usaiem/ZprimeFullHad/'
process_list_ttbar=['BackgroundCycle.MC.TTbar.root']#['BackgroundCycle.MC.TTbarHad.root','BackgroundCycle.MC.TTbarLept.root','BackgroundCycle.MC.TTbarSemi.root']
process_list_qcd=['BackgroundCycle.MC.QCD_HT-250To500.root','BackgroundCycle.MC.QCD_HT-1000ToInf.root','BackgroundCycle.MC.QCD_HT-500To1000.root']#,'BackgroundCycle.MC.QCD_HT-100To250.root']'BackgroundCycle.MC.QCD_HT-250To500.root',
process_data='BackgroundCycle.DATA.DATA'
process_list_data=['BackgroundCycle.DATA.MJDATAB.root','BackgroundCycle.DATA.MJDATAC.root','BackgroundCycle.DATA.MJDATAD.root']
#process_list_data=['BackgroundCycle.DATA.MJDATAD.root']
den_name="den_mistag"
num_name="num_mistag"
folder_old="BaseHistos"
folder_ht="HTDatasetHistos"
folder_quad="QuadJetDatasetHistos"
histo_massshape_name="mass_shape"
#csvcat=['CSV0-L','CSVL-M','CSVM-1']
postfix=''
if len(argv[1:])>0:
postfix='_'+argv[1]
outfile=TFile("bkg"+postfix+".root","RECREATE")
for i in [folder_ht,folder_quad]:
for j in ['MassShape','Mistag']:
outfile.mkdir(i+'/'+j)
print 'hadd plots'
def hadd(inputlist,outputname):
command_list='hadd -f '+path_base+outputname+'.root'# -v 0
for i in inputlist:
command_list+=' '+path_base+i
system(command_list)
return path_base+outputname+'.root'
ttbar_filename= path_base + process_list_ttbar[0]#hadd(process_list_ttbar,'bkg_ttbar')
qcd_filename=hadd(process_list_qcd,'bkg_qcd')
mjdata_filename=hadd(process_list_data,'bkg_mjdata')
ttbar_file=TFile(ttbar_filename,'READ')
qcd_file=TFile(qcd_filename,'READ')
mjdata_file=TFile(mjdata_filename,'READ')
data_file=TFile(path_base+process_data+'.root','READ')
qcdflat_file=TFile(path_base+'BackgroundCycle.MC.QCDFlat.root','READ')
qcd250to500_file=TFile(path_base+'BackgroundCycle.MC.QCD_HT-250To500.root','READ')
qcd500to1000_file=TFile(path_base+'BackgroundCycle.MC.QCD_HT-500To1000.root','READ')
qcd1000toInf_file=TFile(path_base+'BackgroundCycle.MC.QCD_HT-1000ToInf.root','READ')
print 'produce hist'
def slice_and_save(sample,histo):
#histo_stack=0
if '2' in histo.ClassName():
#histoXY=histo.Project3D('yx')
#histoXY.Write(histoXY.GetName())
#histo_stack=THStack(histoXY,'x','Stack_'+histo.GetName(),'')
#else:
histo_stack=THStack(histo,'x','Stack_'+histo.GetName(),'')
histo_1d=histo_stack.GetHists()
#outfile.cd()
histo_stack.Write(histo_stack.GetName()+'_'+sample)
nextinlist=TIter(histo_1d)
obj=nextinlist()
while obj:
obj.Write(obj.GetName()+'_'+sample)
obj=nextinlist()
histo.Write(histo.GetName()+'_'+sample)
def getMistag(sample,sample_file,useCMS=False,isData=False,namepostfix='',folder=folder_old):
num_histo=sample_file.Get(folder+"/"+num_name+namepostfix).Clone('Numerator')#.Clone('Numerator_'+sample)
den_histo=sample_file.Get(folder+"/"+den_name+namepostfix).Clone('Denominator')#.Clone('Denominator_'+sample)
if isData:
num_histo.Add(ttbar_file.Get(folder+"/"+num_name),-1.0)
den_histo.Add(ttbar_file.Get(folder+"/"+den_name),-1.0)
mistag_histo=num_histo.Clone('Mistag')#.Clone('Mistag_'+sample)
mistag_histo.Divide(num_histo,den_histo,1,1,'B')
for ibin in range(1,mistag_histo.GetNbinsX()+1):
for jbin in range(1,mistag_histo.GetNbinsY()+1):
if mistag_histo.GetBinContent(ibin,jbin)<0:
mistag_histo.SetBinContent(ibin,jbin,0.0)
outfolder=folder+'/Mistag'
if useCMS:
outfolder='CMSTagger/Mistag'
outfile.mkdir(outfolder+'/'+sample)
outfile.cd(outfolder+'/'+sample)
slice_and_save(sample,num_histo)
slice_and_save(sample,den_histo)
slice_and_save(sample,mistag_histo)
if '3' in mistag_histo.ClassName():
numXY=num_histo.Project3D('yx')
denXY=den_histo.Project3D('yx')
mistagXY=numXY.Clone('Mistag')
mistagXY.Divide(numXY,denXY,1,1,'B')
slice_and_save(sample,mistagXY)
#mass_shape2D.Write(folder+mass_shape2D.GetName()+'_'+sample)
#mass_shapeStack.Write(folder+mass_shapeStack.GetName()+'_'+sample)
#nextinlist=TList(GetListOfPrimitives())
#obj=nextinlist()
#while obj:
#obj.Write(folder+obj.GetName()+'_'+sample)
#obj=nextinlist()
#def getMistag2(sample,sample_file,useCMS=False,isData=False):
#num_histo=sample_file.Get(folder+"/num_mistagMeasured").Clone('Numerator')#.Clone('Numerator_'+sample)
#den_histo=sample_file.Get(folder+"/den_mistagMeasured").Clone('Denominator')#.Clone('Denominator_'+sample)
#if isData:
#num_histo.Add(ttbar_file.Get(folder+"/"+num_name),-1.0)
#den_histo.Add(ttbar_file.Get(folder+"/"+den_name),-1.0)
#mistag_histo=num_histo.Clone('Mistag')#.Clone('Mistag_'+sample)
#mistag_histo.Divide(num_histo,den_histo,1,1,'B')
#outfolder='HEPTagger/Mistag'
#if useCMS:
#outfolder='CMSTagger/Mistag'
#outfile.mkdir(outfolder+'/'+sample)
#outfile.cd(outfolder+'/'+sample)
#slice_and_save(sample,num_histo)
#slice_and_save(sample,den_histo)
#slice_and_save(sample,mistag_histo)
##mass_shape2D.Write(folder+mass_shape2D.GetName()+'_'+sample)
##mass_shapeStack.Write(folder+mass_shapeStack.GetName()+'_'+sample)
##nextinlist=TList(GetListOfPrimitives())
##obj=nextinlist()
##while obj:
##obj.Write(folder+obj.GetName()+'_'+sample)
##obj=nextinlist()
def getMassShape(sample,sample_file,useCMS=False,isData=False,histoname='', folder=folder_old,rebin=0):
histomassshapename=histo_massshape_name
if histoname!='':
histomassshapename=histoname
mass_shape=sample_file.Get(folder+'/'+histomassshapename).Clone('mass_shape')
mass_shape_high=sample_file.Get(folder+'/'+histomassshapename+'_high').Clone('mass_shape_high')
mass_shape_low =sample_file.Get(folder+'/'+histomassshapename+'_low' ).Clone('mass_shape_low' )
if isData:
mass_shape.Add(ttbar_file.Get(folder+"/"+histomassshapename),-1.0)
mass_shape_high.Add(ttbar_file.Get(folder+"/"+histomassshapename+'_high'),-1.0)
mass_shape_low.Add(ttbar_file.Get(folder+"/"+histomassshapename+'_low'),-1.0)
#mass_shapeStack=THStack(mass_shape2D,'x','Stack','')
#mass_shape1D=mass_shapeStack.GetHists()
outfolder=folder+'/MassShape'
if useCMS:
outfolder='CMSTagger/MassShape'
outfile.mkdir(outfolder+'/'+sample)
outfile.cd(outfolder+'/'+sample)
#slice_and_save(sample,mass_shape2D,folder)
if rebin!=0:
mass_shape_high.Rebin(rebin)
mass_shape_low.Rebin(rebin)
mass_shape.Rebin(rebin)
mass_shape.Write(mass_shape.GetName())
mass_shape_high.Write(mass_shape_high.GetName())#+'_'+sample)
mass_shape_low.Write(mass_shape_low.GetName())
#mass_shape2D.Write(folder+mass_shape2D.GetName()+'_'+sample)
#mass_shapeStack.Write(folder+mass_shapeStack.GetName()+'_'+sample)
#nextinlist=TList(GetListOfPrimitives())
#obj=nextinlist()
#while obj:
#obj.Write(folder+obj.GetName()+'_'+sample)
#obj=nextinlist()
def getMassShape2D(sample,sample_file,useCMS=False,isData=False):
mass_shape=sample_file.Get(folder+'/'+histo_massshape_name+'2D')
if isData:
mass_shape.Add(ttbar_file.Get(folder+"/"+histo_massshape_name+'2D'),-1.0)
#mass_shapeStack=THStack(mass_shape2D,'x','Stack','')
#mass_shape1D=mass_shapeStack.GetHists()
outfolder='HEPTagger/MassShape2D'
if useCMS:
outfolder='CMSTagger/MassShape2D'
outfile.mkdir(outfolder+'/'+sample)
outfile.cd(outfolder+'/'+sample)
#slice_and_save(sample,mass_shape2D,folder)
slice_and_save(sample,mass_shape)
#mass_shape2D.Write(folder+mass_shape2D.GetName()+'_'+sample)
#mass_shapeStack.Write(folder+mass_shapeStack.GetName()+'_'+sample)
#nextinlist=TList(GetListOfPrimitives())
#obj=nextinlist()
#while obj:
#obj.Write(folder+obj.GetName()+'_'+sample)
#obj=nextinlist()
use_cms=True
use_htt=False
ttbar_subtraction=True
no_ttbar_subtraction=False
getMistag('ttbar_htt', ttbar_file, use_htt, no_ttbar_subtraction, '', folder_ht)
getMistag('qcd_htt', qcd_file, use_htt, no_ttbar_subtraction, '', folder_ht)
getMistag('data_htt', data_file, use_htt, ttbar_subtraction, '', folder_ht)
getMistag('data_htt_nosub', data_file, use_htt, no_ttbar_subtraction, '', folder_ht)
getMistag('ttbar_htt', ttbar_file, use_htt, no_ttbar_subtraction, '', folder_quad)
getMistag('qcd_htt', qcd_file, use_htt, no_ttbar_subtraction, '', folder_quad)
getMistag('data_htt', mjdata_file, use_htt, ttbar_subtraction, '', folder_quad)
getMistag('data_htt_nosub', mjdata_file, use_htt, no_ttbar_subtraction, '', folder_quad)
#getMistag('qcdflat', qcdflat_file, use_htt, no_ttbar_subtraction)
#getMistag('qcd250to500', qcd250to500_file, use_htt, no_ttbar_subtraction, '', folder_quad)
getMistag('qcd500to1000', qcd500to1000_file, use_htt, no_ttbar_subtraction, '', folder_quad)
#getMistag('qcd500to1000Measured', qcd500to1000_file, use_htt, no_ttbar_subtraction,'Measured')
#getMistag('qcd500to1000Tag', qcd500to1000_file, use_htt, no_ttbar_subtraction,'Tag')
#getMistag('qcd500to1000R', qcd500to1000_file, use_htt, no_ttbar_subtraction,'R')
#getMistag('qcd500to1000Rmw', qcd500to1000_file, use_htt, no_ttbar_subtraction,'Rmw')
#getMistag('qcd500to1000Invmw', qcd500to1000_file, use_htt, no_ttbar_subtraction,'Invmw')
#getMistag('qcd500to1000Nomw', qcd500to1000_file, use_htt, no_ttbar_subtraction,'Nomw')
#getMistag('qcd1000toInf', qcd1000toInf_file, use_htt, no_ttbar_subtraction)
#getMistag('qcd1000toInfMeasured', qcd1000toInf_file, use_htt, no_ttbar_subtraction,'Measured')
#getMistag('qcd1000toInfTag', qcd1000toInf_file, use_htt, no_ttbar_subtraction,'Tag')
#getMistag('ttbar_cms', ttbar_file, use_cms, no_ttbar_subtraction)
#getMistag('qcd_cms', qcd_file, use_cms, no_ttbar_subtraction)
#getMistag('data_cms', data_file, use_cms, ttbar_subtraction)
#getMistag('data_cms_nosub', data_file, use_cms, no_ttbar_subtraction)
#getMassShape('ttbar_htt', ttbar_file, use_htt, no_ttbar_subtraction, '', folder_ht)
getMassShape('qcd_htt', qcd_file, use_htt, no_ttbar_subtraction, '', folder_ht,4)
#getMassShape('data_htt', data_file, use_htt, ttbar_subtraction, '', folder_ht)
#getMassShape('data_htt_nosub', data_file, use_htt, no_ttbar_subtraction, '', folder_ht)
#getMassShape('ttbar_htt', ttbar_file, use_htt, no_ttbar_subtraction, '', folder_quad)
getMassShape('qcd_htt', qcd_file, use_htt, no_ttbar_subtraction, '', folder_quad,4)
#getMassShape('data_htt', mjdata_file, use_htt, ttbar_subtraction, '', folder_quad)
#getMassShape('data_htt_nosub', mjdata_file, use_htt, no_ttbar_subtraction, '', folder_quad)
#getMassShape('qcd500to1000', qcd500to1000_file, use_htt, no_ttbar_subtraction, '', folder_quad)
#getMassShape('qcd1000toInf', qcd1000toInf_file, use_htt, no_ttbar_subtraction)
#getMassShape('qcd_htt_jec', qcd_file, use_htt, no_ttbar_subtraction, "mass_shapeJEC")
#getMassShape('ttbar_cms', ttbar_file, use_cms, no_ttbar_subtraction)
#getMassShape('qcd_cms', qcd_file, use_cms, no_ttbar_subtraction)
#getMassShape('data_cms', data_file, use_cms, ttbar_subtraction)
#getMassShape('data_cms_nosub', data_file, use_cms, no_ttbar_subtraction)
#getMassShape2D('ttbar_htt', ttbar_file, use_htt, no_ttbar_subtraction)
#getMassShape2D('qcd_htt', qcd_file, use_htt, no_ttbar_subtraction)
#getMassShape2D('data_htt', data_file, use_htt, ttbar_subtraction)
#getMassShape2D('data_htt_nosub', data_file, use_htt, no_ttbar_subtraction)
#nmsfile=TFile('/nfs/dust/cms/user/usaiem/ZprimeFullHad/binmass/theta_qcd.root','READ')
#nmslow=nmsfile.Get('QuadJetDatasetHistos/mass_shape_low').Clone('mass_shape_low')
#nmshigh=nmsfile.Get('QuadJetDatasetHistos/mass_shape_high').Clone('mass_shape_high')
#nmslow2=nmsfile.Get('QuadJetDatasetHistos/mass_shape_lowpt').Clone('mass_shape_lowpt')
#nmshigh2=nmsfile.Get('QuadJetDatasetHistos/mass_shape_highpt').Clone('mass_shape_highpt')
##outfile.mkdir('HEPTagger/MassShape')
#outfile.cd('QuadJetDatasetHistos/MassShape')
#nmslow.Write()
#nmshigh.Write()
##nmslow2.Write()
##nmshigh2.Write()
outfile.cd()
#ckhisto=qcd500to1000_file.Get('BaseHistos/mistag_crosscheck')
#slice_and_save('ck',ckhisto)
print 'make plots'
def make_plot(name,name_list,legend_list,normalize=False,thefile=0):
histo_list=[]
c=TCanvas(name,'',600,600)
c.SetLeftMargin(0.15)#
c.SetRightMargin(0.05)#
c.SetTopMargin(0.10)#
c.SetBottomMargin(0.10)
#legend=TLegend(0.7,0.4,0.95,0.7)
legend=TLegend(0.35,0.2,0.85,0.5)
legend.SetHeader('HEP Tagger, H_{T}>800 GeV')
if folder_quad in name:
legend=TLegend(0.35,0.2,0.85,0.5)
legend.SetHeader('HEP Tagger, H_{T}<800 GeV')
legend.SetTextSize(0.03)
legend.SetBorderSize(0)
legend.SetTextFont(42)
legend.SetLineColor(1)
legend.SetLineStyle(1)
legend.SetLineWidth(1)
legend.SetFillColor(0)
legend.SetFillStyle(0)
for i in range(len(name_list)):
if thefile==0:
histo_list.append(outfile.Get(name_list[i]))
else:
histo_list.append(thefile.Get(name_list[i]))
if normalize:
histo_list[-1].Scale(1.0/histo_list[-1].Integral())
histo_list[-1].SetStats(0)
histo_list[-1].SetLineWidth(3)
histo_list[-1].SetLineColor(i+1)
histo_list[-1].GetXaxis().SetTitle('p_{T} (GeV)')
histo_list[-1].GetYaxis().SetTitle('Misidentification probability')
histo_list[-1].GetYaxis().SetLabelSize(0.043)
histo_list[-1].GetYaxis().SetTitleSize(0.043)
histo_list[-1].GetYaxis().SetTitleOffset(1.5)
histo_list[-1].GetXaxis().SetLabelSize(0.043)
histo_list[-1].GetXaxis().SetTitleSize(0.043)
histo_list[-1].GetXaxis().SetTitleOffset(1.0)
histo_list[-1].GetXaxis().SetMoreLogLabels()
histo_list[-1].GetXaxis().SetNoExponent()
if folder_quad in name:
histo_list[-1].GetXaxis().SetRangeUser(150,400)
else:
histo_list[-1].GetXaxis().SetRangeUser(150,800)
# histo_list[-1].GetYaxis().SetMoreLogLabels()
# histo_list[-1].GetYaxis().SetNoExponent()
c.SetLogy()
c.SetLogx()
legend.AddEntry(histo_list[-1],legend_list[i],'l')
if i==0:
histo_list[-1].Draw()
else:
histo_list[-1].Draw('SAME')
histo_list[0]
legend.Draw()
outfile.cd()
if folder_ht in name:
ROOT.CMS_lumi( c, 2, 33 )
if folder_quad in name:
ROOT.CMS_lumi( c, 1, 33 )
c.Write(name)
c.SaveAs('pdf/'+name+'.pdf')
#thisfile=TFile('dust/BackgroundCycle.MC.ZP500W5.root','READ')
#make_plot('comp_sgn_500',['BaseHistos/MeasuredMtt012','BaseHistos/MeasuredJECMtt012','BaseHistos/MeasuredHTTMtt012','BaseHistos/MeasuredHTTJECMtt012','BaseHistos/MeasuredAK5Mtt012'],['CA15Filtered no JEC','CA15Filtered JEC','HTT CA15 no JEC','HTT CA15 JEC','3xAK5 in R=1.5'],False,thisfile)
#thisfile.Close()
#thisfile=TFile('dust/BackgroundCycle.MC.ZP750W7p5.root','READ')
#make_plot('comp_sgn_750',['BaseHistos/MeasuredMtt012','BaseHistos/MeasuredJECMtt012','BaseHistos/MeasuredHTTMtt012','BaseHistos/MeasuredHTTJECMtt012','BaseHistos/MeasuredAK5Mtt012'],['CA15Filtered no JEC','CA15Filtered JEC','HTT CA15 no JEC','HTT CA15 JEC','3xAK5 in R=1.5'],False,thisfile)
#thisfile.Close()
#thisfile=TFile('dust/BackgroundCycle.MC.ZP1000W10.root','READ')
#make_plot('comp_sgn_1000',['BaseHistos/MeasuredMtt012','BaseHistos/MeasuredJECMtt012','BaseHistos/MeasuredHTTMtt012','BaseHistos/MeasuredHTTJECMtt012','BaseHistos/MeasuredAK5Mtt012'],['CA15Filtered no JEC','CA15Filtered JEC','HTT CA15 no JEC','HTT CA15 JEC','3xAK5 in R=1.5'],False,thisfile)
#thisfile.Close()
#thisfile=TFile('dust/BackgroundCycle.MC.ZP1250W12p5.root','READ')
#make_plot('comp_sgn_1250',['BaseHistos/MeasuredMtt012','BaseHistos/MeasuredJECMtt012','BaseHistos/MeasuredHTTMtt012','BaseHistos/MeasuredHTTJECMtt012','BaseHistos/MeasuredAK5Mtt012'],['CA15Filtered no JEC','CA15Filtered JEC','HTT CA15 no JEC','HTT CA15 JEC','3xAK5 in R=1.5'],False,thisfile)
#thisfile.Close()
#thisfile=TFile('dust/BackgroundCycle.MC.ZP1500W15.root','READ')
#make_plot('comp_sgn_1500',['BaseHistos/MeasuredMtt012','BaseHistos/MeasuredJECMtt012','BaseHistos/MeasuredHTTMtt012','BaseHistos/MeasuredHTTJECMtt012','BaseHistos/MeasuredAK5Mtt012'],['CA15Filtered no JEC','CA15Filtered JEC','HTT CA15 no JEC','HTT CA15 JEC','3xAK5 in R=1.5'],False,thisfile)
#thisfile.Close()
#thisfile=TFile('dust/BackgroundCycle.MC.ZP2000W20.root','READ')
#make_plot('comp_sgn_2000',['BaseHistos/MeasuredMtt012','BaseHistos/MeasuredJECMtt012','BaseHistos/MeasuredHTTMtt012','BaseHistos/MeasuredHTTJECMtt012','BaseHistos/MeasuredAK5Mtt012'],['CA15Filtered no JEC','CA15Filtered JEC','HTT CA15 no JEC','HTT CA15 JEC','3xAK5 in R=1.5'],False,thisfile)
#thisfile.Close()
#thisfile=TFile('dust/BackgroundCycle.MC.ZP3000W30.root','READ')
#make_plot('comp_sgn_3000',['BaseHistos/MeasuredMtt012','BaseHistos/MeasuredJECMtt012','BaseHistos/MeasuredHTTMtt012','BaseHistos/MeasuredHTTJECMtt012','BaseHistos/MeasuredAK5Mtt012'],['CA15Filtered no JEC','CA15Filtered JEC','HTT CA15 no JEC','HTT CA15 JEC','3xAK5 in R=1.5'],False,thisfile)
#thisfile.Close()
#thisfile=TFile('dust/BackgroundCycle.MC.ZP4000W40.root','READ')
#make_plot('comp_sgn_4000',['BaseHistos/MeasuredMtt012','BaseHistos/MeasuredJECMtt012','BaseHistos/MeasuredHTTMtt012','BaseHistos/MeasuredHTTJECMtt012','BaseHistos/MeasuredAK5Mtt012'],['CA15Filtered no JEC','CA15Filtered JEC','HTT CA15 no JEC','HTT CA15 JEC','3xAK5 in R=1.5'],False,thisfile)
#thisfile.Close()
#########################
for i in [folder_ht,folder_quad]:
#make_plot('comp_mistag_'+i,[i+'/Mistag/qcd_htt/Mistag_px4_qcd_htt',i+'/Mistag/qcd_htt/Mistag_px3_qcd_htt',i+'/Mistag/qcd_htt/Mistag_px2_qcd_htt',i+'/Mistag/qcd_htt/Mistag_px1_qcd_htt'],['[CSVM,1]','[CSVL,CSVM]','[0,CSVL]','[-1,0]'])
#make_plot('comp_mistag_data_'+i,[i+'/Mistag/data_htt/Mistag_px4_data_htt',i+'/Mistag/data_htt/Mistag_px3_data_htt',i+'/Mistag/data_htt/Mistag_px2_data_htt',i+'/Mistag/data_htt/Mistag_px1_data_htt'],['[CSVM,1]','[CSVL,CSVM]','[0,CSVL]','[-1,0]'])
#make_plot('comp_mistag_data_nosub_'+i,[i+'/Mistag/data_htt_nosub/Mistag_px4_data_htt_nosub',i+'/Mistag/data_htt_nosub/Mistag_px3_data_htt_nosub',i+'/Mistag/data_htt_nosub/Mistag_px2_data_htt_nosub',i+'/Mistag/data_htt_nosub/Mistag_px1_data_htt_nosub'],['[CSVM,1]','[CSVL,CSVM]','[0,CSVL]','[-1,0]'])
#make_plot('comp_mistag_'+i,[i+'/Mistag/qcd_htt/Mistag_px4_qcd_htt',i+'/Mistag/qcd_htt/Mistag_px3_qcd_htt',i+'/Mistag/qcd_htt/Mistag_px2_qcd_htt'],['[CSVM,1]','[CSVL,CSVM]','[0,CSVL]','[-1,0]'])
make_plot('comp_mistag_data_'+i,[i+'/Mistag/data_htt/Mistag_px4_data_htt',i+'/Mistag/data_htt/Mistag_px3_data_htt',i+'/Mistag/data_htt/Mistag_px2_data_htt'],['[CSVM,1]','[CSVL,CSVM]','[0,CSVL]','[-1,0]'])
make_plot('comp_mistag_data_nosub_'+i,[i+'/Mistag/data_htt_nosub/Mistag_px4_data_htt_nosub',i+'/Mistag/data_htt_nosub/Mistag_px3_data_htt_nosub',i+'/Mistag/data_htt_nosub/Mistag_px2_data_htt_nosub'],['[CSVM,1]','[CSVL,CSVM]','[0,CSVL]','[-1,0]'])
#make_plot('comp_mistag500To1000',['HEPTagger/Mistag/qcd500to1000/Mistag_px4_qcd500to1000','HEPTagger/Mistag/qcd500to1000/Mistag_px3_qcd500to1000','HEPTagger/Mistag/qcd500to1000/Mistag_px2_qcd500to1000','HEPTagger/Mistag/qcd500to1000/Mistag_px1_qcd500to1000'],['[CSVM,1]','[CSVL,CSVM]','[0,CSVL]','[-1,0]'])
#make_plot('comp_mistag500To1000Measured',['HEPTagger/Mistag/qcd500to1000Measured/Mistag_px4_qcd500to1000Measured','HEPTagger/Mistag/qcd500to1000Measured/Mistag_px3_qcd500to1000Measured','HEPTagger/Mistag/qcd500to1000Measured/Mistag_px2_qcd500to1000Measured','HEPTagger/Mistag/qcd500to1000Measured/Mistag_px1_qcd500to1000Measured'],['[CSVM,1]','[CSVL,CSVM]','[0,CSVL]','[-1,0]'])
#make_plot('comp_mistag500To1000Tag',['HEPTagger/Mistag/qcd500to1000Tag/Mistag_px4_qcd500to1000Tag','HEPTagger/Mistag/qcd500to1000Tag/Mistag_px3_qcd500to1000Tag','HEPTagger/Mistag/qcd500to1000Tag/Mistag_px2_qcd500to1000Tag','HEPTagger/Mistag/qcd500to1000Tag/Mistag_px1_qcd500to1000Tag'],['[CSVM,1]','[CSVL,CSVM]','[0,CSVL]','[-1,0]'])
#make_plot('comp_mistag1000ToInf',['HEPTagger/Mistag/qcd1000toInf/Mistag_px4_qcd1000toInf','HEPTagger/Mistag/qcd1000toInf/Mistag_px3_qcd1000toInf','HEPTagger/Mistag/qcd1000toInf/Mistag_px2_qcd1000toInf','HEPTagger/Mistag/qcd1000toInf/Mistag_px1_qcd1000toInf'],['[CSVM,1]','[CSVL,CSVM]','[0,CSVL]','[-1,0]'])
#make_plot('comp_mistag1000ToInfMeasured',['HEPTagger/Mistag/qcd1000toInfMeasured/Mistag_px4_qcd1000toInfMeasured','HEPTagger/Mistag/qcd1000toInfMeasured/Mistag_px3_qcd1000toInfMeasured','HEPTagger/Mistag/qcd1000toInfMeasured/Mistag_px2_qcd1000toInfMeasured','HEPTagger/Mistag/qcd1000toInfMeasured/Mistag_px1_qcd1000toInfMeasured'],['[CSVM,1]','[CSVL,CSVM]','[0,CSVL]','[-1,0]'])
#make_plot('comp_mistag1000ToInfTag',['HEPTagger/Mistag/qcd1000toInfTag/Mistag_px4_qcd1000toInfTag','HEPTagger/Mistag/qcd1000toInfTag/Mistag_px3_qcd1000toInfTag','HEPTagger/Mistag/qcd1000toInfTag/Mistag_px2_qcd1000toInfTag','HEPTagger/Mistag/qcd1000toInfTag/Mistag_px1_qcd1000toInfTag'],['[CSVM,1]','[CSVL,CSVM]','[0,CSVL]','[-1,0]'])
#make_plot('comp_mistag_1000ToInfTag_px4',['HEPTagger/Mistag/qcd1000toInf/Mistag_px4_qcd1000toInf','HEPTagger/Mistag/qcd1000toInfMeasured/Mistag_px4_qcd1000toInfMeasured','HEPTagger/Mistag/qcd1000toInfTag/Mistag_px4_qcd1000toInfTag'],['antitag','no requirement','tag'])
#make_plot('comp_mistag_1000ToInfTag_px3',['HEPTagger/Mistag/qcd1000toInf/Mistag_px3_qcd1000toInf','HEPTagger/Mistag/qcd1000toInfMeasured/Mistag_px3_qcd1000toInfMeasured','HEPTagger/Mistag/qcd1000toInfTag/Mistag_px3_qcd1000toInfTag'],['antitag','no requirement','tag'])
#make_plot('comp_mistag_1000ToInfTag_px2',['HEPTagger/Mistag/qcd1000toInf/Mistag_px2_qcd1000toInf','HEPTagger/Mistag/qcd1000toInfMeasured/Mistag_px2_qcd1000toInfMeasured','HEPTagger/Mistag/qcd1000toInfTag/Mistag_px2_qcd1000toInfTag'],['antitag','no requirement','tag'])
#make_plot('comp_mistag_1000ToInfTag_px1',['HEPTagger/Mistag/qcd1000toInf/Mistag_px1_qcd1000toInf','HEPTagger/Mistag/qcd1000toInfMeasured/Mistag_px1_qcd1000toInfMeasured','HEPTagger/Mistag/qcd1000toInfTag/Mistag_px1_qcd1000toInfTag'],['antitag','no requirement','tag'])
#make_plot('comp_mistag_500To1000Tag_px4',['HEPTagger/Mistag/qcd500to1000/Mistag_px4_qcd500to1000','HEPTagger/Mistag/qcd500to1000Measured/Mistag_px4_qcd500to1000Measured','HEPTagger/Mistag/qcd500to1000Tag/Mistag_px4_qcd500to1000Tag'],['antitag','no requirement','tag'])
#make_plot('comp_mistag_500To1000Tag_px3',['HEPTagger/Mistag/qcd500to1000/Mistag_px3_qcd500to1000','HEPTagger/Mistag/qcd500to1000Measured/Mistag_px3_qcd500to1000Measured','HEPTagger/Mistag/qcd500to1000Tag/Mistag_px3_qcd500to1000Tag'],['antitag','no requirement','tag'])
#make_plot('comp_mistag_500To1000Tag_px2',['HEPTagger/Mistag/qcd500to1000/Mistag_px2_qcd500to1000','HEPTagger/Mistag/qcd500to1000Measured/Mistag_px2_qcd500to1000Measured','HEPTagger/Mistag/qcd500to1000Tag/Mistag_px2_qcd500to1000Tag'],['antitag','no requirement','tag'])
#make_plot('comp_mistag_500To1000Tag_px1',['HEPTagger/Mistag/qcd500to1000/Mistag_px1_qcd500to1000','HEPTagger/Mistag/qcd500to1000Measured/Mistag_px1_qcd500to1000Measured','HEPTagger/Mistag/qcd500to1000Tag/Mistag_px1_qcd500to1000Tag'],['antitag','no requirement','tag'])
#make_plot('comp_mistag_500To1000_px4',['HEPTagger/Mistag/qcd500to1000/Mistag_px4_qcd500to1000','HEPTagger/Mistag/qcd500to1000Tag/Mistag_px4_qcd500to1000Tag','HEPTagger/Mistag/qcd500to1000R/Mistag_px4_qcd500to1000R','HEPTagger/Mistag/qcd500to1000Rmw/Mistag_px4_qcd500to1000Rmw','HEPTagger/Mistag/qcd500to1000Invmw/Mistag_px4_qcd500to1000Invmw','HEPTagger/Mistag/qcd500to1000Nomw/Mistag_px4_qcd500to1000Nomw'],['antitag','tag','R','Rmw','Invmw','Nomw'])
#make_plot('comp_mistag_500To1000_px3',['HEPTagger/Mistag/qcd500to1000/Mistag_px3_qcd500to1000','HEPTagger/Mistag/qcd500to1000Tag/Mistag_px3_qcd500to1000Tag','HEPTagger/Mistag/qcd500to1000R/Mistag_px3_qcd500to1000R','HEPTagger/Mistag/qcd500to1000Rmw/Mistag_px3_qcd500to1000Rmw','HEPTagger/Mistag/qcd500to1000Invmw/Mistag_px3_qcd500to1000Invmw','HEPTagger/Mistag/qcd500to1000Nomw/Mistag_px3_qcd500to1000Nomw'],['antitag','tag','R','Rmw','Invmw','Nomw'])
#make_plot('comp_mistag_500To1000_px2',['HEPTagger/Mistag/qcd500to1000/Mistag_px2_qcd500to1000','HEPTagger/Mistag/qcd500to1000Tag/Mistag_px2_qcd500to1000Tag','HEPTagger/Mistag/qcd500to1000R/Mistag_px2_qcd500to1000R','HEPTagger/Mistag/qcd500to1000Rmw/Mistag_px2_qcd500to1000Rmw','HEPTagger/Mistag/qcd500to1000Invmw/Mistag_px2_qcd500to1000Invmw','HEPTagger/Mistag/qcd500to1000Nomw/Mistag_px2_qcd500to1000Nomw'],['antitag','tag','R','Rmw','Invmw','Nomw'])
#make_plot('comp_mistag_500To1000_px1',['HEPTagger/Mistag/qcd500to1000/Mistag_px1_qcd500to1000','HEPTagger/Mistag/qcd500to1000Tag/Mistag_px1_qcd500to1000Tag','HEPTagger/Mistag/qcd500to1000R/Mistag_px1_qcd500to1000R','HEPTagger/Mistag/qcd500to1000Rmw/Mistag_px1_qcd500to1000Rmw','HEPTagger/Mistag/qcd500to1000Invmw/Mistag_px1_qcd500to1000Invmw','HEPTagger/Mistag/qcd500to1000Nomw/Mistag_px1_qcd500to1000Nomw'],['antitag','tag','R','Rmw','Invmw','Nomw'])
#make_plot('comp_data_nosub_mistag',['HEPTagger/Mistag/data_htt_nosub/Mistag_px4_data_htt_nosub','HEPTagger/Mistag/data_htt_nosub/Mistag_px3_data_htt_nosub','HEPTagger/Mistag/data_htt_nosub/Mistag_px2_data_htt_nosub','HEPTagger/Mistag/data_htt_nosub/Mistag_px1_data_htt_nosub'],['[CSVM,1]','[CSVL,CSVM]','[0,CSVL]','[-1,0]'])
#make_plot('comp_data_mistag',['HEPTagger/Mistag/data_htt/Mistag_px4_data_htt','HEPTagger/Mistag/data_htt/Mistag_px3_data_htt','HEPTagger/Mistag/data_htt/Mistag_px2_data_htt','HEPTagger/Mistag/data_htt/Mistag_px1_data_htt'],['[CSVM,1]','[CSVL,CSVM]','[0,CSVL]','[-1,0]'])
#make_plot('comp_datamc1_mistag',['HEPTagger/Mistag/data_htt/Mistag_px4_data_htt','HEPTagger/Mistag/qcd_htt/Mistag_px4_qcd_htt'],['Data [CSVM,1]','MC [CSVM,1]'])
#make_plot('comp_datamc2_mistag',['HEPTagger/Mistag/data_htt/Mistag_px3_data_htt','HEPTagger/Mistag/qcd_htt/Mistag_px3_qcd_htt'],['Data [CSVL,CSVM]','MC [CSVL,CSVM]'])
#make_plot('comp_datamc3_mistag',['HEPTagger/Mistag/data_htt/Mistag_px2_data_htt','HEPTagger/Mistag/qcd_htt/Mistag_px2_qcd_htt'],['Data [0,CSVL]','MC [0,CSVL]'])
#make_plot('comp_shape',['HEPTagger/MassShape/qcd_htt/mass_shape_qcd_htt','HEPTagger/MassShape/data_htt/mass_shape_data_htt'],['QCD','DATA'],True)
#make_plot('comp_shape_jec',['HEPTagger/MassShape/qcd_htt/mass_shape_qcd_htt','HEPTagger/MassShape/qcd_htt_jec/mass_shapeJEC_qcd_htt_jec'],['NO JEC','JEC'],True)
#make_plot('comp_qcd_shape',['HEPTagger/MassShape2D/qcd_htt/mass_shape2D_px4_qcd_htt','HEPTagger/MassShape2D/qcd_htt/mass_shape2D_px3_qcd_htt','HEPTagger/MassShape2D/qcd_htt/mass_shape2D_px2_qcd_htt'],['[CSVM,1]','[CSVL,CSVM]','[0,CSVL]','[-1,0]'],True)
#make_plot('comp_data_shape',['HEPTagger/MassShape2D/data_htt/mass_shape2D_px4_data_htt','HEPTagger/MassShape2D/data_htt/mass_shape2D_px3_data_htt','HEPTagger/MassShape2D/data_htt/mass_shape2D_px2_data_htt'],['[CSVM,1]','[CSVL,CSVM]','[0,CSVL]','[-1,0]'],True)
print 'close files'
ttbar_file.Close()
qcd_file.Close()
data_file.Close()
outfile.Close()