
(* description = "pp -> ttVV+jets|pp -> ttZZ+jets|pp -> ttZZ|Top-quark pair and two Z-bosons production." *)

SelectInterference = {eQED -> 4};

(* ttxddxzz *) AddProcess[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {F[3,{3}], -F[3,{3}], V[2], V[2]}];
(* uuxttxzz *) AddProcess[FeynArtsProcess -> {F[3,{1}], -F[3,{1}]} -> {F[3,{3}], -F[3,{3}], V[2], V[2]}];
(* ttxzzgg  *) AddProcess[FeynArtsProcess -> {V[5]    ,  V[5]    } -> {F[3,{3}], -F[3,{3}], V[2], V[2]}];
(* ttxbbxzz *) AddProcess[FeynArtsProcess -> {F[4,{3}], -F[4,{3}]} -> {F[3,{3}], -F[3,{3}], V[2], V[2]}];

(* ttxddxzzg *) AddReal[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {F[3,{3}], -F[3,{3}], V[2], V[2], V[5]}];
(* uuxttxzzg *) AddReal[FeynArtsProcess -> {F[3,{1}], -F[3,{1}]} -> {F[3,{3}], -F[3,{3}], V[2], V[2], V[5]}];
(* ttxzzggg  *) AddReal[FeynArtsProcess -> {V[5]    ,  V[5]    } -> {F[3,{3}], -F[3,{3}], V[2], V[2], V[5]}];
(* ttxbbxzzg *) AddReal[FeynArtsProcess -> {F[4,{3}], -F[4,{3}]} -> {F[3,{3}], -F[3,{3}], V[2], V[2], V[5]}];

ChannelMap = {
  "ttxbbxzz" -> "ttxddxzz MB=0",
  (* real *)
  "ttxbbxzzg" -> "ttxddxzzg MB=0"
};
