
(* description = "pp -> VV+jets|pp -> AA+jets|pp -> AA|Photon pair production" *)

SelectInterference = {eQED -> 4};

(* ddxaa *) AddProcess[FeynArtsProcess -> { F[4,{1}], -F[4,{1}]} -> { V[1], V[1]}];
(* uuxaa *) AddProcess[FeynArtsProcess -> { F[3,{1}], -F[3,{1}]} -> { V[1], V[1]}];
(* bbxaa *) AddProcess[FeynArtsProcess -> { F[4,{3}], -F[4,{3}]} -> { V[1], V[1]}];

(* ddxaag *) AddReal[FeynArtsProcess -> { F[4,{1}], -F[4,{1}]} -> { V[1], V[1], V[5]}];
(* uuxaag *) AddReal[FeynArtsProcess -> { F[3,{1}], -F[3,{1}]} -> { V[1], V[1], V[5]}];
(* bbxaag *) AddReal[FeynArtsProcess -> { F[4,{3}], -F[4,{3}]} -> { V[1], V[1], V[5]}];

ChannelMap = {
  "bbxaa" -> "ddxaa MB=0",
  (* real *)
  "bbxaag" -> "ddxaag MB=0"
};
