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

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

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

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

ChannelMap = {
  "uuxttxh" -> "ttxddxh[3,4,1,2,5]",
  "ttxbbxh" -> "ttxddxh MB=0 YB=0",
  (* real *)
  "uuxttxhg" -> "ttxddxhg[3,4,1,2,5,6]",
  "ttxbbxhg" -> "ttxddxhg MB=0 YB=0"
};
