
(* description = "pp -> VV+jets|pp -> llvv+jets|pp -> llvvjj| W+W+ and W-W- vector boson scattering with different flavour leptonic decays at all orders." *)


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


(* as^3*a^4 contribution *)

(* W+W+ + uc channel *)
AddProcess[FeynArtsProcess -> {F[3,{1}], F[3,{2}]} -> {F[1,{1}], -F[2,{1}], F[1,{2}], -F[2,{2}], F[4,{1}], F[4,{2}]},
           SelectInterference -> {eQED -> 8},
           SelectCoupling -> (Exponent[#1, gQCD] === 2 + 2*#2 &),
           UnitaryGauge -> True];

(* W+W+ + uu channel *)
AddProcess[FeynArtsProcess -> {F[3,{1}], F[3,{1}]} -> {F[1,{1}], -F[2,{1}], F[1,{2}], -F[2,{2}], F[4,{1}], F[4,{1}]},
           SelectInterference -> {eQED -> 8},
           SelectCoupling -> (Exponent[#1, gQCD] === 2 + 2*#2 &),
           UnitaryGauge -> True];



(* as^2*a^5 contribution *)

(* W+W+ + uc channel *)
AddProcess[FeynArtsProcess -> {F[3,{1}], F[3,{2}]} -> {F[1,{1}], -F[2,{1}], F[1,{2}], -F[2,{2}], F[4,{1}], F[4,{2}]},
           SelectInterference -> {gQCD -> 4},
           SelectCoupling -> (Exponent[#1, eQED] <= 6 &),
           UnitaryGauge -> False];

(* W+W+ + uu channel *)
AddProcess[FeynArtsProcess -> {F[3,{1}], F[3,{1}]} -> {F[1,{1}], -F[2,{1}], F[1,{2}], -F[2,{2}], F[4,{1}], F[4,{1}]},
           SelectInterference -> {gQCD -> 4},
           SelectCoupling -> (Exponent[#1, eQED] <= 6 &),
           UnitaryGauge -> False];


(* as^1*a^6 contribution *)

(* W+W+ + uc channel *)
AddProcess[FeynArtsProcess -> {F[3,{1}], F[3,{2}]} -> {F[1,{1}], -F[2,{1}], F[1,{2}], -F[2,{2}], F[4,{1}], F[4,{2}]},
           SelectInterference -> {eQED -> 12},
           SelectCoupling -> (Exponent[#1, eQED] === 6 &),
           UnitaryGauge -> False];

(* W+W+ + uu channel *)
AddProcess[FeynArtsProcess -> {F[3,{1}], F[3,{1}]} -> {F[1,{1}], -F[2,{1}], F[1,{2}], -F[2,{2}], F[4,{1}], F[4,{1}]},
           SelectInterference -> {eQED -> 12},
           SelectCoupling -> (Exponent[#1, gQCD] <= 2 &),
           UnitaryGauge -> False];
(*
(* W-W- + uc channel *)
AddProcess[FeynArtsProcess -> {-F[3,{1}], -F[3,{2}]} -> {-F[1,{1}], F[2,{1}], -F[1,{2}], F[2,{2}], -F[4,{1}], -F[4,{2}]},
           SelectInterference -> {eQED -> 12},
           SelectCoupling -> (Exponent[#1, eQED] === 6 &),
           UnitaryGauge -> False];

(* W-W- + uu channel *)
AddProcess[FeynArtsProcess -> {-F[3,{1}], -F[3,{1}]} -> {-F[1,{1}], F[2,{1}], -F[1,{2}], F[2,{2}], -F[4,{1}], -F[4,{1}]},
           SelectInterference -> {eQED -> 12},
           SelectCoupling -> (Exponent[#1, gQCD] <= 2 &),
           UnitaryGauge -> False];
*)


(* a^7 contribution *)

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

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

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

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