Beispiel #1
0
def get_wceff(q2, wc, par, B, M, lep, scale):
    r"""Get a dictionary with the effective $\Delta F=1$ Wilson coefficients
    in the convention appropriate for the generalized angular distributions.
    """
    xi_u = ckm.xi('u',meson_quark[(B,M)])(par)
    xi_t = ckm.xi('t',meson_quark[(B,M)])(par)
    qiqj=meson_quark[(B,M)]
    Yq2 = matrixelements.Y(q2, wc, par, scale, qiqj) + (xi_u/xi_t)*matrixelements.Yu(q2, wc, par, scale, qiqj)
        #   b) NNLO Q1,2
    delta_C7 = matrixelements.delta_C7(par=par, wc=wc, q2=q2, scale=scale, qiqj=qiqj)
    delta_C9 = matrixelements.delta_C9(par=par, wc=wc, q2=q2, scale=scale, qiqj=qiqj)
    mb = running.get_mb(par, scale)
    ll = lep + lep
    c = {}
    c['7']  = wc['C7eff_'+qiqj]      + delta_C7
    c['7p'] = wc['C7effp_'+qiqj]
    c['v']  = wc['C9_'+qiqj+ll]      + delta_C9 + Yq2
    c['vp'] = wc['C9p_'+qiqj+ll]
    c['a']  = wc['C10_'+qiqj+ll]
    c['ap'] = wc['C10p_'+qiqj+ll]
    c['s']  = mb * wc['CS_'+qiqj+ll]
    c['sp'] = mb * wc['CSp_'+qiqj+ll]
    c['p']  = mb * wc['CP_'+qiqj+ll]
    c['pp'] = mb * wc['CPp_'+qiqj+ll]
    c['t']  = 0
    c['tp'] = 0
    return c
Beispiel #2
0
def get_wceff(q2, wc, par, B, M, lep, scale):
    r"""Get a dictionary with the effective $\Delta F=1$ Wilson coefficients
    in the convention appropriate for the generalized angular distributions.
    """
    xi_u = ckm.xi('u', meson_quark[(B, M)])(par)
    xi_t = ckm.xi('t', meson_quark[(B, M)])(par)
    qiqj = meson_quark[(B, M)]
    Yq2 = matrixelements.Y(
        q2, wc, par, scale,
        qiqj) + (xi_u / xi_t) * matrixelements.Yu(q2, wc, par, scale, qiqj)
    #   b) NNLO Q1,2
    delta_C7 = matrixelements.delta_C7(par=par,
                                       wc=wc,
                                       q2=q2,
                                       scale=scale,
                                       qiqj=qiqj)
    delta_C9 = matrixelements.delta_C9(par=par,
                                       wc=wc,
                                       q2=q2,
                                       scale=scale,
                                       qiqj=qiqj)
    mb = running.get_mb(par, scale)
    ll = lep + lep
    c = {}
    c['7'] = wc['C7eff_' + qiqj] + delta_C7
    c['7p'] = wc['C7effp_' + qiqj]
    c['v'] = wc['C9_' + qiqj + ll] + delta_C9 + Yq2
    c['vp'] = wc['C9p_' + qiqj + ll]
    c['a'] = wc['C10_' + qiqj + ll]
    c['ap'] = wc['C10p_' + qiqj + ll]
    c['s'] = mb * wc['CS_' + qiqj + ll]
    c['sp'] = mb * wc['CSp_' + qiqj + ll]
    c['p'] = mb * wc['CP_' + qiqj + ll]
    c['pp'] = mb * wc['CPp_' + qiqj + ll]
    c['t'] = 0
    c['tp'] = 0
    return c