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


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

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

(* ttxaagg *)  AddReal[FeynArtsProcess -> {F[3,{3}], -F[3,{3}]} -> {V[5],V[5],V[1], V[1]}, SelectInterference -> {gQCD -> 4}];
(* uuxttxaa *) AddReal[FeynArtsProcess -> {F[3,{1}], -F[3,{1}]}  -> {F[3,{3}],-F[3,{3}],V[1], V[1]}, SelectInterference -> {gQCD -> 4}];
(* ddxttxaa *) AddReal[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]}  -> {F[3,{3}],-F[3,{3}],V[1], V[1]}, SelectInterference -> {gQCD -> 4}];
(* bbxttxaa *) AddReal[FeynArtsProcess -> {F[3,{3}], -F[4,{3}]}  -> {F[3,{3}],-F[3,{3}],V[1], V[1]}, SelectInterference -> {gQCD -> 4}];

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


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

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


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

(* uuxttxa *) AddTree[FeynArtsProcess -> {F[3,{1}], -F[3,{1}]} -> {F[3,{3}],-F[3,{3}],V[1]}, SelectInterference -> {gQCD -> 2}, SelectCoupling -> (Exponent[#1, gQCD] == 2 - 2 * #2 || Exponent[#1, gQCD] == 0 + 2 * #2 &),  UnitaryGauge -> False];
(* ttxddxa *) AddTree[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {F[3,{3}],-F[3,{3}],V[1]}, SelectInterference -> {gQCD -> 2}, SelectCoupling -> (Exponent[#1, gQCD] == 2 - 2 * #2 || Exponent[#1, gQCD] == 0 + 2 * #2 &),  UnitaryGauge -> False];
(* bbxttxa *) AddTree[FeynArtsProcess -> {F[4,{3}], -F[4,{3}]} -> {F[3,{3}],-F[3,{3}],V[1]}, 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) *)

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


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

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

(* uuxttxaa *) AddReal[FeynArtsProcess -> {F[3,{1}], -F[3,{1}]} -> {F[3,{3}],-F[3,{3}],V[1],V[1]}, SelectInterference -> {eQED -> 6}, SelectCoupling -> (Exponent[#1, gQCD] == 0 || Exponent[#1, gQCD] == 2 &)];
(* ttxddxaa *) AddReal[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {F[3,{3}],-F[3,{3}],V[1],V[1]}, SelectInterference -> {eQED -> 6}, SelectCoupling -> (Exponent[#1, gQCD] == 0 || Exponent[#1, gQCD] == 2 &)];
(* bbxttxaa *) AddReal[FeynArtsProcess -> {F[4,{3}], -F[4,{3}]} -> {F[3,{3}],-F[3,{3}],V[1],V[1]}, SelectInterference -> {eQED -> 6}, SelectCoupling -> (Exponent[#1, gQCD] == 0 || Exponent[#1, gQCD] == 2 &)];
(* ttxaaaa *) AddReal[FeynArtsProcess -> {V[1], V[1]} -> {F[3,{3}],-F[3,{3}],V[1]}, SelectInterference -> {eQED -> 6}, SelectCoupling -> (Exponent[#1, gQCD] == 0 &)];

(* uuxttxag *) AddReal[FeynArtsProcess -> {F[3,{1}], -F[3,{1}]} -> {F[3,{3}],-F[3,{3}],V[1],V[5]}, SelectInterference -> {eQED -> 6}, SelectCoupling -> (Exponent[#1, gQCD] == 1 &)];
(* ttxddxag *) AddReal[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {F[3,{3}],-F[3,{3}],V[1],V[5]}, SelectInterference -> {eQED -> 6}, SelectCoupling -> (Exponent[#1, gQCD] == 1 &)];
(* bbxttxag *) AddReal[FeynArtsProcess -> {F[4,{3}], -F[4,{3}]} -> {F[3,{3}],-F[3,{3}],V[1],V[5]}, SelectInterference -> {eQED -> 6}, SelectCoupling -> (Exponent[#1, gQCD] == 1 &)];
(* ttxaaag *) AddReal[FeynArtsProcess -> {V[1], V[1]} -> {F[3,{3}],-F[3,{3}],V[1],V[5]}, SelectInterference -> {eQED -> 6}, SelectCoupling -> (Exponent[#1, gQCD] == 1 &)];


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

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

(* uuxttxaa *) AddReal[FeynArtsProcess -> {F[3,{1}], -F[3,{1}]} -> {F[3,{3}],-F[3,{3}],V[1],V[1]}, SelectInterference -> {gQCD -> 0}];
(* ttxddxaa *) AddReal[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {F[3,{3}],-F[3,{3}],V[1],V[1]}, SelectInterference -> {gQCD -> 0}];
(* bbxttxaa *) AddReal[FeynArtsProcess -> {F[4,{3}], -F[4,{3}]} -> {F[3,{3}],-F[3,{3}],V[1],V[1]}, SelectInterference -> {gQCD -> 0}];
(* ttxaaaa *) AddReal[FeynArtsProcess -> {V[1], V[1]} -> {F[3,{3}],-F[3,{3}],V[1],V[1]}, SelectInterference -> {gQCD -> 0}];









