
(* description = "pp -> H+jets|pp -> Hbb+jets|pp -> Hbb|Higgs production plus two jets with b/c quarks in the initial and/or final state. 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, MC -> MC, YB -> YB, YC -> YC}];


(* 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}]}];

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

(* ccxbbxh *) AddProcess[FeynArtsProcess -> {F[3,{2}], -F[3,{2}]} -> {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]}];

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

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


ChannelMap = {
  "uuxbbxh" -> "ddxbbxh",
  "uuxccxh" -> "ccxddxh[3,4,1,2,5]",
  (* real *)
  "uuxbbxhg" -> "ddxbbxhg",
  "uuxccxhg" -> "ccxddxhg[3,4,1,2,5,6]",
};


