
(* description = "pp -> HHtt+jets|pp -> HHtt+jets|pp -> HHtt|Di-Higgs plus t t-bar production." *)

OpenLoopsModel = "SM_yuksel";
SelectInterference = {gQCD -> 4};
SetParameters = JoinOptions[{YukB -> 1, YukT -> 1}];


SelectCoupling = ((Exponent[#1,lambdaHHHH] == 0 + 1*#2) && (Exponent[#1, YT] == 1 + 1 * #2 || Exponent[#1, YT] == 2) &);

(* ttxhhgg  *) AddProcess[FeynArtsProcess -> {V[5]    ,  V[5]    } -> {S[1], S[1], F[3,{3}], -F[3,{3}]}];
(* ttxddxhh *) AddProcess[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {S[1], S[1], F[3,{3}], -F[3,{3}]}];


ChannelMap = {
  "uuxttxhh" -> "ttxddxhh[3,4,1,2,5]",
  "ttxbbxhh" -> "ttxddxhh MB=0 YB=0",
  (* real *)
  "uuxttxhhg" -> "ttxddxhhg[3,4,1,2,5,6]",
  "ttxbbxhhg" -> "ttxddxhhg MB=0 YB=0"
};
