예제 #1
0
\hline
Total Background & {incTotal} & {boostTotal} & {vbfTotal} \\
\hline
H$\to\tau\tau$ & {incHTT} & {boostHTT} & {vbfHTT} \\
\hline
Data & {incData} & {boostData} & {vbfData} \\
'''

# Define which sets of bins corresponds to table categories
mt_0jet_bins = ['htt_mt_0_*', 'htt_mt_1_*']
mt_boost_bins = ['htt_mt_2_*', 'htt_mt_3_*']
mt_btag_bins = ['htt_mt_6_*', 'htt_mt_7_*']

tables['mt'] = {
    'incZTT':
    dc.get_rate(mt_0jet_bins, 'ZTT'),
    'boostZTT':
    dc.get_rate(mt_boost_bins, 'ZTT'),
    'vbfZTT':
    dc.get_rate(mt_btag_bins, 'ZTT'),
    'incQCD':
    dc.get_rate(mt_0jet_bins, 'QCD'),
    'boostQCD':
    dc.get_rate(mt_boost_bins, 'QCD'),
    'vbfQCD':
    dc.get_rate(mt_btag_bins, 'QCD'),
    'incWjets':
    dc.get_rate(mt_0jet_bins, 'W'),
    'boostWjets':
    dc.get_rate(mt_boost_bins, 'W'),
    'vbfWjets':
예제 #2
0
\hline
Total Background & {incTotal} & {boostTotal} & {vbfTotal} \\
\hline
H$\to\tau\tau$ & {incHTT} & {boostHTT} & {vbfHTT} \\
\hline
Data & {incData} & {boostData} & {vbfData} \\
'''

# Define which sets of bins corresponds to table categories
mt_0jet_bins = ['htt_mt_0_*', 'htt_mt_1_*']
mt_boost_bins = ['htt_mt_2_*', 'htt_mt_3_*']
mt_vbf_bins = ['htt_mt_5_*']

tables['mt'] = {
    'incZTT':
    dc.get_rate(mt_0jet_bins, 'ZTT'),
    'boostZTT':
    dc.get_rate(mt_boost_bins, 'ZTT'),
    'vbfZTT':
    dc.get_rate(mt_vbf_bins, 'ZTT'),
    'incQCD':
    dc.get_rate(mt_0jet_bins, 'QCD'),
    'boostQCD':
    dc.get_rate(mt_boost_bins, 'QCD'),
    'vbfQCD':
    dc.get_rate(mt_vbf_bins, 'QCD'),
    'incWjets':
    dc.get_rate(mt_0jet_bins, 'W'),
    'boostWjets':
    dc.get_rate(mt_boost_bins, 'W'),
    'vbfWjets':
예제 #3
0
\hline
\hline
Total Background & {incTotal} & {boostTotal} & {vbfTotal} \\
\hline
H$\to\tau\tau$ & {incHTT} & {boostHTT} & {vbfHTT} \\
\hline
Data & {incData} & {boostData} & {vbfData} \\
"""

# Define which sets of bins corresponds to table categories
mt_0jet_bins = ["htt_mt_0_*", "htt_mt_1_*"]
mt_boost_bins = ["htt_mt_2_*", "htt_mt_3_*"]
mt_vbf_bins = ["htt_mt_5_*"]

tables["mt"] = {
    "incZTT": dc.get_rate(mt_0jet_bins, "ZTT"),
    "boostZTT": dc.get_rate(mt_boost_bins, "ZTT"),
    "vbfZTT": dc.get_rate(mt_vbf_bins, "ZTT"),
    "incQCD": dc.get_rate(mt_0jet_bins, "QCD"),
    "boostQCD": dc.get_rate(mt_boost_bins, "QCD"),
    "vbfQCD": dc.get_rate(mt_vbf_bins, "QCD"),
    "incWjets": dc.get_rate(mt_0jet_bins, "W"),
    "boostWjets": dc.get_rate(mt_boost_bins, "W"),
    "vbfWjets": dc.get_rate(mt_vbf_bins, "W"),
    "incZJ": dc.get_rate(mt_0jet_bins, "ZJ") + dc.get_rate(mt_0jet_bins, "ZL"),
    "boostZJ": dc.get_rate(mt_boost_bins, "ZJ") + dc.get_rate(mt_boost_bins, "ZL"),
    "vbfZJ": dc.get_rate(mt_vbf_bins, "ZLL"),
    "incTT": dc.get_rate(mt_0jet_bins, "TT"),
    "boostTT": dc.get_rate(mt_boost_bins, "TT"),
    "vbfTT": dc.get_rate(mt_vbf_bins, "TT"),
    "incVV": dc.get_rate(mt_0jet_bins, "VV"),
예제 #4
0
for channel in channels:

    bins[channel] = [bin for bin in allbins if '_'+channel in bin]
    chBgs = []

    tables[channel] = {}
    uncs[channel] = {}
    sOverSplusBs[channel] = {}
    sigmas[channel] = {}
    widths[channel] = {}

    for cat in bins[channel]:
        for proc in procs:
            try:
                tables[channel][cat+proc] = dc.get_rate(cat, proc).nominal_value
                uncs[channel][cat+proc] = dc.get_rate(cat, proc).std_dev()
            except KeyError:
                pass
            else:
                if proc in backgrounds:
                    chBgs.append(proc)
        if not POSTFIT and channel not in ['ee', 'mm']:
            sOverSplusBs[channel][cat] = getSOverSplusB(catToFile(cat, channel))
            sigma, width = getStdDeviation(catToFile(cat, channel))
            # sigma, width = getStdDeviationDC(catToFileDC(dc, cat, channel))
            sigmas[channel][cat] = sigma
            widths[channel][cat] = width
        tables[channel][cat+'obs'] = dc.get_obs(cat)

    chBgs = set(chBgs)
\hline
\hline
Total Background & {incTotal} & {boostTotal} & {vbfTotal} \\
\hline
H$\to\tau\tau$ & {incHTT} & {boostHTT} & {vbfHTT} \\
\hline
Data & {incData} & {boostData} & {vbfData} \\
'''

# Define which sets of bins corresponds to table categories
mt_0jet_bins = ['htt_mt_0_*', 'htt_mt_1_*']
mt_boost_bins = ['htt_mt_2_*', 'htt_mt_3_*']
mt_btag_bins = ['htt_mt_6_*', 'htt_mt_7_*']

tables['mt'] = {
    'incZTT' : dc.get_rate(mt_0jet_bins, 'ZTT'),
    'boostZTT' : dc.get_rate(mt_boost_bins, 'ZTT'),
    'vbfZTT' : dc.get_rate(mt_btag_bins, 'ZTT'),

    'incQCD' : dc.get_rate(mt_0jet_bins, 'QCD'),
    'boostQCD' : dc.get_rate(mt_boost_bins, 'QCD'),
    'vbfQCD' : dc.get_rate(mt_btag_bins, 'QCD'),

    'incWjets' : dc.get_rate(mt_0jet_bins, 'W'),
    'boostWjets' : dc.get_rate(mt_boost_bins, 'W'),
    'vbfWjets' : dc.get_rate(mt_btag_bins, 'W'),

    'incZJ' : dc.get_rate(mt_0jet_bins, 'ZJ') + dc.get_rate(mt_0jet_bins, 'ZL'),
    'boostZJ' : dc.get_rate(mt_boost_bins, 'ZJ') + dc.get_rate(mt_boost_bins, 'ZL'),
    'vbfZJ' : dc.get_rate(mt_btag_bins, 'ZLL'),