
(* description = "pp -> VV+jets|pp -> llll+jets|gg -> llll|Approximation 'onlyh'; Loop-induced production of four leptons (all combinations of leptons and neutrinos) in gluon fusion: only diagrams which contain a Higgs Yukawa coupling" *)

SelectCoupling = (Exponent[#1, gQCD] == 2 && Exponent[#1, Yuk] == 1 &);
SelectInterference = {eQED -> 8, Yuk -> 2};
CreateTopologiesOptions = {ExcludeTopologies -> {Snails, WFCorrectionCTs, TadpoleCTs, Loops[5|6]}};
OpenLoopsModel = "SM_yuksel";
Approximation = "onlyh";

(* Z Z, same and opposite flavour leptons*)
(* eeexexgg     *) AddProcess[FeynArtsProcess -> {V[5], V[5]} -> {F[2,{1}],  F[2,{1}], -F[2,{1}], -F[2,{1}]}, SetParameters -> JoinOptions[{ME -> 0, YE -> 0}]];
(* eexmmxgg     *) AddProcess[FeynArtsProcess -> {V[5], V[5]} -> {F[2,{1}],  F[2,{2}], -F[2,{1}], -F[2,{2}]}, SetParameters -> JoinOptions[{ME -> 0, YE -> 0, MM -> 0, YM -> 0}]];

(* Z Z, same and opposite flavour neutrinos *)
(* nenenexnexgg *) AddProcess[FeynArtsProcess -> {V[5], V[5]} -> {F[1,{1}],  F[1,{1}], -F[1,{1}], -F[1,{1}]}, SetParameters -> JoinOptions[{ME -> 0, YE -> 0}]];
(* nenexnmnmxgg *) AddProcess[FeynArtsProcess -> {V[5], V[5]} -> {F[1,{1}],  F[1,{2}], -F[1,{1}], -F[1,{2}]}, SetParameters -> JoinOptions[{ME -> 0, YE -> 0, MM -> 0, YM -> 0}]];

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

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

(* W+ W-, opposite flavour *)
(* nexnmemxgg   *) AddProcess[FeynArtsProcess -> {V[5], V[5]} -> {F[2,{2}],  F[1,{1}], -F[2,{1}], -F[1,{2}]}, SetParameters -> JoinOptions[{ME -> 0, YE -> 0, MM -> 0, YM -> 0}]];
