[f_bc2e3080, xs_bc2e3080, n_bc2e3080, lumi],
[f_bc2e80170, xs_bc2e80170, n_bc2e80170, lumi],
[f_gjet40100, xs_gjet40100, n_gjet40100, lumi],
[f_gjet100200, xs_gjet100200, n_gjet100200,lumi],
[f_gjet200Inf, xs_gjet200Inf, n_gjet200Inf,lumi],
]
#======= add the qcd and singleTop histos and apply weights to all ============
for ivar in range(0, len(names)) :
name = names[ivar]
hist_data = weightHisto( dataLabel, tempDir, name, dataSample, kBlack, verbose)
#hist_tp = weightHisto( tprimeLabel, tempDir, name, tprimeSample, 6, verbose)
hist_top = weightHisto( topLabel, tempDir, name, topSample, 206, verbose)
hist_wjets = weightHisto( wjetsLabel, tempDir, name, wjetSample, 210, verbose)
hist_zjets = weightHisto( zjetsLabel, tempDir, name, zjetSample, 215, verbose)
hist_sTop = addHistos( sTopLabel, tempDir, name, singleTopSamples,95, verbose)
hist_qcd = addHistos( qcdLabel, tempDir, name, qcdSamples, 220, verbose)
print 'data area', hist_data.Integral()
allHists.append(hist_data)
allHists.append(hist_qcd)
allHists.append(hist_sTop)
allHists.append(hist_zjets)
allHists.append(hist_wjets)
allHists.append(hist_top)
#allHists.append(hist_tp)
# ============= PLOT ==================
hs = THStack("nEvents","nEvents")
for ihist in allHists[1:]:
ihist.Rebin(options.rebin)
xtitle = "Number of Primary vertices"
elif var == 'nJets3':
xtitle = "number of jets"
else:
xtitle = ""
#======= add the qcd and singleTop histos and apply weights to all ============
for ivar in range(0, len(names)) :
name = names[ivar]
if not useFitter:
hist_data = weightHisto( dataLabel, name, dataSample, kBlack, verbose)
hist_top = weightHisto( topLabel, name, topSample, 206, verbose)
hist_wjets = weightHisto( wjetsLabel, name, wjetSample, 210, verbose)
hist_zjets = weightHisto( zjetsLabel, name, zjetSample, 215, verbose)
hist_sTop = addHistos( sTopLabel, name, singleTopSamples,95, verbose)
hist_qcd = addHistos( qcdLabel, name, qcdSamples, 220, verbose)
if useLooseQCD:
hist_qcdL = addHistos( qcdLabelLoose, name, qcdSamplesLoose, 220, verbose)
if useLooseQCD and hist_qcdL.Integral()!= 0:
hist_qcdL.Scale(hist_qcd.Integral()/hist_qcdL.Integral())
else:
hist_data = addHistosMulti( dataLabel, name, dataSample, kBlack, verbose, None, minTag, maxTag, minJet, maxJet)
hist_top = addHistosMulti( topLabel, name, topSample, 206, verbose, Top_total, minTag, maxTag, minJet, maxJet)
hist_wjets = addHistosMulti( wjetsLabel, name, wjetSample, 210, verbose, WJet_total, minTag, maxTag, minJet, maxJet)
hist_zjets = addHistosMulti( zjetsLabel, name, zjetSample, 215, verbose, ZJet_total, minTag, maxTag, minJet, maxJet)
hist_sTop = addHistosMulti( sTopLabel, name, singleTopSamples,95, verbose, SingleTop_total, minTag, maxTag, minJet, maxJet)
hist_qcd = addHistosMulti( qcdLabel, name, qcdSamples, 220, verbose, QCD_total, minTag, maxTag, minJet, maxJet)
if useLooseQCD:
