
(* description = "pp -> H+jets|pp -> H+jets|pp -> H(j)|Loop induced Higgs plus zero and one jet production (massive b quarks allowed, but not as external states; Yukawas allowed)." *)

OpenLoopsModel = "SM_yuksel";
SelectInterference = {eQED -> 2};


(* hgg *) AddProcess[FeynArtsProcess -> {V[5], V[5]} -> {S[1]}, SetParameters -> JoinOptions[{YukT -> 1, MB->0, YB->0}]];
(* hggg  *) AddProcess[FeynArtsProcess -> {V[5]    ,  V[5]    } -> {S[1], V[5]}, SetParameters -> JoinOptions[{YukT -> 1, MB->0, YB->0}]];
(* uuxhg *) AddProcess[FeynArtsProcess -> {F[3,{1}], -F[3,{1}]} -> {S[1], V[5]}, SetParameters -> JoinOptions[{YukT -> 1, MB->0, YB->0}]];


(* hgg *) AddProcess[FeynArtsProcess -> {V[5], V[5]} -> {S[1]}, SetParameters -> JoinOptions[{YukB -> 1, YukT -> 1}]];
(* hggg  *) AddProcess[FeynArtsProcess -> {V[5]    ,  V[5]    } -> {S[1], V[5]}, SetParameters -> JoinOptions[{YukB -> 1, YukT -> 1}]];
(* uuxhg *) AddProcess[FeynArtsProcess -> {F[3,{1}], -F[3,{1}]} -> {S[1], V[5]}, SetParameters -> JoinOptions[{YukB -> 1, YukT -> 1}]];

ChannelMap = {
  "ddxhg" -> "uuxhg",
  "bbxhg" -> "uuxhg MB=0 YB=0"
};
