def S(wc_obj, par, meson, amplitude, etaCP): M12, G12 = get_M12_G12(wc_obj, par, meson) qp = common.q_over_p(M12, G12) A = amplitude(par) A_bar = amplitude(conjugate_par(par)) xi = etaCP * qp * A / A_bar return -2*xi.imag / ( 1 + abs(xi)**2 )
def S(wc_obj, par, meson, amplitude, etaCP): M12, G12 = get_M12_G12(wc_obj, par, meson) qp = common.q_over_p(M12, G12) A = amplitude(par) A_bar = amplitude(conjugate_par(par)) xi = etaCP * qp * A / A_bar return -2 * xi.imag / (1 + abs(xi)**2)
def S(wc_obj, par, meson, amplitude, etaCP): M12, G12 = get_M12_G12(wc_obj, par, meson) qp = common.q_over_p(M12, G12) DM = common.DeltaM(M12, G12) if DM < 0: qp = -qp # switch the sign of q/p to keep DeltaM > 0 A = amplitude(par) A_bar = amplitude(conjugate_par(par)) xi = etaCP * qp * A / A_bar return -2*xi.imag / ( 1 + abs(xi)**2 )