

SelectInterference = {qQCD -> 0};

Approximation = "lew";
UnitaryGauge = False;

(* ew corrections, no tau and nu_tau, no quarks in loop, no photon self-energies *)

lepew = {
  InsertFieldsOptions -> {Restrictions -> {ExcludeParticles -> {V[4 | 5],U[5],S[1 | 2 | 3], F[1, {3}], F[2, {3}], F[4, __], F[3, __]}}},
  SelectLoopDiagrams -> ( Not[ContainsFermionLoop[##] && NLegsOnLoop[2][##] && ParticlesOnLoop[V[1], 2][##]]  &), 
  QED -> 0,
  GenericModel -> "Lorentz"
};

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

(* mumuxeex *)    AddProcess[FeynArtsProcess -> {F[2,{2}], -F[2,{2}]} -> {F[2,{1}], -F[2,{1}]}, Sequence @@ lepew];
(* mumuxeexa *)   AddProcess[FeynArtsProcess -> {F[2,{2}], -F[2,{2}]} -> {F[2,{1}], -F[2,{1}], V[1]}, Sequence @@ lepew];
(* mumuxeexaa *)  AddProcess[FeynArtsProcess -> {F[2,{2}], -F[2,{2}]} -> {F[2,{1}], -F[2,{1}], V[1], V[1]}, Sequence @@ lepew];
