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

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

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

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

(* old Sherpa *)
(* udxssxhw  *) AddProcess[FeynArtsProcess -> { F[4,{2}], -F[4,{2}]} -> {S[1], -V[3], F[4,{1}], -F[3,{1}]}];
(* uxdssxhwx *) AddProcess[FeynArtsProcess -> { F[4,{2}], -F[4,{2}]} -> {S[1],  V[3], F[3,{1}], -F[4,{1}]}];

(* uddxdxhwg  *) AddReal[FeynArtsProcess -> { F[4,{1}], -F[4,{1}]} -> {S[1], -V[3], V[5], F[4,{1}], -F[3,{1}]}];
(* cddxsxhwg  *) AddReal[FeynArtsProcess -> { F[4,{1}], -F[4,{1}]} -> {S[1], -V[3], V[5], F[4,{2}], -F[3,{2}]}];
(* uuuxdxhwg  *) AddReal[FeynArtsProcess -> {-F[4,{1}],  F[3,{1}]} -> {S[1], -V[3], V[5], F[3,{1}], -F[3,{1}]}];
(* udxhwggg   *) AddReal[FeynArtsProcess -> {-F[4,{1}],  F[3,{1}]} -> {S[1], -V[3], V[5], V[5]    ,  V[5]    }];
(* udxbbxhwg  *) AddReal[FeynArtsProcess -> {-F[4,{1}],  F[3,{1}]} -> {S[1], -V[3], V[5], F[4,{2}], -F[4,{2}]}];

(* uxdddxhwxg *) AddReal[FeynArtsProcess -> { F[4,{1}],  F[4,{1}]} -> {S[1],  V[3], V[5], F[4,{1}],  F[3,{1}]}];
(* cxddxshwxg *) AddReal[FeynArtsProcess -> { F[4,{1}], -F[4,{1}]} -> {S[1],  V[3], V[5], F[3,{2}], -F[4,{2}]}];
(* uuxuxdhwxg *) AddReal[FeynArtsProcess -> { F[4,{1}],  F[3,{1}]} -> {S[1],  V[3], V[5], F[3,{1}],  F[3,{1}]}];
(* uxdhwxggg  *) AddReal[FeynArtsProcess -> { F[4,{1}], -F[3,{1}]} -> {S[1],  V[3], V[5], V[5]    ,  V[5]    }];
(* uxdbbxhwxg *) AddReal[FeynArtsProcess -> { F[4,{1}], -F[3,{1}]} -> {S[1],  V[3], V[5], F[4,{2}], -F[4,{2}]}];

ChannelMap = {
  "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",
  (* reals *)
  "uccxdxhwg"  -> "cddxsxhwg",
  "udxbbxhwg"  -> "cddxsxhwg[1,4,2,3,5,6] MB=0 YB=0",
  "uxccxdhwxg" -> "cxddxshwxg",
  "uxdbbxhwxg" -> "cxddxshwxg[1,4,2,3,5,6] MB=0 YB=0"
};
