
(* description = "ee -> HV+jets|ee -> HV+jets|ee -> HV|Higgs plus on-shell Z-boson production in e+e-/m+m- 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 *)

(* eexhz *)  AddProcess[FeynArtsProcess -> { F[2,{1}], -F[2,{1}]} -> { S[1], V[2]}, UnitaryGauge -> False, SetParameters -> JoinOptions[{ME -> 0}]];
(* eexhz *)  AddProcess[FeynArtsProcess -> { F[2,{1}], -F[2,{1}]} -> { S[1], V[2]}, UnitaryGauge -> False, SetParameters -> JoinOptions[{ME -> ME, MM -> MM}]];
(* mmxhz *)  AddProcess[FeynArtsProcess -> { F[2,{2}], -F[2,{2}]} -> { S[1], V[2]}, UnitaryGauge -> False, SetParameters -> JoinOptions[{ME -> ME, MM -> MM}]];

(* SM + NLO QED *)

(* eexhz *)  AddProcess[FeynArtsProcess -> { F[2,{1}], -F[2,{1}]} -> { S[1], V[2]}, UnitaryGauge -> False, Sequence @@ qed1];
(* mmxhz *)  AddProcess[FeynArtsProcess -> { F[2,{2}], -F[2,{2}]} -> { S[1], V[2]}, UnitaryGauge -> False, Sequence @@ qed1];


(* eexhza *) AddReal[FeynArtsProcess -> { F[2,{1}], -F[2,{1}]} -> { S[1], V[2], V[1]}];
(* eexhza *) AddReal[FeynArtsProcess -> { F[2,{2}], -F[2,{2}]} -> { S[1], V[2], V[1]}];

(* eexhza *) AddReal[FeynArtsProcess -> { F[2,{1}], -F[2,{1}]} -> { S[1], V[2], V[1]}, SetParameters -> JoinOptions[{ME -> ME, MM -> MM}]];
(* mmxhza *) AddReal[FeynArtsProcess -> { F[2,{2}], -F[2,{2}]} -> { S[1], V[2], V[1]}, SetParameters -> JoinOptions[{ME -> ME, MM -> MM}]];

