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

OpenLoopsModel = "SM_yuksel";
SelectInterference = {eQED -> 4};
SetParameters = JoinOptions[{YukB -> 1, YukT -> 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]    }];


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

(* uddxdxhw  *) AddReal[FeynArtsProcess -> { F[4,{1}], -F[4,{1}]} -> {S[1], -V[3], F[4,{1}], -F[3,{1}]}];
(* cddxsxhw  *) AddReal[FeynArtsProcess -> { F[4,{1}], -F[4,{1}]} -> {S[1], -V[3], F[4,{2}], -F[3,{2}]}];
(* uuuxdxhw  *) AddReal[FeynArtsProcess -> {-F[4,{1}],  F[3,{1}]} -> {S[1], -V[3], F[3,{1}], -F[3,{1}]}];
(* udxhwgg   *) AddReal[FeynArtsProcess -> {-F[4,{1}],  F[3,{1}]} -> {S[1], -V[3], V[5]    ,  V[5]    }];
(* udxbbxhw  *) AddReal[FeynArtsProcess -> {-F[4,{1}],  F[3,{1}]} -> {S[1], -V[3], F[4,{2}], -F[4,{2}]}];

(* uxdddxhwx *) AddReal[FeynArtsProcess -> { F[4,{1}],  F[4,{1}]} -> {S[1],  V[3], F[4,{1}],  F[3,{1}]}];
(* cxddxshwx *) AddReal[FeynArtsProcess -> { F[4,{1}], -F[4,{1}]} -> {S[1],  V[3], F[3,{2}], -F[4,{2}]}];
(* uuxuxdhwx *) AddReal[FeynArtsProcess -> { F[4,{1}],  F[3,{1}]} -> {S[1],  V[3], F[3,{1}],  F[3,{1}]}];
(* uxdhwxgg  *) AddReal[FeynArtsProcess -> { F[4,{1}], -F[3,{1}]} -> {S[1],  V[3], V[5]    ,  V[5]    }];
(* uxdbbxhwx *) AddReal[FeynArtsProcess -> { F[4,{1}], -F[3,{1}]} -> {S[1],  V[3], F[4,{2}], -F[4,{2}]}];

ChannelMap = {
(*Wj*)
(* reals *)
  "uccxdxhw"  -> "cddxsxhw",
  "udxbbxhw"  -> "cddxsxhw[1,4,2,3,5,6] MB=0 YB=0",
  "uxccxdhwx" -> "cxddxshwx",
  "uxdbbxhwx" -> "cxddxshwx[1,4,2,3,5,6] MB=0 YB=0",
(*Zj*)
  "bbxhzg" -> "ddxhzg MB=0 YB=0",
  (* reals *)
  "ccxddxhz" -> "uuxddxhz",
  "ddxbbxhz" -> "ddxssxhz MB=0 YB=0",
  "uuxbbxhz" -> "uuxddxhz MB=0 YB=0",
  "bbbxbxhz" -> "dddxdxhz MB=0 YB=0",
  "bbxhzgg"  -> "ddxhzgg MB=0 YB=0"
};
