(* description = "pp -> VV+jets|pp -> VV+jets|pp -> aajj|Photon pair plus two jets production at NLO EW (including all subleading one-loop orders)." *)

OTFColourCorrelations = False;

(* pp -> aajj *)

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

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

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

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

(* dddxdxaa *) AddReal[FeynArtsProcess -> {F[4,{1}], F[4,{1}]}  -> {V[1], V[1], F[4,{1}],F[4,{1}],V[1]}, SelectInterference -> {gQCD -> 4}, SelectCoupling -> (Exponent[#1, eQED] == 3 &)];
(* uuxddxaa *) AddReal[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {V[1], V[1], F[3,{1}],-F[3,{1}],V[1]}, SelectInterference -> {gQCD -> 4}, SelectCoupling -> (Exponent[#1, eQED] == 3 &)];
(* ddxssxaa *) AddReal[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {V[1], V[1], F[4,{2}],-F[4,{2}],V[1]}, SelectInterference -> {gQCD -> 4}, SelectCoupling -> (Exponent[#1, eQED] == 3 &)];
(* ccxddxaa *) AddReal[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {V[1], V[1], F[3,{2}],-F[3,{2}],V[1]}, SelectInterference -> {gQCD -> 4}, SelectCoupling -> (Exponent[#1, eQED] == 3 &)];
(* uuuxuxaa *) AddReal[FeynArtsProcess -> {F[3,{1}], F[3,{1}]}  -> {V[1], V[1], F[3,{1}],F[3,{1}],V[1]}, SelectInterference -> {gQCD -> 4}, SelectCoupling -> (Exponent[#1, eQED] == 3 &)];
(* uuxccxaa *) AddReal[FeynArtsProcess -> {F[3,{1}], -F[3,{1}]} -> {V[1], V[1], F[3,{2}],-F[3,{2}],V[1]}, SelectInterference -> {gQCD -> 4}, SelectCoupling -> (Exponent[#1, eQED] == 3 &)];
(* ddxbbxaa *) AddReal[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {V[1], V[1], F[4,{3}],-F[4,{3}],V[1]}, SelectInterference -> {gQCD -> 4}, SelectCoupling -> (Exponent[#1, eQED] == 3 &)];
(* uuxbbxaa *) AddReal[FeynArtsProcess -> {F[3,{1}], -F[3,{1}]} -> {V[1], V[1], F[4,{3}],-F[4,{3}],V[1]}, SelectInterference -> {gQCD -> 4}, SelectCoupling -> (Exponent[#1, eQED] == 3 &)];
(* bbbxbxaa *) AddReal[FeynArtsProcess -> {F[4,{3}], F[4,{3}]}  -> {V[1], V[1], F[4,{3}],F[4,{3}],V[1]}, SelectInterference -> {gQCD -> 4}, SelectCoupling -> (Exponent[#1, eQED] == 3 &)];

(* dddxdxag *) AddReal[FeynArtsProcess -> {F[4,{1}], F[4,{1}]}  -> {V[1], V[1], F[4,{1}],F[4,{1}],V[5]}, SelectInterference -> {gQCD -> 4}, SelectCoupling -> (Exponent[#1, eQED] == 2 || Exponent[#1, eQED] == 4 &)];
(* uuxddxag *) AddReal[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {V[1], V[1], F[3,{1}],-F[3,{1}],V[5]}, SelectInterference -> {gQCD -> 4}, SelectCoupling -> (Exponent[#1, eQED] == 2 || Exponent[#1, eQED] == 4 &)];
(* ddxssxag *) AddReal[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {V[1], V[1], F[4,{2}],-F[4,{2}],V[5]}, SelectInterference -> {gQCD -> 4}, SelectCoupling -> (Exponent[#1, eQED] == 2 || Exponent[#1, eQED] == 4 &)];
(* ccxddxag *) AddReal[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {V[1], V[1], F[3,{2}],-F[3,{2}],V[5]}, SelectInterference -> {gQCD -> 4}, SelectCoupling -> (Exponent[#1, eQED] == 2 || Exponent[#1, eQED] == 4 &)];
(* uxcdsxag *) AddReal[FeynArtsProcess -> {F[4,{1}], -F[3,{1}]} -> {V[1], V[1], F[4,{2}],-F[3,{2}],V[5]}, SelectInterference -> {gQCD -> 4}, SelectCoupling -> (Exponent[#1, eQED] == 2 || Exponent[#1, eQED] == 4 &)];
(* uuuxuxag *) AddReal[FeynArtsProcess -> {F[3,{1}], F[3,{1}]}  -> {V[1], V[1], F[3,{1}],F[3,{1}],V[5]}, SelectInterference -> {gQCD -> 4}, SelectCoupling -> (Exponent[#1, eQED] == 2 || Exponent[#1, eQED] == 4 &)];
(* uuxccxag *) AddReal[FeynArtsProcess -> {F[3,{1}], -F[3,{1}]} -> {V[1], V[1], F[3,{2}],-F[3,{2}],V[5]}, SelectInterference -> {gQCD -> 4}, SelectCoupling -> (Exponent[#1, eQED] == 2 || Exponent[#1, eQED] == 4 &)];
(* ddxbbxag *) AddReal[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {V[1], V[1], F[4,{3}],-F[4,{3}],V[5]}, SelectInterference -> {gQCD -> 4}, SelectCoupling -> (Exponent[#1, eQED] == 2 || Exponent[#1, eQED] == 4 &)];
(* uuxbbxag *) AddReal[FeynArtsProcess -> {F[3,{1}], -F[3,{1}]} -> {V[1], V[1], F[4,{3}],-F[4,{3}],V[5]}, SelectInterference -> {gQCD -> 4}, SelectCoupling -> (Exponent[#1, eQED] == 2 || Exponent[#1, eQED] == 4 &)];
(* bbbxbxag *) AddReal[FeynArtsProcess -> {F[4,{3}], F[4,{3}]}  -> {V[1], V[1], F[4,{3}],F[4,{3}],V[5]}, SelectInterference -> {gQCD -> 4}, SelectCoupling -> (Exponent[#1, eQED] == 2 || Exponent[#1, eQED] == 4 &)];


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

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

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

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


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

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

(* dddxdxaa *) AddReal[FeynArtsProcess -> {F[4,{1}], F[4,{1}]}  -> {V[1], V[1], F[4,{1}],F[4,{1}],V[1]}, SelectInterference -> {eQED -> 8}, SelectCoupling -> (Exponent[#1, gQCD] == 0 || Exponent[#1, gQCD] == 2 &)];
(* uuxddxaa *) AddReal[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {V[1], V[1], F[3,{1}],-F[3,{1}],V[1]}, SelectInterference -> {eQED -> 8}, SelectCoupling -> (Exponent[#1, gQCD] == 0 || Exponent[#1, gQCD] == 2 &)];
(* ddxssxaa *) AddReal[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {V[1], V[1], F[4,{2}],-F[4,{2}],V[1]}, SelectInterference -> {eQED -> 8}, SelectCoupling -> (Exponent[#1, gQCD] == 0 || Exponent[#1, gQCD] == 2 &)];
(* ccxddxaa *) AddReal[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {V[1], V[1], F[3,{2}],-F[3,{2}],V[1]}, SelectInterference -> {eQED -> 8}, SelectCoupling -> (Exponent[#1, gQCD] == 0 || Exponent[#1, gQCD] == 2 &)];
(* uxcdsxaa *) AddReal[FeynArtsProcess -> {F[4,{1}], -F[3,{1}]} -> {V[1], V[1], F[4,{2}],-F[3,{2}],V[1]}, SelectInterference -> {eQED -> 8}, SelectCoupling -> (Exponent[#1, gQCD] == 0 || Exponent[#1, gQCD] == 2 &)];
(* uuuxuxaa *) AddReal[FeynArtsProcess -> {F[3,{1}], F[3,{1}]}  -> {V[1], V[1], F[3,{1}],F[3,{1}],V[1]}, SelectInterference -> {eQED -> 8}, SelectCoupling -> (Exponent[#1, gQCD] == 0 || Exponent[#1, gQCD] == 2 &)];
(* uuxccxaa *) AddReal[FeynArtsProcess -> {F[3,{1}], -F[3,{1}]} -> {V[1], V[1], F[3,{2}],-F[3,{2}],V[1]}, SelectInterference -> {eQED -> 8}, SelectCoupling -> (Exponent[#1, gQCD] == 0 || Exponent[#1, gQCD] == 2 &)];
(* ddxbbxaa *) AddReal[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {V[1], V[1], F[4,{3}],-F[4,{3}],V[1]}, SelectInterference -> {eQED -> 8}, SelectCoupling -> (Exponent[#1, gQCD] == 0 || Exponent[#1, gQCD] == 2 &)];
(* uuxbbxaa *) AddReal[FeynArtsProcess -> {F[3,{1}], -F[3,{1}]} -> {V[1], V[1], F[4,{3}],-F[4,{3}],V[1]}, SelectInterference -> {eQED -> 8}, SelectCoupling -> (Exponent[#1, gQCD] == 0 || Exponent[#1, gQCD] == 2 &)];
(* bbbxbxaa *) AddReal[FeynArtsProcess -> {F[4,{3}], F[4,{3}]}  -> {V[1], V[1], F[4,{3}],F[4,{3}],V[1]}, SelectInterference -> {eQED -> 8}, SelectCoupling -> (Exponent[#1, gQCD] == 0 || Exponent[#1, gQCD] == 2 &)];

(* dddxdxag *) AddReal[FeynArtsProcess -> {F[4,{1}], F[4,{1}]}  -> {V[1], V[1], F[4,{1}],F[4,{1}],V[5]}, SelectInterference -> {eQED -> 8}, SelectCoupling -> (Exponent[#1, gQCD] == 1 &)];
(* uuxddxag *) AddReal[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {V[1], V[1], F[3,{1}],-F[3,{1}],V[5]}, SelectInterference -> {eQED -> 8}, SelectCoupling -> (Exponent[#1, gQCD] == 1 &)];
(* ddxssxag *) AddReal[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {V[1], V[1], F[4,{2}],-F[4,{2}],V[5]}, SelectInterference -> {eQED -> 8}, SelectCoupling -> (Exponent[#1, gQCD] == 1 &)];
(* ccxddxag *) AddReal[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {V[1], V[1], F[3,{2}],-F[3,{2}],V[5]}, SelectInterference -> {eQED -> 8}, SelectCoupling -> (Exponent[#1, gQCD] == 1 &)];
(* uxcdsxag *) AddReal[FeynArtsProcess -> {F[4,{1}], -F[3,{1}]} -> {V[1], V[1], F[4,{2}],-F[3,{2}],V[5]}, SelectInterference -> {eQED -> 8}, SelectCoupling -> (Exponent[#1, gQCD] == 1 &)];
(* uuuxuxag *) AddReal[FeynArtsProcess -> {F[3,{1}], F[3,{1}]}  -> {V[1], V[1], F[3,{1}],F[3,{1}],V[5]}, SelectInterference -> {eQED -> 8}, SelectCoupling -> (Exponent[#1, gQCD] == 1 &)];
(* uuxccxag *) AddReal[FeynArtsProcess -> {F[3,{1}], -F[3,{1}]} -> {V[1], V[1], F[3,{2}],-F[3,{2}],V[5]}, SelectInterference -> {eQED -> 8}, SelectCoupling -> (Exponent[#1, gQCD] == 1 &)];
(* ddxbbxag *) AddReal[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {V[1], V[1], F[4,{3}],-F[4,{3}],V[5]}, SelectInterference -> {eQED -> 8}, SelectCoupling -> (Exponent[#1, gQCD] == 1 &)];
(* uuxbbxag *) AddReal[FeynArtsProcess -> {F[3,{1}], -F[3,{1}]} -> {V[1], V[1], F[4,{3}],-F[4,{3}],V[5]}, SelectInterference -> {eQED -> 8}, SelectCoupling -> (Exponent[#1, gQCD] == 1 &)];
(* bbbxbxag *) AddReal[FeynArtsProcess -> {F[4,{3}], F[4,{3}]}  -> {V[1], V[1], F[4,{3}],F[4,{3}],V[5]}, SelectInterference -> {eQED -> 8}, SelectCoupling -> (Exponent[#1, gQCD] == 1 &)];

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

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

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


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

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


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

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

(* dddxdxaa *) AddReal[FeynArtsProcess -> {F[4,{1}], F[4,{1}]}  -> {V[1], V[1], F[4,{1}],F[4,{1}],V[1]}, SelectInterference -> {gQCD -> 0}, SelectCoupling -> (Exponent[#1, eQED] == 5 &)];
(* uuxddxaa *) AddReal[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {V[1], V[1], F[3,{1}],-F[3,{1}],V[1]}, SelectInterference -> {gQCD -> 0}, SelectCoupling -> (Exponent[#1, eQED] == 5 &)];
(* ddxssxaa *) AddReal[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {V[1], V[1], F[4,{2}],-F[4,{2}],V[1]}, SelectInterference -> {gQCD -> 0}, SelectCoupling -> (Exponent[#1, eQED] == 5 &)];
(* ccxddxaa *) AddReal[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {V[1], V[1], F[3,{2}],-F[3,{2}],V[1]}, SelectInterference -> {gQCD -> 0}, SelectCoupling -> (Exponent[#1, eQED] == 5 &)];
(* uxcdsxaa *) AddReal[FeynArtsProcess -> {F[4,{1}], -F[3,{1}]} -> {V[1], V[1], F[4,{2}],-F[3,{2}],V[1]}, SelectInterference -> {gQCD -> 0}, SelectCoupling -> (Exponent[#1, eQED] == 5 &)];
(* uuuxuxaa *) AddReal[FeynArtsProcess -> {F[3,{1}], F[3,{1}]}  -> {V[1], V[1], F[3,{1}],F[3,{1}],V[1]}, SelectInterference -> {gQCD -> 0}, SelectCoupling -> (Exponent[#1, eQED] == 5 &)];
(* uuxccxaa *) AddReal[FeynArtsProcess -> {F[3,{1}], -F[3,{1}]} -> {V[1], V[1], F[3,{2}],-F[3,{2}],V[1]}, SelectInterference -> {gQCD -> 0}, SelectCoupling -> (Exponent[#1, eQED] == 5 &)];
(* ddxbbxaa *) AddReal[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {V[1], V[1], F[4,{3}],-F[4,{3}],V[1]}, SelectInterference -> {gQCD -> 0}, SelectCoupling -> (Exponent[#1, eQED] == 5 &)];
(* uuxbbxaa *) AddReal[FeynArtsProcess -> {F[3,{1}], -F[3,{1}]} -> {V[1], V[1], F[4,{3}],-F[4,{3}],V[1]}, SelectInterference -> {gQCD -> 0}, SelectCoupling -> (Exponent[#1, eQED] == 5 &)];
(* bbbxbxaa *) AddReal[FeynArtsProcess -> {F[4,{3}], F[4,{3}]}  -> {V[1], V[1], F[4,{3}],F[4,{3}],V[1]}, SelectInterference -> {gQCD -> 0}, SelectCoupling -> (Exponent[#1, eQED] == 5 &)];

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


(* ddxaaaa *)  AddReal[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {V[1], V[1], V[1],V[1],V[1]}, SelectInterference -> {gQCD -> 0}, SelectCoupling -> (Exponent[#1, eQED] == 5 &)];
(* uuxaaaa *)  AddReal[FeynArtsProcess -> {F[3,{1}], -F[3,{1}]} -> {V[1], V[1], V[1],V[1],V[1]}, SelectInterference -> {gQCD -> 0}, SelectCoupling -> (Exponent[#1, eQED] == 5 &)];
(* bbxaaaa *)  AddReal[FeynArtsProcess -> {F[4,{3}], -F[4,{3}]} -> {V[1], V[1], V[1],V[1],V[1]}, SelectInterference -> {gQCD -> 0}, SelectCoupling -> (Exponent[#1, eQED] == 5 &)];





