
(* description = "pp -> HV+jets|pp -> HZ/A(->ll)+jets|pp -> Hlljj|Higgs plus off-shell Z/A boson and two jets production with leptonic decays (l+l- and nn) where the Higgs couples to a top loop." *)


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

SetParameters -> JoinOptions[{YukT -> 1, ME -> 0, YE -> 0, nf ->5, MB -> 0, YB -> 0}]

SelectCoupling = (Exponent[#1, gQCD] == 2+2*#2 && Exponent[#1, YT] == 0 + 1*#2 && Exponent[#1, YE] == 0 &);

Approximation="top";


(* eexdddxdxh   *) AddProcess[FeynArtsProcess -> { F[4,{1}],  F[4,{1}]} -> {S[1], F[2,{1}], -F[2,{1}], F[4,{1}],  F[4,{1}]}];
(* eexuuxddxh   *) AddProcess[FeynArtsProcess -> { F[4,{1}], -F[4,{1}]} -> {S[1], F[2,{1}], -F[2,{1}], F[3,{1}], -F[3,{1}]}];
(* eexddxssxh   *) AddProcess[FeynArtsProcess -> { F[4,{1}], -F[4,{1}]} -> {S[1], F[2,{1}], -F[2,{1}], F[4,{2}], -F[4,{2}]}];
(* eexddxhgg    *) AddProcess[FeynArtsProcess -> { F[4,{1}], -F[4,{1}]} -> {S[1], F[2,{1}], -F[2,{1}], V[5]    ,  V[5]    }];
(* eexuuuxuxh   *) AddProcess[FeynArtsProcess -> { F[3,{1}],  F[3,{1}]} -> {S[1], F[2,{1}], -F[2,{1}], F[3,{1}],  F[3,{1}]}];
(* eexuuxccxh   *) AddProcess[FeynArtsProcess -> { F[3,{1}], -F[3,{1}]} -> {S[1], F[2,{1}], -F[2,{1}], F[3,{2}], -F[3,{2}]}];
(* eexuuxhgg    *) AddProcess[FeynArtsProcess -> { F[3,{1}], -F[3,{1}]} -> {S[1], F[2,{1}], -F[2,{1}], V[5]    ,  V[5]    }];
(* nenexdddxdxh *) AddProcess[FeynArtsProcess -> { F[4,{1}],  F[4,{1}]} -> {S[1], F[1,{1}], -F[1,{1}], F[4,{1}],  F[4,{1}]}];
(* nenexuuxddxh *) AddProcess[FeynArtsProcess -> { F[4,{1}], -F[4,{1}]} -> {S[1], F[1,{1}], -F[1,{1}], F[3,{1}], -F[3,{1}]}];
(* nenexddxssxh *) AddProcess[FeynArtsProcess -> { F[4,{1}], -F[4,{1}]} -> {S[1], F[1,{1}], -F[1,{1}], F[4,{2}], -F[4,{2}]}];
(* nenexddxhgg  *) AddProcess[FeynArtsProcess -> { F[4,{1}], -F[4,{1}]} -> {S[1], F[1,{1}], -F[1,{1}], V[5]    ,  V[5]    }];
(* nenexuuuxuxh *) AddProcess[FeynArtsProcess -> { F[3,{1}],  F[3,{1}]} -> {S[1], F[1,{1}], -F[1,{1}], F[3,{1}],  F[3,{1}]}];
(* nenexuuxccxh *) AddProcess[FeynArtsProcess -> { F[3,{1}], -F[3,{1}]} -> {S[1], F[1,{1}], -F[1,{1}], F[3,{2}], -F[3,{2}]}];
(* nenexuuxhgg  *) AddProcess[FeynArtsProcess -> { F[3,{1}], -F[3,{1}]} -> {S[1], F[1,{1}], -F[1,{1}], V[5]    ,  V[5]    }];


ChannelMap = {
  "eexccxddxh"   -> "eexuuxddxh",
  "eexddxbbxh"   -> "eexddxssxh MB=0 YB=0",
  "eexuuxbbxh"   -> "eexuuxddxh MB=0 YB=0",
  "eexbbbxbxh"   -> "eexdddxdxh MB=0 YB=0",
  "eexbbxhgg"    -> "eexddxhgg MB=0 YB=0",
  "nenexccxddxh" -> "nenexuuxddxh",
  "nenexddxbbxh" -> "nenexddxssxh MB=0 YB=0",
  "nenexuuxbbxh" -> "nenexuuxddxh MB=0 YB=0",
  "nenexbbbxbxh" -> "nenexdddxdxh MB=0 YB=0",
  "nenexbbxhgg"  -> "nenexddxhgg MB=0 YB=0"
};
