(* description = "pp -> H+jets|pp -> H+jets|pp -> Hjj|EW production of Higgs plus two jets (no VBF approximation)." *)

OpenLoopsModel = "SM_yuksel";
SelectCoupling = (Exponent[#1, eQED] === 3 &);
SelectInterference = {eQED -> 6};
SetParameters = JoinOptions[{MB -> 0, YB -> 0}];
OTFColourCorrelations = False;

(* dddxdxh *) AddProcess[FeynArtsProcess -> {F[4,{1}], F[4,{1}]} -> {S[1], F[4,{1}], F[4,{1}]}, UnitaryGauge -> False];
(* uuuxuxh *) AddProcess[FeynArtsProcess -> {F[3,{1}], F[3,{1}]} -> {S[1], F[3,{1}], F[3,{1}]}, UnitaryGauge -> False];
(* uuxddxh *) AddProcess[FeynArtsProcess -> {F[4,{1}],-F[4,{1}]} -> {S[1], F[3,{1}],-F[3,{1}]}, UnitaryGauge -> False];
(* ccxddxh *) AddProcess[FeynArtsProcess -> {F[4,{1}], F[3,{2}]} -> {S[1], F[4,{1}], F[3,{2}]}, SelectLoopDiagrams -> SameQuarkLineGluon];
(* uuxccxh *) AddProcess[FeynArtsProcess -> {F[3,{1}], F[3,{2}]} -> {S[1], F[3,{1}], F[3,{2}]}, SelectLoopDiagrams -> SameQuarkLineGluon];
(* ddxssxh *) AddProcess[FeynArtsProcess -> {F[4,{1}], F[4,{2}]} -> {S[1], F[4,{1}], F[4,{2}]}, SelectLoopDiagrams -> SameQuarkLineGluon];
(* uxcdsxh *) AddProcess[FeynArtsProcess -> {F[4,{1}],-F[4,{2}]} -> {S[1], F[3,{1}],-F[3,{2}]}, SelectLoopDiagrams -> SameQuarkLineGluon];

(* dddxdxhg *) AddReal[FeynArtsProcess -> {F[4,{1}], F[4,{1}]} -> {S[1], F[4,{1}], F[4,{1}], V[5]}];
(* uuuxuxhg *) AddReal[FeynArtsProcess -> {F[3,{1}], F[3,{1}]} -> {S[1], F[3,{1}], F[3,{1}], V[5]}];
(* uuxddxhg *) AddReal[FeynArtsProcess -> {F[4,{1}],-F[4,{1}]} -> {S[1], F[3,{1}],-F[3,{1}], V[5]}];
(* ccxddxhg *) AddReal[FeynArtsProcess -> {F[4,{1}], F[3,{2}]} -> {S[1], F[4,{1}], F[3,{2}], V[5]}];
(* uuxccxhg *) AddReal[FeynArtsProcess -> {F[3,{1}], F[3,{2}]} -> {S[1], F[3,{1}], F[3,{2}], V[5]}];
(* ddxssxhg *) AddReal[FeynArtsProcess -> {F[4,{1}], F[4,{2}]} -> {S[1], F[4,{1}], F[4,{2}], V[5]}];
(* uxcdsxhg *) AddReal[FeynArtsProcess -> {F[4,{1}],-F[4,{2}]} -> {S[1], F[3,{1}],-F[3,{2}], V[5]}];

ChannelMap = {
  "bbbxbxh" -> "dddxdxh MB=0",
  "uuxbbxh" -> "ccxddxh MB=0",
  "ddxbbxh" -> "ddxssxh MB=0"
};