def _get_angularcoeff(q2, wc_obj, par, B, V, lep, nu): scale = config['renormalization scale']['bvll'] mb = running.get_mb(par, scale) wc = get_wceff_fccc2(wc_obj, par, meson_quark[(B, V)], lep, nu, mb, scale, nf=5) ml = par['m_' + lep] mB = par['m_' + B] mV = par['m_' + V] qi_qj = meson_quark[(B, V)] if qi_qj == 'bu': mlight = 0. # neglecting the up quark mass if qi_qj == 'bc': mlight = running.get_mc( par, scale) # this is needed for scalar contributions N = prefactor(q2, par, B, V, lep) ff = get_ff(q2, par, B, V) h = angular.helicity_amps_v(q2, mB, mV, mb, mlight, ml, 0, ff, wc, N) J = angular.angularcoeffs_general_v(h, q2, mB, mV, mb, mlight, ml, 0) return J
def _get_angularcoeff(q2, wc_obj, par, B, V, lep, nu): scale = config['renormalization scale']['bvll'] mb = running.get_mb(par, scale) wc = get_wceff_fccc(wc_obj, par, meson_quark[(B, V)], lep, nu, mb, scale, nf=5) if lep != nu and all(C == 0 for C in wc.values()): # if all WCs vanish, so does the AC! return { k: 0 for k in ['1s', '1c', '2s', '2c', '6s', '6c', 3, 4, 5, 7, 8, 9] } ml = par['m_' + lep] mB = par['m_' + B] mV = par['m_' + V] qi_qj = meson_quark[(B, V)] if qi_qj == 'bu': mlight = 0. # neglecting the up quark mass if qi_qj == 'bc': mlight = running.get_mc( par, scale) # this is needed for scalar contributions N = prefactor(q2, par, B, V, lep) ff = get_ff(q2, par, B, V) h = angular.helicity_amps_v(q2, mB, mV, mb, mlight, ml, 0, ff, wc, N) J = angular.angularcoeffs_general_v(h, q2, mB, mV, mb, mlight, ml, 0) return J
def get_angularcoeff(q2, wc_obj, par, B, V, lep): scale = config['renormalization scale']['bvll'] mb = running.get_mb(par, scale) wc = get_wceff_fccc(wc_obj, par, meson_quark[(B,V)], lep, mb, scale, nf=5) ml = par['m_'+lep] mB = par['m_'+B] mV = par['m_'+V] qi_qj = meson_quark[(B, V)] if qi_qj == 'bu': mlight = 0. # neglecting the up quark mass if qi_qj == 'bc': mlight = running.get_mc(par, scale) # this is needed for scalar contributions N = prefactor(q2, par, B, V, lep) ff = get_ff(q2, par, B, V) h = angular.helicity_amps_v(q2, mB, mV, mb, mlight, ml, 0, ff, wc, N) J = angular.angularcoeffs_general_v(h, q2, mB, mV, mb, mlight, ml, 0) return J
def helicity_amps_ff(q2, ff, wc_obj, par_dict, B, V, lep, cp_conjugate): par = par_dict.copy() if cp_conjugate: par = conjugate_par(par) scale = config['renormalization scale']['bvll'] label = meson_quark[(B,V)] + lep + lep # e.g. bsmumu, bdtautau wc = wctot_dict(wc_obj, label, scale, par) if cp_conjugate: wc = conjugate_wc(wc) wc_eff = get_wceff(q2, wc, par, B, V, lep, scale) ml = par['m_'+lep] mB = par['m_'+B] mV = par['m_'+V] mb = running.get_mb(par, scale) N = prefactor(q2, par, B, V) h = angular.helicity_amps_v(q2, mB, mV, mb, 0, ml, ml, ff, wc_eff, N) return h
def helicity_amps_ff(q2, ff, wc_obj, par_dict, B, V, lep, cp_conjugate): par = par_dict.copy() if cp_conjugate: par = conjugate_par(par) scale = config['renormalization scale']['bvll'] label = meson_quark[(B, V)] + lep + lep # e.g. bsmumu, bdtautau wc = wctot_dict(wc_obj, label, scale, par) if cp_conjugate: wc = conjugate_wc(wc) wc_eff = get_wceff(q2, wc, par, B, V, lep, scale) ml = par['m_' + lep] mB = par['m_' + B] mV = par['m_' + V] mb = running.get_mb(par, scale) N = prefactor(q2, par, B, V) h = angular.helicity_amps_v(q2, mB, mV, mb, 0, ml, ml, ff, wc_eff, N) return h
def _get_angularcoeff(q2, wc_obj, par, B, V, lep, nu): scale = config['renormalization scale']['bvll'] mb = running.get_mb(par, scale) wc = get_wceff_fccc(wc_obj, par, meson_quark[(B,V)], lep, nu, mb, scale, nf=5) if lep != nu and all(C == 0 for C in wc.values()): # if all WCs vanish, so does the AC! return {k: 0 for k in ['1s', '1c', '2s', '2c', '6s', '6c', 3, 4, 5, 7, 8, 9]} ml = par['m_'+lep] mB = par['m_'+B] mV = par['m_'+V] qi_qj = meson_quark[(B, V)] if qi_qj == 'bu': mlight = 0. # neglecting the up quark mass if qi_qj == 'bc': mlight = running.get_mc(par, scale) # this is needed for scalar contributions N = prefactor(q2, par, B, V, lep) ff = get_ff(q2, par, B, V) h = angular.helicity_amps_v(q2, mB, mV, mb, mlight, ml, 0, ff, wc, N) J = angular.angularcoeffs_general_v(h, q2, mB, mV, mb, mlight, ml, 0) return J