
(* description = "ee -> VV+jets|ee -> VV+jets|ee -> VV|On-shell boson pair production (Photon, W, Z in all combinations) in e+e- at NLO EW and NLO QED." *)

SelectInterference = {gQCD -> 0};

qed1 = {
  InsertFieldsOptions -> {Restrictions -> {ExcludeParticles -> {V[3],U[3],S[3]}, NoQuarkMixing}},
  SelectLoopDiagrams -> ( Not[ContainsFermionLoop[##] && NLegsOnLoop[2][##]] && ParticlesInLoop[V[1]][##]  &),
  QED -> 1,
  SetParameters -> JoinOptions[{ME -> ME, MM -> MM}]
  };


(* NLO EW *)

(* ddxzz *)  AddProcess[FeynArtsProcess -> { F[2,{1}], -F[2,{1}]} -> { V[2], V[2]}, UnitaryGauge -> False, SetParameters -> JoinOptions[{ME -> 0, MM -> 0}]];
(* ddxwwx *) AddProcess[FeynArtsProcess -> { F[2,{1}], -F[2,{1}]} -> {-V[3], V[3]}, UnitaryGauge -> False, SetParameters -> JoinOptions[{ME -> 0, MM -> 0}]];
(* ddxaz *)  AddProcess[FeynArtsProcess -> { F[2,{1}], -F[2,{1}]} -> { V[1], V[2]}, UnitaryGauge -> False, SetParameters -> JoinOptions[{ME -> 0, MM -> 0}]];
(* ddxaa *)  AddProcess[FeynArtsProcess -> { F[2,{1}], -F[2,{1}]} -> { V[1], V[1]}, UnitaryGauge -> False, SetParameters -> JoinOptions[{ME -> 0, MM -> 0}]];

(* NLO QED *)

(* ddxzz *)  AddProcess[FeynArtsProcess -> { F[2,{1}], -F[2,{1}]} -> { V[2], V[2]}, UnitaryGauge -> False, Sequence @@ qed1];
(* ddxaz *)  AddProcess[FeynArtsProcess -> { F[2,{1}], -F[2,{1}]} -> { V[1], V[2]}, UnitaryGauge -> False, Sequence @@ qed1];
(* ddxaa *)  AddProcess[FeynArtsProcess -> { F[2,{1}], -F[2,{1}]} -> { V[1], V[1]}, UnitaryGauge -> False, Sequence @@ qed1];


(* ddxzza *) AddReal[FeynArtsProcess -> { F[2,{1}], -F[2,{1}]} -> { V[2], V[2], V[1]}];
(* ddxwwxa *) AddReal[FeynArtsProcess -> { F[2,{1}], -F[2,{1}]} -> {-V[3], V[3], V[1]}];
(* ddxaaa *) AddReal[FeynArtsProcess -> { F[2,{1}], -F[2,{1}]} -> { V[1], V[1], V[1]}];
(* ddxaza *) AddReal[FeynArtsProcess -> { F[2,{1}], -F[2,{1}]} -> { V[1], V[2], V[1]}];

