
(* description = "pp -> HH+jets|pp -> HH+jets|pp -> HH|Loop-induced Higgs pair production." *)

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

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

(* hhgg *) AddProcess[FeynArtsProcess -> {V[5], V[5]} -> {S[1], S[1]}];

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

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

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