
(* description = "pp -> H+jets|pp -> H+jets|pp -> H(j)|Loop induced Higgs plus zero and one jet production (massive b quarks allowed, but not as external states; Yukawas allowed)." *)


(* hgg *) AddProcess[FeynArtsProcess -> {V[5], V[5]} -> {S[1]}, SetParameters -> JoinOptions[{YukT -> 1, MB->0, YB->0}], SelectInterference -> {eQED -> 2}, OpenLoopsModel -> "SM_yuksel"];
(* hggg  *) AddProcess[FeynArtsProcess -> {V[5]    ,  V[5]    } -> {S[1], V[5]}, SetParameters -> JoinOptions[{YukT -> 1, MB->0, YB->0}], SelectInterference -> {eQED -> 2}, OpenLoopsModel -> "SM_yuksel"];
(* uuxhg *) AddProcess[FeynArtsProcess -> {F[3,{1}], -F[3,{1}]} -> {S[1], V[5]}, SetParameters -> JoinOptions[{YukT -> 1, MB->0, YB->0}], SelectInterference -> {eQED -> 2}, OpenLoopsModel -> "SM_yuksel"];


(* hgg *) AddProcess[FeynArtsProcess -> {V[5], V[5]} -> {S[1]}, SetParameters -> JoinOptions[{YukB -> 1, YukT -> 1}], SelectInterference -> {eQED -> 2}, OpenLoopsModel -> "SM_yuksel"];
(* hggg  *) AddProcess[FeynArtsProcess -> {V[5]    ,  V[5]    } -> {S[1], V[5]}, SetParameters -> JoinOptions[{YukB -> 1, YukT -> 1}], SelectInterference -> {eQED -> 2}, OpenLoopsModel -> "SM_yuksel"];
(* uuxhg *) AddProcess[FeynArtsProcess -> {F[3,{1}], -F[3,{1}]} -> {S[1], V[5]}, SetParameters -> JoinOptions[{YukB -> 1, YukT -> 1}], SelectInterference -> {eQED -> 2}, OpenLoopsModel -> "SM_yuksel"];

(*
higgsapprox = {
	SelectTreeDiagrams -> (False &),
	SelectLoopDiagrams -> ( NFieldPropagators[S[1],1,"Type"->Internal][##]  &),
	SelectCTDiagrams   -> ( NFieldPropagators[S[1],1,"Type"->Internal][##]  &)
};

(* mmxgg *) AddProcess[FeynArtsProcess -> {V[5], V[5]} -> {F[2,{2}], -F[2,{2}]}, SetParameters -> JoinOptions[{YukT -> 1, MB->0, YB->0, MM -> MM, YukM->1}], SelectInterference -> {eQED -> 4}, Sequence @@ higgsapprox, OpenLoopsModel -> "SM_yuksel"];

(* hggg  *) AddProcess[FeynArtsProcess -> {V[5]    ,  V[5]    } -> {F[2,{2}], -F[2,{2}], V[5]}, SetParameters -> JoinOptions[{YukT -> 1, MB->0, YB->0, MM -> MM, YukM->1}], SelectInterference -> {eQED -> 4}, Sequence @@ higgsapprox, OpenLoopsModel -> "SM_yuksel"];
(* uuxhg *) AddProcess[FeynArtsProcess -> {F[3,{1}], -F[3,{1}]} -> {F[2,{2}], -F[2,{2}], V[5]}, SetParameters -> JoinOptions[{YukT -> 1, MB->0, YB->0, MM -> MM, YukM->1}], SelectInterference -> {eQED -> 4}, Sequence @@ higgsapprox, OpenLoopsModel -> "SM_yuksel"];
	


(* aagg *) AddProcess[FeynArtsProcess -> {V[5], V[5]} -> {V[1],V[1]}, SelectInterference -> {DOI ->0}, SelectCoupling -> (Exponent[#1, gQCD] == 2 && Exponent[#1, eQED] == 1 &), OpenLoopsModel -> "HEFTAAH", Sequence @@ higgsapprox, SetParameters -> JoinOptions[{YukT -> 1, MB->0, YB->0, YC -> 0, YU -> 0, YD -> 0, YS -> 0}]];


(* aaggg *) AddProcess[FeynArtsProcess -> {V[5], V[5]} -> {V[1],V[1], V[5]}, SelectInterference -> {DOI ->0}, SelectCoupling -> (Exponent[#1, gQCD] == 3 && Exponent[#1, eQED] == 1 &), OpenLoopsModel -> "HEFTAAH", Sequence @@ higgsapprox, SetParameters -> JoinOptions[{YukT -> 1, MB->0, YB->0, YC -> 0, YU -> 0, YD -> 0, YS -> 0}]];
(* uuxaag *) AddProcess[FeynArtsProcess -> {F[3,{1}], -F[3,{1}]} -> {V[1],V[1], V[5]}, SelectInterference -> {DOI ->0}, SelectCoupling -> (Exponent[#1, gQCD] == 3 && Exponent[#1, eQED] == 1 &), OpenLoopsModel -> "HEFTAAH", Sequence @@ higgsapprox, SetParameters -> JoinOptions[{YukT -> 1, MB->0, YB->0, YC -> 0, YU -> 0, YD -> 0, YS -> 0}]];

*)


ChannelMap = {
  "ddxhg" -> "uuxhg",
  "bbxhg" -> "uuxhg MB=0 YB=0",
  "mmxddxhg" -> "mmxuuxg",
  "mmxbbxhg" -> "mmxuuxg MB=0 YB=0",
  "ddxaag" -> "uuxaag",
  "bbxaag" -> "uuxaag MB=0 YB=0"
};
