
(* description = "pp -> H+jets|pp -> Hbb+jets|pp -> Hbb|Higgs production with b quarks in the initial and/or final state. Can be used with massless b quarks with finite Yukawa coupling." *)

OpenLoopsModel = "SM_yuksel";
SelectInterference = {eQED -> 2};
SetParameters = JoinOptions[{YukB -> 1}];

(* bbxhgg  *) AddProcess[FeynArtsProcess -> {F[4,{3}], -F[4,{3}]} -> {S[1], V[5]    ,  V[5]    }];
(* ddxbbxh *) AddProcess[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {S[1], F[4,{3}], -F[4,{3}]}];
(* bbbxbxh *) AddProcess[FeynArtsProcess -> {F[4,{3}], -F[4,{3}]} -> {S[1], F[4,{3}], -F[4,{3}]}];

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

ChannelMap = {
  "uuxbbxh" -> "ddxbbxh",
  (* real *)
  "uuxbbxhg" -> "ddxbbxhg"
};
