
(* description = "H -> X|pp -> leptons|H -> llll|Higgs decay into four fermions (all combinations) at NLO EW." *)

SelectInterference = {gQCD -> 0, DOI -> 0};
SelectCoupling = (Exponent[#1, eQED] == 3 + 2 * #2 &);
OnTheFlyMode = 0;
UnitaryGauge = False;

nf5 = {
	OpenLoopsModel -> "SM",
	SetParameters -> JoinOptions[{ME -> 0, MB -> 0, nf -> 5}],
	InsertFieldsOptions -> {Restrictions -> {ExcludeParticles -> {F[3,{3}]}, NoQuarkMixing}}
};


(* eeexexh     *) AddProcess[FeynArtsProcess -> {S[1]} -> {F[2,{1}],  F[2,{1}], -F[2,{1}], -F[2,{1}]}, SetParameters -> JoinOptions[{ME -> 0, MB -> 0, YE -> 0, YB -> 0}]];


(* eeexexh     *) AddProcess[FeynArtsProcess -> {S[1]} -> {F[2,{1}],  -F[2,{1}], F[2,{2}], -F[2,{2}]}, SetParameters -> JoinOptions[{MB -> 0, YB -> 0, ME -> 0, YE -> 0, MM -> 0, YM -> 0}]];

(* Z Z, opposite flavour leptons + neutrinos *)
(* nmnmxeexh   *) AddProcess[FeynArtsProcess -> {S[1]} -> {F[2,{1}],  F[1,{2}], -F[2,{1}], -F[1,{2}]}, SetParameters -> JoinOptions[{MB -> 0, YB -> 0, ME -> 0, YE -> 0, MM -> 0, YM -> 0}]];

(* W+ W- / Z Z, same flavour *)
(* nenexeexh   *) AddProcess[FeynArtsProcess -> {S[1]} -> {F[2,{1}],  F[1,{1}], -F[2,{1}], -F[1,{1}]}, SetParameters -> JoinOptions[{MB -> 0, YB -> 0, ME -> 0, YE -> 0, MM -> 0, YM -> 0}]];

(* W+ W-, opposite flavour *)
(* nexnmemxh   *) AddProcess[FeynArtsProcess -> {S[1]} -> {F[2,{2}],  F[1,{1}], -F[2,{1}], -F[1,{2}]}, SetParameters -> JoinOptions[{MB -> 0, YB -> 0, ME -> 0, YE -> 0, MM -> 0, YM -> 0}]];
