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

SelectInterference = {eQED -> 4};

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

(* dddxdxaa *) AddReal[FeynArtsProcess -> { F[4,{1}],  F[4,{1}]} -> { V[1], V[1], F[4,{1}],  F[4,{1}]}];
(* uuuxuxaa *) AddReal[FeynArtsProcess -> { F[3,{1}],  F[3,{1}]} -> { V[1], V[1], F[3,{1}],  F[3,{1}]}];
(* uuxddxaa *) AddReal[FeynArtsProcess -> { F[4,{1}], -F[4,{1}]} -> { V[1], V[1], F[3,{1}], -F[3,{1}]}];
(* ddxssxaa *) AddReal[FeynArtsProcess -> { F[4,{1}], -F[4,{1}]} -> { V[1], V[1], F[4,{2}], -F[4,{2}]}];
(* uuxccxaa *) AddReal[FeynArtsProcess -> { F[3,{1}], -F[3,{1}]} -> { V[1], V[1], F[3,{2}], -F[3,{2}]}];
(* ddxaagg  *) AddReal[FeynArtsProcess -> { F[4,{1}], -F[4,{1}]} -> { V[1], V[1], V[5]    ,  V[5]    }];
(* uuxaagg  *) AddReal[FeynArtsProcess -> { F[3,{1}], -F[3,{1}]} -> { V[1], V[1], V[5]    ,  V[5]    }];

(* bbbxbxaa *) AddReal[FeynArtsProcess -> { F[4,{3}], -F[4,{3}]} -> { V[1], V[1], F[4,{3}], -F[4,{3}]}];
(* uuxbbxaa *) AddReal[FeynArtsProcess -> { F[3,{1}], -F[3,{1}]} -> { V[1], V[1], F[4,{3}], -F[4,{3}]}];
(* ddxbbxaa *) AddReal[FeynArtsProcess -> { F[4,{1}], -F[4,{1}]} -> { V[1], V[1], F[4,{3}], -F[4,{3}]}];
(* bbxaagg  *) AddReal[FeynArtsProcess -> { F[4,{3}], -F[4,{3}]} -> { V[1], V[1], V[5]    ,  V[5]    }];

ChannelMap = {
  "bbxaag" -> "ddxaag MB=0",
  (* reals *)
  "ccxddxaa" -> "uuxddxaa",
  "ddxbbxaa" -> "ddxssxaa MB=0",
  "uuxbbxaa" -> "uuxddxaa MB=0",
  "bbbxbxaa" -> "dddxdxaa MB=0",
  "bbxaagg"  -> "ddxaagg MB=0"
};
