Name: L.X. Xu Date: 11/03/14-06:06:55 AM Z
hi,
the result of (quark,quarkbar->quark,quarkbar) obtained by FeynCalc is
not correct, one term of the result differ by a minus sign from the
correct. We can just find the correct at page 196 of Langacker’s
book<the standard model and beyond> or page 571 of Peskin EQ(17,70).
is there an problem of my code or a bug?
I’m appreciate the help!
Here is my code:
Quit;
«HighEnergyPhysics`FeynCalc`
ClearScalarProducts;
{ScalarProduct[p1, p1] =
ScalarProduct[p2, p2] =
ScalarProduct[p3, p3] = ScalarProduct[p4, p4] = 0,
ScalarProduct[p1, p2] = ScalarProduct[p3, p4] = s/2,
ScalarProduct[p1, p3] = ScalarProduct[p2, p4] =
-t/2,
ScalarProduct[p1, p4] = ScalarProduct[p2, p3] = -u/2
};
ScPr[p_, m_] := -I/(ScalarProduct[p] - m^2) //
ExpandScalarProduct;
ftrace = {DiracTrace -> Tr2, D -> 4};
SUNN = 3;
SetOptions[SUNSimplify, SUNNToCACF -> False];
qav = 6;
f1 = SpinorVBar[p2, 0].QGV[\[Alpha],
k].SpinorU[p1, 0] ScPr[p1 + p2,
0] SpinorUBar[p3, 0].QGV[\[Alpha],
k].SpinorV[p4, 0] // Explicit;
f2 = SpinorUBar[p3, 0].QGV[\[Alpha],
k].SpinorU[p1, 0] ScPr[p1 - p3,
0] SpinorVBar[p2, 0].QGV[\[Alpha],
k].SpinorV[p4, 0] // Explicit;
f = f1 + f2
f1s = SpinorUBar[p1, 0].QGV[\[Beta],
l].SpinorV[p2,
0] (-ScPr[p1 + p2, 0]) SpinorVBar[p4,
0].QGV[\[Beta],
l].SpinorU[p3, 0] // Explicit;
f2s = SpinorUBar[p1, 0].QGV[\[Beta],
l].SpinorU[p3,
0] (-ScPr[p1 - p3, 0]) SpinorVBar[p4,
0].QGV[\[Beta],
l].SpinorV[p2, 0] // Explicit;
fstar = f1s + f2s
Msq = FermionSpinSum[
f fstar // Explicit // Expand]/(qav^2 Gstrong^4) /. ftrace
//
Contract // Simplify // SUNSimplify // Expand
standard = 4/9 ((s^2 + u^2)/t^2 + (u^2 + t^2)/s^2 - 2/3 u^2/(s t))
TrickMandelstam[Msq - standard, {s, t, u, 0}]