
(* description = "pp -> H+jets|QQ -> H+jets|QQ -> H|Higgs production from initial state b/c quarks. Can be used with massless b/c quarks with finite Yukawa coupling." *)

OpenLoopsModel = "SM_yuksel";
SelectInterference = {eQED -> 2};
SetParameters = JoinOptions[{YukB -> 1, YukC -> 1, MB -> MB, YB -> YB, MC -> MC, YC -> YC}];
OnTheFlyMode=0;

(* bbxh *) AddProcess[FeynArtsProcess -> {F[4,{3}], -F[4,{3}]} -> {S[1]}];
(* ccxh *) AddProcess[FeynArtsProcess -> {F[3,{2}], -F[3,{2}]} -> {S[1]}];

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