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

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

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

(* ttxhggg  *) AddReal[FeynArtsProcess -> {V[5]    ,  V[5]    } -> {S[1], S[1], F[3,{3}], -F[3,{3}], V[5]}];
(* ttxddxhg *) AddReal[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {S[1], S[1], F[3,{3}], -F[3,{3}], V[5]}];
(* ttxbbxhg *) AddReal[FeynArtsProcess -> {F[4,{3}], -F[4,{3}]} -> {S[1], S[1], F[3,{3}], -F[3,{3}], V[5]}];

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