added higgs processes and first di-higgs processes

uhh-cms · Feb 20, 2024 · 1146c66 · 1146c66
1 parent 04a894c
commit 1146c66
Showing 1 changed file with 48 additions and 20 deletions.
diff --git 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",