
(* description = "pp -> tt+jets|pp -> tt+jets|pp -> ttbb|Top-quark pair production with a b-bbar pair." *)

SelectInterference = {eQED -> 0};

(* ttxbbbxbx *) AddProcess[FeynArtsProcess -> {F[4,{3}],  F[4,{3}]} -> {F[3,{3}], -F[3,{3}], F[4,{3}],  F[4,{3}]}];
(* ttxddxbbx *) AddProcess[FeynArtsProcess -> {F[4,{3}], -F[4,{3}]} -> {F[3,{3}], -F[3,{3}], F[4,{1}], -F[4,{1}]}];
(* ttxbbxgg  *) AddProcess[FeynArtsProcess -> {F[4,{3}], -F[4,{3}]} -> {F[3,{3}], -F[3,{3}], V[5]    ,  V[5]    }];


(* ttxbbbxbxg *) AddReal[FeynArtsProcess -> {F[4,{3}],  F[4,{3}]} -> {F[3,{3}], -F[3,{3}], F[4,{3}],  F[4,{3}], V[5]}];
(* ttxddxbbxg *) AddReal[FeynArtsProcess -> {F[4,{3}], -F[4,{3}]} -> {F[3,{3}], -F[3,{3}], F[4,{1}], -F[4,{1}], V[5]}];
(* ttxbbxggg  *) AddReal[FeynArtsProcess -> {F[4,{3}], -F[4,{3}]} -> {F[3,{3}], -F[3,{3}], V[5]    ,  V[5]    , V[5]}];

ChannelMap = {
  "uuxttxbbx" -> "ttxddxbbx[3,4,1,2,5,6]",
  (* reals *)
  "uuxttxbbxg" -> "ttxddxbbxg[3,4,1,2,5,6,7]"
};
