def bpll_obs(function, q2, wc_obj, par, B, P, l1, l2): ml1 = par['m_'+l1] ml2 = par['m_'+l2] mB = par['m_'+B] mP = par['m_'+P] if q2 <= (ml1+ml2)**2 or q2 > (mB-mP)**2: return 0 scale = config['renormalization scale']['bpll'] mb = running.get_mb(par, scale) h = helicity_amps(q2, wc_obj, par, B, P, l1, l2) h_bar = helicity_amps_bar(q2, wc_obj, par, B, P, l1, l2) J = angular.angularcoeffs_general_p(h, q2, mB, mP, mb, 0, ml1, ml2) J_bar = angular.angularcoeffs_general_p(h_bar, q2, mB, mP, mb, 0, ml1, ml2) return function(J, J_bar)
def _get_angularcoeff(q2, wc_obj, par, B, P, lep, nu): ml = par['m_' + lep] mB = par['m_' + B] mP = par['m_' + P] scale = config['renormalization scale']['bpll'] mb = running.get_mb(par, scale) wc = get_wceff_fccc2(wc_obj, par, meson_quark[(B, P)], lep, nu, mb, scale, nf=5) N = prefactor(q2, par, B, P, lep) ff = get_ff(q2, par, B, P) qi_qj = meson_quark[(B, P)] 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 h = angular.helicity_amps_p(q2, mB, mP, mb, mlight, ml, 0, ff, wc, N) J = angular.angularcoeffs_general_p(h, q2, mB, mP, mb, mlight, ml, 0) return J
def bpll_obs(function, q2, wc_obj, par, B, P, lep): ml = par['m_'+lep] mB = par['m_'+B] mP = par['m_'+P] if q2 <= (ml+ml)**2 or q2 > (mB-mP)**2: return 0 scale = config['renormalization scale']['bpll'] mb = running.get_mb(par, scale) h = helicity_amps(q2, wc_obj, par, B, P, lep) J = angular.angularcoeffs_general_p(h, q2, mB, mP, mb, 0, ml, ml) if lep == lep: h_bar = helicity_amps_bar(q2, wc_obj, par, B, P, lep) J_bar = angular.angularcoeffs_general_p(h_bar, q2, mB, mP, mb, 0, ml, ml) else: # for LFV decays, don't bother about the CP average. There is no strong phase. J_bar = J return function(J, J_bar)
def bpll_obs(function, q2, wc_obj, par, B, P, lep): ml = par['m_' + lep] mB = par['m_' + B] mP = par['m_' + P] if q2 <= (ml + ml)**2 or q2 > (mB - mP)**2: return 0 scale = config['renormalization scale']['bpll'] mb = running.get_mb(par, scale) h = helicity_amps(q2, wc_obj, par, B, P, lep) J = angular.angularcoeffs_general_p(h, q2, mB, mP, mb, 0, ml, ml) if lep == lep: h_bar = helicity_amps_bar(q2, wc_obj, par, B, P, lep) J_bar = angular.angularcoeffs_general_p(h_bar, q2, mB, mP, mb, 0, ml, ml) else: # for LFV decays, don't bother about the CP average. There is no strong phase. J_bar = J return function(J, J_bar)
def get_angularcoeff(q2, wc_obj, par, B, P, lep): ml = par['m_'+lep] mB = par['m_'+B] mP = par['m_'+P] scale = config['renormalization scale']['bpll'] mb = running.get_mb(par, scale) wc = get_wceff_fccc(wc_obj, par, meson_quark[(B,P)], lep, mb, scale, nf=5) N = prefactor(q2, par, B, P, lep) ff = get_ff(q2, par, B, P) qi_qj = meson_quark[(B, P)] 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 h = angular.helicity_amps_p(q2, mB, mP, mb, mlight, ml, 0, ff, wc, N) J = angular.angularcoeffs_general_p(h, q2, mB, mP, mb, mlight, ml, 0) return J
def _get_angularcoeff(q2, wc_obj, par, B, P, lep, nu): scale = config['renormalization scale']['bpll'] mb = running.get_mb(par, scale) wc = get_wceff_fccc(wc_obj, par, meson_quark[(B,P)], lep, nu, mb, scale, nf=5) if lep != nu and all(C == 0 for C in wc.values()): return {'a': 0, 'b': 0, 'c': 0} # if all WCs vanish, so does the AC! ml = par['m_'+lep] mB = par['m_'+B] mP = par['m_'+P] N = prefactor(q2, par, B, P, lep) ff = get_ff(q2, par, B, P) qi_qj = meson_quark[(B, P)] 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 h = angular.helicity_amps_p(q2, mB, mP, mb, mlight, ml, 0, ff, wc, N) J = angular.angularcoeffs_general_p(h, q2, mB, mP, mb, mlight, ml, 0) return J