hi，recently I've calculated some QCD processes using FeynCalc.Here comes my questions: firstly, the result of the process "quark,quarkbar to quark,quarkbar" obtained by FeynCalc is wrong.Secondly,when calculating the process of quark gluon scattering, the commmand "Contract" seems not working right. There is still Lorentz indexes after using it.
Here are my Codes for two processes:
Q,Qbar to Q,Qbar:
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}]
Quark Gluon Scattering:
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
};
QPr[p_, m_] :=
I (DiracSlash[p] + m)/(ScalarProduct[p] - m^2) //
ExpandScalarProduct;
SUNN = 3;
SetOptions[SUNSimplify, SUNNToCACF -> False];
{qav = 6, gav = 16, eav = 2};
ftrace := {DiracTrace -> Tr2, D -> 4};
f1munu = SpinorUBar[p3, 0].QGV[\[Nu], j].QPr[p1 + p2, 0].QGV[\[Mu],
i].SpinorU[p1, 0] // Explicit
f2munu = SpinorUBar[p3, 0].QGV[\[Mu], i].QPr[p1 - p4, 0].QGV[v,
j].SpinorU[p1, 0] // Explicit
fmunu = f1munu + f2munu;
f1s = -SpinorUBar[p1, 0].QGV[\[Rho], i].QPr[p1 + p2, 0].QGV[\[Sigma],
j].SpinorU[p3, 0] // Explicit
f2s = -SpinorUBar[p1, 0].QGV[\[Sigma], j].QPr[p1 - p4, 0].QGV[\[Rho],
i].SpinorU[p3, 0] // Explicit
fmunus = f1s + f2s;
Pols = PolarizationSum[\[Mu], \[Rho], p2,
p1] PolarizationSum[\[Nu], \[Sigma], p4, p3]
Msqmunu = FermionSpinSum[fmunu fmunus /(gav qav) // Expand]
SUNSimplify[SUNTrace[%], Explicit -> True]
% /. DiracTrace -> TR;
Msq = Pols % // Contract // Simplify
I'm appreciate for the help,cheers!
This archive was generated by hypermail 2b29 : 09/04/20-12:55:05 AM Z CEST