diff --git a/cmsdb/processes/higgs.py b/cmsdb/processes/higgs.py index 0a4c9ff2..77fa4630 100644 --- a/cmsdb/processes/higgs.py +++ b/cmsdb/processes/higgs.py @@ -23,10 +23,14 @@ import cmsdb.constants as const + + # # Single Higgs # +# https://twiki.cern.ch/twiki/bin/view/LHCPhysics/CERNYellowReportPageAt13TeV?rev=24 + h = Process( name="h", id=10000, @@ -37,28 +41,34 @@ name="h_ggf", id=11000, label=r"$H_{ggf}$", - xsecs={13: Number(0.1)}, # TODO -) + xsecs={13: Number(4.858E+01,{ + "pdf" : 3.2j, + "th": (4.6j, 6.7j), + "th_gaussian":3.9j, + "th_gaussian" : 3.2j, + },)}) + h_ggf_tautau = h_ggf.add_process( name="h_ggf_tautau", id=11100, - xsecs={13: Number(0.1)}, # TODO + xsecs={13: h_ggf.xsecs["13"] * const.br_h.tt,}, # TODO ) h_vbf = h.add_process( name="h_vbf", id=12000, - label=r"$H_{vbf}$", - xsecs={13: Number(0.1)}, # TODO -) + label=r"$H_{vbf}$", xsecs={13: Number(3.782E+00),{ + "pdf":2.1j, + "scale":(0.4j,0.3j)}}) h_vbf_tautau = h_vbf.add_process( name="h_vbf_tautau", id=12100, - xsecs={13: Number(0.1)}, # TODO + xsecs={13: h_vbf.xsecs["13"] * const.br_h.tt}, # TODO ) +# empty, no value given vh = h.add_process( name="vh", id=13000, @@ -70,7 +80,10 @@ name="zh", id=13100, label="ZH", - xsecs={13: Number(0.1)}, # TODO + xsecs={13: Number(8.839E-01),{ + "scale":(3.8j,3.1j), + "pdf":1.6j, + }}, # TODO ) zh_tautau = zh.add_process( @@ -97,44 +110,51 @@ wph = vh.add_process( name="wph", id=13200, - xsecs={13: Number(0.1)}, # TODO + xsecs={13: Number(8.400E-01),{ + "pdf":1.9j, + "scale":(0.5j,0.7j) + }}, ) wph_tautau = wph.add_process( name="wph_tautau", id=13210, - xsecs={13: Number(0.1)}, # TODO + xsecs={13: wph.get_xsec(13) * const.br_h.tt}, ) wmh = vh.add_process( name="wmh", id=13300, - xsecs={13: Number(0.1)}, # TODO + xsecs={13: Number(5.328E-01 ),{ + "pdf":1.9j, + "scale":(0.5j,0.7j) + }}, ) + wmh_tautau = wmh.add_process( name="wmh_tautau", id=13310, - xsecs={13: Number(0.1)}, # TODO + xsecs={13: wmh.get_xsec(13) * const.br_h.tt}, # TODO ) ggzh = vh.add_process( name="ggzh", id=14000, - xsecs={13: Number(0.1)}, # TODO + xsecs={13: Number(0.1227 ), {"scale":(25.1j,18.9j),"pdf":1.9j}}, # TODO ) ggzh_llbb = ggzh.add_process( name="ggzh_llbb", id=14100, - xsecs={13: Number(0.1)}, # TODO + xsecs={13: Number(ggzh.get_xsec(13) * const.br_h.bb * const.br_z.clep)}, # TODO ) tth = h.add_process( name="tth", id=15000, label=r"$t\bar{t}H$", - xsecs={13: Number(0.1)}, # TODO + xsecs={13: Number( 5.071E-01 ), {"scale":(5.8j,9.2j),"pdf":3.6j}}, # TODO ) tth_tautau = tth.add_process( @@ -163,18 +183,25 @@ # Basic HH processes # + +# https://twiki.cern.ch/twiki/bin/view/LHCPhysics/LHCHWGHH?rev=90#Current_recommendations_for_HH_c +# scale is according to recommendation: scale + mtop unc +# pdf is according to recommendation: pdf + aS + hh = Process( name="hh", id=20000, label="HH", - xsecs={13: Number(0.1)}, # TODO -) + xsecs={13: Number(0.1)}) # TODO hh_ggf = hh.add_process( name="hh_ggf", id=21000, label=r"$HH_{ggf}$", - xsecs={13: Number(0.1)}, # TODO + xsecs={13: Number(31.05, + {"pdf":3.0j, "scale":(6j,23j)}), + 14: Number(34.43, + {"pdf":3.0j, "scale":(6j,23j)})}, # fb ) # Naming conventions, cross sections and uncertainties are based on: @@ -230,12 +257,13 @@ }, ) +#https://twiki.cern.ch/twiki/bin/view/LHCPhysics/LHCHWGHH?rev=90#HHjj_VBF + hh_vbf = hh.add_process( name="hh_vbf", id=22000, label=r"$HH_{vbf}$", - xsecs={13: Number(0.1)}, # TODO -) + xsecs={13: Number(1.726, {"scale":(0.03j, 0.04j), "pdf":2.1j})}) qqHH_CV_1_C2V_1_kl_1 = hh_vbf.add_process( name="qqHH_CV_1_C2V_1_kl_1",