2HDM: 2HDM.nlo

File 2HDM.nlo, 121.8 KB (added by degrande, 2 years ago)
Line 
1(****************************************************************************************************************)
2(* Model automalically generated by NLOCT-1.02                                                                  *)
3(* date and time of generation : Thu 15 Oct 2015 14:05:33                                                       *)
4(* FeynRules model information :                                                                                *)
5(*   ModelName : 2HDM                                                                                           *)
6(*   Authors : C. Duhr, M. Herquet, C. Degrande                                                                 *)
7(*   Version : 1.2                                                                                              *)
8(*   Date : 12. 10. 2015                                                                                        *)
9(*   Institutions : Universite catholique de Louvain (CP3), NIKHEF Theory Dept., IPPP                           *)
10(*   URLs : http://feynrules.phys.ucl.ac.be/view/Main/2HDM                                                      *)
11(*   Emails : claude.duhr@uclouvain.be, mherquet@nikhef.nl, celine.degrande@durham.ac.uk                        *)
12(****************************************************************************************************************)
13
14(* FeynArts model files used to generated the output *)
15
16CT$Model = "2HDM/2HDM";
17CT$GenericModel = "Lorentz";
18
19(* Those assumptions where made : *)
20
21NLOCT$assumptions = {};
22
23(* Perturbed orders *)
24
25FR$InteractionOrderPerturbativeExpansion = {{QCD, 1}, {QED, 0}, {YB, 0}, {YT, 0}};
26
27(*R2 vertices*)
28
29R2$vertlist = {
30{{{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},
31{{{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},
32{{{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},
33{{{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},
34{{{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},
35{{{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},
36{{{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},
37{{{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)},
38{{{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)},
39{{{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},
40{{{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},
41{{{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},
42{{{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},
43{{{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},
44{{{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},
45{{{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},
46{{{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},
47{{{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)},
48{{{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)},
49{{{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)},
50{{{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)},
51{{{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)},
52{{{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)},
53{{{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)},
54{{{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)},
55{{{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)},
56{{{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)},
57{{{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)},
58{{{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)},
59{{{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)},
60{{{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)},
61{{{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)},
62{{{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)},
63{{{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)},
64{{{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)},
65{{{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},
66{{{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},
67{{{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},
68{{{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},
69{{{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},
70{{{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},
71{{{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},
72{{{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},
73{{{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},
74{{{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},
75{{{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},
76{{{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},
77{{{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)},
78{{{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)},
79{{{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)},
80{{{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)},
81{{{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)},
82{{{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)},
83{{{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},
84{{{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},
85{{{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},
86{{{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},
87{{{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},
88{{{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},
89{{{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)},
90{{{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)},
91{{{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)},
92{{{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)},
93{{{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},
94{{{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},
95{{{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},
96{{{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},
97{{{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},
98{{{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},
99{{{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},
100{{{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},
101{{{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},
102{{{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},
103{{{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},
104{{{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},
105{{{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)},
106{{{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},
107{{{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)},
108{{{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},
109{{{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},
110{{{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},
111{{{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)},
112{{{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)},
113{{{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}};
114
115(*UV vertices*)
116
117UV$vertlist = {
118{{{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))},
119{{{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))},
120{{{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]]])},
121{{{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))},
122{{{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))},
123{{{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]]])},
124{{{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)},
125{{{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)},
126{{{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)},
127{{{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)},
128{{{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)},
129{{{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]]])},
130{{{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]]])},
131{{{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]]])},
132{{{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]]])},
133{{{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]]])},
134{{{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]]])},
135{{{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)},
136{{{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)},
137{{{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)},
138{{{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)},
139{{{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)},
140{{{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)},
141{{{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)},
142{{{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)},
143{{{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)},
144{{{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)},
145{{{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)},
146{{{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)},
147{{{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)},
148{{{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)},
149{{{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])},
150{{{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)},
151{{{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)},
152{{{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])},
153{{{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))},
154{{{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))},
155{{{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]]])},
156{{{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))},
157{{{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))},
158{{{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]]])},
159{{{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))},
160{{{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))},
161{{{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]]])},
162{{{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))},
163{{{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))},
164{{{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]]])},
165{{{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)},
166{{{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)},
167{{{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]},
168{{{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)},
169{{{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)},
170{{{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]},
171{{{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},
172{{{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},
173{{{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},
174{{{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},
175{{{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},
176{{{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},
177{{{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)},
178{{{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)}};
179
180(*CT parameters*)
181
182FR$CTparam = {
183};