
(* description = "pp -> HH+jets|pp -> HH+jets|pp -> HHj|Higgs-pair plus jet production in the Higgs effective field theory." *)

OpenLoopsModel = "HEFT";

SelectInterference = {eQED -> 4};
SetParameters = JoinOptions[{MB -> 0, YB -> 0}];
OnTheFlyMode = 0;

(* hhggg  *) AddProcess[FeynArtsProcess -> {V[5]    ,  V[5]    } -> {S[1], S[1], V[5]}, SelectCoupling -> ( Exponent[#1, gQCD] == 3 + 2*#2 && Exponent[#1, eQED] == 2 &)];
(* ddxhhg *) AddProcess[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {S[1], S[1], V[5]}, SelectCoupling -> ( Exponent[#1, gQCD] == 3 + 2*#2 && Exponent[#1, eQED] == 2 &)];

(* dddxdxhh *) AddReal[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {S[1], S[1], F[4,{1}], -F[4,{1}]}, SelectCoupling -> ( Exponent[#1, gQCD] == 4 && Exponent[#1, eQED] == 2 &)];
(* uuxddxhh *) AddReal[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {S[1], S[1], F[3,{1}], -F[3,{1}]}, SelectCoupling -> ( Exponent[#1, gQCD] == 4 && Exponent[#1, eQED] == 2 &)];
(* ddxhhgg  *) AddReal[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {S[1], S[1], V[5]    ,  V[5]    }, SelectCoupling -> ( Exponent[#1, gQCD] == 4 && Exponent[#1, eQED] == 2 &)];
(* hhgggg   *) AddReal[FeynArtsProcess -> {V[5]    ,  V[5]    } -> {S[1], S[1], V[5]    ,  V[5]    }, SelectCoupling -> ( Exponent[#1, gQCD] == 4 && Exponent[#1, eQED] == 2 &)];

ChannelMap = {
  "uuxhhg" -> "ddxhhg",
  "bbxhhg" -> "ddxhhg MB=0",
  (* reals *)
  "ddxbbxhh" -> "uuxddxhh MB=0",
  "uuxbbxhh" -> "uuxddxhh MB=0",
  "bbbxbxhh" -> "dddxdxhh MB=0",
  "bbxhhgg"  -> "ddxhhgg MB=0",
  "ddxssxhh" -> "uuxddxhh",
  "uuuxuxhh" -> "dddxdxhh",
  "uuxccxhh" -> "uuxddxhh",
  "ccxddxhh" -> "uuxddxhh",
  "uuxhhgg"  -> "ddxhhgg"
};
