
(* description = "pp -> HV+jets|pp -> HW(->ln)+jets|pp -> Hlnjj|Higgs plus off-shell W+/W- boson and two jets production with leptonic decay" *)

OpenLoopsModel = "SM_yuksel";
SelectInterference = {eQED -> 6};
CreateTopologiesOptions = {ExcludeTopologies -> {Snails, WFCorrectionCTs, TadpoleCTs, Loops[7]}};
SetParameters = JoinOptions[{ME -> 0, YE -> 0, YukB -> 1, YukT -> 1}];

(* nexeuddxdxh *) AddProcess[FeynArtsProcess -> { F[4,{1}], -F[4,{1}]} -> {S[1], F[1,{1}], -F[2,{1}], F[4,{1}], -F[3,{1}]}];
(* nexecddxsxh *) AddProcess[FeynArtsProcess -> { F[4,{1}], -F[4,{1}]} -> {S[1], F[1,{1}], -F[2,{1}], F[4,{2}], -F[3,{2}]}];
(* nexeuuuxdxh *) AddProcess[FeynArtsProcess -> {-F[4,{1}],  F[3,{1}]} -> {S[1], F[1,{1}], -F[2,{1}], F[3,{1}], -F[3,{1}]}];
(* nexeudxhgg  *) AddProcess[FeynArtsProcess -> {-F[4,{1}],  F[3,{1}]} -> {S[1], F[1,{1}], -F[2,{1}], V[5]    ,  V[5]    }];
(* nexeudxbbxh *) AddProcess[FeynArtsProcess -> {-F[4,{1}],  F[3,{1}]} -> {S[1], F[1,{1}], -F[2,{1}], F[4,{2}], -F[4,{2}]}];

(* neexuxdddxh *) AddProcess[FeynArtsProcess -> { F[4,{1}],  F[4,{1}]} -> {S[1], F[2,{1}], -F[1,{1}], F[4,{1}],  F[3,{1}]}];
(* neexcxddxsh *) AddProcess[FeynArtsProcess -> { F[4,{1}], -F[4,{1}]} -> {S[1], F[2,{1}], -F[1,{1}], F[3,{2}], -F[4,{2}]}];
(* neexuuxuxdh *) AddProcess[FeynArtsProcess -> { F[4,{1}],  F[3,{1}]} -> {S[1], F[2,{1}], -F[1,{1}], F[3,{1}],  F[3,{1}]}];
(* neexuxdhgg  *) AddProcess[FeynArtsProcess -> { F[4,{1}], -F[3,{1}]} -> {S[1], F[2,{1}], -F[1,{1}], V[5]    ,  V[5]    }];
(* neexuxdbbxh *) AddProcess[FeynArtsProcess -> { F[4,{1}], -F[3,{1}]} -> {S[1], F[2,{1}], -F[1,{1}], F[4,{2}], -F[4,{2}]}];

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

(* nexeuddxdxhg *) AddReal[FeynArtsProcess -> { F[4,{1}], -F[4,{1}]} -> {S[1], F[1,{1}], -F[2,{1}], V[5], F[4,{1}], -F[3,{1}]}];
(* nexecddxsxhg *) AddReal[FeynArtsProcess -> { F[4,{1}], -F[4,{1}]} -> {S[1], F[1,{1}], -F[2,{1}], V[5], F[4,{2}], -F[3,{2}]}];
(* nexeuuuxdxhg *) AddReal[FeynArtsProcess -> {-F[4,{1}],  F[3,{1}]} -> {S[1], F[1,{1}], -F[2,{1}], V[5], F[3,{1}], -F[3,{1}]}];
(* nexeudxhggg  *) AddReal[FeynArtsProcess -> {-F[4,{1}],  F[3,{1}]} -> {S[1], F[1,{1}], -F[2,{1}], V[5], V[5]    ,  V[5]    }];
(* nexeudxbbxhg *) AddReal[FeynArtsProcess -> {-F[4,{1}],  F[3,{1}]} -> {S[1], F[1,{1}], -F[2,{1}], V[5], F[4,{2}], -F[4,{2}]}];

(* neexuxdddxhg *) AddReal[FeynArtsProcess -> { F[4,{1}],  F[4,{1}]} -> {S[1], F[2,{1}], -F[1,{1}], V[5], F[4,{1}],  F[3,{1}]}];
(* neexcxddxshg *) AddReal[FeynArtsProcess -> { F[4,{1}], -F[4,{1}]} -> {S[1], F[2,{1}], -F[1,{1}], V[5], F[3,{2}], -F[4,{2}]}];
(* neexuuxuxdhg *) AddReal[FeynArtsProcess -> { F[4,{1}],  F[3,{1}]} -> {S[1], F[2,{1}], -F[1,{1}], V[5], F[3,{1}],  F[3,{1}]}];
(* neexuxdhggg  *) AddReal[FeynArtsProcess -> { F[4,{1}], -F[3,{1}]} -> {S[1], F[2,{1}], -F[1,{1}], V[5], V[5]    ,  V[5]    }];
(* neexuxdbbxhg *) AddReal[FeynArtsProcess -> { F[4,{1}], -F[3,{1}]} -> {S[1], F[2,{1}], -F[1,{1}], V[5], F[4,{2}], -F[4,{2}]}];

ChannelMap = {
  "nexeuccxdxh" -> "nexecddxsxh",
  "nexeudxbbxh" -> "nexecddxsxh[1,2,3,6,4,5,7] MB=0 YB=0",
  "neexuxccxdh" -> "neexcxddxsh",
  "neexuxdbbxh" -> "neexcxddxsh[1,2,3,6,4,5,7] MB=0 YB=0",
  (* reals *)
  "nexeuccxdxhg" -> "nexecddxsxhg",
  "nexeudxbbxhg" -> "nexecddxsxhg[1,2,3,6,4,5,7] MB=0 YB=0",
  "neexuxccxdhg" -> "neexcxddxshg",
  "neexuxdbbxhg" -> "neexcxddxshg[1,2,3,6,4,5,7] MB=0 YB=0"
};
