
(* description = "pp -> Vtt+jets|pp -> Vtt+jets|pp -> Ztt|On-shell Z plus t t-bar production at NLO EW (including all subleading one-loop orders)." *)


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

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

(* ttxzagg *)  AddReal[FeynArtsProcess -> {F[3,{3}], -F[3,{3}]} -> {V[5],V[5],V[2], V[1]}, SelectInterference -> {gQCD -> 4}];
(* uuxttxza *) AddReal[FeynArtsProcess -> {F[3,{1}], -F[3,{1}]}  -> {F[3,{3}],-F[3,{3}],V[2], V[1]}, SelectInterference -> {gQCD -> 4}];
(* ddxttxza *) AddReal[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]}  -> {F[3,{3}],-F[3,{3}],V[2], V[1]}, SelectInterference -> {gQCD -> 4}];
(* bbxttxza *) AddReal[FeynArtsProcess -> {F[4,{3}], -F[4,{3}]}  -> {F[3,{3}],-F[3,{3}],V[2], V[1]}, SelectInterference -> {gQCD -> 4}];

(* uuxttxzg *) AddReal[FeynArtsProcess -> {F[3,{3}], -F[3,{3}]} -> {F[3,{1}],-F[3,{1}],V[2],V[5]}, SelectInterference -> {gQCD -> 4}, SelectCoupling -> (Exponent[#1, eQED] == 1 || Exponent[#1, eQED] == 3 &)];
(* ddxttxzg *) AddReal[FeynArtsProcess -> {F[3,{3}], -F[3,{3}]} -> {F[4,{1}],-F[4,{1}],V[2],V[5]}, SelectInterference -> {gQCD -> 4}, SelectCoupling -> (Exponent[#1, eQED] == 1 || Exponent[#1, eQED] == 3 &)];
(* bbxttxzg *) AddReal[FeynArtsProcess -> {F[3,{3}], -F[3,{3}]} -> {F[4,{3}],-F[4,{3}],V[2],V[5]}, SelectInterference -> {gQCD -> 4}, SelectCoupling -> (Exponent[#1, eQED] == 1 || Exponent[#1, eQED] == 3 &)];


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

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


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

(* uuxttxz *) AddTree[FeynArtsProcess -> {F[3,{1}], -F[3,{1}]} -> {F[3,{3}],-F[3,{3}],V[2]}, SelectInterference -> {gQCD -> 2}, SelectCoupling -> (Exponent[#1, gQCD] == 2 - 2 * #2 || Exponent[#1, gQCD] == 0 + 2 * #2 &),  UnitaryGauge -> False];
(* ttxddxz *) AddTree[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {F[3,{3}],-F[3,{3}],V[2]}, SelectInterference -> {gQCD -> 2}, SelectCoupling -> (Exponent[#1, gQCD] == 2 - 2 * #2 || Exponent[#1, gQCD] == 0 + 2 * #2 &),  UnitaryGauge -> False];
(* bbxttxz *) AddTree[FeynArtsProcess -> {F[4,{3}], -F[4,{3}]} -> {F[3,{3}],-F[3,{3}],V[2]}, SelectInterference -> {gQCD -> 2}, SelectCoupling -> (Exponent[#1, gQCD] == 2 - 2 * #2 || Exponent[#1, gQCD] == 0 + 2 * #2 &),  UnitaryGauge -> False];


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

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


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

(* uuxttxz *) AddProcess[FeynArtsProcess -> {F[3,{1}], -F[3,{1}]} -> {F[3,{3}],-F[3,{3}],V[2]}, SelectInterference -> {eQED -> 6}, SelectCoupling -> (Exponent[#1, gQCD] == 2 - 2 * #2 || Exponent[#1, gQCD] == 0 + 2 * #2 &),  UnitaryGauge -> False];
(* ttxddxz *) AddProcess[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {F[3,{3}],-F[3,{3}],V[2]}, SelectInterference -> {eQED -> 6}, SelectCoupling -> (Exponent[#1, gQCD] == 2 - 2 * #2 || Exponent[#1, gQCD] == 0 + 2 * #2 &),  UnitaryGauge -> False];
(* bbxttxz *) AddProcess[FeynArtsProcess -> {F[4,{3}], -F[4,{3}]} -> {F[3,{3}],-F[3,{3}],V[2]}, SelectInterference -> {eQED -> 6}, SelectCoupling -> (Exponent[#1, gQCD] == 2 - 2 * #2 || Exponent[#1, gQCD] == 0 + 2 * #2 &),  UnitaryGauge -> False];
(* ttxzaa *) AddProcess[FeynArtsProcess -> {V[1], V[1]} -> {F[3,{3}],-F[3,{3}],V[2]}, SelectInterference -> {eQED -> 6}, SelectCoupling -> (Exponent[#1, gQCD] == 0 + 2 * #2 &),  UnitaryGauge -> False];

(* uuxttxza *) AddReal[FeynArtsProcess -> {F[3,{1}], -F[3,{1}]} -> {F[3,{3}],-F[3,{3}],V[2],V[1]}, SelectInterference -> {eQED -> 6}, SelectCoupling -> (Exponent[#1, gQCD] == 0 || Exponent[#1, gQCD] == 2 &)];
(* ttxddxza *) AddReal[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {F[3,{3}],-F[3,{3}],V[2],V[1]}, SelectInterference -> {eQED -> 6}, SelectCoupling -> (Exponent[#1, gQCD] == 0 || Exponent[#1, gQCD] == 2 &)];
(* bbxttxza *) AddReal[FeynArtsProcess -> {F[4,{3}], -F[4,{3}]} -> {F[3,{3}],-F[3,{3}],V[2],V[1]}, SelectInterference -> {eQED -> 6}, SelectCoupling -> (Exponent[#1, gQCD] == 0 || Exponent[#1, gQCD] == 2 &)];
(* ttxzaaa *) AddReal[FeynArtsProcess -> {V[1], V[1]} -> {F[3,{3}],-F[3,{3}],V[1]}, SelectInterference -> {eQED -> 6}, SelectCoupling -> (Exponent[#1, gQCD] == 0 &)];

(* uuxttxzg *) AddReal[FeynArtsProcess -> {F[3,{1}], -F[3,{1}]} -> {F[3,{3}],-F[3,{3}],V[2],V[5]}, SelectInterference -> {eQED -> 6}, SelectCoupling -> (Exponent[#1, gQCD] == 1 &)];
(* ttxddxzg *) AddReal[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {F[3,{3}],-F[3,{3}],V[2],V[5]}, SelectInterference -> {eQED -> 6}, SelectCoupling -> (Exponent[#1, gQCD] == 1 &)];
(* bbxttxzg *) AddReal[FeynArtsProcess -> {F[4,{3}], -F[4,{3}]} -> {F[3,{3}],-F[3,{3}],V[2],V[5]}, SelectInterference -> {eQED -> 6}, SelectCoupling -> (Exponent[#1, gQCD] == 1 &)];
(* ttxzaag *) AddReal[FeynArtsProcess -> {V[1], V[1]} -> {F[3,{3}],-F[3,{3}],V[2],V[5]}, SelectInterference -> {eQED -> 6}, SelectCoupling -> (Exponent[#1, gQCD] == 1 &)];


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

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

(* uuxttxza *) AddReal[FeynArtsProcess -> {F[3,{1}], -F[3,{1}]} -> {F[3,{3}],-F[3,{3}],V[1],V[2]}, SelectInterference -> {gQCD -> 0}];
(* ttxddxza *) AddReal[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {F[3,{3}],-F[3,{3}],V[1],V[2]}, SelectInterference -> {gQCD -> 0}];
(* bbxttxza *) AddReal[FeynArtsProcess -> {F[4,{3}], -F[4,{3}]} -> {F[3,{3}],-F[3,{3}],V[1],V[2]}, SelectInterference -> {gQCD -> 0}];
(* ttxzaaa *) AddReal[FeynArtsProcess -> {V[1], V[1]} -> {F[3,{3}],-F[3,{3}],V[1],V[2]}, SelectInterference -> {gQCD -> 0}];









