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

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

(* as^3 a^2 + as^3 a^3 *)

(* eexttxggg   *) AddProcess[FeynArtsProcess -> {V[5]	  , V[5]} -> {F[2,{1}], -F[2,{1}], F[3,{3}], -F[3,{3}], V[5]},UnitaryGauge -> False,SelectInterference -> {gQCD -> 6},SelectCoupling -> (Exponent[#1, eQED] == 2+2*#2 &)];
(* eexttxddxg   *) AddProcess[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {F[2,{1}], -F[2,{1}], F[3,{3}], -F[3,{3}], V[5]},UnitaryGauge -> False,SelectInterference -> {gQCD -> 6},SelectCoupling -> (Exponent[#1, eQED] == 4-2*#2 || Exponent[#1, eQED] == 2 + 2 * #2 &)];
(* eexuuxttxg   *) AddProcess[FeynArtsProcess -> {F[3,{1}], -F[3,{1}]} -> {F[2,{1}], -F[2,{1}], F[3,{3}], -F[3,{3}], V[5]},UnitaryGauge -> False,SelectInterference -> {gQCD -> 6},SelectCoupling -> (Exponent[#1, eQED] == 4 - 2 * #2 || Exponent[#1, eQED] == 2 + 2 * #2 &)];
(* eexttxbbxg   *) AddProcess[FeynArtsProcess -> {F[4,{3}], -F[4,{3}]} -> {F[2,{1}], -F[2,{1}], F[3,{3}], -F[3,{3}], V[5]},UnitaryGauge -> False,SelectInterference -> {gQCD -> 6},SelectCoupling -> (Exponent[#1, eQED] == 4 - 2 * #2 || Exponent[#1, eQED] == 2 + 2 * #2 &)];

(* nenexttxggg  *) AddProcess[FeynArtsProcess -> {V[5]    ,  V[5]    } -> {F[1,{1}], -F[1,{1}], F[3,{3}], -F[3,{3}], V[5]},UnitaryGauge -> False,SelectInterference -> {gQCD -> 6},SelectCoupling -> (Exponent[#1, eQED] == 2 + 2 * #2 &)];
(* nenexttxddxg *) AddProcess[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {F[1,{1}], -F[1,{1}], F[3,{3}], -F[3,{3}], V[5]},UnitaryGauge -> False,SelectInterference -> {gQCD -> 6},SelectCoupling -> (Exponent[#1, eQED] == 4 - 2 * #2 || Exponent[#1, eQED] == 2 + 2 * #2 &)];
(* nenexuuxttxg *) AddProcess[FeynArtsProcess -> {F[3,{1}], -F[3,{1}]} -> {F[1,{1}], -F[1,{1}], F[3,{3}], -F[3,{3}], V[5]},UnitaryGauge -> False,SelectInterference -> {gQCD -> 6},SelectCoupling -> (Exponent[#1, eQED] == 4 - 2 * #2 || Exponent[#1, eQED] == 2 + 2 * #2 &)];
(* nenexttxbbxg *) AddProcess[FeynArtsProcess -> {F[4,{3}], -F[4,{3}]} -> {F[1,{1}], -F[1,{1}], F[3,{3}], -F[3,{3}], V[5]},UnitaryGauge -> False,SelectInterference -> {gQCD -> 6},SelectCoupling -> (Exponent[#1, eQED] == 4 - 2 * #2 || Exponent[#1, eQED] == 2 + 2 * #2 &)];


(* eexttxagg    *) AddReal[FeynArtsProcess -> {V[5]    ,  V[5]    } -> {F[2,{1}], -F[2,{1}], F[3,{3}], -F[3,{3}], V[5], V[1]},SelectInterference -> {gQCD -> 6}];
(* eexttxddxa   *) AddReal[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {F[2,{1}], -F[2,{1}], F[3,{3}], -F[3,{3}], V[5], V[1]},SelectInterference -> {gQCD -> 6}];
(* eexuuxttxa   *) AddReal[FeynArtsProcess -> {F[3,{1}], -F[3,{1}]} -> {F[2,{1}], -F[2,{1}], F[3,{3}], -F[3,{3}], V[5], V[1]},SelectInterference -> {gQCD -> 6}];
(* eexttxbbxa   *) AddReal[FeynArtsProcess -> {F[4,{3}], -F[4,{3}]} -> {F[2,{1}], -F[2,{1}], F[3,{3}], -F[3,{3}], V[5], V[1]},SelectInterference -> {gQCD -> 6}];
(* eexttxddxg   *) AddReal[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {F[2,{1}], -F[2,{1}], F[3,{3}], -F[3,{3}], V[5], V[5]},SelectInterference -> {gQCD -> 6},SelectCoupling -> (Exponent[#1, eQED] == 2 || Exponent[#1, eQED] == 4 &)];
(* eexuuxttxg   *) AddReal[FeynArtsProcess -> {F[3,{1}], -F[3,{1}]} -> {F[2,{1}], -F[2,{1}], F[3,{3}], -F[3,{3}], V[5], V[5]},SelectInterference -> {gQCD -> 6},SelectCoupling -> (Exponent[#1, eQED] == 2 || Exponent[#1, eQED] == 4 &)];
(* eexttxbbxg   *) AddReal[FeynArtsProcess -> {F[4,{3}], -F[4,{3}]} -> {F[2,{1}], -F[2,{1}], F[3,{3}], -F[3,{3}], V[5], V[5]},SelectInterference -> {gQCD -> 6},SelectCoupling -> (Exponent[#1, eQED] == 2 || Exponent[#1, eQED] == 4 &)];

(* nenexttxagg  *) AddReal[FeynArtsProcess -> {V[5]    ,  V[5]    } -> {F[1,{1}], -F[1,{1}], F[3,{3}], -F[3,{3}], V[5], V[1]},SelectInterference -> {gQCD -> 6}];
(* nenexttxddxa *) AddReal[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {F[1,{1}], -F[1,{1}], F[3,{3}], -F[3,{3}], V[5], V[1]},SelectInterference -> {gQCD -> 6}];
(* nenexuuxttxa *) AddReal[FeynArtsProcess -> {F[3,{1}], -F[3,{1}]} -> {F[1,{1}], -F[1,{1}], F[3,{3}], -F[3,{3}], V[5], V[1]},SelectInterference -> {gQCD -> 6}];
(* nenexttxbbxa *) AddReal[FeynArtsProcess -> {F[4,{3}], -F[4,{3}]} -> {F[1,{1}], -F[1,{1}], F[3,{3}], -F[3,{3}], V[5], V[1]},SelectInterference -> {gQCD -> 6}];
(* nenexttxddxg *) AddReal[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {F[1,{1}], -F[1,{1}], F[3,{3}], -F[3,{3}], V[5], V[5]},SelectInterference -> {gQCD -> 6},SelectCoupling -> (Exponent[#1, eQED] == 2 || Exponent[#1, eQED] == 4 &)];
(* nenexuuxttxg *) AddReal[FeynArtsProcess -> {F[3,{1}], -F[3,{1}]} -> {F[1,{1}], -F[1,{1}], F[3,{3}], -F[3,{3}], V[5], V[5]},SelectInterference -> {gQCD -> 6},SelectCoupling -> (Exponent[#1, eQED] == 2 || Exponent[#1, eQED] == 4 &)];
(* nenexttxbbxg *) AddReal[FeynArtsProcess -> {F[4,{3}], -F[4,{3}]} -> {F[1,{1}], -F[1,{1}], F[3,{3}], -F[3,{3}], V[5], V[5]},SelectInterference -> {gQCD -> 6},SelectCoupling -> (Exponent[#1, eQED] == 2 || Exponent[#1, eQED] == 4 &)];


(* as^2 a^2 *)

(* eexttxddx   *) AddTree[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {F[2,{1}], -F[2,{1}], F[3,{3}], -F[3,{3}], V[5]},UnitaryGauge -> False,SelectInterference -> {gQCD -> 4},SelectCoupling -> (Exponent[#1, eQED] == 4-2*#2 || Exponent[#1, eQED] == 2 + 2 * #2 &)];
(* eexuuxttx   *) AddTree[FeynArtsProcess -> {F[3,{1}], -F[3,{1}]} -> {F[2,{1}], -F[2,{1}], F[3,{3}], -F[3,{3}], V[5]},UnitaryGauge -> False,SelectInterference -> {gQCD -> 4},SelectCoupling -> (Exponent[#1, eQED] == 4 - 2 * #2 || Exponent[#1, eQED] == 2 + 2 * #2 &)];
(* eexttxbbx   *) AddTree[FeynArtsProcess -> {F[4,{3}], -F[4,{3}]} -> {F[2,{1}], -F[2,{1}], F[3,{3}], -F[3,{3}], V[5]},UnitaryGauge -> False,SelectInterference -> {gQCD -> 4},SelectCoupling -> (Exponent[#1, eQED] == 4 - 2 * #2 || Exponent[#1, eQED] == 2 + 2 * #2 &)];

(* nenexttxddx *) AddTree[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {F[1,{1}], -F[1,{1}], F[3,{3}], -F[3,{3}], V[5]},UnitaryGauge -> False,SelectInterference -> {gQCD -> 4},SelectCoupling -> (Exponent[#1, eQED] == 4 - 2 * #2 || Exponent[#1, eQED] == 2 + 2 * #2 &)];
(* nenexuuxttx *) AddTree[FeynArtsProcess -> {F[3,{1}], -F[3,{1}]} -> {F[1,{1}], -F[1,{1}], F[3,{3}], -F[3,{3}], V[5]},UnitaryGauge -> False,SelectInterference -> {gQCD -> 4},SelectCoupling -> (Exponent[#1, eQED] == 4 - 2 * #2 || Exponent[#1, eQED] == 2 + 2 * #2 &)];
(* nenexttxbbx *) AddTree[FeynArtsProcess -> {F[4,{3}], -F[4,{3}]} -> {F[1,{1}], -F[1,{1}], F[3,{3}], -F[3,{3}], V[5]},UnitaryGauge -> False,SelectInterference -> {gQCD -> 4},SelectCoupling -> (Exponent[#1, eQED] == 4 - 2 * #2 || Exponent[#1, eQED] == 2 + 2 * #2 &)];


(*
(* as a^4 + as^2 a^4 *)

(* eexttxddx   *) AddProcess[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {F[2,{1}], -F[2,{1}], F[3,{3}], -F[3,{3}], V[5]},UnitaryGauge -> False,SelectInterference -> {eQED -> 8},SelectCoupling -> (Exponent[#1, gQCD] == 3 - 2*#2 || Exponent[#1, gQCD] == 1 + 2 * #2 &)];
(* eexuuxttx   *) AddProcess[FeynArtsProcess -> {F[3,{1}], -F[3,{1}]} -> {F[2,{1}], -F[2,{1}], F[3,{3}], -F[3,{3}], V[5]},UnitaryGauge -> False,SelectInterference -> {eQED -> 8},SelectCoupling -> (Exponent[#1, gQCD] == 3 - 2 * #2 || Exponent[#1, gQCD] == 1 + 2 * #2 &)];
(* eexttxbbx   *) AddProcess[FeynArtsProcess -> {F[4,{3}], -F[4,{3}]} -> {F[2,{1}], -F[2,{1}], F[3,{3}], -F[3,{3}], V[5]},UnitaryGauge -> False,SelectInterference -> {eQED -> 8},SelectCoupling -> (Exponent[#1, gQCD] == 3 - 2 * #2 || Exponent[#1, gQCD] == 1 + 2 * #2 &)];

(* nenexttxddx *) AddProcess[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {F[1,{1}], -F[1,{1}], F[3,{3}], -F[3,{3}], V[5]},UnitaryGauge -> False,SelectInterference -> {eQED -> 8},SelectCoupling -> (Exponent[#1, gQCD] == 3 - 2 * #2 || Exponent[#1, gQCD] == 1 + 2 * #2 &)];
(* nenexuuxttx *) AddProcess[FeynArtsProcess -> {F[3,{1}], -F[3,{1}]} -> {F[1,{1}], -F[1,{1}], F[3,{3}], -F[3,{3}], V[5]},UnitaryGauge -> False,SelectInterference -> {eQED -> 8},SelectCoupling -> (Exponent[#1, gQCD] == 3 - 2 * #2 || Exponent[#1, gQCD] == 1 + 2 * #2 &)];
(* nenexttxbbx *) AddProcess[FeynArtsProcess -> {F[4,{3}], -F[4,{3}]} -> {F[1,{1}], -F[1,{1}], F[3,{3}], -F[3,{3}], V[5]},UnitaryGauge -> False,SelectInterference -> {eQED -> 8},SelectCoupling -> (Exponent[#1, gQCD] == 3 - 2 * #2 || Exponent[#1, gQCD] == 1 + 2 * #2 &)];


(* as a^4 + as a^5 *)

(* eexttxddx   *) AddProcess[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {F[2,{1}], -F[2,{1}], F[3,{3}], -F[3,{3}], V[5]},UnitaryGauge -> False,SelectInterference -> {gQCD -> 2},SelectCoupling -> (Exponent[#1, eQED] == 4 + 2*#2 &)];
(* eexuuxttx   *) AddProcess[FeynArtsProcess -> {F[3,{1}], -F[3,{1}]} -> {F[2,{1}], -F[2,{1}], F[3,{3}], -F[3,{3}], V[5]},UnitaryGauge -> False,SelectInterference -> {gQCD -> 2},SelectCoupling -> (Exponent[#1, eQED] == 4 + 2 * #2 &)];
(* eexttxbbx   *) AddProcess[FeynArtsProcess -> {F[4,{3}], -F[4,{3}]} -> {F[2,{1}], -F[2,{1}], F[3,{3}], -F[3,{3}], V[5]},UnitaryGauge -> False,SelectInterference -> {gQCD -> 2},SelectCoupling -> (Exponent[#1, eQED] == 4 + 2 * #2 &)];

(* nenexttxddx *) AddProcess[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {F[1,{1}], -F[1,{1}], F[3,{3}], -F[3,{3}], V[5]},UnitaryGauge -> False,SelectInterference -> {gQCD -> 2},SelectCoupling -> (Exponent[#1, eQED] == 4 + 2 * #2 &)];
(* nenexuuxttx *) AddProcess[FeynArtsProcess -> {F[3,{1}], -F[3,{1}]} -> {F[1,{1}], -F[1,{1}], F[3,{3}], -F[3,{3}], V[5]},UnitaryGauge -> False,SelectInterference -> {gQCD -> 2},SelectCoupling -> (Exponent[#1, eQED] == 4 + 2 * #2 &)];
(* nenexttxbbx *) AddProcess[FeynArtsProcess -> {F[4,{3}], -F[4,{3}]} -> {F[1,{1}], -F[1,{1}], F[3,{3}], -F[3,{3}], V[5]},UnitaryGauge -> False,SelectInterference -> {gQCD -> 2},SelectCoupling -> (Exponent[#1, eQED] == 4 + 2 * #2 &)];

*)

(* eexttxddx   *) AddTree[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {F[2,{1}], -F[2,{1}], F[3,{3}], -F[3,{3}], V[5]},UnitaryGauge -> False,SelectInterference -> {eQED -> 8},SelectCoupling -> (Exponent[#1, gQCD] == 3 - 2*#2 || Exponent[#1, gQCD] == 1 + 2 * #2 &)];
(* eexuuxttx   *) AddTree[FeynArtsProcess -> {F[3,{1}], -F[3,{1}]} -> {F[2,{1}], -F[2,{1}], F[3,{3}], -F[3,{3}], V[5]},UnitaryGauge -> False,SelectInterference -> {eQED -> 8},SelectCoupling -> (Exponent[#1, gQCD] == 3 - 2 * #2 || Exponent[#1, gQCD] == 1 + 2 * #2 &)];
(* eexttxbbx   *) AddTree[FeynArtsProcess -> {F[4,{3}], -F[4,{3}]} -> {F[2,{1}], -F[2,{1}], F[3,{3}], -F[3,{3}], V[5]},UnitaryGauge -> False,SelectInterference -> {eQED -> 8},SelectCoupling -> (Exponent[#1, gQCD] == 3 - 2 * #2 || Exponent[#1, gQCD] == 1 + 2 * #2 &)];

(* nenexttxddx *) AddTree[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {F[1,{1}], -F[1,{1}], F[3,{3}], -F[3,{3}], V[5]},UnitaryGauge -> False,SelectInterference -> {eQED -> 8},SelectCoupling -> (Exponent[#1, gQCD] == 3 - 2 * #2 || Exponent[#1, gQCD] == 1 + 2 * #2 &)];
(* nenexuuxttx *) AddTree[FeynArtsProcess -> {F[3,{1}], -F[3,{1}]} -> {F[1,{1}], -F[1,{1}], F[3,{3}], -F[3,{3}], V[5]},UnitaryGauge -> False,SelectInterference -> {eQED -> 8},SelectCoupling -> (Exponent[#1, gQCD] == 3 - 2 * #2 || Exponent[#1, gQCD] == 1 + 2 * #2 &)];
(* nenexttxbbx *) AddTree[FeynArtsProcess -> {F[4,{3}], -F[4,{3}]} -> {F[1,{1}], -F[1,{1}], F[3,{3}], -F[3,{3}], V[5]},UnitaryGauge -> False,SelectInterference -> {eQED -> 8},SelectCoupling -> (Exponent[#1, gQCD] == 3 - 2 * #2 || Exponent[#1, gQCD] == 1 + 2 * #2 &)];


(* as a^4 + as a^5 *)

(* eexttxddx   *) AddTree[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {F[2,{1}], -F[2,{1}], F[3,{3}], -F[3,{3}], V[5]},UnitaryGauge -> False,SelectInterference -> {gQCD -> 2},SelectCoupling -> (Exponent[#1, eQED] == 4 + 2*#2 &)];
(* eexuuxttx   *) AddTree[FeynArtsProcess -> {F[3,{1}], -F[3,{1}]} -> {F[2,{1}], -F[2,{1}], F[3,{3}], -F[3,{3}], V[5]},UnitaryGauge -> False,SelectInterference -> {gQCD -> 2},SelectCoupling -> (Exponent[#1, eQED] == 4 + 2 * #2 &)];
(* eexttxbbx   *) AddTree[FeynArtsProcess -> {F[4,{3}], -F[4,{3}]} -> {F[2,{1}], -F[2,{1}], F[3,{3}], -F[3,{3}], V[5]},UnitaryGauge -> False,SelectInterference -> {gQCD -> 2},SelectCoupling -> (Exponent[#1, eQED] == 4 + 2 * #2 &)];

(* nenexttxddx *) AddTree[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {F[1,{1}], -F[1,{1}], F[3,{3}], -F[3,{3}], V[5]},UnitaryGauge -> False,SelectInterference -> {gQCD -> 2},SelectCoupling -> (Exponent[#1, eQED] == 4 + 2 * #2 &)];
(* nenexuuxttx *) AddTree[FeynArtsProcess -> {F[3,{1}], -F[3,{1}]} -> {F[1,{1}], -F[1,{1}], F[3,{3}], -F[3,{3}], V[5]},UnitaryGauge -> False,SelectInterference -> {gQCD -> 2},SelectCoupling -> (Exponent[#1, eQED] == 4 + 2 * #2 &)];
(* nenexttxbbx *) AddTree[FeynArtsProcess -> {F[4,{3}], -F[4,{3}]} -> {F[1,{1}], -F[1,{1}], F[3,{3}], -F[3,{3}], V[5]},UnitaryGauge -> False,SelectInterference -> {gQCD -> 2},SelectCoupling -> (Exponent[#1, eQED] == 4 + 2 * #2 &)];


