
(* description = "pp -> tt+jets|pp -> tt|pp -> tt|Top-quark pair production including all subleading one-loop orders." *)


(* O(as^2) + O(as^2 a) *)

(* ttxggg *)  AddProcess[FeynArtsProcess -> {F[3,{3}], -F[3,{3}]} -> {V[5],V[5]}, UnitaryGauge -> False, SelectInterference -> {gQCD -> 4}, SelectCoupling -> (Exponent[#1, eQED] == 0 + 2 * #2 &)];

(* uuxttxg *) AddProcess[FeynArtsProcess -> {F[3,{1}], -F[3,{1}]}  -> {F[3,{3}],-F[3,{3}]}, UnitaryGauge -> False, SelectInterference -> {gQCD -> 4}, SelectCoupling -> (Exponent[#1, eQED] == 2 - 2 * #2 || Exponent[#1, eQED] == 0 + 2 * #2 &)];
(* ttxddxg *) AddProcess[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {F[3,{3}],-F[3,{3}]}, UnitaryGauge -> False, SelectInterference -> {gQCD -> 4}, SelectCoupling -> (Exponent[#1, eQED] == 2 - 2 * #2 || Exponent[#1, eQED] == 0 + 2 * #2 &)];
(* ttxbbxg *) AddProcess[FeynArtsProcess -> {F[4,{3}], -F[4,{3}]} -> {F[3,{3}],-F[3,{3}]}, UnitaryGauge -> False, SelectInterference -> {gQCD -> 4}, SelectCoupling -> (Exponent[#1, eQED] == 2 - 2 * #2 || Exponent[#1, eQED] == 0 + 2 * #2 &)];



(* O(as*a) + O(as a^2) *)

(* uuxttxg *) AddProcess[FeynArtsProcess -> {F[3,{1}], -F[3,{1}]} -> {F[3,{3}],-F[3,{3}]}, SelectInterference -> {gQCD -> 2}, SelectCoupling -> (Exponent[#1, gQCD] == 2 - 2 * #2 || Exponent[#1, gQCD] == 0 + 2 * #2 &),  UnitaryGauge -> False];
(* ttxddxg *) AddProcess[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {F[3,{3}],-F[3,{3}]}, SelectInterference -> {gQCD -> 2}, SelectCoupling -> (Exponent[#1, gQCD] == 2 - 2 * #2 || Exponent[#1, gQCD] == 0 + 2 * #2 &),  UnitaryGauge -> False];
(* bbxttxg *) AddProcess[FeynArtsProcess -> {F[4,{3}], -F[4,{3}]} -> {F[3,{3}],-F[3,{3}]}, SelectInterference -> {gQCD -> 2}, SelectCoupling -> (Exponent[#1, gQCD] == 2 - 2 * #2 || Exponent[#1, gQCD] == 0 + 2 * #2 &),  UnitaryGauge -> False];


(* O(a^2) + O(a^3) *)

(* uuxttxg *) AddProcess[FeynArtsProcess -> {F[3,{1}], -F[3,{1}]} -> {F[3,{3}],-F[3,{3}]}, SelectInterference -> {gQCD -> 0}, SelectCoupling -> (Exponent[#1, eQED] == 2 + 2 * #2 &),  UnitaryGauge -> False];
(* ttxddxg *) AddProcess[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {F[3,{3}],-F[3,{3}]}, SelectInterference -> {gQCD -> 0}, SelectCoupling -> (Exponent[#1, eQED] == 2 + 2 * #2 &),  UnitaryGauge -> False];
(* bbxttxg *) AddProcess[FeynArtsProcess -> {F[4,{3}], -F[4,{3}]} -> {F[3,{3}],-F[3,{3}]}, SelectInterference -> {gQCD -> 0}, SelectCoupling -> (Exponent[#1, eQED] == 2 + 2 * #2 &),  UnitaryGauge -> False];









