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


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

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

(* ttxhagg *)  AddReal[FeynArtsProcess -> {F[3,{3}], -F[3,{3}]} -> {V[5],V[5],S[1], V[1]}, SelectInterference -> {gQCD -> 4}];
(* uuxttxha *) AddReal[FeynArtsProcess -> {F[3,{1}], -F[3,{1}]}  -> {F[3,{3}],-F[3,{3}],S[1], V[1]}, SelectInterference -> {gQCD -> 4}];
(* ddxttxha *) AddReal[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]}  -> {F[3,{3}],-F[3,{3}],S[1], V[1]}, SelectInterference -> {gQCD -> 4}];
(* bbxttxha *) AddReal[FeynArtsProcess -> {F[3,{3}], -F[4,{3}]}  -> {F[3,{3}],-F[3,{3}],S[1], V[1]}, SelectInterference -> {gQCD -> 4}];

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


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

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


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

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

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


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

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

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

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


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

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

(* uuxttxha *) AddReal[FeynArtsProcess -> {F[3,{1}], -F[3,{1}]} -> {F[3,{3}],-F[3,{3}],V[1],S[1]}, SelectInterference -> {gQCD -> 0}];
(* ttxddxha *) AddReal[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {F[3,{3}],-F[3,{3}],V[1],S[1]}, SelectInterference -> {gQCD -> 0}];
(* bbxttxha *) AddReal[FeynArtsProcess -> {F[4,{3}], -F[4,{3}]} -> {F[3,{3}],-F[3,{3}],V[1],S[1]}, SelectInterference -> {gQCD -> 0}];
(* ttxhaaa *) AddReal[FeynArtsProcess -> {V[1], V[1]} -> {F[3,{3}],-F[3,{3}],V[1],S[1]}, SelectInterference -> {gQCD -> 0}];









