Python DataCard.get_obsの例

コード例 #1

ファイル: mssm_h2tau_tables.py プロジェクト: EricBAdamsUMDCP/cmssw-1

    'vbfTT':
    dc.get_rate(mt_btag_bins, 'TT'),
    'incVV':
    dc.get_rate(mt_0jet_bins, 'VV'),
    'boostVV':
    dc.get_rate(mt_boost_bins, 'VV'),
    'vbfVV':
    dc.get_rate(mt_btag_bins, 'VV'),
    'incHTT':
    dc.get_rate(mt_0jet_bins, 'ggH') + dc.get_rate(mt_0jet_bins, 'bbH'),
    'boostHTT':
    dc.get_rate(mt_boost_bins, 'ggH') + dc.get_rate(mt_boost_bins, 'bbH'),
    'vbfHTT':
    dc.get_rate(mt_btag_bins, 'ggH') + dc.get_rate(mt_btag_bins, 'bbH'),
    'incData':
    dc.get_obs(mt_0jet_bins),
    'boostData':
    dc.get_obs(mt_boost_bins),
    'vbfData':
    dc.get_obs(mt_btag_bins),
}
add_up_backgrounds(tables['mt'])

dc.print_systematics(mt_0jet_bins, 'QCD'),

with open('mt_yields_table.tex', 'w') as file:
    file.write(et_mt_template.format(**render_values(tables['mt'])))

###### ET Channel

et_0jet_bins = ['htt_et_0_*', 'htt_et_1_*']

コード例 #2

ファイル: h2tau_tables.py プロジェクト: EricBAdamsUMDCP/cmssw-1

    dc.get_rate(mt_0jet_bins, 'VV'),
    'boostVV':
    dc.get_rate(mt_boost_bins, 'VV'),
    'vbfVV':
    dc.get_rate(mt_vbf_bins, 'VV'),
    'incHTT':
    dc.get_rate(mt_0jet_bins, 'ggH') + dc.get_rate(mt_0jet_bins, 'qqH') +
    dc.get_rate(mt_0jet_bins, 'VH'),
    'boostHTT':
    dc.get_rate(mt_boost_bins, 'ggH') + dc.get_rate(mt_boost_bins, 'qqH') +
    dc.get_rate(mt_boost_bins, 'VH'),
    'vbfHTT':
    dc.get_rate(mt_vbf_bins, 'ggH') + dc.get_rate(mt_vbf_bins, 'qqH') +
    dc.get_rate(mt_vbf_bins, 'VH'),
    'incData':
    dc.get_obs(mt_0jet_bins),
    'boostData':
    dc.get_obs(mt_boost_bins),
    'vbfData':
    dc.get_obs(mt_vbf_bins),
}
add_up_backgrounds(tables['mt'])

dc.print_systematics(mt_0jet_bins, 'QCD'),

with open('mt_yields_table.tex', 'w') as file:
    file.write(et_mt_template.format(**render_values(tables['mt'])))

###### ET Channel

et_0jet_bins = ['htt_et_0_*', 'htt_et_1_*']

コード例 #3

ファイル: h2tau_tables.py プロジェクト: nsahoo/cmssw-1

    "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"),
    "boostVV": dc.get_rate(mt_boost_bins, "VV"),
    "vbfVV": dc.get_rate(mt_vbf_bins, "VV"),
    "incHTT": dc.get_rate(mt_0jet_bins, "ggH") + dc.get_rate(mt_0jet_bins, "qqH") + dc.get_rate(mt_0jet_bins, "VH"),
    "boostHTT": dc.get_rate(mt_boost_bins, "ggH")
    + dc.get_rate(mt_boost_bins, "qqH")
    + dc.get_rate(mt_boost_bins, "VH"),
    "vbfHTT": dc.get_rate(mt_vbf_bins, "ggH") + dc.get_rate(mt_vbf_bins, "qqH") + dc.get_rate(mt_vbf_bins, "VH"),
    "incData": dc.get_obs(mt_0jet_bins),
    "boostData": dc.get_obs(mt_boost_bins),
    "vbfData": dc.get_obs(mt_vbf_bins),
}
add_up_backgrounds(tables["mt"])

dc.print_systematics(mt_0jet_bins, "QCD"),

with open("mt_yields_table.tex", "w") as file:
    file.write(et_mt_template.format(**render_values(tables["mt"])))

###### ET Channel

et_0jet_bins = ["htt_et_0_*", "htt_et_1_*"]
et_boost_bins = ["htt_et_2_*", "htt_et_3_*"]
et_vbf_bins = ["htt_et_5_*"]

コード例 #4

        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)

    # RESCALE LOOP
    if POSTFIT:
        for cat in bins[channel]:
            if 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
            scaleName = cat.replace('htt', 'scales')
            scales = __import__(scaleName).scales
            for proc in procs:

コード例 #5

ファイル: mssm_h2tau_tables.py プロジェクト: JehadMousa/HiggsAnalysis-HiggsToTauTau

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

    'incTT' : dc.get_rate(mt_0jet_bins, 'TT'),
    'boostTT' : dc.get_rate(mt_boost_bins, 'TT'),
    'vbfTT' : dc.get_rate(mt_btag_bins, 'TT'),

    'incVV' : dc.get_rate(mt_0jet_bins, 'VV'),
    'boostVV' : dc.get_rate(mt_boost_bins, 'VV'),
    'vbfVV' : dc.get_rate(mt_btag_bins, 'VV'),

    'incHTT' : dc.get_rate(mt_0jet_bins, 'ggH') + dc.get_rate(mt_0jet_bins, 'bbH'),
    'boostHTT' : dc.get_rate(mt_boost_bins, 'ggH') + dc.get_rate(mt_boost_bins, 'bbH'),
    'vbfHTT' : dc.get_rate(mt_btag_bins, 'ggH') + dc.get_rate(mt_btag_bins, 'bbH'),

    'incData' : dc.get_obs(mt_0jet_bins),
    'boostData' : dc.get_obs(mt_boost_bins),
    'vbfData' : dc.get_obs(mt_btag_bins),
}
add_up_backgrounds(tables['mt'])

dc.print_systematics(mt_0jet_bins, 'QCD'),

with open('mt_yields_table.tex', 'w') as file:
    file.write(et_mt_template.format(**render_values(tables['mt'])))

###### ET Channel

et_0jet_bins = ['htt_et_0_*', 'htt_et_1_*']
et_boost_bins = ['htt_et_2_*', 'htt_et_3_*']
et_btag_bins = ['htt_et_6_*', 'htt_et_7_*']