def T_para_minus_QSS(q2, par, wc, B, V, u, scale):
mB, mb, mc, alpha_s, q, eq, ed, eu, eps_u, qiqj = get_input(par, B, V, scale)
ubar = 1 - u
h_mc = matrixelements.h(ubar*mB**2 + u*q2, mc, scale)
h_mb = matrixelements.h(ubar*mB**2 + u*q2, mb, scale)
h_0 = matrixelements.h(ubar*mB**2 + u*q2, 0, scale)
T_t = (alpha_s/(3*pi)) * eq * 6 * mB/mb*(
h_mc * (-wc['C1_'+qiqj]/6. + wc['C2_'+qiqj] + wc['C4_'+qiqj] + 10*wc['C6_'+qiqj])
+ h_mb * (wc['C3_'+qiqj] + 5*wc['C4_'+qiqj]/6. + 16*wc['C5_'+qiqj] + 22*wc['C6_'+qiqj]/3.)
+ h_0 * (wc['C3_'+qiqj] + 17*wc['C4_'+qiqj]/6. + 16*wc['C5_'+qiqj] + 82*wc['C6_'+qiqj]/3.)
- 8/27. * (-15*wc['C4_'+qiqj]/2. + 12*wc['C5_'+qiqj] - 32*wc['C6_'+qiqj])
)
T_u = ( (alpha_s/(3*pi)) * eq * 6*mB/mb * ( h_mc - h_0 )
* ( wc['C2_'+qiqj] - wc['C1_'+qiqj]/6.) )
return T_t + eps_u * T_u
def T_para_minus_QSS(q2, par, wc, B, V, u, scale):
mB, mb, mc, alpha_s, q, eq, ed, eu, eps_u, qiqj = get_input(par, B, V, scale)
ubar = 1 - u
h_mc = matrixelements.h(ubar*mB**2 + u*q2, mc, scale)
h_mb = matrixelements.h(ubar*mB**2 + u*q2, mb, scale)
h_0 = matrixelements.h(ubar*mB**2 + u*q2, 0, scale)
T_t = (alpha_s/(3*pi)) * eq * 6 * mB/mb*(
h_mc * (-wc['C1_'+qiqj]/6. + wc['C2_'+qiqj] + wc['C4_'+qiqj] + 10*wc['C6_'+qiqj])
+ h_mb * (wc['C3_'+qiqj] + 5*wc['C4_'+qiqj]/6. + 16*wc['C5_'+qiqj] + 22*wc['C6_'+qiqj]/3.)
+ h_0 * (wc['C3_'+qiqj] + 17*wc['C4_'+qiqj]/6. + 16*wc['C5_'+qiqj] + 82*wc['C6_'+qiqj]/3.)
- 8/27. * (-15*wc['C4_'+qiqj]/2. + 12*wc['C5_'+qiqj] - 32*wc['C6_'+qiqj])
)
T_u = ( (alpha_s/(3*pi)) * eq * 6*mB/mb * ( h_mc - h_0 )
* ( wc['C2_'+qiqj] - wc['C1_'+qiqj]/6.) )
return T_t + eps_u * T_u