Add citations for W and Z observables

flav-io · Oct 23, 2020 · 683dd75 · 683dd75
1 parent 8be3df2
commit 683dd75
Show file tree

Hide file tree

Showing 3 changed files with 19 additions and 0 deletions.
diff --git a/flavio/physics/wdecays/mw.py b/flavio/physics/wdecays/mw.py
@@ -10,6 +10,7 @@ def mW_SM(par):
     r"""Two-loop SM prediction for the $W$ pole mass.
 
     Eq. (6) of arXiv:hep-ph/0311148."""
+    flavio.citations.register("Awramik:2003rn")
     dH = log(par['m_h'] / 100)
     dh = (par['m_h'] / 100)**2
     dt = (par['m_t'] / 174.3)**2 - 1
@@ -36,6 +37,7 @@ def dmW_SMEFT(par, C):
     r"""Shift in the $W$ mass due to dimension-6 operators.
 
     Eq. (2) of arXiv:1606.06502."""
+    flavio.citations.register("Bjorn:2016zlr")
     sh2 = smeftew._sinthetahat2(par)
     sh = sqrt(sh2)
     ch2 = 1 - sh2

diff --git a/flavio/physics/zdecays/gammazsm.py b/flavio/physics/zdecays/gammazsm.py
@@ -5,6 +5,7 @@
 
 from math import log
 from scipy import constants
+import flavio
 
 
 # units: GeV=hbar=c=1
@@ -52,6 +53,7 @@ def Zobs(name, m_h, m_t, alpha_s, Dalpha, m_Z):
     r"""Expansion formula for $Z$ partial widths according to eq. (28) of
     arXiv:1401.2447.
     """
+    flavio.citations.register("Freitas:2014hra")
     L_H = log(m_h / 125.7)
     D_t = (m_t / 173.2)**2 - 1
     D_alpha_s = alpha_s / 0.1184 - 1

diff --git a/flavio/physics/zdecays/smeftew.py b/flavio/physics/zdecays/smeftew.py
@@ -4,6 +4,7 @@
 Taken from Appendix A.1 in arXiv:1706.08945."""
 
 from math import sqrt, pi
+import flavio
 
 
 # fermion quantum numbers
@@ -21,14 +22,17 @@
 
 
 def d_GF(par, C):
+    flavio.citations.register("Brivio:2017vri")
     return 1 / (sqrt(2) * par['GF']) * sqrt(2) / 2 * (C['phil3_22'] + C['phil3_11'] - C['ll_1221'] / 2)
 
 
 def _sinthetahat2(par):
     r"""$sin^2(\hat\theta)$"""
+    flavio.citations.register("Brivio:2017vri")
     return 1 / 2 * (1 - sqrt(1 - (4 * pi * par['alpha_e']) / (sqrt(2) * par['GF'] * par['m_Z']**2)))
 
 def d_sh2(par, C):
+    flavio.citations.register("Brivio:2017vri")
     sh2 = _sinthetahat2(par)
     sh = sqrt(sh2)
     ch2 = 1 - sh2
@@ -40,12 +44,14 @@ def d_sh2(par, C):
 
 
 def d_mZ2(par, C):
+    flavio.citations.register("Brivio:2017vri")
     return (1 / (2 * sqrt(2)) * par['m_Z']**2 / par['GF'] * C['phiD']
             + (2**(1 / 4) * sqrt(pi * par['alpha_e']) * par['m_Z'])
             / par['GF']**(3 / 2) * C['phiWB'])
 
 
 def d_gZb(par, C):
+    flavio.citations.register("Brivio:2017vri")
     sh2 = _sinthetahat2(par)
     sh = sqrt(sh2)
     ch = sqrt(1 - sh2)
@@ -70,6 +76,7 @@ def gA_SM(f, par):
 
 
 def d_gVl(f1, f2, par, C):
+    flavio.citations.register("Brivio:2017vri")
     i, j = sorted((_sectors['l'][f1], _sectors['l'][f2]))
     d_g = -(C['phie_{}{}'.format(i, j)] + C['phil1_{}{}'.format(i, j)] + C['phil3_{}{}'.format(i, j)]) / (4 * sqrt(2) * par['GF'])
     if i > j:
@@ -80,6 +87,7 @@ def d_gVl(f1, f2, par, C):
 
 
 def d_gVnu(f1, f2, par, C):
+    flavio.citations.register("Brivio:2017vri")
     i, j = sorted((_sectors['nu'][f1], _sectors['nu'][f2]))
     d_g = -(C['phil1_{}{}'.format(i, j)] - C['phil3_{}{}'.format(i, j)]) / (4 * sqrt(2) * par['GF'])
     if i > j:
@@ -90,6 +98,7 @@ def d_gVnu(f1, f2, par, C):
 
 
 def d_gVu(f1, f2, par, C):
+    flavio.citations.register("Brivio:2017vri")
     i, j = sorted((_sectors['u'][f1], _sectors['u'][f2]))
     d_g = (-C['phiu_{}{}'.format(i, j)] - C['phiq1_{}{}'.format(i, j)] + C['phiq3_{}{}'.format(i, j)]) / (4 * sqrt(2) * par['GF'])
     if i > j:
@@ -100,6 +109,7 @@ def d_gVu(f1, f2, par, C):
 
 
 def d_gVd(f1, f2, par, C):
+    flavio.citations.register("Brivio:2017vri")
     i, j = sorted((_sectors['d'][f1], _sectors['d'][f2]))
     d_g = -(C['phid_{}{}'.format(i, j)] + C['phiq1_{}{}'.format(i, j)] + C['phiq3_{}{}'.format(i, j)]) / (4 * sqrt(2) * par['GF'])
     if i > j:
@@ -110,6 +120,7 @@ def d_gVd(f1, f2, par, C):
 
 
 def d_gAl(f1, f2, par, C):
+    flavio.citations.register("Brivio:2017vri")
     i, j = sorted((_sectors['l'][f1], _sectors['l'][f2]))
     d_g = (C['phie_{}{}'.format(i, j)] - C['phil1_{}{}'.format(i, j)] - C['phil3_{}{}'.format(i, j)]) / (4 * sqrt(2) * par['GF'])
     if i > j:
@@ -124,6 +135,7 @@ def d_gAnu(f1, f2, par, C):
 
 
 def d_gAu(f1, f2, par, C):
+    flavio.citations.register("Brivio:2017vri")
     i, j = sorted((_sectors['u'][f1], _sectors['u'][f2]))
     d_g = -(-C['phiu_{}{}'.format(i, j)] + C['phiq1_{}{}'.format(i, j)] - C['phiq3_{}{}'.format(i, j)]) / (4 * sqrt(2) * par['GF'])
     if i > j:
@@ -134,6 +146,7 @@ def d_gAu(f1, f2, par, C):
 
 
 def d_gAd(f1, f2, par, C):
+    flavio.citations.register("Brivio:2017vri")
     i, j = sorted((_sectors['d'][f1], _sectors['d'][f2]))
     d_g = (C['phid_{}{}'.format(i, j)] - C['phiq1_{}{}'.format(i, j)] - C['phiq3_{}{}'.format(i, j)]) / (4 * sqrt(2) * par['GF'])
     if i > j:
@@ -166,6 +179,7 @@ def d_gA(f1, f2, par, C):
 
 
 def d_gWl(f1, f2, par, C):
+    flavio.citations.register("Brivio:2017vri")
     i, j = sorted((_sectors['l'][f1], _sectors['l'][f2]))
     sh2 = _sinthetahat2(par)
     sh = sqrt(sh2)
@@ -175,6 +189,7 @@ def d_gWl(f1, f2, par, C):
 
 
 def d_gWq(f1, f2, par, C):
+    flavio.citations.register("Brivio:2017vri")
     i, j = sorted((_sectors['u'][f1], _sectors['d'][f2]))
     sh2 = _sinthetahat2(par)
     sh = sqrt(sh2)