print '\n'
if doYields:
### Initialize table maker ###
tableFile = file('yields/.yields_tmp.tex', 'w')
yieldTable = TableMaker('fcncAnalysis/combined_histos/{0}_cut1_{1}_{2}.root'.format(selection, period, batch),
tableFile, scale = LUMIDATA,
delimiter = '&', doSumBG = True)
yieldTable.set_period(period)
yieldTable.add_datasets(samples['all'], Clear = True)
if not doPlots:
#yieldTable.get_scale_factors()
yieldTable.get_scale_factors(['FCNH', 'FCNHUp'])
if do3l:
yieldTable._columnList = ['Rare', 'WZJets3LNu', 'Fakes', 'BG', 'DATA', 'FCNH']
#yieldTable._columnList = ['Rare', 'WZJets3LNu', 'Fakes', 'BG', 'DATA', 'FCNHWW', 'FCNHTauTau', 'FCNHZZ']
#yieldTable._columnList = ['Rare', 'WZJets3LNu', 'Fakes', 'BG', 'DATA', 'FCNHUpWW', 'FCNHUpZZ', 'FCNHUpTauTau']
#yieldTable._columnList = ['Rare', 'WZJets3LNu', 'Fakes', 'BG', 'DATA', 'TTH_M-125']
#yieldTable._columnList = samples['3l_inclusive'] + ['BG', 'DATA', 'FCNH']
#yieldTable._columnList = ['BG', 'DATA', 'FCNHWW', 'FCNHZZ', 'FCNHTauTau']
#yieldTable._columnList = ['BG', 'DATA', 'FCNH']#, 'Significance']
yieldTable.add_datasets(samples['3l_inclusive'], Clear = True)
#yieldTable.add_datasets('TTH_M-125')
yieldTable.add_datasets('FCNH')
#yieldTable.add_datasets('FCNHWW')
#yieldTable.add_datasets('FCNHZZ')
#yieldTable.add_datasets('FCNHTauTau')
yieldTable.add_datasets('DATA')
print '\n'
####################
### MAKE TABLES! ###
####################
if doYields:
doPresel = True
outFile = file('yields/.yields_tmp.tex', 'w')
categoryNames = []
yieldTable = TableMaker('fcncAnalysis/combined_histos/' + selection + '_cut1_' + period + batch + '.root', outFile, scale = LUMIDATA, delimiter = '&', doSumBG = True)
yieldTable.set_period(period)
yieldTable._columnList = ['higgs', 'Triboson', 'ttV', 'Diboson', 'top', 'VJets', 'BG', 'DATA', 'FCNH']#, 'Significance']
#yieldTable._columnList = ['BG', 'DATA', 'FCNH']#, 'Significance']
yieldTable.add_datasets(samples, Clear = True)
yieldTable.add_datasets('FCNH')
yieldTable.add_datasets('DATA')
print '\n\n Printing yields...\n'
if doOS:
categoryNames.extend(catOS)
if doSS:
categoryNames.extend(catSS)
if do3l:
categoryNames.extend(cat3l)
bTable = TableMaker('../HiggsAnalyzer/histos/higgsHistograms_'+dataSet+period+'.root', LUMIDATA)
samples = ['tW', 'ttbar', 'GluGluWW', 'WWJets', 'WZJets', 'ZZJets']
if selection in ['muGamma', 'eGamma']:
samples.append('PhotonJets')
elif selection == 'muEG':
samples.append('DYToTauTau')
elif doZJets:
samples.append('ZJets')
elif not doZJets:
if selection == 'muon':
samples.append('DYToMuMu')
samples.append('DYToTauTau')
if selection == 'electron':
samples.append('DYToEE')
samples.append('DYToTauTau')
bTable.AddDatasets(samples)
bTable.AddDatasets('DATA')
bTable._rowList = ['0 b-jet', '1 b-jet', '2 b-jet', '3 b-jet']
if selection in ['eGamma', 'muGamma']:
bTable._columnList = ['ZZJets', 'WZJets', 'WWJets', 'ttbar', 'tW', 'PhotonJets', 'BG', 'DATA']
if selection in ['electron', 'muon']:
bTable._columnList = ['ZZJets', 'WZJets', 'WWJets', 'ttbar', 'tW', 'ZJets', 'BG', 'DATA']
histDict = bTable.GetHistDict('Jet', 'bJetMultNoVeto')
bTable.PrintTable(histDict, True, True, startBin = 1)
