
(* description = "pp -> H+jets|QQ -> H+jets|QQ -> Hj|Higgs production from (at least one) initial state b/c quark(s). 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;

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

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

(* ccxhgg  *) AddReal[FeynArtsProcess -> {F[3,{2}], -F[3,{2}]} -> {S[1], V[5]    ,  V[5]    }];
(* ddxccxh *) AddReal[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {S[1], F[3,{2}], -F[3,{2}]}];
(* cccxcxh *) AddReal[FeynArtsProcess -> {F[3,{2}], -F[3,{2}]} -> {S[1], F[3,{2}], -F[3,{2}]}];
(* ccxbbxh *) AddReal[FeynArtsProcess -> {F[3,{2}], -F[3,{2}]} -> {S[1], F[4,{3}], -F[4,{3}]}];



ChannelMap = {
  (* real *)
  "uuxbbxh" -> "ddxbbxh",
  "uuxccxh" -> "ccxddxh[3,4,1,2,5]",
  "ccxbbxh" -> "ddxddxh MC=0 YC=0"
};
