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

SelectInterference = {eQED -> 4};

(* ttxddxaa *) AddProcess[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {F[3,{3}], -F[3,{3}], V[1], V[1]}];
(* uuxttxaa *) AddProcess[FeynArtsProcess -> {F[3,{1}], -F[3,{1}]} -> {F[3,{3}], -F[3,{3}], V[1], V[1]}];
(* ttxaagg  *) AddProcess[FeynArtsProcess -> {V[5]    ,  V[5]    } -> {F[3,{3}], -F[3,{3}], V[1], V[1]}];
(* ttxbbxaa *) AddProcess[FeynArtsProcess -> {F[4,{3}], -F[4,{3}]} -> {F[3,{3}], -F[3,{3}], V[1], V[1]}];

(* ttxddxaag *) AddReal[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {F[3,{3}], -F[3,{3}], V[1], V[1], V[5]}];
(* uuxttxaag *) AddReal[FeynArtsProcess -> {F[3,{1}], -F[3,{1}]} -> {F[3,{3}], -F[3,{3}], V[1], V[1], V[5]}];
(* ttxaaggg  *) AddReal[FeynArtsProcess -> {V[5]    ,  V[5]    } -> {F[3,{3}], -F[3,{3}], V[1], V[1], V[5]}];
(* ttxbbxaag *) AddReal[FeynArtsProcess -> {F[4,{3}], -F[4,{3}]} -> {F[3,{3}], -F[3,{3}], V[1], V[1], V[5]}];

ChannelMap = {
  "ttxbbxaa" -> "ttxddxaa MB=0",
  (* real *)
  "ttxbbxaag" -> "ttxddxaag MB=0"
};
