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 )
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 )