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)