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

OpenLoopsModel = "SM_yuksel";
SelectInterference = {eQED -> 4};

(* hhgggg   *) AddProcess[FeynArtsProcess -> {V[5]    ,  V[5]    } -> {S[1], S[1], V[5]    ,  V[5]    }, SetParameters -> JoinOptions[{YukT -> 1, MB -> 0, YB -> 0}]];
(* uuxhhgg  *) AddProcess[FeynArtsProcess -> {F[3,{1}], -F[3,{1}]} -> {S[1], S[1], V[5]    ,  V[5]    }, SetParameters -> JoinOptions[{YukT -> 1, MB -> 0, YB -> 0}]];
(* uuuxuxhh *) AddProcess[FeynArtsProcess -> {F[3,{1}], -F[3,{1}]} -> {S[1], S[1], F[3,{1}], -F[3,{1}]}, SetParameters -> JoinOptions[{YukT -> 1, MB -> 0, YB -> 0}]];
(* uuxddxhh *) AddProcess[FeynArtsProcess -> {F[3,{1}], -F[3,{1}]} -> {S[1], S[1], F[4,{1}], -F[4,{1}]}, SetParameters -> JoinOptions[{YukT -> 1, MB -> 0, YB -> 0}]];


(* hhgggg   *) AddProcess[FeynArtsProcess -> {V[5]    ,  V[5]    } -> {S[1], S[1], V[5]    ,  V[5]    }, SetParameters -> JoinOptions[{YukB -> 1, YukT -> 1}]];
(* uuxhhgg  *) AddProcess[FeynArtsProcess -> {F[3,{1}], -F[3,{1}]} -> {S[1], S[1], V[5]    ,  V[5]    }, SetParameters -> JoinOptions[{YukB -> 1, YukT -> 1}]];
(* uuuxuxhh *) AddProcess[FeynArtsProcess -> {F[3,{1}], -F[3,{1}]} -> {S[1], S[1], F[3,{1}], -F[3,{1}]}, SetParameters -> JoinOptions[{YukB -> 1, YukT -> 1}]];
(* uuxddxhh *) AddProcess[FeynArtsProcess -> {F[3,{1}], -F[3,{1}]} -> {S[1], S[1], F[4,{1}], -F[4,{1}]}, SetParameters -> JoinOptions[{YukB -> 1, YukT -> 1}]];



ChannelMap = {
  "ddxhhgg"  -> "uuxhhgg",
  "dddxdxhh" -> "uuuxuxhh",
  "ddxssxhh" -> "uuxddxhh",
  "uuxccxhh" -> "uuxddxhh",
  "ccxddxhh" -> "uuxddxhh",
  "bbxhhgg"  -> "uuxhhgg MB=0 YB=0",
  "ddxbbxhh" -> "uuxddxhh MB=0 YB=0",
  "uuxbbxhh" -> "uuxddxhh MB=0 YB=0",
  "bbbxbxhh" -> "uuuxuxhh MB=0 YB=0"
}
