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

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

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

(* ttxhgggg   *) AddReal[FeynArtsProcess -> {V[5], V[5]} -> {S[1], F[3,{3}], -F[3,{3}], V[5], V[5]}];
(* ttxddxhgg  *) AddReal[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {S[1], F[3,{3}], -F[3,{3}], V[5], V[5]}];
(* ttxbbxhgg  *) AddReal[FeynArtsProcess -> {F[4,{3}], -F[4,{3}]} -> {S[1], F[3,{3}], -F[3,{3}], V[5], V[5]}];

(* ttxdddxdxh *) AddReal[FeynArtsProcess -> {F[4,{1}],  F[4,{1}]} -> {S[1], F[3,{3}], -F[3,{3}], F[4,{1}],  F[4,{1}]}];
(* uuxttxddxh *) AddReal[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {S[1], F[3,{3}], -F[3,{3}], F[3,{1}], -F[3,{1}]}];
(* ttxbbbxbxh *) AddReal[FeynArtsProcess -> {F[4,{3}],  F[4,{3}]} -> {S[1], F[3,{3}], -F[3,{3}], F[4,{3}],  F[4,{3}]}];
(* ttxddxbbxh *) AddReal[FeynArtsProcess -> {F[4,{3}], -F[4,{3}]} -> {S[1], F[3,{3}], -F[3,{3}], F[4,{1}], -F[4,{1}]}];

ChannelMap = {
  "uuxttxhg" -> "ttxddxhg[3,4,1,2,5,6]",
  "ttxbbxhg" -> "ttxddxhg MB=0",
  (* real *)
  "ttxddxbbxh" -> "uuxttxddxh[3,4,1,2,5,6,7] MB=0",
  "ttxbbbxbxh" -> "ttxdddxdxh MB=0",
  "ttxbbxhgg"  -> "ttxddxhgg MB=0",
  "ttxddxssxh" -> "uuxttxddxh[3,4,1,2,5,6,7]",
  "uuuxuxttxh" -> "ttxdddxdxh[3,4,5,6,1,2,7]",
  "uuxccxttxh" -> "uuxttxddxh[1,2,5,6,3,4,7]",
  "ccxttxddxh" -> "uuxttxddxh",
  "uuxttxhgg"  -> "ttxddxhgg[3,4,1,2,5,6,7]",
  "uuxttxbbxh" -> "ttxddxbbxh[3,4,1,2,5,6,7]"
};
