From 683dd7596e82a9ad1dca7b3f8d5506552508545c Mon Sep 17 00:00:00 2001 From: Matthew Kirk Date: Fri, 23 Oct 2020 17:27:11 +0200 Subject: [PATCH] Add citations for W and Z observables --- flavio/physics/wdecays/mw.py | 2 ++ flavio/physics/zdecays/gammazsm.py | 2 ++ flavio/physics/zdecays/smeftew.py | 15 +++++++++++++++ 3 files changed, 19 insertions(+) diff --git a/flavio/physics/wdecays/mw.py b/flavio/physics/wdecays/mw.py index c6c0f985..16febf7e 100644 --- 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 index ec2a5dc9..23a68b22 100644 --- 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 index 3e2b1712..8e157688 100644 --- 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)