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


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

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

(* uttxdxwa *) AddReal[FeynArtsProcess -> {-F[4,{1}], F[3,{1}]}  -> {F[3,{3}],-F[3,{3}],-V[3], V[1]}, SelectInterference -> {gQCD -> 4}];
(* uxttxdwxa *) AddReal[FeynArtsProcess -> {F[4,{1}], -F[3,{1}]}  -> {F[3,{3}],-F[3,{3}],V[3], V[1]}, SelectInterference -> {gQCD -> 4}];

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


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

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



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

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

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

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


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

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

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









