
(* description = "pp -> HH+jets|pp -> HH+jets|pp -> HHj|Loop-squared Higgs pair production with a jet." *)

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

DefaultProcessOptions = {SelectInterference -> {eQED -> 4}};

(* hhggg  *) AddProcess[FeynArtsProcess -> {V[5]    ,  V[5]    } -> {S[1], S[1], V[5]}];
(* uuxhhg *) AddProcess[FeynArtsProcess -> {F[3,{1}], -F[3,{1}]} -> {S[1], S[1], V[5]}];

(*
(* only diagrams with 3h coupling *)
(* hhggg  *) AddProcess[FeynArtsProcess -> {V[5], V[5]} -> {S[1], S[1], V[5]}, SelectInterference -> {eQED -> 4, Yuk -> 2}, Approximation -> "only3h"];
(* uuxhhg *) AddProcess[FeynArtsProcess -> {F[3,{1}], -F[3,{1}]} -> {S[1], S[1], V[5]}, SelectInterference -> {eQED -> 4, Yuk -> 2}, Approximation -> "only3h"];

(* interference of diagrams with 3h coupling and diagrams without *)
(* hhggg  *) AddProcess[FeynArtsProcess -> {V[5], V[5]} -> {S[1], S[1], V[5]}, SelectInterference -> {eQED -> 4, Yuk -> 3}, Approximation -> "interf3h"];
(* uuxhhg *) AddProcess[FeynArtsProcess -> {F[3,{1}], -F[3,{1}]} -> {S[1], S[1], V[5]}, SelectInterference -> {eQED -> 4, Yuk -> 3}, Approximation -> "interf3h"];

(* only diagrams without 3h coupling *)
(* hhggg  *) AddProcess[FeynArtsProcess -> {V[5], V[5]} -> {S[1], S[1], V[5]}, SelectInterference -> {eQED -> 4, Yuk -> 4}, Approximation -> "no3h"];
(* uuxhhg *) AddProcess[FeynArtsProcess -> {F[3,{1}], -F[3,{1}]} -> {S[1], S[1], V[5]}, SelectInterference -> {eQED -> 4, Yuk -> 4}, Approximation -> "no3h"];
*)

ChannelMap = {
  "ddxhhg" -> "uuxhhg",
  "bbxhhg" -> "uuxhhg MB=0 YB=0"
};
