
(* description = "pp -> tH+jets|pp -> tH+jets|pp -> tHj|Single top plus Higgs boson production and one jet." *)

OpenLoopsModel = "SM_yuksel";
SetParameters = JoinOptions[{YukB -> 1, YukT -> 1}];
SelectInterference = {eQED -> 6};
UnitaryGauge = False;
ColourCorrelations = True;

(* uxtdbxh *) AddProcess[FeynArtsProcess -> { F[4,{1}], -F[4,{3}]} -> {S[1], F[3,{1}], -F[3,{3}]}];
(* utxdxbh *) AddProcess[FeynArtsProcess -> {-F[4,{1}],  F[4,{3}]} -> {S[1], F[3,{3}], -F[3,{1}]}];

(* uxtdbxhg *) AddReal[FeynArtsProcess -> { F[4,{1}], -F[4,{3}]} -> {S[1], F[3,{1}], -F[3,{3}], V[5]}];
(* utxdxbhg *) AddReal[FeynArtsProcess -> {-F[4,{1}],  F[4,{3}]} -> {S[1], F[3,{3}], -F[3,{1}], V[5]}];
