(* description = "pp -> V+jets|pp -> V+jets|pp -> lnjj|Production of an on-shell Z-boson/Photon and two jets." *)

SetParameters = JoinOptions[{ME -> 0}];

(* pp -> z jj *)

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

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

(* dddxdxz *) AddProcess[FeynArtsProcess -> {F[4,{1}], F[4,{1}]}  -> {V[2], F[4,{1}],F[4,{1}]}, UnitaryGauge -> False, SelectInterference -> {gQCD -> 4}, SelectCoupling -> (Exponent[#1, eQED] == 3 - 2 * #2 || Exponent[#1, eQED] == 1 + 2 * #2 &)];
(* uuxddxz *) AddProcess[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {V[2], F[3,{1}],-F[3,{1}]}, UnitaryGauge -> False, SelectInterference -> {gQCD -> 4}, SelectCoupling -> (Exponent[#1, eQED] == 3 - 2 * #2 || Exponent[#1, eQED] == 1 + 2 * #2 &)];
(* ddxssxz *) AddProcess[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {V[2], F[4,{2}],-F[4,{2}]}, UnitaryGauge -> False, SelectInterference -> {gQCD -> 4}, SelectCoupling -> (Exponent[#1, eQED] == 3 - 2 * #2 || Exponent[#1, eQED] == 1 + 2 * #2 &)];
(* ccxddxz *) AddProcess[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {V[2], F[3,{2}],-F[3,{2}]}, UnitaryGauge -> False, SelectInterference -> {gQCD -> 4}, SelectCoupling -> (Exponent[#1, eQED] == 3 - 2 * #2 || Exponent[#1, eQED] == 1 + 2 * #2 &)];
(* uxcdsxz *) AddProcess[FeynArtsProcess -> {F[4,{1}], -F[3,{1}]} -> {V[2], F[4,{2}],-F[3,{2}]}, UnitaryGauge -> False, SelectInterference -> {gQCD -> 4}, SelectCoupling -> (Exponent[#1, eQED] == 3 - 2 * #2 || Exponent[#1, eQED] == 1 + 2 * #2 &)];
(* uuuxuxz *) AddProcess[FeynArtsProcess -> {F[3,{1}], F[3,{1}]}  -> {V[2], F[3,{1}],F[3,{1}]}, UnitaryGauge -> False, SelectInterference -> {gQCD -> 4}, SelectCoupling -> (Exponent[#1, eQED] == 3 - 2 * #2 || Exponent[#1, eQED] == 1 + 2 * #2 &)];
(* uuxccxz *) AddProcess[FeynArtsProcess -> {F[3,{1}], -F[3,{1}]} -> {V[2], F[3,{2}],-F[3,{2}]}, UnitaryGauge -> False, SelectInterference -> {gQCD -> 4}, SelectCoupling -> (Exponent[#1, eQED] == 3 - 2 * #2 || Exponent[#1, eQED] == 1 + 2 * #2 &)];
(* ddxbbxz *) AddProcess[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {V[2], F[4,{3}],-F[4,{3}]}, UnitaryGauge -> False, SelectInterference -> {gQCD -> 4}, SelectCoupling -> (Exponent[#1, eQED] == 3 - 2 * #2 || Exponent[#1, eQED] == 1 + 2 * #2 &)];
(* uuxbbxz *) AddProcess[FeynArtsProcess -> {F[3,{1}], -F[3,{1}]} -> {V[2], F[4,{3}],-F[4,{3}]}, UnitaryGauge -> False, SelectInterference -> {gQCD -> 4}, SelectCoupling -> (Exponent[#1, eQED] == 3 - 2 * #2 || Exponent[#1, eQED] == 1 + 2 * #2 &)];
(* bbbxbxz *) AddProcess[FeynArtsProcess -> {F[4,{3}], F[4,{3}]}  -> {V[2], F[4,{3}],F[4,{3}]}, UnitaryGauge -> False, SelectInterference -> {gQCD -> 4}, SelectCoupling -> (Exponent[#1, eQED] == 3 - 2 * #2 || Exponent[#1, eQED] == 1 + 2 * #2 &)];

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

(* dddxdxaz *) AddReal[FeynArtsProcess -> {F[4,{1}], F[4,{1}]}  -> {V[2], F[4,{1}],F[4,{1}],V[1]}, SelectInterference -> {gQCD -> 4}, SelectCoupling -> (Exponent[#1, eQED] == 2 &)];
(* uuxddxaz *) AddReal[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {V[2], F[3,{1}],-F[3,{1}],V[1]}, SelectInterference -> {gQCD -> 4}, SelectCoupling -> (Exponent[#1, eQED] == 2 &)];
(* ddxssxaz *) AddReal[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {V[2], F[4,{2}],-F[4,{2}],V[1]}, SelectInterference -> {gQCD -> 4}, SelectCoupling -> (Exponent[#1, eQED] == 2 &)];
(* ccxddxaz *) AddReal[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {V[2], F[3,{2}],-F[3,{2}],V[1]}, SelectInterference -> {gQCD -> 4}, SelectCoupling -> (Exponent[#1, eQED] == 2 &)];
(* uuuxuxaz *) AddReal[FeynArtsProcess -> {F[3,{1}], F[3,{1}]}  -> {V[2], F[3,{1}],F[3,{1}],V[1]}, SelectInterference -> {gQCD -> 4}, SelectCoupling -> (Exponent[#1, eQED] == 2 &)];
(* uuxccxaz *) AddReal[FeynArtsProcess -> {F[3,{1}], -F[3,{1}]} -> {V[2], F[3,{2}],-F[3,{2}],V[1]}, SelectInterference -> {gQCD -> 4}, SelectCoupling -> (Exponent[#1, eQED] == 2 &)];
(* ddxbbxaz *) AddReal[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {V[2], F[4,{3}],-F[4,{3}],V[1]}, SelectInterference -> {gQCD -> 4}, SelectCoupling -> (Exponent[#1, eQED] == 2 &)];
(* uuxbbxaz *) AddReal[FeynArtsProcess -> {F[3,{1}], -F[3,{1}]} -> {V[2], F[4,{3}],-F[4,{3}],V[1]}, SelectInterference -> {gQCD -> 4}, SelectCoupling -> (Exponent[#1, eQED] == 2 &)];
(* bbbxbxaz *) AddReal[FeynArtsProcess -> {F[4,{3}], F[4,{3}]}  -> {V[2], F[4,{3}],F[4,{3}],V[1]}, SelectInterference -> {gQCD -> 4}, SelectCoupling -> (Exponent[#1, eQED] == 2 &)];

(* dddxdxzg *) AddReal[FeynArtsProcess -> {F[4,{1}], F[4,{1}]}  -> {V[2], F[4,{1}],F[4,{1}],V[5]}, SelectInterference -> {gQCD -> 4}, SelectCoupling -> (Exponent[#1, eQED] == 1 || Exponent[#1, eQED] == 3 &)];
(* uuxddxzg *) AddReal[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {V[2], F[3,{1}],-F[3,{1}],V[5]}, SelectInterference -> {gQCD -> 4}, SelectCoupling -> (Exponent[#1, eQED] == 1 || Exponent[#1, eQED] == 3 &)];
(* ddxssxzg *) AddReal[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {V[2], F[4,{2}],-F[4,{2}],V[5]}, SelectInterference -> {gQCD -> 4}, SelectCoupling -> (Exponent[#1, eQED] == 1 || Exponent[#1, eQED] == 3 &)];
(* ccxddxzg *) AddReal[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {V[2], F[3,{2}],-F[3,{2}],V[5]}, SelectInterference -> {gQCD -> 4}, SelectCoupling -> (Exponent[#1, eQED] == 1 || Exponent[#1, eQED] == 3 &)];
(* uxcdsxzg *) AddReal[FeynArtsProcess -> {F[4,{1}], -F[3,{1}]} -> {V[2], F[4,{2}],-F[3,{2}],V[5]}, SelectInterference -> {gQCD -> 4}, SelectCoupling -> (Exponent[#1, eQED] == 1 || Exponent[#1, eQED] == 3 &)];
(* uuuxuxzg *) AddReal[FeynArtsProcess -> {F[3,{1}], F[3,{1}]}  -> {V[2], F[3,{1}],F[3,{1}],V[5]}, SelectInterference -> {gQCD -> 4}, SelectCoupling -> (Exponent[#1, eQED] == 1 || Exponent[#1, eQED] == 3 &)];
(* uuxccxzg *) AddReal[FeynArtsProcess -> {F[3,{1}], -F[3,{1}]} -> {V[2], F[3,{2}],-F[3,{2}],V[5]}, SelectInterference -> {gQCD -> 4}, SelectCoupling -> (Exponent[#1, eQED] == 1 || Exponent[#1, eQED] == 3 &)];
(* ddxbbxzg *) AddReal[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {V[2], F[4,{3}],-F[4,{3}],V[5]}, SelectInterference -> {gQCD -> 4}, SelectCoupling -> (Exponent[#1, eQED] == 1 || Exponent[#1, eQED] == 3 &)];
(* uuxbbxzg *) AddReal[FeynArtsProcess -> {F[3,{1}], -F[3,{1}]} -> {V[2], F[4,{3}],-F[4,{3}],V[5]}, SelectInterference -> {gQCD -> 4}, SelectCoupling -> (Exponent[#1, eQED] == 1 || Exponent[#1, eQED] == 3 &)];
(* bbbxbxzg *) AddReal[FeynArtsProcess -> {F[4,{3}], F[4,{3}]}  -> {V[2], F[4,{3}],F[4,{3}],V[5]}, SelectInterference -> {gQCD -> 4}, SelectCoupling -> (Exponent[#1, eQED] == 1 || Exponent[#1, eQED] == 3 &)];


(* pp -> a jj *)

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

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

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

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

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

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


