
(* description = "pp -> HV+jets|pp -> HZ+jets|pp -> HZj|Higgs plus on-shell Z boson and one jet production" *)

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

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

(* dddxdxhz *) AddReal[FeynArtsProcess -> {F[4,{1}],  F[4,{1}]} -> {S[1], V[2], F[4,{1}],  F[4,{1}]}];
(* uuxddxhz *) AddReal[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {S[1], V[2], F[3,{1}], -F[3,{1}]}];
(* ddxssxhz *) AddReal[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {S[1], V[2], F[4,{2}], -F[4,{2}]}];
(* ddxhzgg  *) AddReal[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {S[1], V[2], V[5]    ,  V[5]    }];
(* uuuxuxhz *) AddReal[FeynArtsProcess -> {F[3,{1}],  F[3,{1}]} -> {S[1], V[2], F[3,{1}],  F[3,{1}]}];
(* uuxccxhz *) AddReal[FeynArtsProcess -> {F[3,{1}], -F[3,{1}]} -> {S[1], V[2], F[3,{2}], -F[3,{2}]}];
(* uuxhzgg  *) AddReal[FeynArtsProcess -> {F[3,{1}], -F[3,{1}]} -> {S[1], V[2], V[5]    ,  V[5]    }];

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

ChannelMap = {
  "bbxhzg" -> "ddxhzg MB=0",
  (* reals *)
  "ccxddxhz" -> "uuxddxhz",
  "ddxbbxhz" -> "ddxssxhz MB=0",
  "uuxbbxhz" -> "uuxddxhz MB=0",
  "bbbxbxhz" -> "dddxdxhz MB=0",
  "bbxhzgg"  -> "ddxhzgg MB=0"
};
