
(* description = "pp -> jets|pp -> jets|pp -> jj|Production of two jets at NLO EW including all subleading one-loop orders." *)


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

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

(* dddxdx *) AddProcess[FeynArtsProcess -> {F[4,{1}], F[4,{1}]}  -> {F[4,{1}],F[4,{1}]}, UnitaryGauge -> False, SelectInterference -> {gQCD -> 4}, SelectCoupling -> (Exponent[#1, eQED] == 2 - 2 * #2 || Exponent[#1, eQED] == 0 + 2 * #2 &)];
(* uuxddx *) AddProcess[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {F[3,{1}],-F[3,{1}]}, UnitaryGauge -> False, SelectInterference -> {gQCD -> 4}, SelectCoupling -> (Exponent[#1, eQED] == 2 - 2 * #2 || Exponent[#1, eQED] == 0 + 2 * #2 &)];
(* ddxssx *) AddProcess[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {F[4,{2}],-F[4,{2}]}, UnitaryGauge -> False, SelectInterference -> {gQCD -> 4}, SelectCoupling -> (Exponent[#1, eQED] == 2 - 2 * #2 || Exponent[#1, eQED] == 0 + 2 * #2 &)];
(* ccxddx *) AddProcess[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {F[3,{2}],-F[3,{2}]}, UnitaryGauge -> False, SelectInterference -> {gQCD -> 4}, SelectCoupling -> (Exponent[#1, eQED] == 2 - 2 * #2 || Exponent[#1, eQED] == 0 + 2 * #2 &)];
(* uxcdsx *) AddProcess[FeynArtsProcess -> {F[4,{1}], -F[3,{1}]} -> {F[4,{2}],-F[3,{2}]}, UnitaryGauge -> False, SelectInterference -> {gQCD -> 4}, SelectCoupling -> (Exponent[#1, eQED] == 2 - 2 * #2 || Exponent[#1, eQED] == 0 + 2 * #2 &)];
(* uuuxux *) AddProcess[FeynArtsProcess -> {F[3,{1}], F[3,{1}]}  -> {F[3,{1}],F[3,{1}]}, UnitaryGauge -> False, SelectInterference -> {gQCD -> 4}, SelectCoupling -> (Exponent[#1, eQED] == 2 - 2 * #2 || Exponent[#1, eQED] == 0 + 2 * #2 &)];
(* uuxccx *) AddProcess[FeynArtsProcess -> {F[3,{1}], -F[3,{1}]} -> {F[3,{2}],-F[3,{2}]}, UnitaryGauge -> False, SelectInterference -> {gQCD -> 4}, SelectCoupling -> (Exponent[#1, eQED] == 2 - 2 * #2 || Exponent[#1, eQED] == 0 + 2 * #2 &)];
(* ddxbbx *) AddProcess[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {F[4,{3}],-F[4,{3}]}, UnitaryGauge -> False, SelectInterference -> {gQCD -> 4}, SelectCoupling -> (Exponent[#1, eQED] == 2 - 2 * #2 || Exponent[#1, eQED] == 0 + 2 * #2 &)];
(* uuxbbx *) AddProcess[FeynArtsProcess -> {F[3,{1}], -F[3,{1}]} -> {F[4,{3}],-F[4,{3}]}, UnitaryGauge -> False, SelectInterference -> {gQCD -> 4}, SelectCoupling -> (Exponent[#1, eQED] == 2 - 2 * #2 || Exponent[#1, eQED] == 0 + 2 * #2 &)];
(* bbbxbx *) AddProcess[FeynArtsProcess -> {F[4,{3}], F[4,{3}]}  -> {F[4,{3}],F[4,{3}]}, UnitaryGauge -> False, SelectInterference -> {gQCD -> 4}, SelectCoupling -> (Exponent[#1, eQED] == 2 - 2 * #2 || Exponent[#1, eQED] == 0 + 2 * #2 &)];

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

(* dddxdxa *) AddReal[FeynArtsProcess -> {F[4,{1}], F[4,{1}]}  -> {F[4,{1}],F[4,{1}],V[1]}, SelectInterference -> {gQCD -> 4}, SelectCoupling -> (Exponent[#1, eQED] == 1 &)];
(* uuxddxa *) AddReal[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {F[3,{1}],-F[3,{1}],V[1]}, SelectInterference -> {gQCD -> 4}, SelectCoupling -> (Exponent[#1, eQED] == 1 &)];
(* ddxssxa *) AddReal[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {F[4,{2}],-F[4,{2}],V[1]}, SelectInterference -> {gQCD -> 4}, SelectCoupling -> (Exponent[#1, eQED] == 1 &)];
(* ccxddxa *) AddReal[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {F[3,{2}],-F[3,{2}],V[1]}, SelectInterference -> {gQCD -> 4}, SelectCoupling -> (Exponent[#1, eQED] == 1 &)];
(* uuuxuxa *) AddReal[FeynArtsProcess -> {F[3,{1}], F[3,{1}]}  -> {F[3,{1}],F[3,{1}],V[1]}, SelectInterference -> {gQCD -> 4}, SelectCoupling -> (Exponent[#1, eQED] == 1 &)];
(* uuxccxa *) AddReal[FeynArtsProcess -> {F[3,{1}], -F[3,{1}]} -> {F[3,{2}],-F[3,{2}],V[1]}, SelectInterference -> {gQCD -> 4}, SelectCoupling -> (Exponent[#1, eQED] == 1 &)];
(* ddxbbxa *) AddReal[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {F[4,{3}],-F[4,{3}],V[1]}, SelectInterference -> {gQCD -> 4}, SelectCoupling -> (Exponent[#1, eQED] == 1 &)];
(* uuxbbxa *) AddReal[FeynArtsProcess -> {F[3,{1}], -F[3,{1}]} -> {F[4,{3}],-F[4,{3}],V[1]}, SelectInterference -> {gQCD -> 4}, SelectCoupling -> (Exponent[#1, eQED] == 1 &)];
(* bbbxbxa *) AddReal[FeynArtsProcess -> {F[4,{3}], F[4,{3}]}  -> {F[4,{3}],F[4,{3}],V[1]}, SelectInterference -> {gQCD -> 4}, SelectCoupling -> (Exponent[#1, eQED] == 1 &)];

(* dddxdxg *) AddReal[FeynArtsProcess -> {F[4,{1}], F[4,{1}]}  -> {F[4,{1}],F[4,{1}],V[5]}, SelectInterference -> {gQCD -> 4}, SelectCoupling -> (Exponent[#1, eQED] == 0 || Exponent[#1, eQED] == 2 &)];
(* uuxddxg *) AddReal[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {F[3,{1}],-F[3,{1}],V[5]}, SelectInterference -> {gQCD -> 4}, SelectCoupling -> (Exponent[#1, eQED] == 0 || Exponent[#1, eQED] == 2 &)];
(* ddxssxg *) AddReal[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {F[4,{2}],-F[4,{2}],V[5]}, SelectInterference -> {gQCD -> 4}, SelectCoupling -> (Exponent[#1, eQED] == 0 || Exponent[#1, eQED] == 2 &)];
(* ccxddxg *) AddReal[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {F[3,{2}],-F[3,{2}],V[5]}, SelectInterference -> {gQCD -> 4}, SelectCoupling -> (Exponent[#1, eQED] == 0 || Exponent[#1, eQED] == 2 &)];
(* uxcdsxg *) AddReal[FeynArtsProcess -> {F[4,{1}], -F[3,{1}]} -> {F[4,{2}],-F[3,{2}],V[5]}, SelectInterference -> {gQCD -> 4}, SelectCoupling -> (Exponent[#1, eQED] == 0 || Exponent[#1, eQED] == 2 &)];
(* uuuxuxg *) AddReal[FeynArtsProcess -> {F[3,{1}], F[3,{1}]}  -> {F[3,{1}],F[3,{1}],V[5]}, SelectInterference -> {gQCD -> 4}, SelectCoupling -> (Exponent[#1, eQED] == 0 || Exponent[#1, eQED] == 2 &)];
(* uuxccxg *) AddReal[FeynArtsProcess -> {F[3,{1}], -F[3,{1}]} -> {F[3,{2}],-F[3,{2}],V[5]}, SelectInterference -> {gQCD -> 4}, SelectCoupling -> (Exponent[#1, eQED] == 0 || Exponent[#1, eQED] == 2 &)];
(* ddxbbxg *) AddReal[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {F[4,{3}],-F[4,{3}],V[5]}, SelectInterference -> {gQCD -> 4}, SelectCoupling -> (Exponent[#1, eQED] == 0 || Exponent[#1, eQED] == 2 &)];
(* uuxbbxg *) AddReal[FeynArtsProcess -> {F[3,{1}], -F[3,{1}]} -> {F[4,{3}],-F[4,{3}],V[5]}, SelectInterference -> {gQCD -> 4}, SelectCoupling -> (Exponent[#1, eQED] == 0 || Exponent[#1, eQED] == 2 &)];
(* bbbxbxg *) AddReal[FeynArtsProcess -> {F[4,{3}], F[4,{3}]}  -> {F[4,{3}],F[4,{3}],V[5]}, SelectInterference -> {gQCD -> 4}, SelectCoupling -> (Exponent[#1, eQED] == 0 || Exponent[#1, eQED] == 2 &)];

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

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

(* O(as*a) *)

(* dddxdx *) AddTree[FeynArtsProcess -> {F[4,{1}], F[4,{1}]}  -> {F[4,{1}],F[4,{1}]}, SelectInterference -> {eQED -> 2}, SelectCoupling -> (Exponent[#1, eQED] == 0  || Exponent[#1, eQED] == 2 &)];
(* uuxddx *) AddTree[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {F[3,{1}],-F[3,{1}]}, SelectInterference -> {eQED -> 2}, SelectCoupling -> (Exponent[#1, eQED] == 0  || Exponent[#1, eQED] == 2 &)];
(* ddxssx *) AddTree[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {F[4,{2}],-F[4,{2}]}, SelectInterference -> {eQED -> 2}, SelectCoupling -> (Exponent[#1, eQED] == 0  || Exponent[#1, eQED] == 2 &)];
(* ccxddx *) AddTree[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {F[3,{2}],-F[3,{2}]}, SelectInterference -> {eQED -> 2}, SelectCoupling -> (Exponent[#1, eQED] == 0  || Exponent[#1, eQED] == 2 &)];
(* uxcdsx *) AddTree[FeynArtsProcess -> {F[4,{1}], -F[3,{1}]} -> {F[4,{2}],-F[3,{2}]}, SelectInterference -> {eQED -> 2}, SelectCoupling -> (Exponent[#1, eQED] == 0  || Exponent[#1, eQED] == 2 &)];
(* uuuxux *) AddTree[FeynArtsProcess -> {F[3,{1}], F[3,{1}]}  -> {F[3,{1}],F[3,{1}]}, SelectInterference -> {eQED -> 2}, SelectCoupling -> (Exponent[#1, eQED] == 0  || Exponent[#1, eQED] == 2 &)];
(* uuxccx *) AddTree[FeynArtsProcess -> {F[3,{1}], -F[3,{1}]} -> {F[3,{2}],-F[3,{2}]}, SelectInterference -> {eQED -> 2}, SelectCoupling -> (Exponent[#1, eQED] == 0  || Exponent[#1, eQED] == 2 &)];
(* ddxbbx *) AddTree[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {F[4,{3}],-F[4,{3}]}, SelectInterference -> {eQED -> 2}, SelectCoupling -> (Exponent[#1, eQED] == 0  || Exponent[#1, eQED] == 2 &)];
(* uuxbbx *) AddTree[FeynArtsProcess -> {F[3,{1}], -F[3,{1}]} -> {F[4,{3}],-F[4,{3}]}, SelectInterference -> {eQED -> 2}, SelectCoupling -> (Exponent[#1, eQED] == 0  || Exponent[#1, eQED] == 2 &)];
(* bbbxbx *) AddTree[FeynArtsProcess -> {F[4,{3}], F[4,{3}]}  -> {F[4,{3}],F[4,{3}]}, SelectInterference -> {eQED -> 2}, SelectCoupling -> (Exponent[#1, eQED] == 0  || Exponent[#1, eQED] == 2 &)];


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

(* dddxdx *) AddProcess[FeynArtsProcess -> {F[4,{1}], F[4,{1}]}  -> {F[4,{1}],F[4,{1}]}, UnitaryGauge -> False, SelectInterference -> {eQED -> 4}, SelectCoupling -> (Exponent[#1, gQCD] == 2 - 2 * #2 || Exponent[#1, gQCD] == 0 + 2 * #2 &)];
(* uuxddx *) AddProcess[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {F[3,{1}],-F[3,{1}]}, UnitaryGauge -> False, SelectInterference -> {eQED -> 4}, SelectCoupling -> (Exponent[#1, gQCD] == 2 - 2 * #2 || Exponent[#1, gQCD] == 0 + 2 * #2 &)];
(* ddxssx *) AddProcess[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {F[4,{2}],-F[4,{2}]}, UnitaryGauge -> False, SelectInterference -> {eQED -> 4}, SelectCoupling -> (Exponent[#1, gQCD] == 2 - 2 * #2 || Exponent[#1, gQCD] == 0 + 2 * #2 &)];
(* ccxddx *) AddProcess[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {F[3,{2}],-F[3,{2}]}, UnitaryGauge -> False, SelectInterference -> {eQED -> 4}, SelectCoupling -> (Exponent[#1, gQCD] == 2 - 2 * #2 || Exponent[#1, gQCD] == 0 + 2 * #2 &)];
(* uxcdsx *) AddProcess[FeynArtsProcess -> {F[4,{1}], -F[3,{1}]} -> {F[4,{2}],-F[3,{2}]}, UnitaryGauge -> False, SelectInterference -> {eQED -> 4}, SelectCoupling -> (Exponent[#1, gQCD] == 2 - 2 * #2 || Exponent[#1, gQCD] == 0 + 2 * #2 &)];
(* uuuxux *) AddProcess[FeynArtsProcess -> {F[3,{1}], F[3,{1}]}  -> {F[3,{1}],F[3,{1}]}, UnitaryGauge -> False, SelectInterference -> {eQED -> 4}, SelectCoupling -> (Exponent[#1, gQCD] == 2 - 2 * #2 || Exponent[#1, gQCD] == 0 + 2 * #2 &)];
(* uuxccx *) AddProcess[FeynArtsProcess -> {F[3,{1}], -F[3,{1}]} -> {F[3,{2}],-F[3,{2}]}, UnitaryGauge -> False, SelectInterference -> {eQED -> 4}, SelectCoupling -> (Exponent[#1, gQCD] == 2 - 2 * #2 || Exponent[#1, gQCD] == 0 + 2 * #2 &)];
(* ddxbbx *) AddProcess[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {F[4,{3}],-F[4,{3}]}, UnitaryGauge -> False, SelectInterference -> {eQED -> 4}, SelectCoupling -> (Exponent[#1, gQCD] == 2 - 2 * #2 || Exponent[#1, gQCD] == 0 + 2 * #2 &)];
(* uuxbbx *) AddProcess[FeynArtsProcess -> {F[3,{1}], -F[3,{1}]} -> {F[4,{3}],-F[4,{3}]}, UnitaryGauge -> False, SelectInterference -> {eQED -> 4}, SelectCoupling -> (Exponent[#1, gQCD] == 2 - 2 * #2 || Exponent[#1, gQCD] == 0 + 2 * #2 &)];
(* bbbxbx *) AddProcess[FeynArtsProcess -> {F[4,{3}], F[4,{3}]}  -> {F[4,{3}],F[4,{3}]}, UnitaryGauge -> False, SelectInterference -> {eQED -> 4}, SelectCoupling -> (Exponent[#1, gQCD] == 2 - 2 * #2 || Exponent[#1, gQCD] == 0 + 2 * #2 &)];

(* dddxdxa *) AddReal[FeynArtsProcess -> {F[4,{1}], F[4,{1}]}  -> {F[4,{1}],F[4,{1}],V[1]}, SelectInterference -> {eQED -> 4}, SelectCoupling -> (Exponent[#1, gQCD] == 0 || Exponent[#1, gQCD] == 2 &)];
(* uuxddxa *) AddReal[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {F[3,{1}],-F[3,{1}],V[1]}, SelectInterference -> {eQED -> 4}, SelectCoupling -> (Exponent[#1, gQCD] == 0 || Exponent[#1, gQCD] == 2 &)];
(* ddxssxa *) AddReal[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {F[4,{2}],-F[4,{2}],V[1]}, SelectInterference -> {eQED -> 4}, SelectCoupling -> (Exponent[#1, gQCD] == 0 || Exponent[#1, gQCD] == 2 &)];
(* ccxddxa *) AddReal[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {F[3,{2}],-F[3,{2}],V[1]}, SelectInterference -> {eQED -> 4}, SelectCoupling -> (Exponent[#1, gQCD] == 0 || Exponent[#1, gQCD] == 2 &)];
(* uxcdsxa *) AddReal[FeynArtsProcess -> {F[4,{1}], -F[3,{1}]} -> {F[4,{2}],-F[3,{2}],V[1]}, SelectInterference -> {eQED -> 4}, SelectCoupling -> (Exponent[#1, gQCD] == 0 || Exponent[#1, gQCD] == 2 &)];
(* uuuxuxa *) AddReal[FeynArtsProcess -> {F[3,{1}], F[3,{1}]}  -> {F[3,{1}],F[3,{1}],V[1]}, SelectInterference -> {eQED -> 4}, SelectCoupling -> (Exponent[#1, gQCD] == 0 || Exponent[#1, gQCD] == 2 &)];
(* uuxccxa *) AddReal[FeynArtsProcess -> {F[3,{1}], -F[3,{1}]} -> {F[3,{2}],-F[3,{2}],V[1]}, SelectInterference -> {eQED -> 4}, SelectCoupling -> (Exponent[#1, gQCD] == 0 || Exponent[#1, gQCD] == 2 &)];
(* ddxbbxa *) AddReal[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {F[4,{3}],-F[4,{3}],V[1]}, SelectInterference -> {eQED -> 4}, SelectCoupling -> (Exponent[#1, gQCD] == 0 || Exponent[#1, gQCD] == 2 &)];
(* uuxbbxa *) AddReal[FeynArtsProcess -> {F[3,{1}], -F[3,{1}]} -> {F[4,{3}],-F[4,{3}],V[1]}, SelectInterference -> {eQED -> 4}, SelectCoupling -> (Exponent[#1, gQCD] == 0 || Exponent[#1, gQCD] == 2 &)];
(* bbbxbxa *) AddReal[FeynArtsProcess -> {F[4,{3}], F[4,{3}]}  -> {F[4,{3}],F[4,{3}],V[1]}, SelectInterference -> {eQED -> 4}, SelectCoupling -> (Exponent[#1, gQCD] == 0 || Exponent[#1, gQCD] == 2 &)];

(* dddxdxg *) AddReal[FeynArtsProcess -> {F[4,{1}], F[4,{1}]}  -> {F[4,{1}],F[4,{1}],V[5]}, SelectInterference -> {eQED -> 4}, SelectCoupling -> (Exponent[#1, gQCD] == 1 &)];
(* uuxddxg *) AddReal[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {F[3,{1}],-F[3,{1}],V[5]}, SelectInterference -> {eQED -> 4}, SelectCoupling -> (Exponent[#1, gQCD] == 1 &)];
(* ddxssxg *) AddReal[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {F[4,{2}],-F[4,{2}],V[5]}, SelectInterference -> {eQED -> 4}, SelectCoupling -> (Exponent[#1, gQCD] == 1 &)];
(* ccxddxg *) AddReal[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {F[3,{2}],-F[3,{2}],V[5]}, SelectInterference -> {eQED -> 4}, SelectCoupling -> (Exponent[#1, gQCD] == 1 &)];
(* uxcdsxg *) AddReal[FeynArtsProcess -> {F[4,{1}], -F[3,{1}]} -> {F[4,{2}],-F[3,{2}],V[5]}, SelectInterference -> {eQED -> 4}, SelectCoupling -> (Exponent[#1, gQCD] == 1 &)];
(* uuuxuxg *) AddReal[FeynArtsProcess -> {F[3,{1}], F[3,{1}]}  -> {F[3,{1}],F[3,{1}],V[5]}, SelectInterference -> {eQED -> 4}, SelectCoupling -> (Exponent[#1, gQCD] == 1 &)];
(* uuxccxg *) AddReal[FeynArtsProcess -> {F[3,{1}], -F[3,{1}]} -> {F[3,{2}],-F[3,{2}],V[5]}, SelectInterference -> {eQED -> 4}, SelectCoupling -> (Exponent[#1, gQCD] == 1 &)];
(* ddxbbxg *) AddReal[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {F[4,{3}],-F[4,{3}],V[5]}, SelectInterference -> {eQED -> 4}, SelectCoupling -> (Exponent[#1, gQCD] == 1 &)];
(* uuxbbxg *) AddReal[FeynArtsProcess -> {F[3,{1}], -F[3,{1}]} -> {F[4,{3}],-F[4,{3}],V[5]}, SelectInterference -> {eQED -> 4}, SelectCoupling -> (Exponent[#1, gQCD] == 1 &)];
(* bbbxbxg *) AddReal[FeynArtsProcess -> {F[4,{3}], F[4,{3}]}  -> {F[4,{3}],F[4,{3}],V[5]}, SelectInterference -> {eQED -> 4}, SelectCoupling -> (Exponent[#1, gQCD] == 1 &)];

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

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

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


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

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


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

(* dddxdx *) AddProcess[FeynArtsProcess -> {F[4,{1}], F[4,{1}]}  -> {F[4,{1}],F[4,{1}]}, UnitaryGauge -> False, SelectInterference -> {gQCD -> 0}, SelectCoupling -> (Exponent[#1, eQED] == 2 + 2 * #2 &)];
(* uuxddx *) AddProcess[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {F[3,{1}],-F[3,{1}]}, UnitaryGauge -> False, SelectInterference -> {gQCD -> 0}, SelectCoupling -> (Exponent[#1, eQED] == 2 + 2 * #2 &)];
(* ddxssx *) AddProcess[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {F[4,{2}],-F[4,{2}]}, UnitaryGauge -> False, SelectInterference -> {gQCD -> 0}, SelectCoupling -> (Exponent[#1, eQED] == 2 + 2 * #2 &)];
(* ccxddx *) AddProcess[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {F[3,{2}],-F[3,{2}]}, UnitaryGauge -> False, SelectInterference -> {gQCD -> 0}, SelectCoupling -> (Exponent[#1, eQED] == 2 + 2 * #2 &)];
(* uxcdsx *) AddProcess[FeynArtsProcess -> {F[4,{1}], -F[3,{1}]} -> {F[4,{2}],-F[3,{2}]}, UnitaryGauge -> False, SelectInterference -> {gQCD -> 0}, SelectCoupling -> (Exponent[#1, eQED] == 2 + 2 * #2 &)];
(* uuuxux *) AddProcess[FeynArtsProcess -> {F[3,{1}], F[3,{1}]}  -> {F[3,{1}],F[3,{1}]}, UnitaryGauge -> False, SelectInterference -> {gQCD -> 0}, SelectCoupling -> (Exponent[#1, eQED] == 2 + 2 * #2 &)];
(* uuxccx *) AddProcess[FeynArtsProcess -> {F[3,{1}], -F[3,{1}]} -> {F[3,{2}],-F[3,{2}]}, UnitaryGauge -> False, SelectInterference -> {gQCD -> 0}, SelectCoupling -> (Exponent[#1, eQED] == 2 + 2 * #2 &)];
(* ddxbbx *) AddProcess[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {F[4,{3}],-F[4,{3}]}, UnitaryGauge -> False, SelectInterference -> {gQCD -> 0}, SelectCoupling -> (Exponent[#1, eQED] == 2 + 2 * #2 &)];
(* uuxbbx *) AddProcess[FeynArtsProcess -> {F[3,{1}], -F[3,{1}]} -> {F[4,{3}],-F[4,{3}]}, UnitaryGauge -> False, SelectInterference -> {gQCD -> 0}, SelectCoupling -> (Exponent[#1, eQED] == 2 + 2 * #2 &)];
(* bbbxbx *) AddProcess[FeynArtsProcess -> {F[4,{3}], F[4,{3}]}  -> {F[4,{3}],F[4,{3}]}, UnitaryGauge -> False, SelectInterference -> {gQCD -> 0}, SelectCoupling -> (Exponent[#1, eQED] == 2 + 2 * #2 &)];

(* dddxdxa *) AddReal[FeynArtsProcess -> {F[4,{1}], F[4,{1}]}  -> {F[4,{1}],F[4,{1}],V[1]}, SelectInterference -> {gQCD -> 0}, SelectCoupling -> (Exponent[#1, eQED] == 3 &)];
(* uuxddxa *) AddReal[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {F[3,{1}],-F[3,{1}],V[1]}, SelectInterference -> {gQCD -> 0}, SelectCoupling -> (Exponent[#1, eQED] == 3 &)];
(* ddxssxa *) AddReal[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {F[4,{2}],-F[4,{2}],V[1]}, SelectInterference -> {gQCD -> 0}, SelectCoupling -> (Exponent[#1, eQED] == 3 &)];
(* ccxddxa *) AddReal[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {F[3,{2}],-F[3,{2}],V[1]}, SelectInterference -> {gQCD -> 0}, SelectCoupling -> (Exponent[#1, eQED] == 3 &)];
(* uxcdsxa *) AddReal[FeynArtsProcess -> {F[4,{1}], -F[3,{1}]} -> {F[4,{2}],-F[3,{2}],V[1]}, SelectInterference -> {gQCD -> 0}, SelectCoupling -> (Exponent[#1, eQED] == 3 &)];
(* uuuxuxa *) AddReal[FeynArtsProcess -> {F[3,{1}], F[3,{1}]}  -> {F[3,{1}],F[3,{1}],V[1]}, SelectInterference -> {gQCD -> 0}, SelectCoupling -> (Exponent[#1, eQED] == 3 &)];
(* uuxccxa *) AddReal[FeynArtsProcess -> {F[3,{1}], -F[3,{1}]} -> {F[3,{2}],-F[3,{2}],V[1]}, SelectInterference -> {gQCD -> 0}, SelectCoupling -> (Exponent[#1, eQED] == 3 &)];
(* ddxbbxa *) AddReal[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {F[4,{3}],-F[4,{3}],V[1]}, SelectInterference -> {gQCD -> 0}, SelectCoupling -> (Exponent[#1, eQED] == 3 &)];
(* uuxbbxa *) AddReal[FeynArtsProcess -> {F[3,{1}], -F[3,{1}]} -> {F[4,{3}],-F[4,{3}],V[1]}, SelectInterference -> {gQCD -> 0}, SelectCoupling -> (Exponent[#1, eQED] == 3 &)];
(* bbbxbxa *) AddReal[FeynArtsProcess -> {F[4,{3}], F[4,{3}]}  -> {F[4,{3}],F[4,{3}],V[1]}, SelectInterference -> {gQCD -> 0}, SelectCoupling -> (Exponent[#1, eQED] == 3 &)];

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

(* ddxaaa *)  AddReal[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {V[1],V[1],V[1]}, SelectInterference -> {gQCD -> 0}, SelectCoupling -> (Exponent[#1, eQED] == 3 &)];
(* uuxaaa *)  AddReal[FeynArtsProcess -> {F[3,{1}], -F[3,{1}]} -> {V[1],V[1],V[1]}, SelectInterference -> {gQCD -> 0}, SelectCoupling -> (Exponent[#1, eQED] == 3 &)];
(* bbxaaa *)  AddReal[FeynArtsProcess -> {F[4,{3}], -F[4,{3}]} -> {V[1],V[1],V[1]}, SelectInterference -> {gQCD -> 0}, SelectCoupling -> (Exponent[#1, eQED] == 3 &)];






