2HDM: 2HDM.nlo

(****************************************************************************************************************)
(* Model automalically generated by NLOCT-1.02                                                                  *)
(* date and time of generation : Thu 15 Oct 2015 14:05:33                                                       *)
(* FeynRules model information :                                                                                *)
(*   ModelName : 2HDM                                                                                           *)
(*   Authors : C. Duhr, M. Herquet, C. Degrande                                                                 *)
(*   Version : 1.2                                                                                              *)
(*   Date : 12. 10. 2015                                                                                        *)
(*   Institutions : Universite catholique de Louvain (CP3), NIKHEF Theory Dept., IPPP                           *)
(*   URLs : http://feynrules.phys.ucl.ac.be/view/Main/2HDM                                                      *)
(*   Emails : claude.duhr@uclouvain.be, mherquet@nikhef.nl, celine.degrande@durham.ac.uk                        *)
(****************************************************************************************************************)

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

CT$Model = "2HDM/2HDM";
CT$GenericModel = "Lorentz";

(* Those assumptions where made : *)

NLOCT$assumptions = {};

(* Perturbed orders *)

FR$InteractionOrderPerturbativeExpansion = {{QCD, 1}, {QED, 0}, {YB, 0}, {YT, 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}]*(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]]]*((2*IPL[{c}] + 2*IPL[{d}] + 9*IPL[{G}] + 2*IPL[{s}] + 2*IPL[{u}])*SP[2, 2] + 2*IPL[{b}]*(-6*MB^2 + SP[2, 2]) + 2*IPL[{t}]*(-6*MT^2 + SP[2, 2]))))/Pi^2},
{{{anti[t], 1}, {t, 2}, {G0, 3}}, -(gs^2*yt*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*Sqrt[2]*Pi^2)},
{{{anti[b], 1}, {b, 2}, {G0, 3}}, (gs^2*yb*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*Sqrt[2]*Pi^2)},
{{{anti[b], 1}, {t, 2}, {anti[GP], 3}}, ((I/3)*gs^2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{t, b, G}]*(yb*ProjM[Index[Spin, Ext[1]], Index[Spin, Ext[2]]] - yt*ProjP[Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/Pi^2},
{{{anti[t], 1}, {b, 2}, {GP, 3}}, ((I/3)*gs^2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{t, b, G}]*(-(yt*ProjM[Index[Spin, Ext[1]], Index[Spin, Ext[2]]]) + yb*ProjP[Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/Pi^2},
{{{anti[d], 1}, {u, 2}, {anti[Hc], 3}}, ((I/3)*gs^2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{u, d, G}]*(Conjugate[GD1x1]*ProjM[Index[Spin, Ext[1]], Index[Spin, Ext[2]]] - GU1x1*ProjP[Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/Pi^2},
{{{anti[s], 1}, {c, 2}, {anti[Hc], 3}}, ((I/3)*gs^2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{c, s, G}]*(Conjugate[GD2x2]*ProjM[Index[Spin, Ext[1]], Index[Spin, Ext[2]]] - GU2x2*ProjP[Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/Pi^2},
{{{anti[b], 1}, {t, 2}, {anti[Hc], 3}}, ((I/3)*gs^2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{t, b, G}]*(Conjugate[GD3x3]*ProjM[Index[Spin, Ext[1]], Index[Spin, Ext[2]]] - GU3x3*ProjP[Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/Pi^2},
{{{anti[u], 1}, {d, 2}, {Hc, 3}}, ((I/3)*gs^2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{u, d, G}]*(-(Conjugate[GU1x1]*ProjM[Index[Spin, Ext[1]], Index[Spin, Ext[2]]]) + GD1x1*ProjP[Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/Pi^2},
{{{anti[c], 1}, {s, 2}, {Hc, 3}}, ((I/3)*gs^2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{c, s, G}]*(-(Conjugate[GU2x2]*ProjM[Index[Spin, Ext[1]], Index[Spin, Ext[2]]]) + GD2x2*ProjP[Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/Pi^2},
{{{anti[t], 1}, {b, 2}, {Hc, 3}}, ((I/3)*gs^2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{t, b, G}]*(-(Conjugate[GU3x3]*ProjM[Index[Spin, Ext[1]], Index[Spin, Ext[2]]]) + GD3x3*ProjP[Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/Pi^2},
{{{anti[u], 1}, {u, 2}, {h1, 3}}, ((I/3)*gs^2*TH2x1*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{u, G}]*(Conjugate[GU1x1]*ProjM[Index[Spin, Ext[1]], Index[Spin, Ext[2]]] + GU1x1*ProjP[Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(Sqrt[2]*Pi^2)},
{{{anti[c], 1}, {c, 2}, {h1, 3}}, ((I/3)*gs^2*TH2x1*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{c, G}]*(Conjugate[GU2x2]*ProjM[Index[Spin, Ext[1]], Index[Spin, Ext[2]]] + GU2x2*ProjP[Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(Sqrt[2]*Pi^2)},
{{{anti[t], 1}, {t, 2}, {h1, 3}}, ((I/3)*gs^2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{t, G}]*((TH1x1*yt + TH2x1*Conjugate[GU3x3])*ProjM[Index[Spin, Ext[1]], Index[Spin, Ext[2]]] + (GU3x3*TH2x1 + TH1x1*yt)*ProjP[Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(Sqrt[2]*Pi^2)},
{{{anti[d], 1}, {d, 2}, {h1, 3}}, ((I/3)*gs^2*TH2x1*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{d, G}]*(Conjugate[GD1x1]*ProjM[Index[Spin, Ext[1]], Index[Spin, Ext[2]]] + GD1x1*ProjP[Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(Sqrt[2]*Pi^2)},
{{{anti[s], 1}, {s, 2}, {h1, 3}}, ((I/3)*gs^2*TH2x1*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{s, G}]*(Conjugate[GD2x2]*ProjM[Index[Spin, Ext[1]], Index[Spin, Ext[2]]] + GD2x2*ProjP[Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(Sqrt[2]*Pi^2)},
{{{anti[b], 1}, {b, 2}, {h1, 3}}, ((I/3)*gs^2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{b, G}]*((TH1x1*yb + TH2x1*Conjugate[GD3x3])*ProjM[Index[Spin, Ext[1]], Index[Spin, Ext[2]]] + (GD3x3*TH2x1 + TH1x1*yb)*ProjP[Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(Sqrt[2]*Pi^2)},
{{{anti[u], 1}, {u, 2}, {h2, 3}}, ((I/3)*gs^2*TH2x2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{u, G}]*(Conjugate[GU1x1]*ProjM[Index[Spin, Ext[1]], Index[Spin, Ext[2]]] + GU1x1*ProjP[Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(Sqrt[2]*Pi^2)},
{{{anti[c], 1}, {c, 2}, {h2, 3}}, ((I/3)*gs^2*TH2x2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{c, G}]*(Conjugate[GU2x2]*ProjM[Index[Spin, Ext[1]], Index[Spin, Ext[2]]] + GU2x2*ProjP[Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(Sqrt[2]*Pi^2)},
{{{anti[t], 1}, {t, 2}, {h2, 3}}, ((I/3)*gs^2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{t, G}]*((TH1x2*yt + TH2x2*Conjugate[GU3x3])*ProjM[Index[Spin, Ext[1]], Index[Spin, Ext[2]]] + (GU3x3*TH2x2 + TH1x2*yt)*ProjP[Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(Sqrt[2]*Pi^2)},
{{{anti[d], 1}, {d, 2}, {h2, 3}}, ((I/3)*gs^2*TH2x2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{d, G}]*(Conjugate[GD1x1]*ProjM[Index[Spin, Ext[1]], Index[Spin, Ext[2]]] + GD1x1*ProjP[Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(Sqrt[2]*Pi^2)},
{{{anti[s], 1}, {s, 2}, {h2, 3}}, ((I/3)*gs^2*TH2x2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{s, G}]*(Conjugate[GD2x2]*ProjM[Index[Spin, Ext[1]], Index[Spin, Ext[2]]] + GD2x2*ProjP[Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(Sqrt[2]*Pi^2)},
{{{anti[b], 1}, {b, 2}, {h2, 3}}, ((I/3)*gs^2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{b, G}]*((TH1x2*yb + TH2x2*Conjugate[GD3x3])*ProjM[Index[Spin, Ext[1]], Index[Spin, Ext[2]]] + (GD3x3*TH2x2 + TH1x2*yb)*ProjP[Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(Sqrt[2]*Pi^2)},
{{{anti[u], 1}, {u, 2}, {h3, 3}}, -(gs^2*TH3x3*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{u, G}]*(Conjugate[GU1x1]*ProjM[Index[Spin, Ext[1]], Index[Spin, Ext[2]]] - GU1x1*ProjP[Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(3*Sqrt[2]*Pi^2)},
{{{anti[c], 1}, {c, 2}, {h3, 3}}, -(gs^2*TH3x3*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{c, G}]*(Conjugate[GU2x2]*ProjM[Index[Spin, Ext[1]], Index[Spin, Ext[2]]] - GU2x2*ProjP[Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(3*Sqrt[2]*Pi^2)},
{{{anti[t], 1}, {t, 2}, {h3, 3}}, -(gs^2*TH3x3*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{t, G}]*(Conjugate[GU3x3]*ProjM[Index[Spin, Ext[1]], Index[Spin, Ext[2]]] - GU3x3*ProjP[Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(3*Sqrt[2]*Pi^2)},
{{{anti[d], 1}, {d, 2}, {h3, 3}}, (gs^2*TH3x3*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{d, G}]*(Conjugate[GD1x1]*ProjM[Index[Spin, Ext[1]], Index[Spin, Ext[2]]] - GD1x1*ProjP[Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(3*Sqrt[2]*Pi^2)},
{{{anti[s], 1}, {s, 2}, {h3, 3}}, (gs^2*TH3x3*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{s, G}]*(Conjugate[GD2x2]*ProjM[Index[Spin, Ext[1]], Index[Spin, Ext[2]]] - GD2x2*ProjP[Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(3*Sqrt[2]*Pi^2)},
{{{anti[b], 1}, {b, 2}, {h3, 3}}, (gs^2*TH3x3*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{b, G}]*(Conjugate[GD3x3]*ProjM[Index[Spin, Ext[1]], Index[Spin, Ext[2]]] - GD3x3*ProjP[Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(3*Sqrt[2]*Pi^2)},
{{{anti[u], 1}, {u, 2}, {A, 3}}, ((-I/36)*gs^2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{u, G}]*((cw*g1 + 3*gw*sw)*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]] + 4*cw*g1*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/36)*gs^2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{c, G}]*((cw*g1 + 3*gw*sw)*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]] + 4*cw*g1*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/36)*gs^2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{t, G}]*((cw*g1 + 3*gw*sw)*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]] + 4*cw*g1*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/36)*gs^2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{d, G}]*((cw*g1 - 3*gw*sw)*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]] - 2*cw*g1*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/36)*gs^2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{s, G}]*((cw*g1 - 3*gw*sw)*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]] - 2*cw*g1*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/36)*gs^2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{b, G}]*((cw*g1 - 3*gw*sw)*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]] - 2*cw*g1*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)*gs^2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{u, G}]*((3*cw*gw - g1*sw)*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]] - 4*g1*sw*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjP][Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/Pi^2},
{{{anti[c], 1}, {c, 2}, {Z, 3}}, ((-I/36)*gs^2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{c, G}]*((3*cw*gw - g1*sw)*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]] - 4*g1*sw*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjP][Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/Pi^2},
{{{anti[t], 1}, {t, 2}, {Z, 3}}, ((-I/36)*gs^2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{t, G}]*((3*cw*gw - g1*sw)*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]] - 4*g1*sw*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjP][Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/Pi^2},
{{{anti[d], 1}, {d, 2}, {Z, 3}}, ((I/36)*gs^2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{d, G}]*((3*cw*gw + g1*sw)*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]] - 2*g1*sw*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjP][Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/Pi^2},
{{{anti[s], 1}, {s, 2}, {Z, 3}}, ((I/36)*gs^2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{s, G}]*((3*cw*gw + g1*sw)*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]] - 2*g1*sw*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjP][Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/Pi^2},
{{{anti[b], 1}, {b, 2}, {Z, 3}}, ((I/36)*gs^2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{b, G}]*((3*cw*gw + g1*sw)*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]] - 2*g1*sw*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjP][Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/Pi^2},
{{{anti[d], 1}, {u, 2}, {anti[W], 3}}, ((-I/6)*gs^2*gw*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)},
{{{anti[s], 1}, {c, 2}, {anti[W], 3}}, ((-I/6)*gs^2*gw*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)},
{{{anti[b], 1}, {t, 2}, {anti[W], 3}}, ((-I/6)*gs^2*gw*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)},
{{{anti[u], 1}, {d, 2}, {W, 3}}, ((-I/6)*gs^2*gw*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)},
{{{anti[c], 1}, {s, 2}, {W, 3}}, ((-I/6)*gs^2*gw*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)},
{{{anti[t], 1}, {b, 2}, {W, 3}}, ((-I/6)*gs^2*gw*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)},
{{{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}, {h1, 3}}, ((-I/16)*gs^2*IndexDelta[Index[Gluon, Ext[1]], Index[Gluon, Ext[2]]]*(MB*(GD3x3*TH2x1 + 2*TH1x1*yb + TH2x1*Conjugate[GD3x3])*IPL[{b}] + MT*(GU3x3*TH2x1 + 2*TH1x1*yt + TH2x1*Conjugate[GU3x3])*IPL[{t}])*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[2]]])/(Sqrt[2]*Pi^2)},
{{{G, 1}, {G, 2}, {h2, 3}}, ((-I/16)*gs^2*IndexDelta[Index[Gluon, Ext[1]], Index[Gluon, Ext[2]]]*(MB*(GD3x3*TH2x2 + 2*TH1x2*yb + TH2x2*Conjugate[GD3x3])*IPL[{b}] + MT*(GU3x3*TH2x2 + 2*TH1x2*yt + TH2x2*Conjugate[GU3x3])*IPL[{t}])*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[2]]])/(Sqrt[2]*Pi^2)},
{{{G, 1}, {G, 2}, {h3, 3}}, -(gs^2*TH3x3*IndexDelta[Index[Gluon, Ext[1]], Index[Gluon, Ext[2]]]*(MB*(-GD3x3 + Conjugate[GD3x3])*IPL[{b}] + MT*(GU3x3 - Conjugate[GU3x3])*IPL[{t}])*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[2]]])/(16*Sqrt[2]*Pi^2)},
{{{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)*gs^2*gw^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},
{{{Z, 1}, {Z, 2}, {G, 3}, {G, 4}}, ((I/1728)*gs^2*IndexDelta[Index[Gluon, Ext[3]], Index[Gluon, Ext[4]]]*(9*cw^2*gw^2*IPL[{b}] + 6*cw*g1*gw*sw*IPL[{b}] + 5*g1^2*sw^2*IPL[{b}] + (9*cw^2*gw^2 - 6*cw*g1*gw*sw + 17*g1^2*sw^2)*IPL[{c}] + 9*cw^2*gw^2*IPL[{d}] + 6*cw*g1*gw*sw*IPL[{d}] + 5*g1^2*sw^2*IPL[{d}] + 9*cw^2*gw^2*IPL[{s}] + 6*cw*g1*gw*sw*IPL[{s}] + 5*g1^2*sw^2*IPL[{s}] + 9*cw^2*gw^2*IPL[{t}] - 6*cw*g1*gw*sw*IPL[{t}] + 17*g1^2*sw^2*IPL[{t}] + (9*cw^2*gw^2 - 6*cw*g1*gw*sw + 17*g1^2*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]]]))/Pi^2},
{{{A, 1}, {Z, 2}, {G, 3}, {G, 4}}, ((I/1728)*gs^2*IndexDelta[Index[Gluon, Ext[3]], Index[Gluon, Ext[4]]]*(-3*cw^2*g1*gw*IPL[{b}] - 5*cw*g1^2*sw*IPL[{b}] + 9*cw*gw^2*sw*IPL[{b}] + 3*g1*gw*sw^2*IPL[{b}] + (3*cw^2*g1*gw - 17*cw*g1^2*sw + 9*cw*gw^2*sw - 3*g1*gw*sw^2)*IPL[{c}] - 3*cw^2*g1*gw*IPL[{d}] - 5*cw*g1^2*sw*IPL[{d}] + 9*cw*gw^2*sw*IPL[{d}] + 3*g1*gw*sw^2*IPL[{d}] - 3*cw^2*g1*gw*IPL[{s}] - 5*cw*g1^2*sw*IPL[{s}] + 9*cw*gw^2*sw*IPL[{s}] + 3*g1*gw*sw^2*IPL[{s}] + 3*cw^2*g1*gw*IPL[{t}] - 17*cw*g1^2*sw*IPL[{t}] + 9*cw*gw^2*sw*IPL[{t}] - 3*g1*gw*sw^2*IPL[{t}] + (3*cw^2*g1*gw - 17*cw*g1^2*sw + 9*cw*gw^2*sw - 3*g1*gw*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]]]))/Pi^2},
{{{A, 1}, {A, 2}, {G, 3}, {G, 4}}, ((I/1728)*gs^2*IndexDelta[Index[Gluon, Ext[3]], Index[Gluon, Ext[4]]]*(5*cw^2*g1^2*IPL[{b}] - 6*cw*g1*gw*sw*IPL[{b}] + 9*gw^2*sw^2*IPL[{b}] + (17*cw^2*g1^2 + 6*cw*g1*gw*sw + 9*gw^2*sw^2)*IPL[{c}] + 5*cw^2*g1^2*IPL[{d}] - 6*cw*g1*gw*sw*IPL[{d}] + 9*gw^2*sw^2*IPL[{d}] + 5*cw^2*g1^2*IPL[{s}] - 6*cw*g1*gw*sw*IPL[{s}] + 9*gw^2*sw^2*IPL[{s}] + 17*cw^2*g1^2*IPL[{t}] + 6*cw*g1*gw*sw*IPL[{t}] + 9*gw^2*sw^2*IPL[{t}] + (17*cw^2*g1^2 + 6*cw*g1*gw*sw + 9*gw^2*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]]]))/Pi^2},
{{{Z, 1}, {G, 2}, {G, 3}, {G, 4}}, ((I/576)*gs^3*(-27*(cw*gw + g1*sw)*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*gw*IPL[{b}] + g1*sw*IPL[{b}] + (3*cw*gw - 5*g1*sw)*IPL[{c}] - 3*cw*gw*IPL[{d}] + g1*sw*IPL[{d}] - 3*cw*gw*IPL[{s}] + g1*sw*IPL[{s}] + 3*cw*gw*IPL[{t}] - 5*g1*sw*IPL[{t}] + (3*cw*gw - 5*g1*sw)*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]]])))/Pi^2},
{{{A, 1}, {G, 2}, {G, 3}, {G, 4}}, ((I/576)*gs^3*(27*(cw*g1 - gw*sw)*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]]] + (-(cw*g1*IPL[{b}]) - 3*gw*sw*IPL[{b}] + (5*cw*g1 + 3*gw*sw)*IPL[{c}] - cw*g1*IPL[{d}] - 3*gw*sw*IPL[{d}] - cw*g1*IPL[{s}] - 3*gw*sw*IPL[{s}] + 5*cw*g1*IPL[{t}] + 3*gw*sw*IPL[{t}] + (5*cw*g1 + 3*gw*sw)*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}, {G0, 3}, {G0, 4}}, ((-I/16)*gs^2*IndexDelta[Index[Gluon, Ext[1]], Index[Gluon, Ext[2]]]*(yb^2*IPL[{b}] + yt^2*IPL[{t}])*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[2]]])/Pi^2},
{{{G, 1}, {G, 2}, {anti[GP], 3}, {GP, 4}}, ((-I/16)*gs^2*(yb^2 + yt^2)*IndexDelta[Index[Gluon, Ext[1]], Index[Gluon, Ext[2]]]*IPL[{t, b}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[2]]])/Pi^2},
{{{G, 1}, {G, 2}, {anti[Hc], 3}, {GP, 4}}, ((-I/16)*gs^2*(GU3x3*yt + yb*Conjugate[GD3x3])*IndexDelta[Index[Gluon, Ext[1]], Index[Gluon, Ext[2]]]*IPL[{t, b}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[2]]])/Pi^2},
{{{G, 1}, {G, 2}, {anti[GP], 3}, {Hc, 4}}, ((-I/16)*gs^2*(GD3x3*yb + yt*Conjugate[GU3x3])*IndexDelta[Index[Gluon, Ext[1]], Index[Gluon, Ext[2]]]*IPL[{t, b}]*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[2]]])/Pi^2},
{{{G, 1}, {G, 2}, {anti[Hc], 3}, {Hc, 4}}, ((-I/16)*gs^2*IndexDelta[Index[Gluon, Ext[1]], Index[Gluon, Ext[2]]]*(GD2x2*Conjugate[GD2x2]*IPL[{c, s}] + GU2x2*Conjugate[GU2x2]*IPL[{c, s}] + GD3x3*Conjugate[GD3x3]*IPL[{t, b}] + GU3x3*Conjugate[GU3x3]*IPL[{t, b}] + GD1x1*Conjugate[GD1x1]*IPL[{u, d}] + GU1x1*Conjugate[GU1x1]*IPL[{u, d}])*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[2]]])/Pi^2},
{{{G, 1}, {G, 2}, {G0, 3}, {h1, 4}}, (gs^2*TH2x1*IndexDelta[Index[Gluon, Ext[1]], Index[Gluon, Ext[2]]]*(yb*(-GD3x3 + Conjugate[GD3x3])*IPL[{b}] + yt*(GU3x3 - Conjugate[GU3x3])*IPL[{t}])*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[2]]])/(32*Pi^2)},
{{{G, 1}, {G, 2}, {h1, 3}, {h1, 4}}, ((-I/16)*gs^2*IndexDelta[Index[Gluon, Ext[1]], Index[Gluon, Ext[2]]]*(GD3x3*TH1x1*TH2x1*yb*IPL[{b}] + TH1x1^2*yb^2*IPL[{b}] + GD3x3*TH2x1^2*Conjugate[GD3x3]*IPL[{b}] + TH1x1*TH2x1*yb*Conjugate[GD3x3]*IPL[{b}] + GU2x2*TH2x1^2*Conjugate[GU2x2]*IPL[{c}] + GD1x1*TH2x1^2*Conjugate[GD1x1]*IPL[{d}] + GD2x2*TH2x1^2*Conjugate[GD2x2]*IPL[{s}] + GU3x3*TH1x1*TH2x1*yt*IPL[{t}] + TH1x1^2*yt^2*IPL[{t}] + GU3x3*TH2x1^2*Conjugate[GU3x3]*IPL[{t}] + TH1x1*TH2x1*yt*Conjugate[GU3x3]*IPL[{t}] + GU1x1*TH2x1^2*Conjugate[GU1x1]*IPL[{u}])*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[2]]])/Pi^2},
{{{G, 1}, {G, 2}, {G0, 3}, {h2, 4}}, (gs^2*TH2x2*IndexDelta[Index[Gluon, Ext[1]], Index[Gluon, Ext[2]]]*(yb*(-GD3x3 + Conjugate[GD3x3])*IPL[{b}] + yt*(GU3x3 - Conjugate[GU3x3])*IPL[{t}])*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[2]]])/(32*Pi^2)},
{{{G, 1}, {G, 2}, {h1, 3}, {h2, 4}}, ((-I/32)*gs^2*IndexDelta[Index[Gluon, Ext[1]], Index[Gluon, Ext[2]]]*(GD3x3*TH1x2*TH2x1*yb*IPL[{b}] + GD3x3*TH1x1*TH2x2*yb*IPL[{b}] + 2*TH1x1*TH1x2*yb^2*IPL[{b}] + 2*GD3x3*TH2x1*TH2x2*Conjugate[GD3x3]*IPL[{b}] + TH1x2*TH2x1*yb*Conjugate[GD3x3]*IPL[{b}] + TH1x1*TH2x2*yb*Conjugate[GD3x3]*IPL[{b}] + 2*GU2x2*TH2x1*TH2x2*Conjugate[GU2x2]*IPL[{c}] + 2*GD1x1*TH2x1*TH2x2*Conjugate[GD1x1]*IPL[{d}] + 2*GD2x2*TH2x1*TH2x2*Conjugate[GD2x2]*IPL[{s}] + GU3x3*TH1x2*TH2x1*yt*IPL[{t}] + GU3x3*TH1x1*TH2x2*yt*IPL[{t}] + 2*TH1x1*TH1x2*yt^2*IPL[{t}] + 2*GU3x3*TH2x1*TH2x2*Conjugate[GU3x3]*IPL[{t}] + TH1x2*TH2x1*yt*Conjugate[GU3x3]*IPL[{t}] + TH1x1*TH2x2*yt*Conjugate[GU3x3]*IPL[{t}] + 2*GU1x1*TH2x1*TH2x2*Conjugate[GU1x1]*IPL[{u}])*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[2]]])/Pi^2},
{{{G, 1}, {G, 2}, {h2, 3}, {h2, 4}}, ((-I/16)*gs^2*IndexDelta[Index[Gluon, Ext[1]], Index[Gluon, Ext[2]]]*(GD3x3*TH1x2*TH2x2*yb*IPL[{b}] + TH1x2^2*yb^2*IPL[{b}] + GD3x3*TH2x2^2*Conjugate[GD3x3]*IPL[{b}] + TH1x2*TH2x2*yb*Conjugate[GD3x3]*IPL[{b}] + GU2x2*TH2x2^2*Conjugate[GU2x2]*IPL[{c}] + GD1x1*TH2x2^2*Conjugate[GD1x1]*IPL[{d}] + GD2x2*TH2x2^2*Conjugate[GD2x2]*IPL[{s}] + GU3x3*TH1x2*TH2x2*yt*IPL[{t}] + TH1x2^2*yt^2*IPL[{t}] + GU3x3*TH2x2^2*Conjugate[GU3x3]*IPL[{t}] + TH1x2*TH2x2*yt*Conjugate[GU3x3]*IPL[{t}] + GU1x1*TH2x2^2*Conjugate[GU1x1]*IPL[{u}])*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[2]]])/Pi^2},
{{{G, 1}, {G, 2}, {G0, 3}, {h3, 4}}, ((-I/32)*gs^2*TH3x3*IndexDelta[Index[Gluon, Ext[1]], Index[Gluon, Ext[2]]]*(yb*(GD3x3 + Conjugate[GD3x3])*IPL[{b}] + yt*(GU3x3 + Conjugate[GU3x3])*IPL[{t}])*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[2]]])/Pi^2},
{{{G, 1}, {G, 2}, {h1, 3}, {h3, 4}}, -(gs^2*TH1x1*TH3x3*IndexDelta[Index[Gluon, Ext[1]], Index[Gluon, Ext[2]]]*(yb*(-GD3x3 + Conjugate[GD3x3])*IPL[{b}] + yt*(GU3x3 - Conjugate[GU3x3])*IPL[{t}])*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[2]]])/(32*Pi^2)},
{{{G, 1}, {G, 2}, {h2, 3}, {h3, 4}}, -(gs^2*TH1x2*TH3x3*IndexDelta[Index[Gluon, Ext[1]], Index[Gluon, Ext[2]]]*(yb*(-GD3x3 + Conjugate[GD3x3])*IPL[{b}] + yt*(GU3x3 - Conjugate[GU3x3])*IPL[{t}])*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[2]]])/(32*Pi^2)},
{{{G, 1}, {G, 2}, {h3, 3}, {h3, 4}}, ((-I/16)*gs^2*TH3x3^2*IndexDelta[Index[Gluon, Ext[1]], Index[Gluon, Ext[2]]]*(GD3x3*Conjugate[GD3x3]*IPL[{b}] + GU2x2*Conjugate[GU2x2]*IPL[{c}] + GD1x1*Conjugate[GD1x1]*IPL[{d}] + GD2x2*Conjugate[GD2x2]*IPL[{s}] + GU3x3*Conjugate[GU3x3]*IPL[{t}] + GU1x1*Conjugate[GU1x1]*IPL[{u}])*ME[Index[Lorentz, Ext[1]], Index[Lorentz, Ext[2]]])/Pi^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}]*(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[{c, G}]*TensDot[SlashedP[2], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]])/(12*FR$Eps*Pi^2) + (gs^2*IPL[{c, G}]*TensDot[SlashedP[2], ProjP][Index[Spin, Ext[1]], Index[Spin, Ext[2]]])/(12*FR$Eps*Pi^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*If[MT == 0, 0, 1]*IPL[{t, G}])/(3*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*If[MT == 0, 0, 1]*IPL[{t, G}])/(3*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*If[MB == 0, 0, 1]*IPL[{b, G}])/(3*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*If[MB == 0, 0, 1]*IPL[{b, G}])/(3*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*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[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[t], 1}, {t, 2}, {G0, 3}}, -(gs^2*yt*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*Sqrt[2]*FR$Eps*Pi^2) + (yt*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*If[MT == 0, 0, 1]*IPL[{t, G}])/(3*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*Sqrt[2]*MT)},
{{{anti[b], 1}, {b, 2}, {G0, 3}}, (gs^2*yb*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*Sqrt[2]*FR$Eps*Pi^2) - (yb*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*If[MB == 0, 0, 1]*IPL[{b, G}])/(3*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*Sqrt[2]*MB)},
{{{anti[b], 1}, {t, 2}, {anti[GP], 3}}, ((I/3)*gs^2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{t, b, G}]*(yb*ProjM[Index[Spin, Ext[1]], Index[Spin, Ext[2]]] - yt*ProjP[Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(FR$Eps*Pi^2) + (I/2)*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*(-((yb*(2*(-(gs^2*MB*If[MB == 0, 0, 1]*IPL[{b, G}])/(3*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]]])/MB) + (yt*(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*If[MT == 0, 0, 1]*IPL[{t, G}])/(3*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]]])/MT)},
{{{anti[t], 1}, {b, 2}, {GP, 3}}, ((I/3)*gs^2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{t, b, G}]*(-(yt*ProjM[Index[Spin, Ext[1]], Index[Spin, Ext[2]]]) + yb*ProjP[Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(FR$Eps*Pi^2) + (I/2)*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*((yt*(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*If[MT == 0, 0, 1]*IPL[{t, G}])/(3*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 - (yb*(2*(-(gs^2*MB*If[MB == 0, 0, 1]*IPL[{b, G}])/(3*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]]])/MB)},
{{{anti[d], 1}, {u, 2}, {anti[Hc], 3}}, ((I/3)*gs^2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{u, d, G}]*(Conjugate[GD1x1]*ProjM[Index[Spin, Ext[1]], Index[Spin, Ext[2]]] - GU1x1*ProjP[Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(FR$Eps*Pi^2) + (I/2)*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*(-(Conjugate[GD1x1]*((gs^2*IPL[{d, G}])/(12*FR$Eps*Pi^2) + (gs^2*IPL[{u, G}])/(12*FR$Eps*Pi^2))*ProjM[Index[Spin, Ext[1]], Index[Spin, Ext[2]]]) + GU1x1*((gs^2*IPL[{d, G}])/(12*FR$Eps*Pi^2) + (gs^2*IPL[{u, G}])/(12*FR$Eps*Pi^2))*ProjP[Index[Spin, Ext[1]], Index[Spin, Ext[2]]])},
{{{anti[s], 1}, {c, 2}, {anti[Hc], 3}}, ((I/3)*gs^2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{c, s, G}]*(Conjugate[GD2x2]*ProjM[Index[Spin, Ext[1]], Index[Spin, Ext[2]]] - GU2x2*ProjP[Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(FR$Eps*Pi^2) + (I/2)*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*(-(Conjugate[GD2x2]*((gs^2*IPL[{c, G}])/(12*FR$Eps*Pi^2) + (gs^2*IPL[{s, G}])/(12*FR$Eps*Pi^2))*ProjM[Index[Spin, Ext[1]], Index[Spin, Ext[2]]]) + GU2x2*((gs^2*IPL[{c, G}])/(12*FR$Eps*Pi^2) + (gs^2*IPL[{s, G}])/(12*FR$Eps*Pi^2))*ProjP[Index[Spin, Ext[1]], Index[Spin, Ext[2]]])},
{{{anti[b], 1}, {t, 2}, {anti[Hc], 3}}, ((I/3)*gs^2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{t, b, G}]*(Conjugate[GD3x3]*ProjM[Index[Spin, Ext[1]], Index[Spin, Ext[2]]] - GU3x3*ProjP[Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(FR$Eps*Pi^2) + (I/2)*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*(-(Conjugate[GD3x3]*((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]]]) + GU3x3*((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]]])},
{{{anti[u], 1}, {d, 2}, {Hc, 3}}, ((I/3)*gs^2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{u, d, G}]*(-(Conjugate[GU1x1]*ProjM[Index[Spin, Ext[1]], Index[Spin, Ext[2]]]) + GD1x1*ProjP[Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(FR$Eps*Pi^2) + (I/2)*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*(Conjugate[GU1x1]*((gs^2*IPL[{d, G}])/(12*FR$Eps*Pi^2) + (gs^2*IPL[{u, G}])/(12*FR$Eps*Pi^2))*ProjM[Index[Spin, Ext[1]], Index[Spin, Ext[2]]] - GD1x1*((gs^2*IPL[{d, G}])/(12*FR$Eps*Pi^2) + (gs^2*IPL[{u, G}])/(12*FR$Eps*Pi^2))*ProjP[Index[Spin, Ext[1]], Index[Spin, Ext[2]]])},
{{{anti[c], 1}, {s, 2}, {Hc, 3}}, ((I/3)*gs^2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{c, s, G}]*(-(Conjugate[GU2x2]*ProjM[Index[Spin, Ext[1]], Index[Spin, Ext[2]]]) + GD2x2*ProjP[Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(FR$Eps*Pi^2) + (I/2)*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*(Conjugate[GU2x2]*((gs^2*IPL[{c, G}])/(12*FR$Eps*Pi^2) + (gs^2*IPL[{s, G}])/(12*FR$Eps*Pi^2))*ProjM[Index[Spin, Ext[1]], Index[Spin, Ext[2]]] - GD2x2*((gs^2*IPL[{c, G}])/(12*FR$Eps*Pi^2) + (gs^2*IPL[{s, G}])/(12*FR$Eps*Pi^2))*ProjP[Index[Spin, Ext[1]], Index[Spin, Ext[2]]])},
{{{anti[t], 1}, {b, 2}, {Hc, 3}}, ((I/3)*gs^2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{t, b, G}]*(-(Conjugate[GU3x3]*ProjM[Index[Spin, Ext[1]], Index[Spin, Ext[2]]]) + GD3x3*ProjP[Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(FR$Eps*Pi^2) + (I/2)*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*(Conjugate[GU3x3]*((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]]] - GD3x3*((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]]])},
{{{anti[u], 1}, {u, 2}, {h1, 3}}, ((I/4)*gs^2*TH2x1*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{u, G}]*(Conjugate[GU1x1]*ProjM[Index[Spin, Ext[1]], Index[Spin, Ext[2]]] + GU1x1*ProjP[Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(Sqrt[2]*FR$Eps*Pi^2)},
{{{anti[c], 1}, {c, 2}, {h1, 3}}, ((I/4)*gs^2*TH2x1*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{c, G}]*(Conjugate[GU2x2]*ProjM[Index[Spin, Ext[1]], Index[Spin, Ext[2]]] + GU2x2*ProjP[Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(Sqrt[2]*FR$Eps*Pi^2)},
{{{anti[t], 1}, {t, 2}, {h1, 3}}, ((I/3)*gs^2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{t, G}]*((TH1x1*yt + TH2x1*Conjugate[GU3x3])*ProjM[Index[Spin, Ext[1]], Index[Spin, Ext[2]]] + (GU3x3*TH2x1 + TH1x1*yt)*ProjP[Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(Sqrt[2]*FR$Eps*Pi^2) - ((I/2)*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*(2*TH1x1*yt*(-(gs^2*MT*If[MT == 0, 0, 1]*IPL[{t, G}])/(3*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]]]) + 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)*((TH1x1*yt + TH2x1*Conjugate[GU3x3])*ProjM[Index[Spin, Ext[1]], Index[Spin, Ext[2]]] + (GU3x3*TH2x1 + TH1x1*yt)*ProjP[Index[Spin, Ext[1]], Index[Spin, Ext[2]]])))/(Sqrt[2]*MT)},
{{{anti[d], 1}, {d, 2}, {h1, 3}}, ((I/4)*gs^2*TH2x1*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{d, G}]*(Conjugate[GD1x1]*ProjM[Index[Spin, Ext[1]], Index[Spin, Ext[2]]] + GD1x1*ProjP[Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(Sqrt[2]*FR$Eps*Pi^2)},
{{{anti[s], 1}, {s, 2}, {h1, 3}}, ((I/4)*gs^2*TH2x1*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{s, G}]*(Conjugate[GD2x2]*ProjM[Index[Spin, Ext[1]], Index[Spin, Ext[2]]] + GD2x2*ProjP[Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(Sqrt[2]*FR$Eps*Pi^2)},
{{{anti[b], 1}, {b, 2}, {h1, 3}}, ((I/3)*gs^2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{b, G}]*((TH1x1*yb + TH2x1*Conjugate[GD3x3])*ProjM[Index[Spin, Ext[1]], Index[Spin, Ext[2]]] + (GD3x3*TH2x1 + TH1x1*yb)*ProjP[Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(Sqrt[2]*FR$Eps*Pi^2) - ((I/2)*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*(2*TH1x1*yb*(-(gs^2*MB*If[MB == 0, 0, 1]*IPL[{b, G}])/(3*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]]]) + 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)*((TH1x1*yb + TH2x1*Conjugate[GD3x3])*ProjM[Index[Spin, Ext[1]], Index[Spin, Ext[2]]] + (GD3x3*TH2x1 + TH1x1*yb)*ProjP[Index[Spin, Ext[1]], Index[Spin, Ext[2]]])))/(Sqrt[2]*MB)},
{{{anti[u], 1}, {u, 2}, {h2, 3}}, ((I/4)*gs^2*TH2x2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{u, G}]*(Conjugate[GU1x1]*ProjM[Index[Spin, Ext[1]], Index[Spin, Ext[2]]] + GU1x1*ProjP[Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(Sqrt[2]*FR$Eps*Pi^2)},
{{{anti[c], 1}, {c, 2}, {h2, 3}}, ((I/4)*gs^2*TH2x2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{c, G}]*(Conjugate[GU2x2]*ProjM[Index[Spin, Ext[1]], Index[Spin, Ext[2]]] + GU2x2*ProjP[Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(Sqrt[2]*FR$Eps*Pi^2)},
{{{anti[t], 1}, {t, 2}, {h2, 3}}, ((I/3)*gs^2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{t, G}]*((TH1x2*yt + TH2x2*Conjugate[GU3x3])*ProjM[Index[Spin, Ext[1]], Index[Spin, Ext[2]]] + (GU3x3*TH2x2 + TH1x2*yt)*ProjP[Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(Sqrt[2]*FR$Eps*Pi^2) - ((I/2)*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*(2*TH1x2*yt*(-(gs^2*MT*If[MT == 0, 0, 1]*IPL[{t, G}])/(3*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]]]) + 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)*((TH1x2*yt + TH2x2*Conjugate[GU3x3])*ProjM[Index[Spin, Ext[1]], Index[Spin, Ext[2]]] + (GU3x3*TH2x2 + TH1x2*yt)*ProjP[Index[Spin, Ext[1]], Index[Spin, Ext[2]]])))/(Sqrt[2]*MT)},
{{{anti[d], 1}, {d, 2}, {h2, 3}}, ((I/4)*gs^2*TH2x2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{d, G}]*(Conjugate[GD1x1]*ProjM[Index[Spin, Ext[1]], Index[Spin, Ext[2]]] + GD1x1*ProjP[Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(Sqrt[2]*FR$Eps*Pi^2)},
{{{anti[s], 1}, {s, 2}, {h2, 3}}, ((I/4)*gs^2*TH2x2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{s, G}]*(Conjugate[GD2x2]*ProjM[Index[Spin, Ext[1]], Index[Spin, Ext[2]]] + GD2x2*ProjP[Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(Sqrt[2]*FR$Eps*Pi^2)},
{{{anti[b], 1}, {b, 2}, {h2, 3}}, ((I/3)*gs^2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{b, G}]*((TH1x2*yb + TH2x2*Conjugate[GD3x3])*ProjM[Index[Spin, Ext[1]], Index[Spin, Ext[2]]] + (GD3x3*TH2x2 + TH1x2*yb)*ProjP[Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(Sqrt[2]*FR$Eps*Pi^2) - ((I/2)*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*(2*TH1x2*yb*(-(gs^2*MB*If[MB == 0, 0, 1]*IPL[{b, G}])/(3*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]]]) + 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)*((TH1x2*yb + TH2x2*Conjugate[GD3x3])*ProjM[Index[Spin, Ext[1]], Index[Spin, Ext[2]]] + (GD3x3*TH2x2 + TH1x2*yb)*ProjP[Index[Spin, Ext[1]], Index[Spin, Ext[2]]])))/(Sqrt[2]*MB)},
{{{anti[u], 1}, {u, 2}, {h3, 3}}, -(gs^2*TH3x3*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{u, G}]*(Conjugate[GU1x1]*ProjM[Index[Spin, Ext[1]], Index[Spin, Ext[2]]] - GU1x1*ProjP[Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(4*Sqrt[2]*FR$Eps*Pi^2)},
{{{anti[c], 1}, {c, 2}, {h3, 3}}, -(gs^2*TH3x3*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{c, G}]*(Conjugate[GU2x2]*ProjM[Index[Spin, Ext[1]], Index[Spin, Ext[2]]] - GU2x2*ProjP[Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(4*Sqrt[2]*FR$Eps*Pi^2)},
{{{anti[t], 1}, {t, 2}, {h3, 3}}, -(gs^2*TH3x3*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{t, G}]*(Conjugate[GU3x3]*ProjM[Index[Spin, Ext[1]], Index[Spin, Ext[2]]] - GU3x3*ProjP[Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(3*Sqrt[2]*FR$Eps*Pi^2) + (TH3x3*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*((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)*(Conjugate[GU3x3]*ProjM[Index[Spin, Ext[1]], Index[Spin, Ext[2]]] - GU3x3*ProjP[Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(2*Sqrt[2])},
{{{anti[d], 1}, {d, 2}, {h3, 3}}, (gs^2*TH3x3*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{d, G}]*(Conjugate[GD1x1]*ProjM[Index[Spin, Ext[1]], Index[Spin, Ext[2]]] - GD1x1*ProjP[Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(4*Sqrt[2]*FR$Eps*Pi^2)},
{{{anti[s], 1}, {s, 2}, {h3, 3}}, (gs^2*TH3x3*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{s, G}]*(Conjugate[GD2x2]*ProjM[Index[Spin, Ext[1]], Index[Spin, Ext[2]]] - GD2x2*ProjP[Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(4*Sqrt[2]*FR$Eps*Pi^2)},
{{{anti[b], 1}, {b, 2}, {h3, 3}}, (gs^2*TH3x3*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{b, G}]*(Conjugate[GD3x3]*ProjM[Index[Spin, Ext[1]], Index[Spin, Ext[2]]] - GD3x3*ProjP[Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(3*Sqrt[2]*FR$Eps*Pi^2) - (TH3x3*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*((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)*(Conjugate[GD3x3]*ProjM[Index[Spin, Ext[1]], Index[Spin, Ext[2]]] - GD3x3*ProjP[Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(2*Sqrt[2])},
{{{anti[u], 1}, {u, 2}, {A, 3}}, ((-I/72)*gs^2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{u, G}]*((cw*g1 + 3*gw*sw)*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]] + 4*cw*g1*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjP][Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(FR$Eps*Pi^2) + (I/6)*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*((gs^2*(cw*g1 + 3*gw*sw)*IPL[{u, G}]*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]])/(12*FR$Eps*Pi^2) + (cw*g1*gs^2*IPL[{u, G}]*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjP][Index[Spin, Ext[1]], Index[Spin, Ext[2]]])/(3*FR$Eps*Pi^2))},
{{{anti[c], 1}, {c, 2}, {A, 3}}, ((-I/72)*gs^2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{c, G}]*((cw*g1 + 3*gw*sw)*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]] + 4*cw*g1*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjP][Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(FR$Eps*Pi^2) + (I/6)*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*((gs^2*(cw*g1 + 3*gw*sw)*IPL[{c, G}]*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]])/(12*FR$Eps*Pi^2) + (cw*g1*gs^2*IPL[{c, G}]*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjP][Index[Spin, Ext[1]], Index[Spin, Ext[2]]])/(3*FR$Eps*Pi^2))},
{{{anti[t], 1}, {t, 2}, {A, 3}}, ((-I/72)*gs^2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{t, G}]*((cw*g1 + 3*gw*sw)*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]] + 4*cw*g1*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjP][Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(FR$Eps*Pi^2) + (I/6)*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*((cw*g1 + 3*gw*sw)*((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*cw*g1*((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/72)*gs^2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{d, G}]*((cw*g1 - 3*gw*sw)*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]] - 2*cw*g1*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjP][Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(FR$Eps*Pi^2) + (I/6)*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*((gs^2*(cw*g1 - 3*gw*sw)*IPL[{d, G}]*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]])/(12*FR$Eps*Pi^2) - (cw*g1*gs^2*IPL[{d, G}]*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjP][Index[Spin, Ext[1]], Index[Spin, Ext[2]]])/(6*FR$Eps*Pi^2))},
{{{anti[s], 1}, {s, 2}, {A, 3}}, ((-I/72)*gs^2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{s, G}]*((cw*g1 - 3*gw*sw)*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]] - 2*cw*g1*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjP][Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(FR$Eps*Pi^2) + (I/6)*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*((gs^2*(cw*g1 - 3*gw*sw)*IPL[{s, G}]*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]])/(12*FR$Eps*Pi^2) - (cw*g1*gs^2*IPL[{s, G}]*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjP][Index[Spin, Ext[1]], Index[Spin, Ext[2]]])/(6*FR$Eps*Pi^2))},
{{{anti[b], 1}, {b, 2}, {A, 3}}, ((-I/72)*gs^2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{b, G}]*((cw*g1 - 3*gw*sw)*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]] - 2*cw*g1*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjP][Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(FR$Eps*Pi^2) + (I/6)*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*((cw*g1 - 3*gw*sw)*((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*cw*g1*((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)*gs^2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{u, G}]*((3*cw*gw - g1*sw)*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]] - 4*g1*sw*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjP][Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(FR$Eps*Pi^2) + (I/6)*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*((gs^2*(3*cw*gw - g1*sw)*IPL[{u, G}]*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]])/(12*FR$Eps*Pi^2) - (g1*gs^2*sw*IPL[{u, G}]*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjP][Index[Spin, Ext[1]], Index[Spin, Ext[2]]])/(3*FR$Eps*Pi^2))},
{{{anti[c], 1}, {c, 2}, {Z, 3}}, ((-I/72)*gs^2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{c, G}]*((3*cw*gw - g1*sw)*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]] - 4*g1*sw*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjP][Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(FR$Eps*Pi^2) + (I/6)*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*((gs^2*(3*cw*gw - g1*sw)*IPL[{c, G}]*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]])/(12*FR$Eps*Pi^2) - (g1*gs^2*sw*IPL[{c, G}]*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjP][Index[Spin, Ext[1]], Index[Spin, Ext[2]]])/(3*FR$Eps*Pi^2))},
{{{anti[t], 1}, {t, 2}, {Z, 3}}, ((-I/72)*gs^2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{t, G}]*((3*cw*gw - g1*sw)*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]] - 4*g1*sw*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjP][Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(FR$Eps*Pi^2) + (I/6)*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*((3*cw*gw - g1*sw)*((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*g1*sw*((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}, {Z, 3}}, ((I/72)*gs^2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{d, G}]*((3*cw*gw + g1*sw)*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]] - 2*g1*sw*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjP][Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(FR$Eps*Pi^2) + (I/6)*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*(-(gs^2*(3*cw*gw + g1*sw)*IPL[{d, G}]*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]])/(12*FR$Eps*Pi^2) + (g1*gs^2*sw*IPL[{d, G}]*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjP][Index[Spin, Ext[1]], Index[Spin, Ext[2]]])/(6*FR$Eps*Pi^2))},
{{{anti[s], 1}, {s, 2}, {Z, 3}}, ((I/72)*gs^2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{s, G}]*((3*cw*gw + g1*sw)*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]] - 2*g1*sw*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjP][Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(FR$Eps*Pi^2) + (I/6)*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*(-(gs^2*(3*cw*gw + g1*sw)*IPL[{s, G}]*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]])/(12*FR$Eps*Pi^2) + (g1*gs^2*sw*IPL[{s, G}]*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjP][Index[Spin, Ext[1]], Index[Spin, Ext[2]]])/(6*FR$Eps*Pi^2))},
{{{anti[b], 1}, {b, 2}, {Z, 3}}, ((I/72)*gs^2*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*IPL[{b, G}]*((3*cw*gw + g1*sw)*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]] - 2*g1*sw*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjP][Index[Spin, Ext[1]], Index[Spin, Ext[2]]]))/(FR$Eps*Pi^2) + (I/6)*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*(-((3*cw*gw + g1*sw)*((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*g1*sw*((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[d], 1}, {u, 2}, {anti[W], 3}}, ((I/2)*gw*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] - ((I/12)*gs^2*gw*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)},
{{{anti[s], 1}, {c, 2}, {anti[W], 3}}, ((I/2)*gw*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*((gs^2*IPL[{c, G}])/(12*FR$Eps*Pi^2) + (gs^2*IPL[{s, G}])/(12*FR$Eps*Pi^2))*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]])/Sqrt[2] - ((I/12)*gs^2*gw*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)},
{{{anti[b], 1}, {t, 2}, {anti[W], 3}}, ((-I/12)*gs^2*gw*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) + ((I/2)*gw*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]},
{{{anti[u], 1}, {d, 2}, {W, 3}}, ((I/2)*gw*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] - ((I/12)*gs^2*gw*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)},
{{{anti[c], 1}, {s, 2}, {W, 3}}, ((I/2)*gw*IndexDelta[Index[Colour, Ext[1]], Index[Colour, Ext[2]]]*((gs^2*IPL[{c, G}])/(12*FR$Eps*Pi^2) + (gs^2*IPL[{s, G}])/(12*FR$Eps*Pi^2))*TensDot[Ga[Index[Lorentz, Ext[3]]], ProjM][Index[Spin, Ext[1]], Index[Spin, Ext[2]]])/Sqrt[2] - ((I/12)*gs^2*gw*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)},
{{{anti[t], 1}, {b, 2}, {W, 3}}, ((-I/12)*gs^2*gw*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) + ((I/2)*gw*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]},
{{{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[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*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[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[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*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[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[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*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[{c, 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[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[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*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[{c, 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[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[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*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[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[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*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[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[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*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[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[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*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[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[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*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[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[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*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[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[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*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[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[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*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[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[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*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[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[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*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[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 = {
};