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

OpenLoopsModel = "HEFT";

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

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

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

ChannelMap = {
  (* reals *)
  "uuxhhg" -> "ddxhhg",
  "bbxhhg" -> "ddxhhg MB=0"
};
