Module flavio.physics.zdecays.gammazsm

$Z$ partial widths in the SM.

Based on arXiv:1401.2447

Functions

def GammaZ_SM(par, f)

Expand source code

def GammaZ_SM(par, f):
    if f in ['e', 'mu']:
        name = 'Gammae,mu'
    elif f in ['d', 's']:
        name = 'Gammad,s'
    elif 'nu' in f:
        name = 'Gammanu'
    else:
        name = 'Gamma' + f
    GSM = Zobs(name, par['m_h'], par['m_t'], par['alpha_s'], 0.059, par['m_Z'])
    return GSM + par['delta_' + name]

def Zobs(name, m_h, m_t, alpha_s, Dalpha, m_Z)

Expand source code

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
    D_alpha = Dalpha / 0.059 - 1
    D_Z = m_Z / 91.1876 - 1
    c = cdict[name]
    return (c[0] + c[1] * L_H + c[2] * D_t + c[3] * D_alpha_s
            + c[4] * D_alpha_s**2 + c[5] * D_alpha_s * D_t
            + c[6] * D_alpha + c[7] * D_Z) * units[name]

Expansion formula for $Z$ partial widths according to eq. (28) of arXiv:1401.2447.