def _br_klnu(wc_obj, par, P, lep, nu):
# CKM element
if P == 'K+':
Vij = flavio.physics.ckm.get_ckm(par)[0, 1]
qiqj = 'su'
elif P == 'pi+':
Vij = flavio.physics.ckm.get_ckm(par)[0, 0]
qiqj = 'du'
# renormalization scale is m_rho
scale = par['m_rho0']
# Wilson coefficients
wc = wc_obj.get_wc(qiqj + lep + 'nu' + nu, scale, par, nf_out=3)
# add SM contribution to Wilson coefficient
if lep == nu:
wc['CVL_' + qiqj + lep + 'nu' +
nu] += flavio.physics.bdecays.wilsoncoefficients.get_CVLSM(par,
scale,
nf=3)
ms = flavio.physics.running.running.get_ms(par, scale)
mu = flavio.physics.running.running.get_mu(par, scale)
return br_plnu_general(wc,
par,
Vij,
P,
qiqj,
lep,
nu,
ms,
mu,
delta=delta_Plnu(par, P, lep))
def _br_plnu(wc_obj, par, P, lep, nu):
# CKM element
if P == 'K+':
Vij = flavio.physics.ckm.get_ckm(par)[0, 1]
qiqj = 'su'
elif P == 'pi+':
Vij = flavio.physics.ckm.get_ckm(par)[0, 0]
qiqj = 'du'
# renormalization scale is m_rho
scale = par['m_rho0']
# Wilson coefficients
wc = wc_obj.get_wc(qiqj + lep + 'nu' + nu, scale, par, nf_out=3)
# add SM contribution to Wilson coefficient
if lep == nu:
# for the SM contribution, need the Fermi constant with possible
# NP effects in mu->enunu subtracted, not the measured one
r_GF = GFeff(wc_obj, par) / par['GF']
wc['CVL_' + qiqj + lep + 'nu' +
nu] += flavio.physics.bdecays.wilsoncoefficients.get_CVLSM(
par, scale, nf=3) * r_GF
if P == 'K+':
mq = flavio.physics.running.running.get_ms(par, scale)
elif P == 'pi+':
mq = flavio.physics.running.running.get_md(par, scale)
mu = flavio.physics.running.running.get_mu(par, scale)
return br_plnu_general(wc,
par,
Vij,
P,
qiqj,
lep,
nu,
mq,
mu,
delta=delta_Plnu(par, P, lep))
def br_plnu(wc_obj, par, P, lep):
r"""Branching ratio of $P^+\to\ell^+\nu_\ell$."""
# CKM element
if P == 'K+':
Vij = flavio.physics.ckm.get_ckm(par)[0, 1]
qiqj = 'su'
elif P == 'pi+':
Vij = flavio.physics.ckm.get_ckm(par)[0, 0]
qiqj = 'du'
# renormalization scale is m_rho
scale = par['m_rho0']
# Wilson coefficients
wc = wc_obj.get_wc(qiqj + lep + 'nu', scale, par)
# add SM contribution to Wilson coefficient
wc['CV_' + qiqj + lep +
'nu'] += flavio.physics.bdecays.wilsoncoefficients.get_CVSM(par,
scale,
nf=3)
return br_plnu_general(wc,
par,
Vij,
P,
qiqj,
lep,
delta=delta_Plnu(par, P, lep))
def br_Dlnu(wc_obj, par, P, lep):
r"""Branching ratio of $D^+\to\ell^+\nu_\ell$."""
# CKM element
if P=='D+':
Vij = flavio.physics.ckm.get_ckm(par)[1,0].conj() # Vcd*
qiqj = 'dc'
elif P=='Ds':
Vij = flavio.physics.ckm.get_ckm(par)[1,1].conj() # Vcs*
qiqj = 'sc'
scale = flavio.config['renormalization scale']['dll']
# Wilson coefficients
wc = wc_obj.get_wc(qiqj + lep + 'nu', scale, par, nf_out=4)
# add SM contribution to Wilson coefficient
wc['CV_'+qiqj+lep+'nu'] += flavio.physics.bdecays.wilsoncoefficients.get_CVSM(par, scale, nf=4)
return br_plnu_general(wc, par, Vij, P, qiqj, lep, delta=0)
def br_plnu(wc_obj, par, P, lep):
r"""Branching ratio of $P^+\to\ell^+\nu_\ell$."""
# CKM element
if P=='K+':
Vij = flavio.physics.ckm.get_ckm(par)[0,1]
qiqj = 'su'
elif P=='pi+':
Vij = flavio.physics.ckm.get_ckm(par)[0,0]
qiqj = 'du'
# renormalization scale is m_rho
scale = par['m_rho0']
# Wilson coefficients
wc = wc_obj.get_wc(qiqj + lep + 'nu', scale, par)
# add SM contribution to Wilson coefficient
wc['CV_'+qiqj+lep+'nu'] += flavio.physics.bdecays.wilsoncoefficients.get_CVSM(par, scale, nf=3)
return br_plnu_general(wc, par, Vij, P, qiqj, lep, delta=delta_Plnu(par, P, lep))
def _br_Dlnu(wc_obj, par, P, lep, nu):
# CKM element
if P=='D+':
Vij = flavio.physics.ckm.get_ckm(par)[1,0].conj() # Vcd*
qiqj = 'dc'
elif P=='Ds':
Vij = flavio.physics.ckm.get_ckm(par)[1,1].conj() # Vcs*
qiqj = 'sc'
scale = flavio.config['renormalization scale']['dll']
# Wilson coefficients
wc = wc_obj.get_wc(qiqj + lep + 'nu' + nu, scale, par, nf_out=4)
# add SM contribution to Wilson coefficient
if lep == nu:
wc['CVL_'+qiqj+lep+'nu' + nu] += flavio.physics.bdecays.wilsoncoefficients.get_CVLSM(par, scale, nf=4)
mb = flavio.physics.running.running.get_mb(par, scale)
mc = flavio.physics.running.running.get_mc(par, scale)
return br_plnu_general(wc, par, Vij, P, qiqj, lep, nu, mb, mc, delta=0)
def _br_Dlnu(wc_obj, par, P, lep):
r"""Branching ratio of $D^+\to\ell^+\nu_\ell$."""
# CKM element
if P == 'D+':
Vij = flavio.physics.ckm.get_ckm(par)[1, 0].conj() # Vcd*
qiqj = 'dc'
elif P == 'Ds':
Vij = flavio.physics.ckm.get_ckm(par)[1, 1].conj() # Vcs*
qiqj = 'sc'
scale = flavio.config['renormalization scale']['dll']
# Wilson coefficients
wc = wc_obj.get_wc(qiqj + lep + 'nu', scale, par, nf_out=4)
# add SM contribution to Wilson coefficient
wc['CV_' + qiqj + lep +
'nu'] += flavio.physics.bdecays.wilsoncoefficients.get_CVSM(par,
scale,
nf=4)
return br_plnu_general(wc, par, Vij, P, qiqj, lep, delta=0)
def br_Dlnu2(wc_obj, par, P, lep, nu):
if lep == nu:
return _br_Dlnu(wc_obj, par, P, lep)
if P == 'D+':
Vij = flavio.physics.ckm.get_ckm(par)[1, 0].conj() # Vcd*
qiqj = 'dc'
elif P == 'Ds':
Vij = flavio.physics.ckm.get_ckm(par)[1, 1].conj() # Vcs*
qiqj = 'sc'
scale = flavio.config['renormalization scale']['dll']
# Wilson coefficients
wc = wc_obj.get_wc(qiqj + lep + 'nu' + nu, scale, par, nf_out=4)
qqlnu = qiqj + lep + 'nu'
wc['CV_' + qqlnu] = wc['CV_' + qqlnu + nu]
wc['CVp_' + qqlnu] = wc['CVp_' + qqlnu + nu]
wc['CS_' + qqlnu] = wc['CS_' + qqlnu + nu]
wc['CSp_' + qqlnu] = wc['CSp_' + qqlnu + nu]
return br_plnu_general(wc, par, Vij, P, qiqj, lep, delta=0)
def _br_klnu(wc_obj, par, P, lep, nu):
# CKM element
if P=='K+':
Vij = flavio.physics.ckm.get_ckm(par)[0,1]
qiqj = 'su'
elif P=='pi+':
Vij = flavio.physics.ckm.get_ckm(par)[0,0]
qiqj = 'du'
# renormalization scale is m_rho
scale = par['m_rho0']
# Wilson coefficients
wc = wc_obj.get_wc(qiqj + lep + 'nu' + nu, scale, par, nf_out=3)
# add SM contribution to Wilson coefficient
if lep == nu:
wc['CVL_'+qiqj+lep+'nu'+nu] += flavio.physics.bdecays.wilsoncoefficients.get_CVLSM(par, scale, nf=3)
ms = flavio.physics.running.running.get_ms(par, scale)
mu = flavio.physics.running.running.get_mu(par, scale)
return br_plnu_general(wc, par, Vij, P, qiqj, lep, nu, ms, mu, delta=delta_Plnu(par, P, lep))
def _br_Dlnu(wc_obj, par, P, lep, nu):
# CKM element
if P == 'D+':
Vij = flavio.physics.ckm.get_ckm(par)[1, 0].conj() # Vcd*
qiqj = 'dc'
elif P == 'Ds':
Vij = flavio.physics.ckm.get_ckm(par)[1, 1].conj() # Vcs*
qiqj = 'sc'
scale = flavio.config['renormalization scale']['dll']
# Wilson coefficients
wc = wc_obj.get_wc(qiqj + lep + 'nu' + nu, scale, par, nf_out=4)
# add SM contribution to Wilson coefficient
if lep == nu:
wc['CVL_' + qiqj + lep + 'nu' +
nu] += flavio.physics.bdecays.wilsoncoefficients.get_CVLSM(par,
scale,
nf=4)
mb = flavio.physics.running.running.get_mb(par, scale)
mc = flavio.physics.running.running.get_mc(par, scale)
return br_plnu_general(wc, par, Vij, P, qiqj, lep, nu, mb, mc, delta=0)