StandardModel: SMQCD.nlo

(****************************************************************************************************************)
(* Model automalically generated by NLOCT-0.3                                                                   *)
(* date and time of generation : Tue 4 Mar 2014 14:20:53                                                        *)
(* FeynRules model information :                                                                                *)
(*   ModelName : Standard Model                                                                                 *)
(*   Authors : N. Christensen, C. Duhr, B. Fuks                                                                 *)
(*   Version : 1.4.5                                                                                            *)
(*   Date : 21. 11. 2012                                                                                        *)
(*   Institutions : Michigan State University, Universite catholique de Louvain (CP3), IPHC Strasbourg / Univer *)
(*   Emails : neil@pa.msu.edu, claude.duhr@uclouvain.be, benjamin.fuks@cnrs.in2p3.fr                            *)
(*   URLs : http://feynrules.phys.ucl.ac.be/view/Main/StandardModel                                             *)
(****************************************************************************************************************)

(* FeynArts model files used to generated the output *)

CT$Model = "SMQCDrenoL/SMQCDrenoL";
CT$GenericModel = "Lorentz";

(* Those assumptions where made : *)

NLOCT$assumptions = {MH < 2*MW, MH < 2*MZ, MW < MZ, MZ < 2*MW, MH > 0, MZ > 0, MW > 0};

(* Perturbed orders *)

FR$InteractionOrderPerturbativeExpansion = {{QCD, 1}, {QED, 0}};

(*R2 vertices*)

R2$vertlist = {
{{{anti[u], 1}, {u, 2}}, ((-I/12)*gs^2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{u, G}]*(TensDot[SlashedP[2], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]] + TensDot[SlashedP[2], ProjP][Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/Pi^2},
{{{anti[c], 1}, {c, 2}}, ((I/12)*gs^2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{c, G}]*(2*MC*ProjM[Index[Spin, Ext[1]], Index[Spin, Ext[2]]] + 2*MC*ProjP[Index[Spin, Ext[1]], Index[Spin, Ext[2]]] - TensDot[SlashedP[2], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]] - TensDot[SlashedP[2], ProjP][Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/Pi^2},
{{{anti[t], 1}, {t, 2}}, ((I/12)*gs^2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{t, G}]*(2*MT*ProjM[Index[Spin, Ext[1]], Index[Spin, Ext[2]]] + 2*MT*ProjP[Index[Spin, Ext[1]], Index[Spin, Ext[2]]] - TensDot[SlashedP[2], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]] - TensDot[SlashedP[2], ProjP][Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/Pi^2},
{{{anti[d], 1}, {d, 2}}, ((-I/12)*gs^2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{d, G}]*(TensDot[SlashedP[2], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]] + TensDot[SlashedP[2], ProjP][Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/Pi^2},
{{{anti[s], 1}, {s, 2}}, ((-I/12)*gs^2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{s, G}]*(TensDot[SlashedP[2], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]] + TensDot[SlashedP[2], ProjP][Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/Pi^2},
{{{anti[b], 1}, {b, 2}}, ((I/12)*gs^2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{b, G}]*(2*MB*ProjM[Index[Spin, Ext[1]], Index[Spin, Ext[2]]] + 2*MB*ProjP[Index[Spin, Ext[1]], Index[Spin, Ext[2]]] - TensDot[SlashedP[2], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]] - TensDot[SlashedP[2], ProjP][Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/Pi^2},
{{{G, 1}, {G, 2}}, ((I/96)*gs^2*IndexDelta[Index[Gluon, Ext[1]], Index[Gluon, Ext[2]]]*(-6*FV[2, Index[Lorentz, Ext[1]]]*FV[2, Index[Lorentz, Ext[2]]]*IPL[{G}] + ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[2]]]*(-12*MB^2*IPL[{b}] + 2*IPL[{b}]*SP[2, 2] + 2*IPL[{d}]*SP[2, 2] + 9*IPL[{G}]*SP[2, 2] + 2*IPL[{s}]*SP[2, 2] + 2*IPL[{u}]*SP[2, 2] + 2*IPL[{c}]*(-6*MC^2 + SP[2, 2]) + 2*IPL[{t}]*(-6*MT^2 + SP[2, 2]))))/Pi^2},
{{{anti[c], 1}, {c, 2}, {H, 3}}, ((I/3)*gs^2*MC*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{c, G}]*(ProjM[Index[Spin, Ext[1]], Index[Spin, Ext[2]]] + ProjP[Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(Pi^2*vev)},
{{{anti[t], 1}, {t, 2}, {H, 3}}, ((I/3)*gs^2*MT*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{t, G}]*(ProjM[Index[Spin, Ext[1]], Index[Spin, Ext[2]]] + ProjP[Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(Pi^2*vev)},
{{{anti[b], 1}, {b, 2}, {H, 3}}, ((I/3)*gs^2*MB*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{b, G}]*(ProjM[Index[Spin, Ext[1]], Index[Spin, Ext[2]]] + ProjP[Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(Pi^2*vev)},
{{{anti[c], 1}, {c, 2}, {G0, 3}}, -(gs^2*MC*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{c, G}]*(ProjM[Index[Spin, Ext[1]], Index[Spin, Ext[2]]] - ProjP[Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(3*Pi^2*vev)},
{{{anti[t], 1}, {t, 2}, {G0, 3}}, -(gs^2*MT*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{t, G}]*(ProjM[Index[Spin, Ext[1]], Index[Spin, Ext[2]]] - ProjP[Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(3*Pi^2*vev)},
{{{anti[b], 1}, {b, 2}, {G0, 3}}, (gs^2*MB*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{b, G}]*(ProjM[Index[Spin, Ext[1]], Index[Spin, Ext[2]]] - ProjP[Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(3*Pi^2*vev)},
{{{anti[s], 1}, {c, 2}, {anti[GP], 3}}, (Sqrt[2]*gs^2*MC*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{c, s, G}]*ProjP[Index[Spin, Ext[1]], Index[Spin, Ext[2]]])/(3*Pi^2*vev)},
{{{anti[b], 1}, {t, 2}, {anti[GP], 3}}, -(Sqrt[2]*gs^2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{t, b, G}]*(MB*ProjM[Index[Spin, Ext[1]], Index[Spin, Ext[2]]] - MT*ProjP[Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(3*Pi^2*vev)},
{{{anti[c], 1}, {s, 2}, {GP, 3}}, -(Sqrt[2]*gs^2*MC*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{c, s, G}]*ProjM[Index[Spin, Ext[1]], Index[Spin, Ext[2]]])/(3*Pi^2*vev)},
{{{anti[t], 1}, {b, 2}, {GP, 3}}, -(Sqrt[2]*gs^2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{t, b, G}]*(MT*ProjM[Index[Spin, Ext[1]], Index[Spin, Ext[2]]] - MB*ProjP[Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(3*Pi^2*vev)},
{{{anti[u], 1}, {u, 2}, {A, 3}}, ((-I/9)*ee*gs^2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{u, G}]*(TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]] + TensDot[Ga[Index[Lorentz, Ext[3]]], ProjP][Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/Pi^2},
{{{anti[c], 1}, {c, 2}, {A, 3}}, ((-I/9)*ee*gs^2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{c, G}]*(TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]] + TensDot[Ga[Index[Lorentz, Ext[3]]], ProjP][Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/Pi^2},
{{{anti[t], 1}, {t, 2}, {A, 3}}, ((-I/9)*ee*gs^2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{t, G}]*(TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]] + TensDot[Ga[Index[Lorentz, Ext[3]]], ProjP][Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/Pi^2},
{{{anti[d], 1}, {d, 2}, {A, 3}}, ((I/18)*ee*gs^2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{d, G}]*(TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]] + TensDot[Ga[Index[Lorentz, Ext[3]]], ProjP][Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/Pi^2},
{{{anti[s], 1}, {s, 2}, {A, 3}}, ((I/18)*ee*gs^2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{s, G}]*(TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]] + TensDot[Ga[Index[Lorentz, Ext[3]]], ProjP][Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/Pi^2},
{{{anti[b], 1}, {b, 2}, {A, 3}}, ((I/18)*ee*gs^2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{b, G}]*(TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]] + TensDot[Ga[Index[Lorentz, Ext[3]]], ProjP][Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/Pi^2},
{{{anti[u], 1}, {u, 2}, {Z, 3}}, ((-I/36)*ee*gs^2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{u, G}]*((3*cw^2 - sw^2)*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]] - 4*sw^2*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjP][Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(cw*Pi^2*sw)},
{{{anti[c], 1}, {c, 2}, {Z, 3}}, ((-I/36)*ee*gs^2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{c, G}]*((3*cw^2 - sw^2)*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]] - 4*sw^2*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjP][Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(cw*Pi^2*sw)},
{{{anti[t], 1}, {t, 2}, {Z, 3}}, ((-I/36)*ee*gs^2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{t, G}]*((3*cw^2 - sw^2)*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]] - 4*sw^2*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjP][Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(cw*Pi^2*sw)},
{{{anti[d], 1}, {d, 2}, {Z, 3}}, ((I/36)*ee*gs^2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{d, G}]*((3*cw^2 + sw^2)*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]] - 2*sw^2*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjP][Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(cw*Pi^2*sw)},
{{{anti[s], 1}, {s, 2}, {Z, 3}}, ((I/36)*ee*gs^2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{s, G}]*((3*cw^2 + sw^2)*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]] - 2*sw^2*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjP][Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(cw*Pi^2*sw)},
{{{anti[b], 1}, {b, 2}, {Z, 3}}, ((I/36)*ee*gs^2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{b, G}]*((3*cw^2 + sw^2)*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]] - 2*sw^2*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjP][Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(cw*Pi^2*sw)},
{{{anti[d], 1}, {u, 2}, {anti[W], 3}}, ((-I/6)*ee*gs^2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{u, d, G}]*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]])/(Sqrt[2]*Pi^2*sw)},
{{{anti[s], 1}, {c, 2}, {anti[W], 3}}, ((-I/6)*ee*gs^2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{c, s, G}]*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]])/(Sqrt[2]*Pi^2*sw)},
{{{anti[b], 1}, {t, 2}, {anti[W], 3}}, ((-I/6)*ee*gs^2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{t, b, G}]*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]])/(Sqrt[2]*Pi^2*sw)},
{{{anti[u], 1}, {d, 2}, {W, 3}}, ((-I/6)*ee*gs^2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{u, d, G}]*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]])/(Sqrt[2]*Pi^2*sw)},
{{{anti[c], 1}, {s, 2}, {W, 3}}, ((-I/6)*ee*gs^2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{c, s, G}]*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]])/(Sqrt[2]*Pi^2*sw)},
{{{anti[t], 1}, {b, 2}, {W, 3}}, ((-I/6)*ee*gs^2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{t, b, G}]*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]])/(Sqrt[2]*Pi^2*sw)},
{{{anti[u], 1}, {u, 2}, {G, 3}}, ((-I/6)*gs^3*IPL[{u, G}]*SUNT[Index[Gluon, Ext[3]], Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*(TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]] + TensDot[Ga[Index[Lorentz, Ext[3]]], ProjP][Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/Pi^2},
{{{anti[c], 1}, {c, 2}, {G, 3}}, ((-I/6)*gs^3*IPL[{c, G}]*SUNT[Index[Gluon, Ext[3]], Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*(TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]] + TensDot[Ga[Index[Lorentz, Ext[3]]], ProjP][Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/Pi^2},
{{{anti[t], 1}, {t, 2}, {G, 3}}, ((-I/6)*gs^3*IPL[{t, G}]*SUNT[Index[Gluon, Ext[3]], Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*(TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]] + TensDot[Ga[Index[Lorentz, Ext[3]]], ProjP][Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/Pi^2},
{{{anti[d], 1}, {d, 2}, {G, 3}}, ((-I/6)*gs^3*IPL[{d, G}]*SUNT[Index[Gluon, Ext[3]], Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*(TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]] + TensDot[Ga[Index[Lorentz, Ext[3]]], ProjP][Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/Pi^2},
{{{anti[s], 1}, {s, 2}, {G, 3}}, ((-I/6)*gs^3*IPL[{s, G}]*SUNT[Index[Gluon, Ext[3]], Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*(TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]] + TensDot[Ga[Index[Lorentz, Ext[3]]], ProjP][Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/Pi^2},
{{{anti[b], 1}, {b, 2}, {G, 3}}, ((-I/6)*gs^3*IPL[{b, G}]*SUNT[Index[Gluon, Ext[3]], Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*(TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]] + TensDot[Ga[Index[Lorentz, Ext[3]]], ProjP][Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/Pi^2},
{{{G, 1}, {G, 2}, {G, 3}}, -(gs^3*(8*IPL[{b}] + 8*IPL[{c}] + 8*IPL[{d}] + 33*IPL[{G}] + 8*IPL[{s}] + 8*IPL[{t}] + 8*IPL[{u}])*(FV[1, Index[Lorentz, Ext[3]]]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[2]]] - FV[2, Index[Lorentz, Ext[3]]]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[2]]] - FV[1, Index[Lorentz, Ext[2]]]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[3]]] + FV[3, Index[Lorentz, Ext[2]]]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[3]]] + FV[2, Index[Lorentz, Ext[1]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[3]]] - FV[3, Index[Lorentz, Ext[1]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[3]]])*SUNF[Index[Gluon, Ext[1]], Index[Gluon, Ext[2]], Index[Gluon, Ext[3]]])/(192*Pi^2)},
{{{G, 1}, {G, 2}, {H, 3}}, ((-I/8)*gs^2*IndexDelta[Index[Gluon, Ext[1]], Index[Gluon, Ext[2]]]*(MB^2*IPL[{b}] + MC^2*IPL[{c}] + MT^2*IPL[{t}])*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[2]]])/(Pi^2*vev)},
{{{G, 1}, {G, 2}, {G, 3}, {G, 4}}, ((I/96)*gs^4*(-2*IndexDelta[Index[Gluon, Ext[1]], Index[Gluon, Ext[2]]]*IndexDelta[Index[Gluon, Ext[3]], Index[Gluon, Ext[4]]]*IPL[{G}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[4]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[3]]] - 2*IndexDelta[Index[Gluon, Ext[1]], Index[Gluon, Ext[2]]]*IndexDelta[Index[Gluon, Ext[3]], Index[Gluon, Ext[4]]]*IPL[{G}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[3]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[4]]] - 2*IndexDelta[Index[Gluon, Ext[1]], Index[Gluon, Ext[2]]]*IndexDelta[Index[Gluon, Ext[3]], Index[Gluon, Ext[4]]]*IPL[{G}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[2]]]*ME[Index[Lorentz, Ext[3]], Index[Lorentz, Ext[4]]] - 2*IndexDelta[Index[Gluon, Ext[1]], Index[Gluon, Ext[4]]]*IndexDelta[Index[Gluon, Ext[2]], Index[Gluon, Ext[3]]]*IPL[{G}]*(ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[4]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[3]]] + ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[3]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[4]]] + ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[2]]]*ME[Index[Lorentz, Ext[3]], Index[Lorentz, Ext[4]]]) - 2*IndexDelta[Index[Gluon, Ext[1]], Index[Gluon, Ext[3]]]*IndexDelta[Index[Gluon, Ext[2]], Index[Gluon, Ext[4]]]*IPL[{G}]*(ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[4]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[3]]] + ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[3]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[4]]] + ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[2]]]*ME[Index[Lorentz, Ext[3]], Index[Lorentz, Ext[4]]]) - 6*IPL[{b}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[4]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[3]]]*SUNF[Index[Gluon, Ext[1]], Index[Gluon, Ext[2]], Index[Gluon, Ext[3]], Index[Gluon, Ext[4]]] - 6*IPL[{c}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[4]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[3]]]*SUNF[Index[Gluon, Ext[1]], Index[Gluon, Ext[2]], Index[Gluon, Ext[3]], Index[Gluon, Ext[4]]] - 6*IPL[{d}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[4]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[3]]]*SUNF[Index[Gluon, Ext[1]], Index[Gluon, Ext[2]], Index[Gluon, Ext[3]], Index[Gluon, Ext[4]]] - 21*IPL[{G}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[4]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[3]]]*SUNF[Index[Gluon, Ext[1]], Index[Gluon, Ext[2]], Index[Gluon, Ext[3]], Index[Gluon, Ext[4]]] - 6*IPL[{s}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[4]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[3]]]*SUNF[Index[Gluon, Ext[1]], Index[Gluon, Ext[2]], Index[Gluon, Ext[3]], Index[Gluon, Ext[4]]] - 6*IPL[{t}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[4]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[3]]]*SUNF[Index[Gluon, Ext[1]], Index[Gluon, Ext[2]], Index[Gluon, Ext[3]], Index[Gluon, Ext[4]]] - 6*IPL[{u}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[4]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[3]]]*SUNF[Index[Gluon, Ext[1]], Index[Gluon, Ext[2]], Index[Gluon, Ext[3]], Index[Gluon, Ext[4]]] + 10*IPL[{b}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[3]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[4]]]*SUNF[Index[Gluon, Ext[1]], Index[Gluon, Ext[2]], Index[Gluon, Ext[3]], Index[Gluon, Ext[4]]] + 10*IPL[{c}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[3]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[4]]]*SUNF[Index[Gluon, Ext[1]], Index[Gluon, Ext[2]], Index[Gluon, Ext[3]], Index[Gluon, Ext[4]]] + 10*IPL[{d}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[3]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[4]]]*SUNF[Index[Gluon, Ext[1]], Index[Gluon, Ext[2]], Index[Gluon, Ext[3]], Index[Gluon, Ext[4]]] + 57*IPL[{G}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[3]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[4]]]*SUNF[Index[Gluon, Ext[1]], Index[Gluon, Ext[2]], Index[Gluon, Ext[3]], Index[Gluon, Ext[4]]] + 10*IPL[{s}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[3]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[4]]]*SUNF[Index[Gluon, Ext[1]], Index[Gluon, Ext[2]], Index[Gluon, Ext[3]], Index[Gluon, Ext[4]]] + 10*IPL[{t}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[3]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[4]]]*SUNF[Index[Gluon, Ext[1]], Index[Gluon, Ext[2]], Index[Gluon, Ext[3]], Index[Gluon, Ext[4]]] + 10*IPL[{u}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[3]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[4]]]*SUNF[Index[Gluon, Ext[1]], Index[Gluon, Ext[2]], Index[Gluon, Ext[3]], Index[Gluon, Ext[4]]] - 6*IPL[{b}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[2]]]*ME[Index[Lorentz, Ext[3]], Index[Lorentz, Ext[4]]]*SUNF[Index[Gluon, Ext[1]], Index[Gluon, Ext[2]], Index[Gluon, Ext[3]], Index[Gluon, Ext[4]]] - 6*IPL[{c}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[2]]]*ME[Index[Lorentz, Ext[3]], Index[Lorentz, Ext[4]]]*SUNF[Index[Gluon, Ext[1]], Index[Gluon, Ext[2]], Index[Gluon, Ext[3]], Index[Gluon, Ext[4]]] - 6*IPL[{d}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[2]]]*ME[Index[Lorentz, Ext[3]], Index[Lorentz, Ext[4]]]*SUNF[Index[Gluon, Ext[1]], Index[Gluon, Ext[2]], Index[Gluon, Ext[3]], Index[Gluon, Ext[4]]] + 24*IPL[{G}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[2]]]*ME[Index[Lorentz, Ext[3]], Index[Lorentz, Ext[4]]]*SUNF[Index[Gluon, Ext[1]], Index[Gluon, Ext[2]], Index[Gluon, Ext[3]], Index[Gluon, Ext[4]]] - 6*IPL[{s}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[2]]]*ME[Index[Lorentz, Ext[3]], Index[Lorentz, Ext[4]]]*SUNF[Index[Gluon, Ext[1]], Index[Gluon, Ext[2]], Index[Gluon, Ext[3]], Index[Gluon, Ext[4]]] - 6*IPL[{t}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[2]]]*ME[Index[Lorentz, Ext[3]], Index[Lorentz, Ext[4]]]*SUNF[Index[Gluon, Ext[1]], Index[Gluon, Ext[2]], Index[Gluon, Ext[3]], Index[Gluon, Ext[4]]] - 6*IPL[{u}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[2]]]*ME[Index[Lorentz, Ext[3]], Index[Lorentz, Ext[4]]]*SUNF[Index[Gluon, Ext[1]], Index[Gluon, Ext[2]], Index[Gluon, Ext[3]], Index[Gluon, Ext[4]]] - 4*IPL[{b}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[4]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[3]]]*SUNF[Index[Gluon, Ext[1]], Index[Gluon, Ext[3]], Index[Gluon, Ext[2]], Index[Gluon, Ext[4]]] - 4*IPL[{c}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[4]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[3]]]*SUNF[Index[Gluon, Ext[1]], Index[Gluon, Ext[3]], Index[Gluon, Ext[2]], Index[Gluon, Ext[4]]] - 4*IPL[{d}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[4]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[3]]]*SUNF[Index[Gluon, Ext[1]], Index[Gluon, Ext[3]], Index[Gluon, Ext[2]], Index[Gluon, Ext[4]]] - 27*IPL[{G}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[4]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[3]]]*SUNF[Index[Gluon, Ext[1]], Index[Gluon, Ext[3]], Index[Gluon, Ext[2]], Index[Gluon, Ext[4]]] - 4*IPL[{s}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[4]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[3]]]*SUNF[Index[Gluon, Ext[1]], Index[Gluon, Ext[3]], Index[Gluon, Ext[2]], Index[Gluon, Ext[4]]] - 4*IPL[{t}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[4]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[3]]]*SUNF[Index[Gluon, Ext[1]], Index[Gluon, Ext[3]], Index[Gluon, Ext[2]], Index[Gluon, Ext[4]]] - 4*IPL[{u}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[4]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[3]]]*SUNF[Index[Gluon, Ext[1]], Index[Gluon, Ext[3]], Index[Gluon, Ext[2]], Index[Gluon, Ext[4]]] - 4*IPL[{b}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[3]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[4]]]*SUNF[Index[Gluon, Ext[1]], Index[Gluon, Ext[3]], Index[Gluon, Ext[2]], Index[Gluon, Ext[4]]] - 4*IPL[{c}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[3]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[4]]]*SUNF[Index[Gluon, Ext[1]], Index[Gluon, Ext[3]], Index[Gluon, Ext[2]], Index[Gluon, Ext[4]]] - 4*IPL[{d}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[3]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[4]]]*SUNF[Index[Gluon, Ext[1]], Index[Gluon, Ext[3]], Index[Gluon, Ext[2]], Index[Gluon, Ext[4]]] - 36*IPL[{G}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[3]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[4]]]*SUNF[Index[Gluon, Ext[1]], Index[Gluon, Ext[3]], Index[Gluon, Ext[2]], Index[Gluon, Ext[4]]] - 4*IPL[{s}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[3]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[4]]]*SUNF[Index[Gluon, Ext[1]], Index[Gluon, Ext[3]], Index[Gluon, Ext[2]], Index[Gluon, Ext[4]]] - 4*IPL[{t}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[3]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[4]]]*SUNF[Index[Gluon, Ext[1]], Index[Gluon, Ext[3]], Index[Gluon, Ext[2]], Index[Gluon, Ext[4]]] - 4*IPL[{u}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[3]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[4]]]*SUNF[Index[Gluon, Ext[1]], Index[Gluon, Ext[3]], Index[Gluon, Ext[2]], Index[Gluon, Ext[4]]] + 12*IPL[{b}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[2]]]*ME[Index[Lorentz, Ext[3]], Index[Lorentz, Ext[4]]]*SUNF[Index[Gluon, Ext[1]], Index[Gluon, Ext[3]], Index[Gluon, Ext[2]], Index[Gluon, Ext[4]]] + 12*IPL[{c}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[2]]]*ME[Index[Lorentz, Ext[3]], Index[Lorentz, Ext[4]]]*SUNF[Index[Gluon, Ext[1]], Index[Gluon, Ext[3]], Index[Gluon, Ext[2]], Index[Gluon, Ext[4]]] + 12*IPL[{d}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[2]]]*ME[Index[Lorentz, Ext[3]], Index[Lorentz, Ext[4]]]*SUNF[Index[Gluon, Ext[1]], Index[Gluon, Ext[3]], Index[Gluon, Ext[2]], Index[Gluon, Ext[4]]] - 3*IPL[{G}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[2]]]*ME[Index[Lorentz, Ext[3]], Index[Lorentz, Ext[4]]]*SUNF[Index[Gluon, Ext[1]], Index[Gluon, Ext[3]], Index[Gluon, Ext[2]], Index[Gluon, Ext[4]]] + 12*IPL[{s}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[2]]]*ME[Index[Lorentz, Ext[3]], Index[Lorentz, Ext[4]]]*SUNF[Index[Gluon, Ext[1]], Index[Gluon, Ext[3]], Index[Gluon, Ext[2]], Index[Gluon, Ext[4]]] + 12*IPL[{t}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[2]]]*ME[Index[Lorentz, Ext[3]], Index[Lorentz, Ext[4]]]*SUNF[Index[Gluon, Ext[1]], Index[Gluon, Ext[3]], Index[Gluon, Ext[2]], Index[Gluon, Ext[4]]] + 12*IPL[{u}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[2]]]*ME[Index[Lorentz, Ext[3]], Index[Lorentz, Ext[4]]]*SUNF[Index[Gluon, Ext[1]], Index[Gluon, Ext[3]], Index[Gluon, Ext[2]], Index[Gluon, Ext[4]]] + 9*IPL[{G}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[3]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[4]]]*SUNF[Index[Gluon, Ext[1]], Index[Gluon, Ext[4]], Index[Gluon, Ext[2]], Index[Gluon, Ext[3]]] + 45*IPL[{G}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[2]]]*ME[Index[Lorentz, Ext[3]], Index[Lorentz, Ext[4]]]*SUNF[Index[Gluon, Ext[1]], Index[Gluon, Ext[4]], Index[Gluon, Ext[2]], Index[Gluon, Ext[3]]] + 4*IPL[{b}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[4]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[3]]]*SUNT[Index[Gluon, Ext[1]], Index[Gluon, Ext[3]], Index[Gluon, Ext[2]], Index[Gluon, Ext[4]]] + 4*IPL[{c}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[4]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[3]]]*SUNT[Index[Gluon, Ext[1]], Index[Gluon, Ext[3]], Index[Gluon, Ext[2]], Index[Gluon, Ext[4]]] + 4*IPL[{d}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[4]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[3]]]*SUNT[Index[Gluon, Ext[1]], Index[Gluon, Ext[3]], Index[Gluon, Ext[2]], Index[Gluon, Ext[4]]] - 12*IPL[{G}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[4]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[3]]]*SUNT[Index[Gluon, Ext[1]], Index[Gluon, Ext[3]], Index[Gluon, Ext[2]], Index[Gluon, Ext[4]]] + 4*IPL[{s}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[4]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[3]]]*SUNT[Index[Gluon, Ext[1]], Index[Gluon, Ext[3]], Index[Gluon, Ext[2]], Index[Gluon, Ext[4]]] + 4*IPL[{t}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[4]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[3]]]*SUNT[Index[Gluon, Ext[1]], Index[Gluon, Ext[3]], Index[Gluon, Ext[2]], Index[Gluon, Ext[4]]] + 4*IPL[{u}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[4]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[3]]]*SUNT[Index[Gluon, Ext[1]], Index[Gluon, Ext[3]], Index[Gluon, Ext[2]], Index[Gluon, Ext[4]]] + 4*IPL[{b}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[3]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[4]]]*SUNT[Index[Gluon, Ext[1]], Index[Gluon, Ext[3]], Index[Gluon, Ext[2]], Index[Gluon, Ext[4]]] + 4*IPL[{c}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[3]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[4]]]*SUNT[Index[Gluon, Ext[1]], Index[Gluon, Ext[3]], Index[Gluon, Ext[2]], Index[Gluon, Ext[4]]] + 4*IPL[{d}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[3]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[4]]]*SUNT[Index[Gluon, Ext[1]], Index[Gluon, Ext[3]], Index[Gluon, Ext[2]], Index[Gluon, Ext[4]]] - 12*IPL[{G}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[3]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[4]]]*SUNT[Index[Gluon, Ext[1]], Index[Gluon, Ext[3]], Index[Gluon, Ext[2]], Index[Gluon, Ext[4]]] + 4*IPL[{s}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[3]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[4]]]*SUNT[Index[Gluon, Ext[1]], Index[Gluon, Ext[3]], Index[Gluon, Ext[2]], Index[Gluon, Ext[4]]] + 4*IPL[{t}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[3]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[4]]]*SUNT[Index[Gluon, Ext[1]], Index[Gluon, Ext[3]], Index[Gluon, Ext[2]], Index[Gluon, Ext[4]]] + 4*IPL[{u}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[3]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[4]]]*SUNT[Index[Gluon, Ext[1]], Index[Gluon, Ext[3]], Index[Gluon, Ext[2]], Index[Gluon, Ext[4]]] + 4*IPL[{b}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[2]]]*ME[Index[Lorentz, Ext[3]], Index[Lorentz, Ext[4]]]*SUNT[Index[Gluon, Ext[1]], Index[Gluon, Ext[3]], Index[Gluon, Ext[2]], Index[Gluon, Ext[4]]] + 4*IPL[{c}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[2]]]*ME[Index[Lorentz, Ext[3]], Index[Lorentz, Ext[4]]]*SUNT[Index[Gluon, Ext[1]], Index[Gluon, Ext[3]], Index[Gluon, Ext[2]], Index[Gluon, Ext[4]]] + 4*IPL[{d}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[2]]]*ME[Index[Lorentz, Ext[3]], Index[Lorentz, Ext[4]]]*SUNT[Index[Gluon, Ext[1]], Index[Gluon, Ext[3]], Index[Gluon, Ext[2]], Index[Gluon, Ext[4]]] - 12*IPL[{G}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[2]]]*ME[Index[Lorentz, Ext[3]], Index[Lorentz, Ext[4]]]*SUNT[Index[Gluon, Ext[1]], Index[Gluon, Ext[3]], Index[Gluon, Ext[2]], Index[Gluon, Ext[4]]] + 4*IPL[{s}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[2]]]*ME[Index[Lorentz, Ext[3]], Index[Lorentz, Ext[4]]]*SUNT[Index[Gluon, Ext[1]], Index[Gluon, Ext[3]], Index[Gluon, Ext[2]], Index[Gluon, Ext[4]]] + 4*IPL[{t}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[2]]]*ME[Index[Lorentz, Ext[3]], Index[Lorentz, Ext[4]]]*SUNT[Index[Gluon, Ext[1]], Index[Gluon, Ext[3]], Index[Gluon, Ext[2]], Index[Gluon, Ext[4]]] + 4*IPL[{u}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[2]]]*ME[Index[Lorentz, Ext[3]], Index[Lorentz, Ext[4]]]*SUNT[Index[Gluon, Ext[1]], Index[Gluon, Ext[3]], Index[Gluon, Ext[2]], Index[Gluon, Ext[4]]] + 4*IPL[{b}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[4]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[3]]]*SUNT[Index[Gluon, Ext[1]], Index[Gluon, Ext[4]], Index[Gluon, Ext[2]], Index[Gluon, Ext[3]]] + 4*IPL[{c}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[4]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[3]]]*SUNT[Index[Gluon, Ext[1]], Index[Gluon, Ext[4]], Index[Gluon, Ext[2]], Index[Gluon, Ext[3]]] + 4*IPL[{d}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[4]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[3]]]*SUNT[Index[Gluon, Ext[1]], Index[Gluon, Ext[4]], Index[Gluon, Ext[2]], Index[Gluon, Ext[3]]] - 12*IPL[{G}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[4]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[3]]]*SUNT[Index[Gluon, Ext[1]], Index[Gluon, Ext[4]], Index[Gluon, Ext[2]], Index[Gluon, Ext[3]]] + 4*IPL[{s}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[4]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[3]]]*SUNT[Index[Gluon, Ext[1]], Index[Gluon, Ext[4]], Index[Gluon, Ext[2]], Index[Gluon, Ext[3]]] + 4*IPL[{t}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[4]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[3]]]*SUNT[Index[Gluon, Ext[1]], Index[Gluon, Ext[4]], Index[Gluon, Ext[2]], Index[Gluon, Ext[3]]] + 4*IPL[{u}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[4]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[3]]]*SUNT[Index[Gluon, Ext[1]], Index[Gluon, Ext[4]], Index[Gluon, Ext[2]], Index[Gluon, Ext[3]]] + 4*IPL[{b}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[3]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[4]]]*SUNT[Index[Gluon, Ext[1]], Index[Gluon, Ext[4]], Index[Gluon, Ext[2]], Index[Gluon, Ext[3]]] + 4*IPL[{c}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[3]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[4]]]*SUNT[Index[Gluon, Ext[1]], Index[Gluon, Ext[4]], Index[Gluon, Ext[2]], Index[Gluon, Ext[3]]] + 4*IPL[{d}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[3]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[4]]]*SUNT[Index[Gluon, Ext[1]], Index[Gluon, Ext[4]], Index[Gluon, Ext[2]], Index[Gluon, Ext[3]]] - 12*IPL[{G}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[3]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[4]]]*SUNT[Index[Gluon, Ext[1]], Index[Gluon, Ext[4]], Index[Gluon, Ext[2]], Index[Gluon, Ext[3]]] + 4*IPL[{s}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[3]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[4]]]*SUNT[Index[Gluon, Ext[1]], Index[Gluon, Ext[4]], Index[Gluon, Ext[2]], Index[Gluon, Ext[3]]] + 4*IPL[{t}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[3]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[4]]]*SUNT[Index[Gluon, Ext[1]], Index[Gluon, Ext[4]], Index[Gluon, Ext[2]], Index[Gluon, Ext[3]]] + 4*IPL[{u}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[3]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[4]]]*SUNT[Index[Gluon, Ext[1]], Index[Gluon, Ext[4]], Index[Gluon, Ext[2]], Index[Gluon, Ext[3]]] + 4*IPL[{b}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[2]]]*ME[Index[Lorentz, Ext[3]], Index[Lorentz, Ext[4]]]*SUNT[Index[Gluon, Ext[1]], Index[Gluon, Ext[4]], Index[Gluon, Ext[2]], Index[Gluon, Ext[3]]] + 4*IPL[{c}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[2]]]*ME[Index[Lorentz, Ext[3]], Index[Lorentz, Ext[4]]]*SUNT[Index[Gluon, Ext[1]], Index[Gluon, Ext[4]], Index[Gluon, Ext[2]], Index[Gluon, Ext[3]]] + 4*IPL[{d}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[2]]]*ME[Index[Lorentz, Ext[3]], Index[Lorentz, Ext[4]]]*SUNT[Index[Gluon, Ext[1]], Index[Gluon, Ext[4]], Index[Gluon, Ext[2]], Index[Gluon, Ext[3]]] - 12*IPL[{G}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[2]]]*ME[Index[Lorentz, Ext[3]], Index[Lorentz, Ext[4]]]*SUNT[Index[Gluon, Ext[1]], Index[Gluon, Ext[4]], Index[Gluon, Ext[2]], Index[Gluon, Ext[3]]] + 4*IPL[{s}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[2]]]*ME[Index[Lorentz, Ext[3]], Index[Lorentz, Ext[4]]]*SUNT[Index[Gluon, Ext[1]], Index[Gluon, Ext[4]], Index[Gluon, Ext[2]], Index[Gluon, Ext[3]]] + 4*IPL[{t}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[2]]]*ME[Index[Lorentz, Ext[3]], Index[Lorentz, Ext[4]]]*SUNT[Index[Gluon, Ext[1]], Index[Gluon, Ext[4]], Index[Gluon, Ext[2]], Index[Gluon, Ext[3]]] + 4*IPL[{u}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[2]]]*ME[Index[Lorentz, Ext[3]], Index[Lorentz, Ext[4]]]*SUNT[Index[Gluon, Ext[1]], Index[Gluon, Ext[4]], Index[Gluon, Ext[2]], Index[Gluon, Ext[3]]]))/Pi^2},
{{{anti[W], 1}, {W, 2}, {G, 3}, {G, 4}}, ((I/96)*ee^2*gs^2*IndexDelta[Index[Gluon, Ext[3]], Index[Gluon, Ext[4]]]*(IPL[{c, s}] + IPL[{t, b}] + IPL[{u, d}])*(ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[4]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[3]]] + ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[3]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[4]]] + ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[2]]]*ME[Index[Lorentz, Ext[3]], Index[Lorentz, Ext[4]]]))/(Pi^2*sw^2)},
{{{A, 1}, {Z, 2}, {G, 3}, {G, 4}}, ((I/864)*ee^2*gs^2*IndexDelta[Index[Gluon, Ext[3]], Index[Gluon, Ext[4]]]*(3*cw^2*IPL[{b}] - sw^2*IPL[{b}] + 2*(3*cw^2 - 5*sw^2)*IPL[{c}] + 3*cw^2*IPL[{d}] - sw^2*IPL[{d}] + 3*cw^2*IPL[{s}] - sw^2*IPL[{s}] + 6*cw^2*IPL[{t}] - 10*sw^2*IPL[{t}] + 2*(3*cw^2 - 5*sw^2)*IPL[{u}])*(ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[4]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[3]]] + ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[3]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[4]]] + ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[2]]]*ME[Index[Lorentz, Ext[3]], Index[Lorentz, Ext[4]]]))/(cw*Pi^2*sw)},
{{{Z, 1}, {Z, 2}, {G, 3}, {G, 4}}, ((I/1728)*ee^2*gs^2*IndexDelta[Index[Gluon, Ext[3]], Index[Gluon, Ext[4]]]*(9*cw^4*IPL[{b}] + 6*cw^2*sw^2*IPL[{b}] + 5*sw^4*IPL[{b}] + (9*cw^4 - 6*cw^2*sw^2 + 17*sw^4)*IPL[{c}] + 9*cw^4*IPL[{d}] + 6*cw^2*sw^2*IPL[{d}] + 5*sw^4*IPL[{d}] + 9*cw^4*IPL[{s}] + 6*cw^2*sw^2*IPL[{s}] + 5*sw^4*IPL[{s}] + 9*cw^4*IPL[{t}] - 6*cw^2*sw^2*IPL[{t}] + 17*sw^4*IPL[{t}] + (9*cw^4 - 6*cw^2*sw^2 + 17*sw^4)*IPL[{u}])*(ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[4]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[3]]] + ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[3]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[4]]] + ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[2]]]*ME[Index[Lorentz, Ext[3]], Index[Lorentz, Ext[4]]]))/(cw^2*Pi^2*sw^2)},
{{{A, 1}, {A, 2}, {G, 3}, {G, 4}}, ((I/216)*ee^2*gs^2*IndexDelta[Index[Gluon, Ext[3]], Index[Gluon, Ext[4]]]*(IPL[{b}] + 4*IPL[{c}] + IPL[{d}] + IPL[{s}] + 4*IPL[{t}] + 4*IPL[{u}])*(ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[4]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[3]]] + ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[3]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[4]]] + ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[2]]]*ME[Index[Lorentz, Ext[3]], Index[Lorentz, Ext[4]]]))/Pi^2},
{{{Z, 1}, {G, 2}, {G, 3}, {G, 4}}, ((I/576)*ee*gs^3*(-27*(cw^2 + sw^2)*Eps[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[2]], Index[Lorentz, Ext[3]], Index[Lorentz, Ext[4]]]*(-IPL[{b}] + IPL[{c}] - IPL[{d}] - IPL[{s}] + IPL[{t}] + IPL[{u}])*SUNF[Index[Gluon, Ext[2]], Index[Gluon, Ext[3]], Index[Gluon, Ext[4]]] - I*(-3*cw^2*IPL[{b}] + sw^2*IPL[{b}] + (3*cw^2 - 5*sw^2)*IPL[{c}] - 3*cw^2*IPL[{d}] + sw^2*IPL[{d}] - 3*cw^2*IPL[{s}] + sw^2*IPL[{s}] + 3*cw^2*IPL[{t}] - 5*sw^2*IPL[{t}] + (3*cw^2 - 5*sw^2)*IPL[{u}])*(ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[4]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[3]]] + ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[3]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[4]]] + ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[2]]]*ME[Index[Lorentz, Ext[3]], Index[Lorentz, Ext[4]]])*(SUNF[Index[Gluon, Ext[2]], Index[Gluon, Ext[3]], Index[Gluon, Ext[4]]] + (4*I)*SUNT[Index[Gluon, Ext[2]], Index[Gluon, Ext[3]], Index[Gluon, Ext[4]]])))/(cw*Pi^2*sw)},
{{{A, 1}, {G, 2}, {G, 3}, {G, 4}}, ((I/144)*ee*gs^3*(-IPL[{b}] + 2*IPL[{c}] - IPL[{d}] - IPL[{s}] + 2*IPL[{t}] + 2*IPL[{u}])*(ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[4]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[3]]] + ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[3]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[4]]] + ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[2]]]*ME[Index[Lorentz, Ext[3]], Index[Lorentz, Ext[4]]])*((-I)*SUNF[Index[Gluon, Ext[2]], Index[Gluon, Ext[3]], Index[Gluon, Ext[4]]] + 4*SUNT[Index[Gluon, Ext[2]], Index[Gluon, Ext[3]], Index[Gluon, Ext[4]]]))/Pi^2},
{{{G, 1}, {G, 2}, {H, 3}, {H, 4}}, ((-I/8)*gs^2*IndexDelta[Index[Gluon, Ext[1]], Index[Gluon, Ext[2]]]*(MB^2*IPL[{b}] + MC^2*IPL[{c}] + MT^2*IPL[{t}])*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[2]]])/(Pi^2*vev^2)},
{{{G, 1}, {G, 2}, {G0, 3}, {G0, 4}}, ((-I/8)*gs^2*IndexDelta[Index[Gluon, Ext[1]], Index[Gluon, Ext[2]]]*(MB^2*IPL[{b}] + MC^2*IPL[{c}] + MT^2*IPL[{t}])*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[2]]])/(Pi^2*vev^2)},
{{{G, 1}, {G, 2}, {anti[GP], 3}, {GP, 4}}, ((-I/8)*gs^2*IndexDelta[Index[Gluon, Ext[1]], Index[Gluon, Ext[2]]]*(MC^2*IPL[{c, s}] + (MB^2 + MT^2)*IPL[{t, b}])*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[2]]])/(Pi^2*vev^2)}};

(*UV vertices*)

UV$vertlist = {
{{{anti[u], 1}, {u, 2}}, ((-I/12)*gs^2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{u, G}]*(TensDot[SlashedP[2], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]] + TensDot[SlashedP[2], ProjP][Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(FR$Eps*Pi^2) + I*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*((gs^2*IPL[{u, G}]*TensDot[SlashedP[2], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]])/(12*FR$Eps*Pi^2) + (gs^2*IPL[{u, G}]*TensDot[SlashedP[2], ProjP][Index[Spin, Ext[1]], Index[Spin, Ext[2]]])/(12*FR$Eps*Pi^2))},
{{{anti[c], 1}, {c, 2}}, ((I/12)*gs^2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{c, G}]*(4*MC*ProjM[Index[Spin, Ext[1]], Index[Spin, Ext[2]]] + 4*MC*ProjP[Index[Spin, Ext[1]], Index[Spin, Ext[2]]] - TensDot[SlashedP[2], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]] - TensDot[SlashedP[2], ProjP][Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(FR$Eps*Pi^2) + (I/2)*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*(-((MC*((gs^2*IPL[{c, G}])/(6*Pi^2) - (5*gs^2*If[MC == 0, 0, 1]*IPL[{c, G}])/(6*Pi^2) - (gs^2*If[MC == 0, 0, 1]*IPL[{c, G}])/(3*FR$Eps*Pi^2) - (gs^2*If[MC == 0, 1, 0]*IPL[{c, G}])/(6*Pi^2) + (gs^2*If[MC == 0, 1, 0]*IPL[{c, G}])/(6*FR$Eps*Pi^2) + (gs^2*If[MC == 0, 0, 1]*IPL[{c, G}]*Log[MC/FR$MU])/Pi^2) + 2*((gs^2*MC*IPL[{c, G}])/(12*Pi^2) - (5*gs^2*MC*If[MC == 0, 0, 1]*IPL[{c, G}])/(12*Pi^2) - (gs^2*MC*If[MC == 0, 1, 0]*IPL[{c, G}])/(12*Pi^2) + (gs^2*MC*If[MC == 0, 1, 0]*IPL[{c, G}])/(4*FR$Eps*Pi^2) + (gs^2*MC*If[MC == 0, 0, 1]*IPL[{c, G}]*Log[MC/FR$MU])/(2*Pi^2)))*ProjM[Index[Spin, Ext[1]], Index[Spin, Ext[2]]]) - (MC*((gs^2*IPL[{c, G}])/(6*Pi^2) - (5*gs^2*If[MC == 0, 0, 1]*IPL[{c, G}])/(6*Pi^2) - (gs^2*If[MC == 0, 0, 1]*IPL[{c, G}])/(3*FR$Eps*Pi^2) - (gs^2*If[MC == 0, 1, 0]*IPL[{c, G}])/(6*Pi^2) + (gs^2*If[MC == 0, 1, 0]*IPL[{c, G}])/(6*FR$Eps*Pi^2) + (gs^2*If[MC == 0, 0, 1]*IPL[{c, G}]*Log[MC/FR$MU])/Pi^2) + 2*((gs^2*MC*IPL[{c, G}])/(12*Pi^2) - (5*gs^2*MC*If[MC == 0, 0, 1]*IPL[{c, G}])/(12*Pi^2) - (gs^2*MC*If[MC == 0, 1, 0]*IPL[{c, G}])/(12*Pi^2) + (gs^2*MC*If[MC == 0, 1, 0]*IPL[{c, G}])/(4*FR$Eps*Pi^2) + (gs^2*MC*If[MC == 0, 0, 1]*IPL[{c, G}]*Log[MC/FR$MU])/(2*Pi^2)))*ProjP[Index[Spin, Ext[1]], Index[Spin, Ext[2]]] + 2*((gs^2*IPL[{c, G}])/(12*Pi^2) - (5*gs^2*If[MC == 0, 0, 1]*IPL[{c, G}])/(12*Pi^2) - (gs^2*If[MC == 0, 0, 1]*IPL[{c, G}])/(6*FR$Eps*Pi^2) - (gs^2*If[MC == 0, 1, 0]*IPL[{c, G}])/(12*Pi^2) + (gs^2*If[MC == 0, 1, 0]*IPL[{c, G}])/(12*FR$Eps*Pi^2) + (gs^2*If[MC == 0, 0, 1]*IPL[{c, G}]*Log[MC/FR$MU])/(2*Pi^2))*TensDot[SlashedP[2], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]] + 2*((gs^2*IPL[{c, G}])/(12*Pi^2) - (5*gs^2*If[MC == 0, 0, 1]*IPL[{c, G}])/(12*Pi^2) - (gs^2*If[MC == 0, 0, 1]*IPL[{c, G}])/(6*FR$Eps*Pi^2) - (gs^2*If[MC == 0, 1, 0]*IPL[{c, G}])/(12*Pi^2) + (gs^2*If[MC == 0, 1, 0]*IPL[{c, G}])/(12*FR$Eps*Pi^2) + (gs^2*If[MC == 0, 0, 1]*IPL[{c, G}]*Log[MC/FR$MU])/(2*Pi^2))*TensDot[SlashedP[2], ProjP][Index[Spin, Ext[1]], Index[Spin, Ext[2]]])},
{{{anti[t], 1}, {t, 2}}, ((I/12)*gs^2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{t, G}]*(4*MT*ProjM[Index[Spin, Ext[1]], Index[Spin, Ext[2]]] + 4*MT*ProjP[Index[Spin, Ext[1]], Index[Spin, Ext[2]]] - TensDot[SlashedP[2], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]] - TensDot[SlashedP[2], ProjP][Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(FR$Eps*Pi^2) + (I/2)*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*(-((MT*((gs^2*IPL[{t, G}])/(6*Pi^2) - (5*gs^2*If[MT == 0, 0, 1]*IPL[{t, G}])/(6*Pi^2) - (gs^2*If[MT == 0, 0, 1]*IPL[{t, G}])/(3*FR$Eps*Pi^2) - (gs^2*If[MT == 0, 1, 0]*IPL[{t, G}])/(6*Pi^2) + (gs^2*If[MT == 0, 1, 0]*IPL[{t, G}])/(6*FR$Eps*Pi^2) + (gs^2*If[MT == 0, 0, 1]*IPL[{t, G}]*Log[MT/FR$MU])/Pi^2) + 2*((gs^2*MT*IPL[{t, G}])/(12*Pi^2) - (5*gs^2*MT*If[MT == 0, 0, 1]*IPL[{t, G}])/(12*Pi^2) - (gs^2*MT*If[MT == 0, 1, 0]*IPL[{t, G}])/(12*Pi^2) + (gs^2*MT*If[MT == 0, 1, 0]*IPL[{t, G}])/(4*FR$Eps*Pi^2) + (gs^2*MT*If[MT == 0, 0, 1]*IPL[{t, G}]*Log[MT/FR$MU])/(2*Pi^2)))*ProjM[Index[Spin, Ext[1]], Index[Spin, Ext[2]]]) - (MT*((gs^2*IPL[{t, G}])/(6*Pi^2) - (5*gs^2*If[MT == 0, 0, 1]*IPL[{t, G}])/(6*Pi^2) - (gs^2*If[MT == 0, 0, 1]*IPL[{t, G}])/(3*FR$Eps*Pi^2) - (gs^2*If[MT == 0, 1, 0]*IPL[{t, G}])/(6*Pi^2) + (gs^2*If[MT == 0, 1, 0]*IPL[{t, G}])/(6*FR$Eps*Pi^2) + (gs^2*If[MT == 0, 0, 1]*IPL[{t, G}]*Log[MT/FR$MU])/Pi^2) + 2*((gs^2*MT*IPL[{t, G}])/(12*Pi^2) - (5*gs^2*MT*If[MT == 0, 0, 1]*IPL[{t, G}])/(12*Pi^2) - (gs^2*MT*If[MT == 0, 1, 0]*IPL[{t, G}])/(12*Pi^2) + (gs^2*MT*If[MT == 0, 1, 0]*IPL[{t, G}])/(4*FR$Eps*Pi^2) + (gs^2*MT*If[MT == 0, 0, 1]*IPL[{t, G}]*Log[MT/FR$MU])/(2*Pi^2)))*ProjP[Index[Spin, Ext[1]], Index[Spin, Ext[2]]] + 2*((gs^2*IPL[{t, G}])/(12*Pi^2) - (5*gs^2*If[MT == 0, 0, 1]*IPL[{t, G}])/(12*Pi^2) - (gs^2*If[MT == 0, 0, 1]*IPL[{t, G}])/(6*FR$Eps*Pi^2) - (gs^2*If[MT == 0, 1, 0]*IPL[{t, G}])/(12*Pi^2) + (gs^2*If[MT == 0, 1, 0]*IPL[{t, G}])/(12*FR$Eps*Pi^2) + (gs^2*If[MT == 0, 0, 1]*IPL[{t, G}]*Log[MT/FR$MU])/(2*Pi^2))*TensDot[SlashedP[2], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]] + 2*((gs^2*IPL[{t, G}])/(12*Pi^2) - (5*gs^2*If[MT == 0, 0, 1]*IPL[{t, G}])/(12*Pi^2) - (gs^2*If[MT == 0, 0, 1]*IPL[{t, G}])/(6*FR$Eps*Pi^2) - (gs^2*If[MT == 0, 1, 0]*IPL[{t, G}])/(12*Pi^2) + (gs^2*If[MT == 0, 1, 0]*IPL[{t, G}])/(12*FR$Eps*Pi^2) + (gs^2*If[MT == 0, 0, 1]*IPL[{t, G}]*Log[MT/FR$MU])/(2*Pi^2))*TensDot[SlashedP[2], ProjP][Index[Spin, Ext[1]], Index[Spin, Ext[2]]])},
{{{anti[d], 1}, {d, 2}}, ((-I/12)*gs^2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{d, G}]*(TensDot[SlashedP[2], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]] + TensDot[SlashedP[2], ProjP][Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(FR$Eps*Pi^2) + I*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*((gs^2*IPL[{d, G}]*TensDot[SlashedP[2], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]])/(12*FR$Eps*Pi^2) + (gs^2*IPL[{d, G}]*TensDot[SlashedP[2], ProjP][Index[Spin, Ext[1]], Index[Spin, Ext[2]]])/(12*FR$Eps*Pi^2))},
{{{anti[s], 1}, {s, 2}}, ((-I/12)*gs^2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{s, G}]*(TensDot[SlashedP[2], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]] + TensDot[SlashedP[2], ProjP][Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(FR$Eps*Pi^2) + I*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*((gs^2*IPL[{s, G}]*TensDot[SlashedP[2], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]])/(12*FR$Eps*Pi^2) + (gs^2*IPL[{s, G}]*TensDot[SlashedP[2], ProjP][Index[Spin, Ext[1]], Index[Spin, Ext[2]]])/(12*FR$Eps*Pi^2))},
{{{anti[b], 1}, {b, 2}}, ((I/12)*gs^2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{b, G}]*(4*MB*ProjM[Index[Spin, Ext[1]], Index[Spin, Ext[2]]] + 4*MB*ProjP[Index[Spin, Ext[1]], Index[Spin, Ext[2]]] - TensDot[SlashedP[2], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]] - TensDot[SlashedP[2], ProjP][Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(FR$Eps*Pi^2) + (I/2)*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*(-((MB*((gs^2*IPL[{b, G}])/(6*Pi^2) - (5*gs^2*If[MB == 0, 0, 1]*IPL[{b, G}])/(6*Pi^2) - (gs^2*If[MB == 0, 0, 1]*IPL[{b, G}])/(3*FR$Eps*Pi^2) - (gs^2*If[MB == 0, 1, 0]*IPL[{b, G}])/(6*Pi^2) + (gs^2*If[MB == 0, 1, 0]*IPL[{b, G}])/(6*FR$Eps*Pi^2) + (gs^2*If[MB == 0, 0, 1]*IPL[{b, G}]*Log[MB/FR$MU])/Pi^2) + 2*((gs^2*MB*IPL[{b, G}])/(12*Pi^2) - (5*gs^2*MB*If[MB == 0, 0, 1]*IPL[{b, G}])/(12*Pi^2) - (gs^2*MB*If[MB == 0, 1, 0]*IPL[{b, G}])/(12*Pi^2) + (gs^2*MB*If[MB == 0, 1, 0]*IPL[{b, G}])/(4*FR$Eps*Pi^2) + (gs^2*MB*If[MB == 0, 0, 1]*IPL[{b, G}]*Log[MB/FR$MU])/(2*Pi^2)))*ProjM[Index[Spin, Ext[1]], Index[Spin, Ext[2]]]) - (MB*((gs^2*IPL[{b, G}])/(6*Pi^2) - (5*gs^2*If[MB == 0, 0, 1]*IPL[{b, G}])/(6*Pi^2) - (gs^2*If[MB == 0, 0, 1]*IPL[{b, G}])/(3*FR$Eps*Pi^2) - (gs^2*If[MB == 0, 1, 0]*IPL[{b, G}])/(6*Pi^2) + (gs^2*If[MB == 0, 1, 0]*IPL[{b, G}])/(6*FR$Eps*Pi^2) + (gs^2*If[MB == 0, 0, 1]*IPL[{b, G}]*Log[MB/FR$MU])/Pi^2) + 2*((gs^2*MB*IPL[{b, G}])/(12*Pi^2) - (5*gs^2*MB*If[MB == 0, 0, 1]*IPL[{b, G}])/(12*Pi^2) - (gs^2*MB*If[MB == 0, 1, 0]*IPL[{b, G}])/(12*Pi^2) + (gs^2*MB*If[MB == 0, 1, 0]*IPL[{b, G}])/(4*FR$Eps*Pi^2) + (gs^2*MB*If[MB == 0, 0, 1]*IPL[{b, G}]*Log[MB/FR$MU])/(2*Pi^2)))*ProjP[Index[Spin, Ext[1]], Index[Spin, Ext[2]]] + 2*((gs^2*IPL[{b, G}])/(12*Pi^2) - (5*gs^2*If[MB == 0, 0, 1]*IPL[{b, G}])/(12*Pi^2) - (gs^2*If[MB == 0, 0, 1]*IPL[{b, G}])/(6*FR$Eps*Pi^2) - (gs^2*If[MB == 0, 1, 0]*IPL[{b, G}])/(12*Pi^2) + (gs^2*If[MB == 0, 1, 0]*IPL[{b, G}])/(12*FR$Eps*Pi^2) + (gs^2*If[MB == 0, 0, 1]*IPL[{b, G}]*Log[MB/FR$MU])/(2*Pi^2))*TensDot[SlashedP[2], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]] + 2*((gs^2*IPL[{b, G}])/(12*Pi^2) - (5*gs^2*If[MB == 0, 0, 1]*IPL[{b, G}])/(12*Pi^2) - (gs^2*If[MB == 0, 0, 1]*IPL[{b, G}])/(6*FR$Eps*Pi^2) - (gs^2*If[MB == 0, 1, 0]*IPL[{b, G}])/(12*Pi^2) + (gs^2*If[MB == 0, 1, 0]*IPL[{b, G}])/(12*FR$Eps*Pi^2) + (gs^2*If[MB == 0, 0, 1]*IPL[{b, G}]*Log[MB/FR$MU])/(2*Pi^2))*TensDot[SlashedP[2], ProjP][Index[Spin, Ext[1]], Index[Spin, Ext[2]]])},
{{{G, 1}, {G, 2}}, I*IndexDelta[Index[Gluon, Ext[1]], Index[Gluon, Ext[2]]]*((gs^2*If[MB == 0, 1, 0]*IPL[{b}])/(24*FR$Eps*Pi^2) + (gs^2*If[MC == 0, 1, 0]*IPL[{c}])/(24*FR$Eps*Pi^2) + (gs^2*IPL[{d}])/(24*FR$Eps*Pi^2) - (19*gs^2*IPL[{G}])/(64*FR$Eps*Pi^2) - (gs^2*IPL[{ghG}])/(64*FR$Eps*Pi^2) + (gs^2*IPL[{s}])/(24*FR$Eps*Pi^2) + (gs^2*If[MT == 0, 1, 0]*IPL[{t}])/(24*FR$Eps*Pi^2) + (gs^2*IPL[{u}])/(24*FR$Eps*Pi^2) + (gs^2*If[MB == 0, 0, 1]*IPL[{b}]*Log[MB/FR$MU])/(12*Pi^2) + (gs^2*If[MC == 0, 0, 1]*IPL[{c}]*Log[MC/FR$MU])/(12*Pi^2) + (gs^2*If[MT == 0, 0, 1]*IPL[{t}]*Log[MT/FR$MU])/(12*Pi^2))*(FV[2, Index[Lorentz, Ext[1]]]*FV[2, Index[Lorentz, Ext[2]]] - ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[2]]]*SP[2, 2]) - ((I/192)*gs^2*IndexDelta[Index[Gluon, Ext[1]], Index[Gluon, Ext[2]]]*(2*FV[2, Index[Lorentz, Ext[1]]]*FV[2, Index[Lorentz, Ext[2]]]*(4*IPL[{b}] + 4*IPL[{c}] + 4*IPL[{d}] - 33*IPL[{G}] + 3*IPL[{ghG}] + 4*IPL[{s}] + 4*IPL[{t}] + 4*IPL[{u}]) - (8*IPL[{b}] + 8*IPL[{c}] + 8*IPL[{d}] - 57*IPL[{G}] - 3*IPL[{ghG}] + 8*IPL[{s}] + 8*IPL[{t}] + 8*IPL[{u}])*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[2]]]*SP[2, 2]))/(FR$Eps*Pi^2)},
{{{anti[c], 1}, {c, 2}, {H, 3}}, ((I/3)*gs^2*MC*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{c, G}]*(ProjM[Index[Spin, Ext[1]], Index[Spin, Ext[2]]] + ProjP[Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(FR$Eps*Pi^2*vev) - ((I/2)*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*(MC*((gs^2*IPL[{c, G}])/(6*Pi^2) - (5*gs^2*If[MC == 0, 0, 1]*IPL[{c, G}])/(6*Pi^2) - (gs^2*If[MC == 0, 0, 1]*IPL[{c, G}])/(3*FR$Eps*Pi^2) - (gs^2*If[MC == 0, 1, 0]*IPL[{c, G}])/(6*Pi^2) + (gs^2*If[MC == 0, 1, 0]*IPL[{c, G}])/(6*FR$Eps*Pi^2) + (gs^2*If[MC == 0, 0, 1]*IPL[{c, G}]*Log[MC/FR$MU])/Pi^2) + 2*((gs^2*MC*IPL[{c, G}])/(12*Pi^2) - (5*gs^2*MC*If[MC == 0, 0, 1]*IPL[{c, G}])/(12*Pi^2) - (gs^2*MC*If[MC == 0, 1, 0]*IPL[{c, G}])/(12*Pi^2) + (gs^2*MC*If[MC == 0, 1, 0]*IPL[{c, G}])/(4*FR$Eps*Pi^2) + (gs^2*MC*If[MC == 0, 0, 1]*IPL[{c, G}]*Log[MC/FR$MU])/(2*Pi^2)))*(ProjM[Index[Spin, Ext[1]], Index[Spin, Ext[2]]] + ProjP[Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/vev},
{{{anti[t], 1}, {t, 2}, {H, 3}}, ((I/3)*gs^2*MT*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{t, G}]*(ProjM[Index[Spin, Ext[1]], Index[Spin, Ext[2]]] + ProjP[Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(FR$Eps*Pi^2*vev) - ((I/2)*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*(MT*((gs^2*IPL[{t, G}])/(6*Pi^2) - (5*gs^2*If[MT == 0, 0, 1]*IPL[{t, G}])/(6*Pi^2) - (gs^2*If[MT == 0, 0, 1]*IPL[{t, G}])/(3*FR$Eps*Pi^2) - (gs^2*If[MT == 0, 1, 0]*IPL[{t, G}])/(6*Pi^2) + (gs^2*If[MT == 0, 1, 0]*IPL[{t, G}])/(6*FR$Eps*Pi^2) + (gs^2*If[MT == 0, 0, 1]*IPL[{t, G}]*Log[MT/FR$MU])/Pi^2) + 2*((gs^2*MT*IPL[{t, G}])/(12*Pi^2) - (5*gs^2*MT*If[MT == 0, 0, 1]*IPL[{t, G}])/(12*Pi^2) - (gs^2*MT*If[MT == 0, 1, 0]*IPL[{t, G}])/(12*Pi^2) + (gs^2*MT*If[MT == 0, 1, 0]*IPL[{t, G}])/(4*FR$Eps*Pi^2) + (gs^2*MT*If[MT == 0, 0, 1]*IPL[{t, G}]*Log[MT/FR$MU])/(2*Pi^2)))*(ProjM[Index[Spin, Ext[1]], Index[Spin, Ext[2]]] + ProjP[Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/vev},
{{{anti[b], 1}, {b, 2}, {H, 3}}, ((I/3)*gs^2*MB*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{b, G}]*(ProjM[Index[Spin, Ext[1]], Index[Spin, Ext[2]]] + ProjP[Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(FR$Eps*Pi^2*vev) - ((I/2)*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*(MB*((gs^2*IPL[{b, G}])/(6*Pi^2) - (5*gs^2*If[MB == 0, 0, 1]*IPL[{b, G}])/(6*Pi^2) - (gs^2*If[MB == 0, 0, 1]*IPL[{b, G}])/(3*FR$Eps*Pi^2) - (gs^2*If[MB == 0, 1, 0]*IPL[{b, G}])/(6*Pi^2) + (gs^2*If[MB == 0, 1, 0]*IPL[{b, G}])/(6*FR$Eps*Pi^2) + (gs^2*If[MB == 0, 0, 1]*IPL[{b, G}]*Log[MB/FR$MU])/Pi^2) + 2*((gs^2*MB*IPL[{b, G}])/(12*Pi^2) - (5*gs^2*MB*If[MB == 0, 0, 1]*IPL[{b, G}])/(12*Pi^2) - (gs^2*MB*If[MB == 0, 1, 0]*IPL[{b, G}])/(12*Pi^2) + (gs^2*MB*If[MB == 0, 1, 0]*IPL[{b, G}])/(4*FR$Eps*Pi^2) + (gs^2*MB*If[MB == 0, 0, 1]*IPL[{b, G}]*Log[MB/FR$MU])/(2*Pi^2)))*(ProjM[Index[Spin, Ext[1]], Index[Spin, Ext[2]]] + ProjP[Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/vev},
{{{anti[c], 1}, {c, 2}, {G0, 3}}, -(gs^2*MC*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{c, G}]*(ProjM[Index[Spin, Ext[1]], Index[Spin, Ext[2]]] - ProjP[Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(3*FR$Eps*Pi^2*vev) + (IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*(MC*((gs^2*IPL[{c, G}])/(6*Pi^2) - (5*gs^2*If[MC == 0, 0, 1]*IPL[{c, G}])/(6*Pi^2) - (gs^2*If[MC == 0, 0, 1]*IPL[{c, G}])/(3*FR$Eps*Pi^2) - (gs^2*If[MC == 0, 1, 0]*IPL[{c, G}])/(6*Pi^2) + (gs^2*If[MC == 0, 1, 0]*IPL[{c, G}])/(6*FR$Eps*Pi^2) + (gs^2*If[MC == 0, 0, 1]*IPL[{c, G}]*Log[MC/FR$MU])/Pi^2) + 2*((gs^2*MC*IPL[{c, G}])/(12*Pi^2) - (5*gs^2*MC*If[MC == 0, 0, 1]*IPL[{c, G}])/(12*Pi^2) - (gs^2*MC*If[MC == 0, 1, 0]*IPL[{c, G}])/(12*Pi^2) + (gs^2*MC*If[MC == 0, 1, 0]*IPL[{c, G}])/(4*FR$Eps*Pi^2) + (gs^2*MC*If[MC == 0, 0, 1]*IPL[{c, G}]*Log[MC/FR$MU])/(2*Pi^2)))*(ProjM[Index[Spin, Ext[1]], Index[Spin, Ext[2]]] - ProjP[Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(2*vev)},
{{{anti[t], 1}, {t, 2}, {G0, 3}}, -(gs^2*MT*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{t, G}]*(ProjM[Index[Spin, Ext[1]], Index[Spin, Ext[2]]] - ProjP[Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(3*FR$Eps*Pi^2*vev) + (IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*(MT*((gs^2*IPL[{t, G}])/(6*Pi^2) - (5*gs^2*If[MT == 0, 0, 1]*IPL[{t, G}])/(6*Pi^2) - (gs^2*If[MT == 0, 0, 1]*IPL[{t, G}])/(3*FR$Eps*Pi^2) - (gs^2*If[MT == 0, 1, 0]*IPL[{t, G}])/(6*Pi^2) + (gs^2*If[MT == 0, 1, 0]*IPL[{t, G}])/(6*FR$Eps*Pi^2) + (gs^2*If[MT == 0, 0, 1]*IPL[{t, G}]*Log[MT/FR$MU])/Pi^2) + 2*((gs^2*MT*IPL[{t, G}])/(12*Pi^2) - (5*gs^2*MT*If[MT == 0, 0, 1]*IPL[{t, G}])/(12*Pi^2) - (gs^2*MT*If[MT == 0, 1, 0]*IPL[{t, G}])/(12*Pi^2) + (gs^2*MT*If[MT == 0, 1, 0]*IPL[{t, G}])/(4*FR$Eps*Pi^2) + (gs^2*MT*If[MT == 0, 0, 1]*IPL[{t, G}]*Log[MT/FR$MU])/(2*Pi^2)))*(ProjM[Index[Spin, Ext[1]], Index[Spin, Ext[2]]] - ProjP[Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(2*vev)},
{{{anti[b], 1}, {b, 2}, {G0, 3}}, (gs^2*MB*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{b, G}]*(ProjM[Index[Spin, Ext[1]], Index[Spin, Ext[2]]] - ProjP[Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(3*FR$Eps*Pi^2*vev) - (IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*(MB*((gs^2*IPL[{b, G}])/(6*Pi^2) - (5*gs^2*If[MB == 0, 0, 1]*IPL[{b, G}])/(6*Pi^2) - (gs^2*If[MB == 0, 0, 1]*IPL[{b, G}])/(3*FR$Eps*Pi^2) - (gs^2*If[MB == 0, 1, 0]*IPL[{b, G}])/(6*Pi^2) + (gs^2*If[MB == 0, 1, 0]*IPL[{b, G}])/(6*FR$Eps*Pi^2) + (gs^2*If[MB == 0, 0, 1]*IPL[{b, G}]*Log[MB/FR$MU])/Pi^2) + 2*((gs^2*MB*IPL[{b, G}])/(12*Pi^2) - (5*gs^2*MB*If[MB == 0, 0, 1]*IPL[{b, G}])/(12*Pi^2) - (gs^2*MB*If[MB == 0, 1, 0]*IPL[{b, G}])/(12*Pi^2) + (gs^2*MB*If[MB == 0, 1, 0]*IPL[{b, G}])/(4*FR$Eps*Pi^2) + (gs^2*MB*If[MB == 0, 0, 1]*IPL[{b, G}]*Log[MB/FR$MU])/(2*Pi^2)))*(ProjM[Index[Spin, Ext[1]], Index[Spin, Ext[2]]] - ProjP[Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(2*vev)},
{{{anti[s], 1}, {c, 2}, {anti[GP], 3}}, (Sqrt[2]*gs^2*MC*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{c, s, G}]*ProjP[Index[Spin, Ext[1]], Index[Spin, Ext[2]]])/(3*FR$Eps*Pi^2*vev) - (IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*(MC*((gs^2*IPL[{c, G}])/(12*Pi^2) - (5*gs^2*If[MC == 0, 0, 1]*IPL[{c, G}])/(12*Pi^2) - (gs^2*If[MC == 0, 0, 1]*IPL[{c, G}])/(6*FR$Eps*Pi^2) - (gs^2*If[MC == 0, 1, 0]*IPL[{c, G}])/(12*Pi^2) + (gs^2*If[MC == 0, 1, 0]*IPL[{c, G}])/(12*FR$Eps*Pi^2) + (gs^2*IPL[{s, G}])/(12*FR$Eps*Pi^2) + (gs^2*If[MC == 0, 0, 1]*IPL[{c, G}]*Log[MC/FR$MU])/(2*Pi^2)) + 2*((gs^2*MC*IPL[{c, G}])/(12*Pi^2) - (5*gs^2*MC*If[MC == 0, 0, 1]*IPL[{c, G}])/(12*Pi^2) - (gs^2*MC*If[MC == 0, 1, 0]*IPL[{c, G}])/(12*Pi^2) + (gs^2*MC*If[MC == 0, 1, 0]*IPL[{c, G}])/(4*FR$Eps*Pi^2) + (gs^2*MC*If[MC == 0, 0, 1]*IPL[{c, G}]*Log[MC/FR$MU])/(2*Pi^2)))*ProjP[Index[Spin, Ext[1]], Index[Spin, Ext[2]]])/(Sqrt[2]*vev)},
{{{anti[b], 1}, {t, 2}, {anti[GP], 3}}, -(Sqrt[2]*gs^2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{t, b, G}]*(MB*ProjM[Index[Spin, Ext[1]], Index[Spin, Ext[2]]] - MT*ProjP[Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(3*FR$Eps*Pi^2*vev) + (IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*((2*((gs^2*MB*IPL[{b, G}])/(12*Pi^2) - (5*gs^2*MB*If[MB == 0, 0, 1]*IPL[{b, G}])/(12*Pi^2) - (gs^2*MB*If[MB == 0, 1, 0]*IPL[{b, G}])/(12*Pi^2) + (gs^2*MB*If[MB == 0, 1, 0]*IPL[{b, G}])/(4*FR$Eps*Pi^2) + (gs^2*MB*If[MB == 0, 0, 1]*IPL[{b, G}]*Log[MB/FR$MU])/(2*Pi^2)) + MB*((gs^2*IPL[{b, G}])/(12*Pi^2) - (5*gs^2*If[MB == 0, 0, 1]*IPL[{b, G}])/(12*Pi^2) - (gs^2*If[MB == 0, 0, 1]*IPL[{b, G}])/(6*FR$Eps*Pi^2) - (gs^2*If[MB == 0, 1, 0]*IPL[{b, G}])/(12*Pi^2) + (gs^2*If[MB == 0, 1, 0]*IPL[{b, G}])/(12*FR$Eps*Pi^2) + (gs^2*IPL[{t, G}])/(12*Pi^2) - (5*gs^2*If[MT == 0, 0, 1]*IPL[{t, G}])/(12*Pi^2) - (gs^2*If[MT == 0, 0, 1]*IPL[{t, G}])/(6*FR$Eps*Pi^2) - (gs^2*If[MT == 0, 1, 0]*IPL[{t, G}])/(12*Pi^2) + (gs^2*If[MT == 0, 1, 0]*IPL[{t, G}])/(12*FR$Eps*Pi^2) + (gs^2*If[MB == 0, 0, 1]*IPL[{b, G}]*Log[MB/FR$MU])/(2*Pi^2) + (gs^2*If[MT == 0, 0, 1]*IPL[{t, G}]*Log[MT/FR$MU])/(2*Pi^2)))*ProjM[Index[Spin, Ext[1]], Index[Spin, Ext[2]]] - (MT*((gs^2*IPL[{b, G}])/(12*Pi^2) - (5*gs^2*If[MB == 0, 0, 1]*IPL[{b, G}])/(12*Pi^2) - (gs^2*If[MB == 0, 0, 1]*IPL[{b, G}])/(6*FR$Eps*Pi^2) - (gs^2*If[MB == 0, 1, 0]*IPL[{b, G}])/(12*Pi^2) + (gs^2*If[MB == 0, 1, 0]*IPL[{b, G}])/(12*FR$Eps*Pi^2) + (gs^2*IPL[{t, G}])/(12*Pi^2) - (5*gs^2*If[MT == 0, 0, 1]*IPL[{t, G}])/(12*Pi^2) - (gs^2*If[MT == 0, 0, 1]*IPL[{t, G}])/(6*FR$Eps*Pi^2) - (gs^2*If[MT == 0, 1, 0]*IPL[{t, G}])/(12*Pi^2) + (gs^2*If[MT == 0, 1, 0]*IPL[{t, G}])/(12*FR$Eps*Pi^2) + (gs^2*If[MB == 0, 0, 1]*IPL[{b, G}]*Log[MB/FR$MU])/(2*Pi^2) + (gs^2*If[MT == 0, 0, 1]*IPL[{t, G}]*Log[MT/FR$MU])/(2*Pi^2)) + 2*((gs^2*MT*IPL[{t, G}])/(12*Pi^2) - (5*gs^2*MT*If[MT == 0, 0, 1]*IPL[{t, G}])/(12*Pi^2) - (gs^2*MT*If[MT == 0, 1, 0]*IPL[{t, G}])/(12*Pi^2) + (gs^2*MT*If[MT == 0, 1, 0]*IPL[{t, G}])/(4*FR$Eps*Pi^2) + (gs^2*MT*If[MT == 0, 0, 1]*IPL[{t, G}]*Log[MT/FR$MU])/(2*Pi^2)))*ProjP[Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(Sqrt[2]*vev)},
{{{anti[c], 1}, {s, 2}, {GP, 3}}, -(Sqrt[2]*gs^2*MC*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{c, s, G}]*ProjM[Index[Spin, Ext[1]], Index[Spin, Ext[2]]])/(3*FR$Eps*Pi^2*vev) + (IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*(MC*((gs^2*IPL[{c, G}])/(12*Pi^2) - (5*gs^2*If[MC == 0, 0, 1]*IPL[{c, G}])/(12*Pi^2) - (gs^2*If[MC == 0, 0, 1]*IPL[{c, G}])/(6*FR$Eps*Pi^2) - (gs^2*If[MC == 0, 1, 0]*IPL[{c, G}])/(12*Pi^2) + (gs^2*If[MC == 0, 1, 0]*IPL[{c, G}])/(12*FR$Eps*Pi^2) + (gs^2*IPL[{s, G}])/(12*FR$Eps*Pi^2) + (gs^2*If[MC == 0, 0, 1]*IPL[{c, G}]*Log[MC/FR$MU])/(2*Pi^2)) + 2*((gs^2*MC*IPL[{c, G}])/(12*Pi^2) - (5*gs^2*MC*If[MC == 0, 0, 1]*IPL[{c, G}])/(12*Pi^2) - (gs^2*MC*If[MC == 0, 1, 0]*IPL[{c, G}])/(12*Pi^2) + (gs^2*MC*If[MC == 0, 1, 0]*IPL[{c, G}])/(4*FR$Eps*Pi^2) + (gs^2*MC*If[MC == 0, 0, 1]*IPL[{c, G}]*Log[MC/FR$MU])/(2*Pi^2)))*ProjM[Index[Spin, Ext[1]], Index[Spin, Ext[2]]])/(Sqrt[2]*vev)},
{{{anti[t], 1}, {b, 2}, {GP, 3}}, -(Sqrt[2]*gs^2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{t, b, G}]*(MT*ProjM[Index[Spin, Ext[1]], Index[Spin, Ext[2]]] - MB*ProjP[Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(3*FR$Eps*Pi^2*vev) + (IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*((MT*((gs^2*IPL[{b, G}])/(12*Pi^2) - (5*gs^2*If[MB == 0, 0, 1]*IPL[{b, G}])/(12*Pi^2) - (gs^2*If[MB == 0, 0, 1]*IPL[{b, G}])/(6*FR$Eps*Pi^2) - (gs^2*If[MB == 0, 1, 0]*IPL[{b, G}])/(12*Pi^2) + (gs^2*If[MB == 0, 1, 0]*IPL[{b, G}])/(12*FR$Eps*Pi^2) + (gs^2*IPL[{t, G}])/(12*Pi^2) - (5*gs^2*If[MT == 0, 0, 1]*IPL[{t, G}])/(12*Pi^2) - (gs^2*If[MT == 0, 0, 1]*IPL[{t, G}])/(6*FR$Eps*Pi^2) - (gs^2*If[MT == 0, 1, 0]*IPL[{t, G}])/(12*Pi^2) + (gs^2*If[MT == 0, 1, 0]*IPL[{t, G}])/(12*FR$Eps*Pi^2) + (gs^2*If[MB == 0, 0, 1]*IPL[{b, G}]*Log[MB/FR$MU])/(2*Pi^2) + (gs^2*If[MT == 0, 0, 1]*IPL[{t, G}]*Log[MT/FR$MU])/(2*Pi^2)) + 2*((gs^2*MT*IPL[{t, G}])/(12*Pi^2) - (5*gs^2*MT*If[MT == 0, 0, 1]*IPL[{t, G}])/(12*Pi^2) - (gs^2*MT*If[MT == 0, 1, 0]*IPL[{t, G}])/(12*Pi^2) + (gs^2*MT*If[MT == 0, 1, 0]*IPL[{t, G}])/(4*FR$Eps*Pi^2) + (gs^2*MT*If[MT == 0, 0, 1]*IPL[{t, G}]*Log[MT/FR$MU])/(2*Pi^2)))*ProjM[Index[Spin, Ext[1]], Index[Spin, Ext[2]]] - (2*((gs^2*MB*IPL[{b, G}])/(12*Pi^2) - (5*gs^2*MB*If[MB == 0, 0, 1]*IPL[{b, G}])/(12*Pi^2) - (gs^2*MB*If[MB == 0, 1, 0]*IPL[{b, G}])/(12*Pi^2) + (gs^2*MB*If[MB == 0, 1, 0]*IPL[{b, G}])/(4*FR$Eps*Pi^2) + (gs^2*MB*If[MB == 0, 0, 1]*IPL[{b, G}]*Log[MB/FR$MU])/(2*Pi^2)) + MB*((gs^2*IPL[{b, G}])/(12*Pi^2) - (5*gs^2*If[MB == 0, 0, 1]*IPL[{b, G}])/(12*Pi^2) - (gs^2*If[MB == 0, 0, 1]*IPL[{b, G}])/(6*FR$Eps*Pi^2) - (gs^2*If[MB == 0, 1, 0]*IPL[{b, G}])/(12*Pi^2) + (gs^2*If[MB == 0, 1, 0]*IPL[{b, G}])/(12*FR$Eps*Pi^2) + (gs^2*IPL[{t, G}])/(12*Pi^2) - (5*gs^2*If[MT == 0, 0, 1]*IPL[{t, G}])/(12*Pi^2) - (gs^2*If[MT == 0, 0, 1]*IPL[{t, G}])/(6*FR$Eps*Pi^2) - (gs^2*If[MT == 0, 1, 0]*IPL[{t, G}])/(12*Pi^2) + (gs^2*If[MT == 0, 1, 0]*IPL[{t, G}])/(12*FR$Eps*Pi^2) + (gs^2*If[MB == 0, 0, 1]*IPL[{b, G}]*Log[MB/FR$MU])/(2*Pi^2) + (gs^2*If[MT == 0, 0, 1]*IPL[{t, G}]*Log[MT/FR$MU])/(2*Pi^2)))*ProjP[Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(Sqrt[2]*vev)},
{{{anti[u], 1}, {u, 2}, {A, 3}}, ((-I/18)*ee*gs^2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{u, G}]*(TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]] + TensDot[Ga[Index[Lorentz, Ext[3]]], ProjP][Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(FR$Eps*Pi^2) + ((2*I)/3)*ee*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*((gs^2*IPL[{u, G}]*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]])/(12*FR$Eps*Pi^2) + (gs^2*IPL[{u, G}]*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjP][Index[Spin, Ext[1]], Index[Spin, Ext[2]]])/(12*FR$Eps*Pi^2))},
{{{anti[c], 1}, {c, 2}, {A, 3}}, ((-I/18)*ee*gs^2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{c, G}]*(TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]] + TensDot[Ga[Index[Lorentz, Ext[3]]], ProjP][Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(FR$Eps*Pi^2) + ((2*I)/3)*ee*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*(((gs^2*IPL[{c, G}])/(12*Pi^2) - (5*gs^2*If[MC == 0, 0, 1]*IPL[{c, G}])/(12*Pi^2) - (gs^2*If[MC == 0, 0, 1]*IPL[{c, G}])/(6*FR$Eps*Pi^2) - (gs^2*If[MC == 0, 1, 0]*IPL[{c, G}])/(12*Pi^2) + (gs^2*If[MC == 0, 1, 0]*IPL[{c, G}])/(12*FR$Eps*Pi^2) + (gs^2*If[MC == 0, 0, 1]*IPL[{c, G}]*Log[MC/FR$MU])/(2*Pi^2))*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]] + ((gs^2*IPL[{c, G}])/(12*Pi^2) - (5*gs^2*If[MC == 0, 0, 1]*IPL[{c, G}])/(12*Pi^2) - (gs^2*If[MC == 0, 0, 1]*IPL[{c, G}])/(6*FR$Eps*Pi^2) - (gs^2*If[MC == 0, 1, 0]*IPL[{c, G}])/(12*Pi^2) + (gs^2*If[MC == 0, 1, 0]*IPL[{c, G}])/(12*FR$Eps*Pi^2) + (gs^2*If[MC == 0, 0, 1]*IPL[{c, G}]*Log[MC/FR$MU])/(2*Pi^2))*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjP][Index[Spin, Ext[1]], Index[Spin, Ext[2]]])},
{{{anti[t], 1}, {t, 2}, {A, 3}}, ((-I/18)*ee*gs^2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{t, G}]*(TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]] + TensDot[Ga[Index[Lorentz, Ext[3]]], ProjP][Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(FR$Eps*Pi^2) + ((2*I)/3)*ee*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*(((gs^2*IPL[{t, G}])/(12*Pi^2) - (5*gs^2*If[MT == 0, 0, 1]*IPL[{t, G}])/(12*Pi^2) - (gs^2*If[MT == 0, 0, 1]*IPL[{t, G}])/(6*FR$Eps*Pi^2) - (gs^2*If[MT == 0, 1, 0]*IPL[{t, G}])/(12*Pi^2) + (gs^2*If[MT == 0, 1, 0]*IPL[{t, G}])/(12*FR$Eps*Pi^2) + (gs^2*If[MT == 0, 0, 1]*IPL[{t, G}]*Log[MT/FR$MU])/(2*Pi^2))*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]] + ((gs^2*IPL[{t, G}])/(12*Pi^2) - (5*gs^2*If[MT == 0, 0, 1]*IPL[{t, G}])/(12*Pi^2) - (gs^2*If[MT == 0, 0, 1]*IPL[{t, G}])/(6*FR$Eps*Pi^2) - (gs^2*If[MT == 0, 1, 0]*IPL[{t, G}])/(12*Pi^2) + (gs^2*If[MT == 0, 1, 0]*IPL[{t, G}])/(12*FR$Eps*Pi^2) + (gs^2*If[MT == 0, 0, 1]*IPL[{t, G}]*Log[MT/FR$MU])/(2*Pi^2))*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjP][Index[Spin, Ext[1]], Index[Spin, Ext[2]]])},
{{{anti[d], 1}, {d, 2}, {A, 3}}, ((I/36)*ee*gs^2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{d, G}]*(TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]] + TensDot[Ga[Index[Lorentz, Ext[3]]], ProjP][Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(FR$Eps*Pi^2) - (I/3)*ee*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*((gs^2*IPL[{d, G}]*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]])/(12*FR$Eps*Pi^2) + (gs^2*IPL[{d, G}]*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjP][Index[Spin, Ext[1]], Index[Spin, Ext[2]]])/(12*FR$Eps*Pi^2))},
{{{anti[s], 1}, {s, 2}, {A, 3}}, ((I/36)*ee*gs^2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{s, G}]*(TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]] + TensDot[Ga[Index[Lorentz, Ext[3]]], ProjP][Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(FR$Eps*Pi^2) - (I/3)*ee*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*((gs^2*IPL[{s, G}]*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]])/(12*FR$Eps*Pi^2) + (gs^2*IPL[{s, G}]*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjP][Index[Spin, Ext[1]], Index[Spin, Ext[2]]])/(12*FR$Eps*Pi^2))},
{{{anti[b], 1}, {b, 2}, {A, 3}}, ((I/36)*ee*gs^2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{b, G}]*(TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]] + TensDot[Ga[Index[Lorentz, Ext[3]]], ProjP][Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(FR$Eps*Pi^2) - (I/3)*ee*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*(((gs^2*IPL[{b, G}])/(12*Pi^2) - (5*gs^2*If[MB == 0, 0, 1]*IPL[{b, G}])/(12*Pi^2) - (gs^2*If[MB == 0, 0, 1]*IPL[{b, G}])/(6*FR$Eps*Pi^2) - (gs^2*If[MB == 0, 1, 0]*IPL[{b, G}])/(12*Pi^2) + (gs^2*If[MB == 0, 1, 0]*IPL[{b, G}])/(12*FR$Eps*Pi^2) + (gs^2*If[MB == 0, 0, 1]*IPL[{b, G}]*Log[MB/FR$MU])/(2*Pi^2))*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]] + ((gs^2*IPL[{b, G}])/(12*Pi^2) - (5*gs^2*If[MB == 0, 0, 1]*IPL[{b, G}])/(12*Pi^2) - (gs^2*If[MB == 0, 0, 1]*IPL[{b, G}])/(6*FR$Eps*Pi^2) - (gs^2*If[MB == 0, 1, 0]*IPL[{b, G}])/(12*Pi^2) + (gs^2*If[MB == 0, 1, 0]*IPL[{b, G}])/(12*FR$Eps*Pi^2) + (gs^2*If[MB == 0, 0, 1]*IPL[{b, G}]*Log[MB/FR$MU])/(2*Pi^2))*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjP][Index[Spin, Ext[1]], Index[Spin, Ext[2]]])},
{{{anti[u], 1}, {u, 2}, {Z, 3}}, ((-I/72)*ee*gs^2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{u, G}]*((3*cw^2 - sw^2)*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]] - 4*sw^2*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjP][Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(cw*FR$Eps*Pi^2*sw) + ((I/6)*ee*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*((gs^2*(3*cw^2 - sw^2)*IPL[{u, G}]*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]])/(12*FR$Eps*Pi^2) - (gs^2*sw^2*IPL[{u, G}]*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjP][Index[Spin, Ext[1]], Index[Spin, Ext[2]]])/(3*FR$Eps*Pi^2)))/(cw*sw)},
{{{anti[c], 1}, {c, 2}, {Z, 3}}, ((-I/72)*ee*gs^2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{c, G}]*((3*cw^2 - sw^2)*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]] - 4*sw^2*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjP][Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(cw*FR$Eps*Pi^2*sw) + ((I/6)*ee*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*((3*cw^2 - sw^2)*((gs^2*IPL[{c, G}])/(12*Pi^2) - (5*gs^2*If[MC == 0, 0, 1]*IPL[{c, G}])/(12*Pi^2) - (gs^2*If[MC == 0, 0, 1]*IPL[{c, G}])/(6*FR$Eps*Pi^2) - (gs^2*If[MC == 0, 1, 0]*IPL[{c, G}])/(12*Pi^2) + (gs^2*If[MC == 0, 1, 0]*IPL[{c, G}])/(12*FR$Eps*Pi^2) + (gs^2*If[MC == 0, 0, 1]*IPL[{c, G}]*Log[MC/FR$MU])/(2*Pi^2))*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]] - 4*sw^2*((gs^2*IPL[{c, G}])/(12*Pi^2) - (5*gs^2*If[MC == 0, 0, 1]*IPL[{c, G}])/(12*Pi^2) - (gs^2*If[MC == 0, 0, 1]*IPL[{c, G}])/(6*FR$Eps*Pi^2) - (gs^2*If[MC == 0, 1, 0]*IPL[{c, G}])/(12*Pi^2) + (gs^2*If[MC == 0, 1, 0]*IPL[{c, G}])/(12*FR$Eps*Pi^2) + (gs^2*If[MC == 0, 0, 1]*IPL[{c, G}]*Log[MC/FR$MU])/(2*Pi^2))*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjP][Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(cw*sw)},
{{{anti[t], 1}, {t, 2}, {Z, 3}}, ((-I/72)*ee*gs^2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{t, G}]*((3*cw^2 - sw^2)*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]] - 4*sw^2*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjP][Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(cw*FR$Eps*Pi^2*sw) + ((I/6)*ee*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*((3*cw^2 - sw^2)*((gs^2*IPL[{t, G}])/(12*Pi^2) - (5*gs^2*If[MT == 0, 0, 1]*IPL[{t, G}])/(12*Pi^2) - (gs^2*If[MT == 0, 0, 1]*IPL[{t, G}])/(6*FR$Eps*Pi^2) - (gs^2*If[MT == 0, 1, 0]*IPL[{t, G}])/(12*Pi^2) + (gs^2*If[MT == 0, 1, 0]*IPL[{t, G}])/(12*FR$Eps*Pi^2) + (gs^2*If[MT == 0, 0, 1]*IPL[{t, G}]*Log[MT/FR$MU])/(2*Pi^2))*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]] - 4*sw^2*((gs^2*IPL[{t, G}])/(12*Pi^2) - (5*gs^2*If[MT == 0, 0, 1]*IPL[{t, G}])/(12*Pi^2) - (gs^2*If[MT == 0, 0, 1]*IPL[{t, G}])/(6*FR$Eps*Pi^2) - (gs^2*If[MT == 0, 1, 0]*IPL[{t, G}])/(12*Pi^2) + (gs^2*If[MT == 0, 1, 0]*IPL[{t, G}])/(12*FR$Eps*Pi^2) + (gs^2*If[MT == 0, 0, 1]*IPL[{t, G}]*Log[MT/FR$MU])/(2*Pi^2))*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjP][Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(cw*sw)},
{{{anti[d], 1}, {d, 2}, {Z, 3}}, ((I/72)*ee*gs^2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{d, G}]*((3*cw^2 + sw^2)*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]] - 2*sw^2*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjP][Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(cw*FR$Eps*Pi^2*sw) - ((I/6)*ee*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*((gs^2*(3*cw^2 + sw^2)*IPL[{d, G}]*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]])/(12*FR$Eps*Pi^2) - (gs^2*sw^2*IPL[{d, G}]*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjP][Index[Spin, Ext[1]], Index[Spin, Ext[2]]])/(6*FR$Eps*Pi^2)))/(cw*sw)},
{{{anti[s], 1}, {s, 2}, {Z, 3}}, ((I/72)*ee*gs^2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{s, G}]*((3*cw^2 + sw^2)*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]] - 2*sw^2*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjP][Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(cw*FR$Eps*Pi^2*sw) - ((I/6)*ee*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*((gs^2*(3*cw^2 + sw^2)*IPL[{s, G}]*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]])/(12*FR$Eps*Pi^2) - (gs^2*sw^2*IPL[{s, G}]*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjP][Index[Spin, Ext[1]], Index[Spin, Ext[2]]])/(6*FR$Eps*Pi^2)))/(cw*sw)},
{{{anti[b], 1}, {b, 2}, {Z, 3}}, ((I/72)*ee*gs^2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{b, G}]*((3*cw^2 + sw^2)*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]] - 2*sw^2*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjP][Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(cw*FR$Eps*Pi^2*sw) - ((I/6)*ee*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*((3*cw^2 + sw^2)*((gs^2*IPL[{b, G}])/(12*Pi^2) - (5*gs^2*If[MB == 0, 0, 1]*IPL[{b, G}])/(12*Pi^2) - (gs^2*If[MB == 0, 0, 1]*IPL[{b, G}])/(6*FR$Eps*Pi^2) - (gs^2*If[MB == 0, 1, 0]*IPL[{b, G}])/(12*Pi^2) + (gs^2*If[MB == 0, 1, 0]*IPL[{b, G}])/(12*FR$Eps*Pi^2) + (gs^2*If[MB == 0, 0, 1]*IPL[{b, G}]*Log[MB/FR$MU])/(2*Pi^2))*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]] - 2*sw^2*((gs^2*IPL[{b, G}])/(12*Pi^2) - (5*gs^2*If[MB == 0, 0, 1]*IPL[{b, G}])/(12*Pi^2) - (gs^2*If[MB == 0, 0, 1]*IPL[{b, G}])/(6*FR$Eps*Pi^2) - (gs^2*If[MB == 0, 1, 0]*IPL[{b, G}])/(12*Pi^2) + (gs^2*If[MB == 0, 1, 0]*IPL[{b, G}])/(12*FR$Eps*Pi^2) + (gs^2*If[MB == 0, 0, 1]*IPL[{b, G}]*Log[MB/FR$MU])/(2*Pi^2))*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjP][Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(cw*sw)},
{{{anti[d], 1}, {u, 2}, {anti[W], 3}}, ((I/2)*ee*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*((gs^2*IPL[{d, G}])/(12*FR$Eps*Pi^2) + (gs^2*IPL[{u, G}])/(12*FR$Eps*Pi^2))*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]])/(Sqrt[2]*sw) - ((I/12)*ee*gs^2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{u, d, G}]*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]])/(Sqrt[2]*FR$Eps*Pi^2*sw)},
{{{anti[s], 1}, {c, 2}, {anti[W], 3}}, ((-I/12)*ee*gs^2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{c, s, G}]*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]])/(Sqrt[2]*FR$Eps*Pi^2*sw) + ((I/2)*ee*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*((gs^2*IPL[{c, G}])/(12*Pi^2) - (5*gs^2*If[MC == 0, 0, 1]*IPL[{c, G}])/(12*Pi^2) - (gs^2*If[MC == 0, 0, 1]*IPL[{c, G}])/(6*FR$Eps*Pi^2) - (gs^2*If[MC == 0, 1, 0]*IPL[{c, G}])/(12*Pi^2) + (gs^2*If[MC == 0, 1, 0]*IPL[{c, G}])/(12*FR$Eps*Pi^2) + (gs^2*IPL[{s, G}])/(12*FR$Eps*Pi^2) + (gs^2*If[MC == 0, 0, 1]*IPL[{c, G}]*Log[MC/FR$MU])/(2*Pi^2))*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]])/(Sqrt[2]*sw)},
{{{anti[b], 1}, {t, 2}, {anti[W], 3}}, ((-I/12)*ee*gs^2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{t, b, G}]*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]])/(Sqrt[2]*FR$Eps*Pi^2*sw) + ((I/2)*ee*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*((gs^2*IPL[{b, G}])/(12*Pi^2) - (5*gs^2*If[MB == 0, 0, 1]*IPL[{b, G}])/(12*Pi^2) - (gs^2*If[MB == 0, 0, 1]*IPL[{b, G}])/(6*FR$Eps*Pi^2) - (gs^2*If[MB == 0, 1, 0]*IPL[{b, G}])/(12*Pi^2) + (gs^2*If[MB == 0, 1, 0]*IPL[{b, G}])/(12*FR$Eps*Pi^2) + (gs^2*IPL[{t, G}])/(12*Pi^2) - (5*gs^2*If[MT == 0, 0, 1]*IPL[{t, G}])/(12*Pi^2) - (gs^2*If[MT == 0, 0, 1]*IPL[{t, G}])/(6*FR$Eps*Pi^2) - (gs^2*If[MT == 0, 1, 0]*IPL[{t, G}])/(12*Pi^2) + (gs^2*If[MT == 0, 1, 0]*IPL[{t, G}])/(12*FR$Eps*Pi^2) + (gs^2*If[MB == 0, 0, 1]*IPL[{b, G}]*Log[MB/FR$MU])/(2*Pi^2) + (gs^2*If[MT == 0, 0, 1]*IPL[{t, G}]*Log[MT/FR$MU])/(2*Pi^2))*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]])/(Sqrt[2]*sw)},
{{{anti[u], 1}, {d, 2}, {W, 3}}, ((I/2)*ee*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*((gs^2*IPL[{d, G}])/(12*FR$Eps*Pi^2) + (gs^2*IPL[{u, G}])/(12*FR$Eps*Pi^2))*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]])/(Sqrt[2]*sw) - ((I/12)*ee*gs^2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{u, d, G}]*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]])/(Sqrt[2]*FR$Eps*Pi^2*sw)},
{{{anti[c], 1}, {s, 2}, {W, 3}}, ((-I/12)*ee*gs^2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{c, s, G}]*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]])/(Sqrt[2]*FR$Eps*Pi^2*sw) + ((I/2)*ee*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*((gs^2*IPL[{c, G}])/(12*Pi^2) - (5*gs^2*If[MC == 0, 0, 1]*IPL[{c, G}])/(12*Pi^2) - (gs^2*If[MC == 0, 0, 1]*IPL[{c, G}])/(6*FR$Eps*Pi^2) - (gs^2*If[MC == 0, 1, 0]*IPL[{c, G}])/(12*Pi^2) + (gs^2*If[MC == 0, 1, 0]*IPL[{c, G}])/(12*FR$Eps*Pi^2) + (gs^2*IPL[{s, G}])/(12*FR$Eps*Pi^2) + (gs^2*If[MC == 0, 0, 1]*IPL[{c, G}]*Log[MC/FR$MU])/(2*Pi^2))*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]])/(Sqrt[2]*sw)},
{{{anti[t], 1}, {b, 2}, {W, 3}}, ((-I/12)*ee*gs^2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{t, b, G}]*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]])/(Sqrt[2]*FR$Eps*Pi^2*sw) + ((I/2)*ee*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*((gs^2*IPL[{b, G}])/(12*Pi^2) - (5*gs^2*If[MB == 0, 0, 1]*IPL[{b, G}])/(12*Pi^2) - (gs^2*If[MB == 0, 0, 1]*IPL[{b, G}])/(6*FR$Eps*Pi^2) - (gs^2*If[MB == 0, 1, 0]*IPL[{b, G}])/(12*Pi^2) + (gs^2*If[MB == 0, 1, 0]*IPL[{b, G}])/(12*FR$Eps*Pi^2) + (gs^2*IPL[{t, G}])/(12*Pi^2) - (5*gs^2*If[MT == 0, 0, 1]*IPL[{t, G}])/(12*Pi^2) - (gs^2*If[MT == 0, 0, 1]*IPL[{t, G}])/(6*FR$Eps*Pi^2) - (gs^2*If[MT == 0, 1, 0]*IPL[{t, G}])/(12*Pi^2) + (gs^2*If[MT == 0, 1, 0]*IPL[{t, G}])/(12*FR$Eps*Pi^2) + (gs^2*If[MB == 0, 0, 1]*IPL[{b, G}]*Log[MB/FR$MU])/(2*Pi^2) + (gs^2*If[MT == 0, 0, 1]*IPL[{t, G}]*Log[MT/FR$MU])/(2*Pi^2))*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]])/(Sqrt[2]*sw)},
{{{anti[u], 1}, {u, 2}, {G, 3}}, (((-13*I)/48)*gs^3*IPL[{u, G}]*SUNT[Index[Gluon, Ext[3]], Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*(TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]] + TensDot[Ga[Index[Lorentz, Ext[3]]], ProjP][Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(FR$Eps*Pi^2) + ((I/2)*gs*SUNT[Index[Gluon, Ext[3]], Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*(((aS*(-8*gs^2*If[MB == 0, 0, 1]*IPL[{b}]*Log[MB/FR$MU] - 8*gs^2*If[MC == 0, 0, 1]*IPL[{c}]*Log[MC/FR$MU] - 8*gs^2*If[MT == 0, 0, 1]*IPL[{t}]*Log[MT/FR$MU]))/(96*Pi^2) + aS*((gs^2*If[MB == 0, 1, 0]*IPL[{b}])/(24*FR$Eps*Pi^2) + (gs^2*If[MC == 0, 1, 0]*IPL[{c}])/(24*FR$Eps*Pi^2) + (gs^2*IPL[{d}])/(24*FR$Eps*Pi^2) - (19*gs^2*IPL[{G}])/(64*FR$Eps*Pi^2) - (gs^2*IPL[{ghG}])/(64*FR$Eps*Pi^2) + (gs^2*IPL[{s}])/(24*FR$Eps*Pi^2) + (gs^2*If[MT == 0, 1, 0]*IPL[{t}])/(24*FR$Eps*Pi^2) + (gs^2*IPL[{u}])/(24*FR$Eps*Pi^2) + (gs^2*IPL[{u, G}])/(6*FR$Eps*Pi^2) + (gs^2*If[MB == 0, 0, 1]*IPL[{b}]*Log[MB/FR$MU])/(12*Pi^2) + (gs^2*If[MC == 0, 0, 1]*IPL[{c}]*Log[MC/FR$MU])/(12*Pi^2) + (gs^2*If[MT == 0, 0, 1]*IPL[{t}]*Log[MT/FR$MU])/(12*Pi^2)))*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]] + ((aS*(-8*gs^2*If[MB == 0, 0, 1]*IPL[{b}]*Log[MB/FR$MU] - 8*gs^2*If[MC == 0, 0, 1]*IPL[{c}]*Log[MC/FR$MU] - 8*gs^2*If[MT == 0, 0, 1]*IPL[{t}]*Log[MT/FR$MU]))/(96*Pi^2) + aS*((gs^2*If[MB == 0, 1, 0]*IPL[{b}])/(24*FR$Eps*Pi^2) + (gs^2*If[MC == 0, 1, 0]*IPL[{c}])/(24*FR$Eps*Pi^2) + (gs^2*IPL[{d}])/(24*FR$Eps*Pi^2) - (19*gs^2*IPL[{G}])/(64*FR$Eps*Pi^2) - (gs^2*IPL[{ghG}])/(64*FR$Eps*Pi^2) + (gs^2*IPL[{s}])/(24*FR$Eps*Pi^2) + (gs^2*If[MT == 0, 1, 0]*IPL[{t}])/(24*FR$Eps*Pi^2) + (gs^2*IPL[{u}])/(24*FR$Eps*Pi^2) + (gs^2*IPL[{u, G}])/(6*FR$Eps*Pi^2) + (gs^2*If[MB == 0, 0, 1]*IPL[{b}]*Log[MB/FR$MU])/(12*Pi^2) + (gs^2*If[MC == 0, 0, 1]*IPL[{c}]*Log[MC/FR$MU])/(12*Pi^2) + (gs^2*If[MT == 0, 0, 1]*IPL[{t}]*Log[MT/FR$MU])/(12*Pi^2)))*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjP][Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/aS},
{{{anti[c], 1}, {c, 2}, {G, 3}}, (((-13*I)/48)*gs^3*IPL[{c, G}]*SUNT[Index[Gluon, Ext[3]], Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*(TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]] + TensDot[Ga[Index[Lorentz, Ext[3]]], ProjP][Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(FR$Eps*Pi^2) + ((I/2)*gs*SUNT[Index[Gluon, Ext[3]], Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*(((aS*(-8*gs^2*If[MB == 0, 0, 1]*IPL[{b}]*Log[MB/FR$MU] - 8*gs^2*If[MC == 0, 0, 1]*IPL[{c}]*Log[MC/FR$MU] - 8*gs^2*If[MT == 0, 0, 1]*IPL[{t}]*Log[MT/FR$MU]))/(96*Pi^2) + aS*((gs^2*If[MB == 0, 1, 0]*IPL[{b}])/(24*FR$Eps*Pi^2) + (gs^2*If[MC == 0, 1, 0]*IPL[{c}])/(24*FR$Eps*Pi^2) + (gs^2*IPL[{d}])/(24*FR$Eps*Pi^2) - (19*gs^2*IPL[{G}])/(64*FR$Eps*Pi^2) - (gs^2*IPL[{ghG}])/(64*FR$Eps*Pi^2) + (gs^2*IPL[{s}])/(24*FR$Eps*Pi^2) + (gs^2*If[MT == 0, 1, 0]*IPL[{t}])/(24*FR$Eps*Pi^2) + (gs^2*IPL[{u}])/(24*FR$Eps*Pi^2) + (gs^2*If[MB == 0, 0, 1]*IPL[{b}]*Log[MB/FR$MU])/(12*Pi^2) + (gs^2*If[MC == 0, 0, 1]*IPL[{c}]*Log[MC/FR$MU])/(12*Pi^2) + 2*((gs^2*IPL[{c, G}])/(12*Pi^2) - (5*gs^2*If[MC == 0, 0, 1]*IPL[{c, G}])/(12*Pi^2) - (gs^2*If[MC == 0, 0, 1]*IPL[{c, G}])/(6*FR$Eps*Pi^2) - (gs^2*If[MC == 0, 1, 0]*IPL[{c, G}])/(12*Pi^2) + (gs^2*If[MC == 0, 1, 0]*IPL[{c, G}])/(12*FR$Eps*Pi^2) + (gs^2*If[MC == 0, 0, 1]*IPL[{c, G}]*Log[MC/FR$MU])/(2*Pi^2)) + (gs^2*If[MT == 0, 0, 1]*IPL[{t}]*Log[MT/FR$MU])/(12*Pi^2)))*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]] + ((aS*(-8*gs^2*If[MB == 0, 0, 1]*IPL[{b}]*Log[MB/FR$MU] - 8*gs^2*If[MC == 0, 0, 1]*IPL[{c}]*Log[MC/FR$MU] - 8*gs^2*If[MT == 0, 0, 1]*IPL[{t}]*Log[MT/FR$MU]))/(96*Pi^2) + aS*((gs^2*If[MB == 0, 1, 0]*IPL[{b}])/(24*FR$Eps*Pi^2) + (gs^2*If[MC == 0, 1, 0]*IPL[{c}])/(24*FR$Eps*Pi^2) + (gs^2*IPL[{d}])/(24*FR$Eps*Pi^2) - (19*gs^2*IPL[{G}])/(64*FR$Eps*Pi^2) - (gs^2*IPL[{ghG}])/(64*FR$Eps*Pi^2) + (gs^2*IPL[{s}])/(24*FR$Eps*Pi^2) + (gs^2*If[MT == 0, 1, 0]*IPL[{t}])/(24*FR$Eps*Pi^2) + (gs^2*IPL[{u}])/(24*FR$Eps*Pi^2) + (gs^2*If[MB == 0, 0, 1]*IPL[{b}]*Log[MB/FR$MU])/(12*Pi^2) + (gs^2*If[MC == 0, 0, 1]*IPL[{c}]*Log[MC/FR$MU])/(12*Pi^2) + 2*((gs^2*IPL[{c, G}])/(12*Pi^2) - (5*gs^2*If[MC == 0, 0, 1]*IPL[{c, G}])/(12*Pi^2) - (gs^2*If[MC == 0, 0, 1]*IPL[{c, G}])/(6*FR$Eps*Pi^2) - (gs^2*If[MC == 0, 1, 0]*IPL[{c, G}])/(12*Pi^2) + (gs^2*If[MC == 0, 1, 0]*IPL[{c, G}])/(12*FR$Eps*Pi^2) + (gs^2*If[MC == 0, 0, 1]*IPL[{c, G}]*Log[MC/FR$MU])/(2*Pi^2)) + (gs^2*If[MT == 0, 0, 1]*IPL[{t}]*Log[MT/FR$MU])/(12*Pi^2)))*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjP][Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/aS},
{{{anti[t], 1}, {t, 2}, {G, 3}}, (((-13*I)/48)*gs^3*IPL[{t, G}]*SUNT[Index[Gluon, Ext[3]], Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*(TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]] + TensDot[Ga[Index[Lorentz, Ext[3]]], ProjP][Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(FR$Eps*Pi^2) + ((I/2)*gs*SUNT[Index[Gluon, Ext[3]], Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*(((aS*(-8*gs^2*If[MB == 0, 0, 1]*IPL[{b}]*Log[MB/FR$MU] - 8*gs^2*If[MC == 0, 0, 1]*IPL[{c}]*Log[MC/FR$MU] - 8*gs^2*If[MT == 0, 0, 1]*IPL[{t}]*Log[MT/FR$MU]))/(96*Pi^2) + aS*((gs^2*If[MB == 0, 1, 0]*IPL[{b}])/(24*FR$Eps*Pi^2) + (gs^2*If[MC == 0, 1, 0]*IPL[{c}])/(24*FR$Eps*Pi^2) + (gs^2*IPL[{d}])/(24*FR$Eps*Pi^2) - (19*gs^2*IPL[{G}])/(64*FR$Eps*Pi^2) - (gs^2*IPL[{ghG}])/(64*FR$Eps*Pi^2) + (gs^2*IPL[{s}])/(24*FR$Eps*Pi^2) + (gs^2*If[MT == 0, 1, 0]*IPL[{t}])/(24*FR$Eps*Pi^2) + (gs^2*IPL[{u}])/(24*FR$Eps*Pi^2) + (gs^2*If[MB == 0, 0, 1]*IPL[{b}]*Log[MB/FR$MU])/(12*Pi^2) + (gs^2*If[MC == 0, 0, 1]*IPL[{c}]*Log[MC/FR$MU])/(12*Pi^2) + (gs^2*If[MT == 0, 0, 1]*IPL[{t}]*Log[MT/FR$MU])/(12*Pi^2) + 2*((gs^2*IPL[{t, G}])/(12*Pi^2) - (5*gs^2*If[MT == 0, 0, 1]*IPL[{t, G}])/(12*Pi^2) - (gs^2*If[MT == 0, 0, 1]*IPL[{t, G}])/(6*FR$Eps*Pi^2) - (gs^2*If[MT == 0, 1, 0]*IPL[{t, G}])/(12*Pi^2) + (gs^2*If[MT == 0, 1, 0]*IPL[{t, G}])/(12*FR$Eps*Pi^2) + (gs^2*If[MT == 0, 0, 1]*IPL[{t, G}]*Log[MT/FR$MU])/(2*Pi^2))))*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]] + ((aS*(-8*gs^2*If[MB == 0, 0, 1]*IPL[{b}]*Log[MB/FR$MU] - 8*gs^2*If[MC == 0, 0, 1]*IPL[{c}]*Log[MC/FR$MU] - 8*gs^2*If[MT == 0, 0, 1]*IPL[{t}]*Log[MT/FR$MU]))/(96*Pi^2) + aS*((gs^2*If[MB == 0, 1, 0]*IPL[{b}])/(24*FR$Eps*Pi^2) + (gs^2*If[MC == 0, 1, 0]*IPL[{c}])/(24*FR$Eps*Pi^2) + (gs^2*IPL[{d}])/(24*FR$Eps*Pi^2) - (19*gs^2*IPL[{G}])/(64*FR$Eps*Pi^2) - (gs^2*IPL[{ghG}])/(64*FR$Eps*Pi^2) + (gs^2*IPL[{s}])/(24*FR$Eps*Pi^2) + (gs^2*If[MT == 0, 1, 0]*IPL[{t}])/(24*FR$Eps*Pi^2) + (gs^2*IPL[{u}])/(24*FR$Eps*Pi^2) + (gs^2*If[MB == 0, 0, 1]*IPL[{b}]*Log[MB/FR$MU])/(12*Pi^2) + (gs^2*If[MC == 0, 0, 1]*IPL[{c}]*Log[MC/FR$MU])/(12*Pi^2) + (gs^2*If[MT == 0, 0, 1]*IPL[{t}]*Log[MT/FR$MU])/(12*Pi^2) + 2*((gs^2*IPL[{t, G}])/(12*Pi^2) - (5*gs^2*If[MT == 0, 0, 1]*IPL[{t, G}])/(12*Pi^2) - (gs^2*If[MT == 0, 0, 1]*IPL[{t, G}])/(6*FR$Eps*Pi^2) - (gs^2*If[MT == 0, 1, 0]*IPL[{t, G}])/(12*Pi^2) + (gs^2*If[MT == 0, 1, 0]*IPL[{t, G}])/(12*FR$Eps*Pi^2) + (gs^2*If[MT == 0, 0, 1]*IPL[{t, G}]*Log[MT/FR$MU])/(2*Pi^2))))*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjP][Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/aS},
{{{anti[d], 1}, {d, 2}, {G, 3}}, (((-13*I)/48)*gs^3*IPL[{d, G}]*SUNT[Index[Gluon, Ext[3]], Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*(TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]] + TensDot[Ga[Index[Lorentz, Ext[3]]], ProjP][Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(FR$Eps*Pi^2) + ((I/2)*gs*SUNT[Index[Gluon, Ext[3]], Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*(((aS*(-8*gs^2*If[MB == 0, 0, 1]*IPL[{b}]*Log[MB/FR$MU] - 8*gs^2*If[MC == 0, 0, 1]*IPL[{c}]*Log[MC/FR$MU] - 8*gs^2*If[MT == 0, 0, 1]*IPL[{t}]*Log[MT/FR$MU]))/(96*Pi^2) + aS*((gs^2*If[MB == 0, 1, 0]*IPL[{b}])/(24*FR$Eps*Pi^2) + (gs^2*If[MC == 0, 1, 0]*IPL[{c}])/(24*FR$Eps*Pi^2) + (gs^2*IPL[{d}])/(24*FR$Eps*Pi^2) - (19*gs^2*IPL[{G}])/(64*FR$Eps*Pi^2) - (gs^2*IPL[{ghG}])/(64*FR$Eps*Pi^2) + (gs^2*IPL[{s}])/(24*FR$Eps*Pi^2) + (gs^2*If[MT == 0, 1, 0]*IPL[{t}])/(24*FR$Eps*Pi^2) + (gs^2*IPL[{u}])/(24*FR$Eps*Pi^2) + (gs^2*IPL[{d, G}])/(6*FR$Eps*Pi^2) + (gs^2*If[MB == 0, 0, 1]*IPL[{b}]*Log[MB/FR$MU])/(12*Pi^2) + (gs^2*If[MC == 0, 0, 1]*IPL[{c}]*Log[MC/FR$MU])/(12*Pi^2) + (gs^2*If[MT == 0, 0, 1]*IPL[{t}]*Log[MT/FR$MU])/(12*Pi^2)))*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]] + ((aS*(-8*gs^2*If[MB == 0, 0, 1]*IPL[{b}]*Log[MB/FR$MU] - 8*gs^2*If[MC == 0, 0, 1]*IPL[{c}]*Log[MC/FR$MU] - 8*gs^2*If[MT == 0, 0, 1]*IPL[{t}]*Log[MT/FR$MU]))/(96*Pi^2) + aS*((gs^2*If[MB == 0, 1, 0]*IPL[{b}])/(24*FR$Eps*Pi^2) + (gs^2*If[MC == 0, 1, 0]*IPL[{c}])/(24*FR$Eps*Pi^2) + (gs^2*IPL[{d}])/(24*FR$Eps*Pi^2) - (19*gs^2*IPL[{G}])/(64*FR$Eps*Pi^2) - (gs^2*IPL[{ghG}])/(64*FR$Eps*Pi^2) + (gs^2*IPL[{s}])/(24*FR$Eps*Pi^2) + (gs^2*If[MT == 0, 1, 0]*IPL[{t}])/(24*FR$Eps*Pi^2) + (gs^2*IPL[{u}])/(24*FR$Eps*Pi^2) + (gs^2*IPL[{d, G}])/(6*FR$Eps*Pi^2) + (gs^2*If[MB == 0, 0, 1]*IPL[{b}]*Log[MB/FR$MU])/(12*Pi^2) + (gs^2*If[MC == 0, 0, 1]*IPL[{c}]*Log[MC/FR$MU])/(12*Pi^2) + (gs^2*If[MT == 0, 0, 1]*IPL[{t}]*Log[MT/FR$MU])/(12*Pi^2)))*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjP][Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/aS},
{{{anti[s], 1}, {s, 2}, {G, 3}}, (((-13*I)/48)*gs^3*IPL[{s, G}]*SUNT[Index[Gluon, Ext[3]], Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*(TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]] + TensDot[Ga[Index[Lorentz, Ext[3]]], ProjP][Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(FR$Eps*Pi^2) + ((I/2)*gs*SUNT[Index[Gluon, Ext[3]], Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*(((aS*(-8*gs^2*If[MB == 0, 0, 1]*IPL[{b}]*Log[MB/FR$MU] - 8*gs^2*If[MC == 0, 0, 1]*IPL[{c}]*Log[MC/FR$MU] - 8*gs^2*If[MT == 0, 0, 1]*IPL[{t}]*Log[MT/FR$MU]))/(96*Pi^2) + aS*((gs^2*If[MB == 0, 1, 0]*IPL[{b}])/(24*FR$Eps*Pi^2) + (gs^2*If[MC == 0, 1, 0]*IPL[{c}])/(24*FR$Eps*Pi^2) + (gs^2*IPL[{d}])/(24*FR$Eps*Pi^2) - (19*gs^2*IPL[{G}])/(64*FR$Eps*Pi^2) - (gs^2*IPL[{ghG}])/(64*FR$Eps*Pi^2) + (gs^2*IPL[{s}])/(24*FR$Eps*Pi^2) + (gs^2*If[MT == 0, 1, 0]*IPL[{t}])/(24*FR$Eps*Pi^2) + (gs^2*IPL[{u}])/(24*FR$Eps*Pi^2) + (gs^2*IPL[{s, G}])/(6*FR$Eps*Pi^2) + (gs^2*If[MB == 0, 0, 1]*IPL[{b}]*Log[MB/FR$MU])/(12*Pi^2) + (gs^2*If[MC == 0, 0, 1]*IPL[{c}]*Log[MC/FR$MU])/(12*Pi^2) + (gs^2*If[MT == 0, 0, 1]*IPL[{t}]*Log[MT/FR$MU])/(12*Pi^2)))*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]] + ((aS*(-8*gs^2*If[MB == 0, 0, 1]*IPL[{b}]*Log[MB/FR$MU] - 8*gs^2*If[MC == 0, 0, 1]*IPL[{c}]*Log[MC/FR$MU] - 8*gs^2*If[MT == 0, 0, 1]*IPL[{t}]*Log[MT/FR$MU]))/(96*Pi^2) + aS*((gs^2*If[MB == 0, 1, 0]*IPL[{b}])/(24*FR$Eps*Pi^2) + (gs^2*If[MC == 0, 1, 0]*IPL[{c}])/(24*FR$Eps*Pi^2) + (gs^2*IPL[{d}])/(24*FR$Eps*Pi^2) - (19*gs^2*IPL[{G}])/(64*FR$Eps*Pi^2) - (gs^2*IPL[{ghG}])/(64*FR$Eps*Pi^2) + (gs^2*IPL[{s}])/(24*FR$Eps*Pi^2) + (gs^2*If[MT == 0, 1, 0]*IPL[{t}])/(24*FR$Eps*Pi^2) + (gs^2*IPL[{u}])/(24*FR$Eps*Pi^2) + (gs^2*IPL[{s, G}])/(6*FR$Eps*Pi^2) + (gs^2*If[MB == 0, 0, 1]*IPL[{b}]*Log[MB/FR$MU])/(12*Pi^2) + (gs^2*If[MC == 0, 0, 1]*IPL[{c}]*Log[MC/FR$MU])/(12*Pi^2) + (gs^2*If[MT == 0, 0, 1]*IPL[{t}]*Log[MT/FR$MU])/(12*Pi^2)))*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjP][Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/aS},
{{{anti[b], 1}, {b, 2}, {G, 3}}, (((-13*I)/48)*gs^3*IPL[{b, G}]*SUNT[Index[Gluon, Ext[3]], Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*(TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]] + TensDot[Ga[Index[Lorentz, Ext[3]]], ProjP][Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(FR$Eps*Pi^2) + ((I/2)*gs*SUNT[Index[Gluon, Ext[3]], Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*(((aS*(-8*gs^2*If[MB == 0, 0, 1]*IPL[{b}]*Log[MB/FR$MU] - 8*gs^2*If[MC == 0, 0, 1]*IPL[{c}]*Log[MC/FR$MU] - 8*gs^2*If[MT == 0, 0, 1]*IPL[{t}]*Log[MT/FR$MU]))/(96*Pi^2) + aS*((gs^2*If[MB == 0, 1, 0]*IPL[{b}])/(24*FR$Eps*Pi^2) + (gs^2*If[MC == 0, 1, 0]*IPL[{c}])/(24*FR$Eps*Pi^2) + (gs^2*IPL[{d}])/(24*FR$Eps*Pi^2) - (19*gs^2*IPL[{G}])/(64*FR$Eps*Pi^2) - (gs^2*IPL[{ghG}])/(64*FR$Eps*Pi^2) + (gs^2*IPL[{s}])/(24*FR$Eps*Pi^2) + (gs^2*If[MT == 0, 1, 0]*IPL[{t}])/(24*FR$Eps*Pi^2) + (gs^2*IPL[{u}])/(24*FR$Eps*Pi^2) + (gs^2*If[MB == 0, 0, 1]*IPL[{b}]*Log[MB/FR$MU])/(12*Pi^2) + 2*((gs^2*IPL[{b, G}])/(12*Pi^2) - (5*gs^2*If[MB == 0, 0, 1]*IPL[{b, G}])/(12*Pi^2) - (gs^2*If[MB == 0, 0, 1]*IPL[{b, G}])/(6*FR$Eps*Pi^2) - (gs^2*If[MB == 0, 1, 0]*IPL[{b, G}])/(12*Pi^2) + (gs^2*If[MB == 0, 1, 0]*IPL[{b, G}])/(12*FR$Eps*Pi^2) + (gs^2*If[MB == 0, 0, 1]*IPL[{b, G}]*Log[MB/FR$MU])/(2*Pi^2)) + (gs^2*If[MC == 0, 0, 1]*IPL[{c}]*Log[MC/FR$MU])/(12*Pi^2) + (gs^2*If[MT == 0, 0, 1]*IPL[{t}]*Log[MT/FR$MU])/(12*Pi^2)))*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]] + ((aS*(-8*gs^2*If[MB == 0, 0, 1]*IPL[{b}]*Log[MB/FR$MU] - 8*gs^2*If[MC == 0, 0, 1]*IPL[{c}]*Log[MC/FR$MU] - 8*gs^2*If[MT == 0, 0, 1]*IPL[{t}]*Log[MT/FR$MU]))/(96*Pi^2) + aS*((gs^2*If[MB == 0, 1, 0]*IPL[{b}])/(24*FR$Eps*Pi^2) + (gs^2*If[MC == 0, 1, 0]*IPL[{c}])/(24*FR$Eps*Pi^2) + (gs^2*IPL[{d}])/(24*FR$Eps*Pi^2) - (19*gs^2*IPL[{G}])/(64*FR$Eps*Pi^2) - (gs^2*IPL[{ghG}])/(64*FR$Eps*Pi^2) + (gs^2*IPL[{s}])/(24*FR$Eps*Pi^2) + (gs^2*If[MT == 0, 1, 0]*IPL[{t}])/(24*FR$Eps*Pi^2) + (gs^2*IPL[{u}])/(24*FR$Eps*Pi^2) + (gs^2*If[MB == 0, 0, 1]*IPL[{b}]*Log[MB/FR$MU])/(12*Pi^2) + 2*((gs^2*IPL[{b, G}])/(12*Pi^2) - (5*gs^2*If[MB == 0, 0, 1]*IPL[{b, G}])/(12*Pi^2) - (gs^2*If[MB == 0, 0, 1]*IPL[{b, G}])/(6*FR$Eps*Pi^2) - (gs^2*If[MB == 0, 1, 0]*IPL[{b, G}])/(12*Pi^2) + (gs^2*If[MB == 0, 1, 0]*IPL[{b, G}])/(12*FR$Eps*Pi^2) + (gs^2*If[MB == 0, 0, 1]*IPL[{b, G}]*Log[MB/FR$MU])/(2*Pi^2)) + (gs^2*If[MC == 0, 0, 1]*IPL[{c}]*Log[MC/FR$MU])/(12*Pi^2) + (gs^2*If[MT == 0, 0, 1]*IPL[{t}]*Log[MT/FR$MU])/(12*Pi^2)))*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjP][Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/aS},
{{{G, 1}, {G, 2}, {G, 3}}, -(gs*((aS*(-8*gs^2*If[MB == 0, 0, 1]*IPL[{b}]*Log[MB/FR$MU] - 8*gs^2*If[MC == 0, 0, 1]*IPL[{c}]*Log[MC/FR$MU] - 8*gs^2*If[MT == 0, 0, 1]*IPL[{t}]*Log[MT/FR$MU]))/(96*Pi^2) + 3*aS*((gs^2*If[MB == 0, 1, 0]*IPL[{b}])/(24*FR$Eps*Pi^2) + (gs^2*If[MC == 0, 1, 0]*IPL[{c}])/(24*FR$Eps*Pi^2) + (gs^2*IPL[{d}])/(24*FR$Eps*Pi^2) - (19*gs^2*IPL[{G}])/(64*FR$Eps*Pi^2) - (gs^2*IPL[{ghG}])/(64*FR$Eps*Pi^2) + (gs^2*IPL[{s}])/(24*FR$Eps*Pi^2) + (gs^2*If[MT == 0, 1, 0]*IPL[{t}])/(24*FR$Eps*Pi^2) + (gs^2*IPL[{u}])/(24*FR$Eps*Pi^2) + (gs^2*If[MB == 0, 0, 1]*IPL[{b}]*Log[MB/FR$MU])/(12*Pi^2) + (gs^2*If[MC == 0, 0, 1]*IPL[{c}]*Log[MC/FR$MU])/(12*Pi^2) + (gs^2*If[MT == 0, 0, 1]*IPL[{t}]*Log[MT/FR$MU])/(12*Pi^2)))*((-FV[1, Index[Lorentz, Ext[3]]] + FV[2, Index[Lorentz, Ext[3]]])*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[2]]] + (FV[1, Index[Lorentz, Ext[2]]] - FV[3, Index[Lorentz, Ext[2]]])*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[3]]] + (-FV[2, Index[Lorentz, Ext[1]]] + FV[3, Index[Lorentz, Ext[1]]])*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[3]]])*SUNF[Index[Gluon, Ext[1]], Index[Gluon, Ext[2]], Index[Gluon, Ext[3]]])/(2*aS) - (gs^3*(27*FV[3, Index[Lorentz, Ext[3]]]*IPL[{G}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[2]]] - 3*FV[3, Index[Lorentz, Ext[3]]]*IPL[{ghG}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[2]]] - 8*FV[2, Index[Lorentz, Ext[3]]]*(2*IPL[{b}] + 2*IPL[{c}] + 2*IPL[{d}] - 9*IPL[{G}] + 2*IPL[{s}] + 2*IPL[{t}] + 2*IPL[{u}])*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[2]]] + 2*FV[1, Index[Lorentz, Ext[3]]]*(8*IPL[{b}] + 8*IPL[{c}] + 8*IPL[{d}] - 9*IPL[{G}] - 3*IPL[{ghG}] + 8*IPL[{s}] + 8*IPL[{t}] + 8*IPL[{u}])*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[2]]] - 16*FV[1, Index[Lorentz, Ext[2]]]*IPL[{b}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[3]]] + 16*FV[3, Index[Lorentz, Ext[2]]]*IPL[{b}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[3]]] - 16*FV[1, Index[Lorentz, Ext[2]]]*IPL[{c}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[3]]] + 16*FV[3, Index[Lorentz, Ext[2]]]*IPL[{c}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[3]]] - 16*FV[1, Index[Lorentz, Ext[2]]]*IPL[{d}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[3]]] + 16*FV[3, Index[Lorentz, Ext[2]]]*IPL[{d}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[3]]] + 18*FV[1, Index[Lorentz, Ext[2]]]*IPL[{G}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[3]]] - 27*FV[2, Index[Lorentz, Ext[2]]]*IPL[{G}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[3]]] - 72*FV[3, Index[Lorentz, Ext[2]]]*IPL[{G}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[3]]] + 6*FV[1, Index[Lorentz, Ext[2]]]*IPL[{ghG}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[3]]] + 3*FV[2, Index[Lorentz, Ext[2]]]*IPL[{ghG}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[3]]] - 16*FV[1, Index[Lorentz, Ext[2]]]*IPL[{s}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[3]]] + 16*FV[3, Index[Lorentz, Ext[2]]]*IPL[{s}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[3]]] - 16*FV[1, Index[Lorentz, Ext[2]]]*IPL[{t}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[3]]] + 16*FV[3, Index[Lorentz, Ext[2]]]*IPL[{t}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[3]]] - 16*FV[1, Index[Lorentz, Ext[2]]]*IPL[{u}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[3]]] + 16*FV[3, Index[Lorentz, Ext[2]]]*IPL[{u}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[3]]] + 16*FV[2, Index[Lorentz, Ext[1]]]*IPL[{b}]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[3]]] - 16*FV[3, Index[Lorentz, Ext[1]]]*IPL[{b}]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[3]]] + 16*FV[2, Index[Lorentz, Ext[1]]]*IPL[{c}]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[3]]] - 16*FV[3, Index[Lorentz, Ext[1]]]*IPL[{c}]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[3]]] + 16*FV[2, Index[Lorentz, Ext[1]]]*IPL[{d}]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[3]]] - 16*FV[3, Index[Lorentz, Ext[1]]]*IPL[{d}]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[3]]] - 45*FV[2, Index[Lorentz, Ext[1]]]*IPL[{G}]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[3]]] + 45*FV[3, Index[Lorentz, Ext[1]]]*IPL[{G}]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[3]]] - 3*FV[2, Index[Lorentz, Ext[1]]]*IPL[{ghG}]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[3]]] + 3*FV[3, Index[Lorentz, Ext[1]]]*IPL[{ghG}]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[3]]] + 16*FV[2, Index[Lorentz, Ext[1]]]*IPL[{s}]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[3]]] - 16*FV[3, Index[Lorentz, Ext[1]]]*IPL[{s}]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[3]]] + 16*FV[2, Index[Lorentz, Ext[1]]]*IPL[{t}]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[3]]] - 16*FV[3, Index[Lorentz, Ext[1]]]*IPL[{t}]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[3]]] + 16*FV[2, Index[Lorentz, Ext[1]]]*IPL[{u}]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[3]]] - 16*FV[3, Index[Lorentz, Ext[1]]]*IPL[{u}]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[3]]])*SUNF[Index[Gluon, Ext[1]], Index[Gluon, Ext[2]], Index[Gluon, Ext[3]]])/(384*FR$Eps*Pi^2)},
{{{G, 1}, {G, 2}, {G, 3}, {G, 4}}, (I*gs^2*((aS*(-8*gs^2*If[MB == 0, 0, 1]*IPL[{b}]*Log[MB/FR$MU] - 8*gs^2*If[MC == 0, 0, 1]*IPL[{c}]*Log[MC/FR$MU] - 8*gs^2*If[MT == 0, 0, 1]*IPL[{t}]*Log[MT/FR$MU]))/(96*Pi^2) + 2*aS*((gs^2*If[MB == 0, 1, 0]*IPL[{b}])/(24*FR$Eps*Pi^2) + (gs^2*If[MC == 0, 1, 0]*IPL[{c}])/(24*FR$Eps*Pi^2) + (gs^2*IPL[{d}])/(24*FR$Eps*Pi^2) - (19*gs^2*IPL[{G}])/(64*FR$Eps*Pi^2) - (gs^2*IPL[{ghG}])/(64*FR$Eps*Pi^2) + (gs^2*IPL[{s}])/(24*FR$Eps*Pi^2) + (gs^2*If[MT == 0, 1, 0]*IPL[{t}])/(24*FR$Eps*Pi^2) + (gs^2*IPL[{u}])/(24*FR$Eps*Pi^2) + (gs^2*If[MB == 0, 0, 1]*IPL[{b}]*Log[MB/FR$MU])/(12*Pi^2) + (gs^2*If[MC == 0, 0, 1]*IPL[{c}]*Log[MC/FR$MU])/(12*Pi^2) + (gs^2*If[MT == 0, 0, 1]*IPL[{t}]*Log[MT/FR$MU])/(12*Pi^2)))*(ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[4]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[3]]]*(SUNF[Index[Gluon, Ext[1]], Index[Gluon, Ext[2]], Index[Gluon, Ext[3]], Index[Gluon, Ext[4]]] + SUNF[Index[Gluon, Ext[1]], Index[Gluon, Ext[3]], Index[Gluon, Ext[2]], Index[Gluon, Ext[4]]]) + ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[3]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[4]]]*(-SUNF[Index[Gluon, Ext[1]], Index[Gluon, Ext[2]], Index[Gluon, Ext[3]], Index[Gluon, Ext[4]]] + SUNF[Index[Gluon, Ext[1]], Index[Gluon, Ext[4]], Index[Gluon, Ext[2]], Index[Gluon, Ext[3]]]) - ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[2]]]*ME[Index[Lorentz, Ext[3]], Index[Lorentz, Ext[4]]]*(SUNF[Index[Gluon, Ext[1]], Index[Gluon, Ext[3]], Index[Gluon, Ext[2]], Index[Gluon, Ext[4]]] + SUNF[Index[Gluon, Ext[1]], Index[Gluon, Ext[4]], Index[Gluon, Ext[2]], Index[Gluon, Ext[3]]])))/aS + ((I/384)*gs^4*(-3*IndexDelta[Index[Gluon, Ext[1]], Index[Gluon, Ext[2]]]*IndexDelta[Index[Gluon, Ext[3]], Index[Gluon, Ext[4]]]*IPL[{G}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[4]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[3]]] + 3*IndexDelta[Index[Gluon, Ext[1]], Index[Gluon, Ext[2]]]*IndexDelta[Index[Gluon, Ext[3]], Index[Gluon, Ext[4]]]*IPL[{ghG}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[4]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[3]]] - 3*IndexDelta[Index[Gluon, Ext[1]], Index[Gluon, Ext[2]]]*IndexDelta[Index[Gluon, Ext[3]], Index[Gluon, Ext[4]]]*IPL[{G}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[3]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[4]]] + 3*IndexDelta[Index[Gluon, Ext[1]], Index[Gluon, Ext[2]]]*IndexDelta[Index[Gluon, Ext[3]], Index[Gluon, Ext[4]]]*IPL[{ghG}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[3]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[4]]] - 3*IndexDelta[Index[Gluon, Ext[1]], Index[Gluon, Ext[2]]]*IndexDelta[Index[Gluon, Ext[3]], Index[Gluon, Ext[4]]]*IPL[{G}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[2]]]*ME[Index[Lorentz, Ext[3]], Index[Lorentz, Ext[4]]] + 3*IndexDelta[Index[Gluon, Ext[1]], Index[Gluon, Ext[2]]]*IndexDelta[Index[Gluon, Ext[3]], Index[Gluon, Ext[4]]]*IPL[{ghG}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[2]]]*ME[Index[Lorentz, Ext[3]], Index[Lorentz, Ext[4]]] + 3*IndexDelta[Index[Gluon, Ext[1]], Index[Gluon, Ext[4]]]*IndexDelta[Index[Gluon, Ext[2]], Index[Gluon, Ext[3]]]*(-IPL[{G}] + IPL[{ghG}])*(ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[4]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[3]]] + ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[3]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[4]]] + ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[2]]]*ME[Index[Lorentz, Ext[3]], Index[Lorentz, Ext[4]]]) + 3*IndexDelta[Index[Gluon, Ext[1]], Index[Gluon, Ext[3]]]*IndexDelta[Index[Gluon, Ext[2]], Index[Gluon, Ext[4]]]*(-IPL[{G}] + IPL[{ghG}])*(ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[4]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[3]]] + ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[3]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[4]]] + ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[2]]]*ME[Index[Lorentz, Ext[3]], Index[Lorentz, Ext[4]]]) - 16*IPL[{b}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[4]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[3]]]*SUNF[Index[Gluon, Ext[1]], Index[Gluon, Ext[2]], Index[Gluon, Ext[3]], Index[Gluon, Ext[4]]] - 16*IPL[{c}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[4]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[3]]]*SUNF[Index[Gluon, Ext[1]], Index[Gluon, Ext[2]], Index[Gluon, Ext[3]], Index[Gluon, Ext[4]]] - 16*IPL[{d}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[4]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[3]]]*SUNF[Index[Gluon, Ext[1]], Index[Gluon, Ext[2]], Index[Gluon, Ext[3]], Index[Gluon, Ext[4]]] - 21*IPL[{G}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[4]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[3]]]*SUNF[Index[Gluon, Ext[1]], Index[Gluon, Ext[2]], Index[Gluon, Ext[3]], Index[Gluon, Ext[4]]] - 3*IPL[{ghG}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[4]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[3]]]*SUNF[Index[Gluon, Ext[1]], Index[Gluon, Ext[2]], Index[Gluon, Ext[3]], Index[Gluon, Ext[4]]] - 16*IPL[{s}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[4]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[3]]]*SUNF[Index[Gluon, Ext[1]], Index[Gluon, Ext[2]], Index[Gluon, Ext[3]], Index[Gluon, Ext[4]]] - 16*IPL[{t}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[4]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[3]]]*SUNF[Index[Gluon, Ext[1]], Index[Gluon, Ext[2]], Index[Gluon, Ext[3]], Index[Gluon, Ext[4]]] - 16*IPL[{u}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[4]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[3]]]*SUNF[Index[Gluon, Ext[1]], Index[Gluon, Ext[2]], Index[Gluon, Ext[3]], Index[Gluon, Ext[4]]] + 32*IPL[{b}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[3]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[4]]]*SUNF[Index[Gluon, Ext[1]], Index[Gluon, Ext[2]], Index[Gluon, Ext[3]], Index[Gluon, Ext[4]]] + 32*IPL[{c}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[3]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[4]]]*SUNF[Index[Gluon, Ext[1]], Index[Gluon, Ext[2]], Index[Gluon, Ext[3]], Index[Gluon, Ext[4]]] + 32*IPL[{d}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[3]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[4]]]*SUNF[Index[Gluon, Ext[1]], Index[Gluon, Ext[2]], Index[Gluon, Ext[3]], Index[Gluon, Ext[4]]] + 213*IPL[{G}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[3]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[4]]]*SUNF[Index[Gluon, Ext[1]], Index[Gluon, Ext[2]], Index[Gluon, Ext[3]], Index[Gluon, Ext[4]]] - 3*IPL[{ghG}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[3]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[4]]]*SUNF[Index[Gluon, Ext[1]], Index[Gluon, Ext[2]], Index[Gluon, Ext[3]], Index[Gluon, Ext[4]]] + 32*IPL[{s}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[3]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[4]]]*SUNF[Index[Gluon, Ext[1]], Index[Gluon, Ext[2]], Index[Gluon, Ext[3]], Index[Gluon, Ext[4]]] + 32*IPL[{t}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[3]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[4]]]*SUNF[Index[Gluon, Ext[1]], Index[Gluon, Ext[2]], Index[Gluon, Ext[3]], Index[Gluon, Ext[4]]] + 32*IPL[{u}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[3]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[4]]]*SUNF[Index[Gluon, Ext[1]], Index[Gluon, Ext[2]], Index[Gluon, Ext[3]], Index[Gluon, Ext[4]]] - 16*IPL[{b}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[2]]]*ME[Index[Lorentz, Ext[3]], Index[Lorentz, Ext[4]]]*SUNF[Index[Gluon, Ext[1]], Index[Gluon, Ext[2]], Index[Gluon, Ext[3]], Index[Gluon, Ext[4]]] - 16*IPL[{c}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[2]]]*ME[Index[Lorentz, Ext[3]], Index[Lorentz, Ext[4]]]*SUNF[Index[Gluon, Ext[1]], Index[Gluon, Ext[2]], Index[Gluon, Ext[3]], Index[Gluon, Ext[4]]] - 16*IPL[{d}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[2]]]*ME[Index[Lorentz, Ext[3]], Index[Lorentz, Ext[4]]]*SUNF[Index[Gluon, Ext[1]], Index[Gluon, Ext[2]], Index[Gluon, Ext[3]], Index[Gluon, Ext[4]]] + 141*IPL[{G}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[2]]]*ME[Index[Lorentz, Ext[3]], Index[Lorentz, Ext[4]]]*SUNF[Index[Gluon, Ext[1]], Index[Gluon, Ext[2]], Index[Gluon, Ext[3]], Index[Gluon, Ext[4]]] - 3*IPL[{ghG}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[2]]]*ME[Index[Lorentz, Ext[3]], Index[Lorentz, Ext[4]]]*SUNF[Index[Gluon, Ext[1]], Index[Gluon, Ext[2]], Index[Gluon, Ext[3]], Index[Gluon, Ext[4]]] - 16*IPL[{s}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[2]]]*ME[Index[Lorentz, Ext[3]], Index[Lorentz, Ext[4]]]*SUNF[Index[Gluon, Ext[1]], Index[Gluon, Ext[2]], Index[Gluon, Ext[3]], Index[Gluon, Ext[4]]] - 16*IPL[{t}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[2]]]*ME[Index[Lorentz, Ext[3]], Index[Lorentz, Ext[4]]]*SUNF[Index[Gluon, Ext[1]], Index[Gluon, Ext[2]], Index[Gluon, Ext[3]], Index[Gluon, Ext[4]]] - 16*IPL[{u}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[2]]]*ME[Index[Lorentz, Ext[3]], Index[Lorentz, Ext[4]]]*SUNF[Index[Gluon, Ext[1]], Index[Gluon, Ext[2]], Index[Gluon, Ext[3]], Index[Gluon, Ext[4]]] - 16*IPL[{b}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[4]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[3]]]*SUNF[Index[Gluon, Ext[1]], Index[Gluon, Ext[3]], Index[Gluon, Ext[2]], Index[Gluon, Ext[4]]] - 16*IPL[{c}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[4]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[3]]]*SUNF[Index[Gluon, Ext[1]], Index[Gluon, Ext[3]], Index[Gluon, Ext[2]], Index[Gluon, Ext[4]]] - 16*IPL[{d}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[4]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[3]]]*SUNF[Index[Gluon, Ext[1]], Index[Gluon, Ext[3]], Index[Gluon, Ext[2]], Index[Gluon, Ext[4]]] - 30*IPL[{G}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[4]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[3]]]*SUNF[Index[Gluon, Ext[1]], Index[Gluon, Ext[3]], Index[Gluon, Ext[2]], Index[Gluon, Ext[4]]] + 6*IPL[{ghG}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[4]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[3]]]*SUNF[Index[Gluon, Ext[1]], Index[Gluon, Ext[3]], Index[Gluon, Ext[2]], Index[Gluon, Ext[4]]] - 16*IPL[{s}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[4]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[3]]]*SUNF[Index[Gluon, Ext[1]], Index[Gluon, Ext[3]], Index[Gluon, Ext[2]], Index[Gluon, Ext[4]]] - 16*IPL[{t}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[4]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[3]]]*SUNF[Index[Gluon, Ext[1]], Index[Gluon, Ext[3]], Index[Gluon, Ext[2]], Index[Gluon, Ext[4]]] - 16*IPL[{u}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[4]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[3]]]*SUNF[Index[Gluon, Ext[1]], Index[Gluon, Ext[3]], Index[Gluon, Ext[2]], Index[Gluon, Ext[4]]] - 16*IPL[{b}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[3]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[4]]]*SUNF[Index[Gluon, Ext[1]], Index[Gluon, Ext[3]], Index[Gluon, Ext[2]], Index[Gluon, Ext[4]]] - 16*IPL[{c}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[3]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[4]]]*SUNF[Index[Gluon, Ext[1]], Index[Gluon, Ext[3]], Index[Gluon, Ext[2]], Index[Gluon, Ext[4]]] - 16*IPL[{d}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[3]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[4]]]*SUNF[Index[Gluon, Ext[1]], Index[Gluon, Ext[3]], Index[Gluon, Ext[2]], Index[Gluon, Ext[4]]] - 192*IPL[{G}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[3]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[4]]]*SUNF[Index[Gluon, Ext[1]], Index[Gluon, Ext[3]], Index[Gluon, Ext[2]], Index[Gluon, Ext[4]]] + 6*IPL[{ghG}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[3]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[4]]]*SUNF[Index[Gluon, Ext[1]], Index[Gluon, Ext[3]], Index[Gluon, Ext[2]], Index[Gluon, Ext[4]]] - 16*IPL[{s}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[3]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[4]]]*SUNF[Index[Gluon, Ext[1]], Index[Gluon, Ext[3]], Index[Gluon, Ext[2]], Index[Gluon, Ext[4]]] - 16*IPL[{t}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[3]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[4]]]*SUNF[Index[Gluon, Ext[1]], Index[Gluon, Ext[3]], Index[Gluon, Ext[2]], Index[Gluon, Ext[4]]] - 16*IPL[{u}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[3]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[4]]]*SUNF[Index[Gluon, Ext[1]], Index[Gluon, Ext[3]], Index[Gluon, Ext[2]], Index[Gluon, Ext[4]]] + 32*IPL[{b}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[2]]]*ME[Index[Lorentz, Ext[3]], Index[Lorentz, Ext[4]]]*SUNF[Index[Gluon, Ext[1]], Index[Gluon, Ext[3]], Index[Gluon, Ext[2]], Index[Gluon, Ext[4]]] + 32*IPL[{c}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[2]]]*ME[Index[Lorentz, Ext[3]], Index[Lorentz, Ext[4]]]*SUNF[Index[Gluon, Ext[1]], Index[Gluon, Ext[3]], Index[Gluon, Ext[2]], Index[Gluon, Ext[4]]] + 32*IPL[{d}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[2]]]*ME[Index[Lorentz, Ext[3]], Index[Lorentz, Ext[4]]]*SUNF[Index[Gluon, Ext[1]], Index[Gluon, Ext[3]], Index[Gluon, Ext[2]], Index[Gluon, Ext[4]]] - 120*IPL[{G}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[2]]]*ME[Index[Lorentz, Ext[3]], Index[Lorentz, Ext[4]]]*SUNF[Index[Gluon, Ext[1]], Index[Gluon, Ext[3]], Index[Gluon, Ext[2]], Index[Gluon, Ext[4]]] + 6*IPL[{ghG}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[2]]]*ME[Index[Lorentz, Ext[3]], Index[Lorentz, Ext[4]]]*SUNF[Index[Gluon, Ext[1]], Index[Gluon, Ext[3]], Index[Gluon, Ext[2]], Index[Gluon, Ext[4]]] + 32*IPL[{s}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[2]]]*ME[Index[Lorentz, Ext[3]], Index[Lorentz, Ext[4]]]*SUNF[Index[Gluon, Ext[1]], Index[Gluon, Ext[3]], Index[Gluon, Ext[2]], Index[Gluon, Ext[4]]] + 32*IPL[{t}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[2]]]*ME[Index[Lorentz, Ext[3]], Index[Lorentz, Ext[4]]]*SUNF[Index[Gluon, Ext[1]], Index[Gluon, Ext[3]], Index[Gluon, Ext[2]], Index[Gluon, Ext[4]]] + 32*IPL[{u}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[2]]]*ME[Index[Lorentz, Ext[3]], Index[Lorentz, Ext[4]]]*SUNF[Index[Gluon, Ext[1]], Index[Gluon, Ext[3]], Index[Gluon, Ext[2]], Index[Gluon, Ext[4]]] + 162*IPL[{G}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[3]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[4]]]*SUNF[Index[Gluon, Ext[1]], Index[Gluon, Ext[4]], Index[Gluon, Ext[2]], Index[Gluon, Ext[3]]] + 162*IPL[{G}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[2]]]*ME[Index[Lorentz, Ext[3]], Index[Lorentz, Ext[4]]]*SUNF[Index[Gluon, Ext[1]], Index[Gluon, Ext[4]], Index[Gluon, Ext[2]], Index[Gluon, Ext[3]]] - 18*IPL[{G}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[4]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[3]]]*SUNT[Index[Gluon, Ext[1]], Index[Gluon, Ext[3]], Index[Gluon, Ext[2]], Index[Gluon, Ext[4]]] + 18*IPL[{ghG}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[4]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[3]]]*SUNT[Index[Gluon, Ext[1]], Index[Gluon, Ext[3]], Index[Gluon, Ext[2]], Index[Gluon, Ext[4]]] - 18*IPL[{G}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[3]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[4]]]*SUNT[Index[Gluon, Ext[1]], Index[Gluon, Ext[3]], Index[Gluon, Ext[2]], Index[Gluon, Ext[4]]] + 18*IPL[{ghG}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[3]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[4]]]*SUNT[Index[Gluon, Ext[1]], Index[Gluon, Ext[3]], Index[Gluon, Ext[2]], Index[Gluon, Ext[4]]] - 18*IPL[{G}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[2]]]*ME[Index[Lorentz, Ext[3]], Index[Lorentz, Ext[4]]]*SUNT[Index[Gluon, Ext[1]], Index[Gluon, Ext[3]], Index[Gluon, Ext[2]], Index[Gluon, Ext[4]]] + 18*IPL[{ghG}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[2]]]*ME[Index[Lorentz, Ext[3]], Index[Lorentz, Ext[4]]]*SUNT[Index[Gluon, Ext[1]], Index[Gluon, Ext[3]], Index[Gluon, Ext[2]], Index[Gluon, Ext[4]]] - 18*IPL[{G}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[4]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[3]]]*SUNT[Index[Gluon, Ext[1]], Index[Gluon, Ext[4]], Index[Gluon, Ext[2]], Index[Gluon, Ext[3]]] + 18*IPL[{ghG}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[4]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[3]]]*SUNT[Index[Gluon, Ext[1]], Index[Gluon, Ext[4]], Index[Gluon, Ext[2]], Index[Gluon, Ext[3]]] - 18*IPL[{G}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[3]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[4]]]*SUNT[Index[Gluon, Ext[1]], Index[Gluon, Ext[4]], Index[Gluon, Ext[2]], Index[Gluon, Ext[3]]] + 18*IPL[{ghG}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[3]]]*ME[Index[Lorentz, Ext[2]], Index[Lorentz, Ext[4]]]*SUNT[Index[Gluon, Ext[1]], Index[Gluon, Ext[4]], Index[Gluon, Ext[2]], Index[Gluon, Ext[3]]] - 18*IPL[{G}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[2]]]*ME[Index[Lorentz, Ext[3]], Index[Lorentz, Ext[4]]]*SUNT[Index[Gluon, Ext[1]], Index[Gluon, Ext[4]], Index[Gluon, Ext[2]], Index[Gluon, Ext[3]]] + 18*IPL[{ghG}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[2]]]*ME[Index[Lorentz, Ext[3]], Index[Lorentz, Ext[4]]]*SUNT[Index[Gluon, Ext[1]], Index[Gluon, Ext[4]], Index[Gluon, Ext[2]], Index[Gluon, Ext[3]]]))/(FR$Eps*Pi^2)}};

(*CT parameters*)

FR$CTparam = {
};