
(* description = "pp -> HV+jets|pp -> HV+jets|pp -> HV|Higgs plus on-shell boson (Z/W) production"*)

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

(* HW *)

(* udxhw  *) AddProcess[FeynArtsProcess -> {-F[4,{1}],  F[3,{1}]} -> {S[1], -V[3]}];
(* uxdhwx *) AddProcess[FeynArtsProcess -> { F[4,{1}], -F[3,{1}]} -> {S[1],  V[3]}];

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


(* HZ *)

(* ddxhz *) AddProcess[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {S[1], V[2]}];
(* uuxhz *) AddProcess[FeynArtsProcess -> {F[3,{1}], -F[3,{1}]} -> {S[1], V[2]}];
(* bbxhz *) AddProcess[FeynArtsProcess -> {F[4,{3}], -F[4,{3}]} -> {S[1], V[2]}];

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


ChannelMap = {
  "bbxhz" -> "ddxhz MB=0 YB=0",
  (* real *)
  "bbxhzg" -> "ddxhzg MB=0 YB=0"
};
