
(* description = "pp -> HV+jets|pp -> HZ/A(->ll)+jets|pp -> Hll|Higgs plus off-shell Z/A boson production with leptonic decays (l+l- and nn)" *)

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

(* eexddxh   *) AddProcess[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {S[1], F[2,{1}], -F[2,{1}]}];
(* eexuuxh   *) AddProcess[FeynArtsProcess -> {F[3,{1}], -F[3,{1}]} -> {S[1], F[2,{1}], -F[2,{1}]}];
(* eexbbxh   *) AddProcess[FeynArtsProcess -> {F[4,{3}], -F[4,{3}]} -> {S[1], F[2,{1}], -F[2,{1}]}];
(* nenexddxh *) AddProcess[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {S[1], F[1,{1}], -F[1,{1}]}];
(* nenexuuxh *) AddProcess[FeynArtsProcess -> {F[3,{1}], -F[3,{1}]} -> {S[1], F[1,{1}], -F[1,{1}]}];
(* nenexbbxh *) AddProcess[FeynArtsProcess -> {F[4,{3}], -F[4,{3}]} -> {S[1], F[1,{1}], -F[1,{1}]}];

(* eexddxhg   *) AddReal[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {S[1], F[2,{1}], -F[2,{1}], V[5]}];
(* eexuuxhg   *) AddReal[FeynArtsProcess -> {F[3,{1}], -F[3,{1}]} -> {S[1], F[2,{1}], -F[2,{1}], V[5]}];
(* eexbbxhg   *) AddReal[FeynArtsProcess -> {F[4,{3}], -F[4,{3}]} -> {S[1], F[2,{1}], -F[2,{1}], V[5]}];
(* nenexddxhg *) AddReal[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {S[1], F[1,{1}], -F[1,{1}], V[5]}];
(* nenexuuxhg *) AddReal[FeynArtsProcess -> {F[3,{1}], -F[3,{1}]} -> {S[1], F[1,{1}], -F[1,{1}], V[5]}];
(* nenexbbxhg *) AddReal[FeynArtsProcess -> {F[4,{3}], -F[4,{3}]} -> {S[1], F[1,{1}], -F[1,{1}], V[5]}];

ChannelMap = {
  "eexbbxh"   -> "eexddxh MB=0 YB=0",
  "nenexbbxh" -> "nenexddxh MB=0 YB=0",
  (* reals *)
  "eexbbxhg"   -> "eexddxhg MB=0 YB=0",
  "nenexbbxhg" -> "nenexddxhg MB=0 YB=0"
};
