(* description = "pp -> HH+jets|pp -> HH+jets|pp -> HHjj|EW production of a Higgs-pair plus two jets (no VBF approximation)." *)

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

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

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

ChannelMap = {
  "bbbxbxhh" -> "dddxdxhh MB=0",
  "uuxbbxhh" -> "ccxddxhh MB=0",
  "ddxbbxhh" -> "ddxssxhh MB=0"
};