
(* description = "pp -> Vtt+jets|pp -> Z/A(->ll)tt+jets|pp -> lltt|Off-shell Z/A boson plus t t-bar production with leptonic decays (l+l- and nn) at NLO EW." *)

SelectInterference = {eQED -> 4};
SetParameters = JoinOptions[{ME -> 0}];
ColourCorrelations = True;

(* eexttxgg    *) AddProcess[FeynArtsProcess -> {V[5]    ,  V[5]    } -> {F[2,{1}], -F[2,{1}], F[3,{3}], -F[3,{3}]}];
(* eexttxddx   *) AddProcess[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {F[2,{1}], -F[2,{1}], F[3,{3}], -F[3,{3}]}];
(* eexuuxttx   *) AddProcess[FeynArtsProcess -> {F[3,{1}], -F[3,{1}]} -> {F[2,{1}], -F[2,{1}], F[3,{3}], -F[3,{3}]}];
(* eexttxbbx   *) AddProcess[FeynArtsProcess -> {F[4,{3}], -F[4,{3}]} -> {F[2,{1}], -F[2,{1}], F[3,{3}], -F[3,{3}]}];

(* nenexttxgg  *) AddProcess[FeynArtsProcess -> {V[5]    ,  V[5]    } -> {F[1,{1}], -F[1,{1}], F[3,{3}], -F[3,{3}]}];
(* nenexttxddx *) AddProcess[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {F[1,{1}], -F[1,{1}], F[3,{3}], -F[3,{3}]}];
(* nenexuuxttx *) AddProcess[FeynArtsProcess -> {F[3,{1}], -F[3,{1}]} -> {F[1,{1}], -F[1,{1}], F[3,{3}], -F[3,{3}]}];
(* nenexttxbbx *) AddProcess[FeynArtsProcess -> {F[4,{3}], -F[4,{3}]} -> {F[1,{1}], -F[1,{1}], F[3,{3}], -F[3,{3}]}];

(* eexttxggg    *) AddReal[FeynArtsProcess -> {V[5]    ,  V[5]    } -> {F[2,{1}], -F[2,{1}], F[3,{3}], -F[3,{3}], V[5]}];
(* eexttxddxg   *) AddReal[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {F[2,{1}], -F[2,{1}], F[3,{3}], -F[3,{3}], V[5]}];
(* eexuuxttxg   *) AddReal[FeynArtsProcess -> {F[3,{1}], -F[3,{1}]} -> {F[2,{1}], -F[2,{1}], F[3,{3}], -F[3,{3}], V[5]}];
(* eexttxbbxg   *) AddReal[FeynArtsProcess -> {F[4,{3}], -F[4,{3}]} -> {F[2,{1}], -F[2,{1}], F[3,{3}], -F[3,{3}], V[5]}];

(* nenexttxggg  *) AddReal[FeynArtsProcess -> {V[5]    ,  V[5]    } -> {F[1,{1}], -F[1,{1}], F[3,{3}], -F[3,{3}], V[5]}];
(* nenexttxddxg *) AddReal[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {F[1,{1}], -F[1,{1}], F[3,{3}], -F[3,{3}], V[5]}];
(* nenexuuxttxg *) AddReal[FeynArtsProcess -> {F[3,{1}], -F[3,{1}]} -> {F[1,{1}], -F[1,{1}], F[3,{3}], -F[3,{3}], V[5]}];
(* nenexttxbbxg *) AddReal[FeynArtsProcess -> {F[4,{3}], -F[4,{3}]} -> {F[1,{1}], -F[1,{1}], F[3,{3}], -F[3,{3}], V[5]}];

ChannelMap = {
  "eexttxbbx"   -> "eexttxddx MB=0",
  "nenexttxbbx" -> "nenexttxddx MB=0",
  (* real *)
  "eexttxbbxg"   -> "eexttxddxg MB=0",
  "nenexttxbbxg" -> "nenexttxddxg MB=0"
};
