
(* description = "pp -> H+jets|pp -> H+jets|pp -> H(j)|Higgs production in the Higgs effective field theory with leptonic (mu+mu- and 4f) decays at NLO EW." *)

OpenLoopsModel = "HEFT";
SelectInterference = {gQCD -> 4, DOI -> 0};
OnTheFlyMode = 0;
UnitaryGauge = False;

SelectTreeDiagrams = (HasResonance[S[1], {1,2}][##] &);
SelectLoopDiagrams = (HasResonance[S[1], {1,2}][##] &);
SelectCTDiagrams = (HasResonance[S[1], {1,2}][##] &);

(* 2f *)

(* mmxgg     *) AddProcess[FeynArtsProcess -> {V[5], V[5]} -> {F[2,{2}],  -F[2,{2}]}, SetParameters -> JoinOptions[{MB -> 0, YB -> 0, MM -> MM, YM -> YM}], SelectCoupling -> (Exponent[#1, eQED] == 2 + 2 * #2 && Exponent[#1, gQCD] == 2 &)];


(* 4f *)

(* eeexexgg     *) AddProcess[FeynArtsProcess -> {V[5], V[5]} -> {F[2,{1}],  F[2,{1}], -F[2,{1}], -F[2,{1}]}, SetParameters -> JoinOptions[{MB -> 0, YB -> 0, ME -> 0, YE -> 0, MM -> 0, YM -> 0}], SelectCoupling -> (Exponent[#1, eQED] == 4 + 2 * #2 && Exponent[#1, gQCD] == 2 &)];
(* eeexexgg     *) AddProcess[FeynArtsProcess -> {V[5], V[5]} -> {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}], SelectCoupling -> (Exponent[#1, eQED] == 4 + 2 * #2 && Exponent[#1, gQCD] == 2 &)];

(* 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[{MB -> 0, YB -> 0, ME -> 0, YE -> 0, MM -> 0, YM -> 0}], SelectCoupling -> (Exponent[#1, eQED] == 4 + 2 * #2 && Exponent[#1, gQCD] == 2 &)];

(* 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[{MB -> 0, YB -> 0, ME -> 0, YE -> 0, MM -> 0, YM -> 0}], SelectCoupling -> (Exponent[#1, eQED] == 4 + 2 * #2 && Exponent[#1, gQCD] == 2 &)];

(* 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[{MB -> 0, YB -> 0, ME -> 0, YE -> 0, MM -> 0, YM -> 0}], SelectCoupling -> (Exponent[#1, eQED] == 4 + 2 * #2 && Exponent[#1, gQCD] == 2 &)];


