(* description = "pp -> HtV+jets|pp -> HtW(->ln)+jets|pp -> Htln|Higgs plus Single top plus on-shell W+/W- production with leptonic decays. Real corrections in DR scheme." *)

SelectInterference = {eQED -> 6};
UnitaryGauge = True;

SetParameters = JoinOptions[{ME -> 0, MB -> 0}];



(*
(* neextxbhgg *) AddReal[FeynArtsProcess -> {F[1,{1}], -F[2,{1}]} -> { -F[4,{3}],  F[3,{3}], S[1], V[5], V[5]}];
(* neextxbhgg *) AddReal[FeynArtsProcess -> {F[1,{1}], -F[2,{1}]} -> { -F[4,{3}],  F[3,{3}], S[1], V[5], V[5]}, SelectTreeDiagrams -> (HasResonance[F[3,{3}], {1,2,3,5}][##] || HasResonance[F[3,{3}], {1,2,3}][##] &)];
*)




(* neextxbhg *) AddProcess[FeynArtsProcess -> { F[4,{3}], V[5]} -> { F[2,{1}], -F[1,{1}],  F[3,{3}], S[1]}, UnitaryGauge -> False];
(* nexetbxhg *) AddProcess[FeynArtsProcess -> {-F[4,{3}], V[5]} -> {-F[2,{1}],  F[1,{1}], -F[3,{3}], S[1]}, UnitaryGauge -> False];

(* neextxbhgg *) AddReal[FeynArtsProcess -> {F[1,{1}], -F[2,{1}]} -> {F[3,{3}],  -F[4,{3}],  S[1], V[5], V[5]}, SelectTreeDiagrams -> (Not[HasResonance[F[3,{3}], {1,2,4}][##] || HasResonance[F[3,{3}], {1,2,4,5}][##]] &), UnitaryGauge -> True];
(* neexuuxtxbh *) AddReal[FeynArtsProcess -> { F[4,{3}], F[3,{1}]} -> { F[2,{1}], -F[1,{1}],  F[3,{3}], S[1], F[3,{1}]}, UnitaryGauge -> True];
(* neextxddxbhgg *) AddReal[FeynArtsProcess -> { F[4,{3}], F[4,{1}]} -> { F[2,{1}], -F[1,{1}],  F[3,{3}], S[1], F[4,{1}]}, UnitaryGauge -> True];

(* nexetbxhgg *) AddReal[FeynArtsProcess -> {-F[1,{1}], F[2,{1}]} -> {-F[3,{3}], F[4,{3}], S[1], V[5], V[5]}, SelectTreeDiagrams -> (Not[HasResonance[F[3,{3}], {1,2,4}][##] || HasResonance[F[3,{3}], {1,2,4,5}][##]] &), UnitaryGauge -> True];
(* nexeuuxtbxh *) AddReal[FeynArtsProcess -> {-F[4,{3}], F[3,{1}]} -> {-F[2,{1}],  F[1,{1}], -F[3,{3}], S[1], F[3,{1}]}, UnitaryGauge -> True];
(* nexetddxbxh *) AddReal[FeynArtsProcess -> {-F[4,{3}], F[4,{1}]} -> {-F[2,{1}],  F[1,{1}], -F[3,{3}], S[1], F[4,{1}]}, UnitaryGauge -> True];