
(* description = "pp -> H+jets|pp -> H+jets|pp -> Hjj|Higgs plus two jets production in the Higgs effective field theory." *)

OpenLoopsModel = "HEFT";
SelectInterference = {eQED -> 2, DOI -> 0};
SetParameters = JoinOptions[{MB -> 0, YB -> 0}];
OnTheFlyMode = 0;

(* dddxdxh *) AddProcess[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {S[1], F[4,{1}], -F[4,{1}]}];
(* uuxddxh *) AddProcess[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {S[1], F[3,{1}], -F[3,{1}]}];
(* ddxhgg  *) AddProcess[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {S[1], V[5]    ,  V[5]    }];
(* hgggg   *) AddProcess[FeynArtsProcess -> {V[5]    ,  V[5]    } -> {S[1], V[5]    ,  V[5]    }];

(* dddxdxhg *) AddReal[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {S[1], F[4,{1}], -F[4,{1}], V[5]}];
(* uuxddxhg *) AddReal[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {S[1], F[3,{1}], -F[3,{1}], V[5]}];
(* ddxhggg  *) AddReal[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {S[1], V[5]    ,  V[5]    , V[5]}];
(* hggggg   *) AddReal[FeynArtsProcess -> {V[5]    ,  V[5]    } -> {S[1], V[5]    ,  V[5]    , V[5]}];

ChannelMap = {
  "ddxbbxh" -> "uuxddxh MB=0",
  "uuxbbxh" -> "uuxddxh MB=0",
  "bbbxbxh" -> "dddxdxh MB=0",
  "bbxhgg"  -> "ddxhgg MB=0",
  "ddxssxh" -> "uuxddxh",
  "uuuxuxh" -> "dddxdxh",
  "uuxccxh" -> "uuxddxh",
  "ccxddxh" -> "uuxddxh",
  "uuxhgg"  -> "ddxhgg",
  (* reals *)
  "ddxbbxhg" -> "uuxddxhg MB=0",
  "uuxbbxhg" -> "uuxddxhg MB=0",
  "bbbxbxhg" -> "dddxdxhg MB=0",
  "bbxhggg"  -> "ddxhggg MB=0",
  "ddxssxhg" -> "uuxddxhg",
  "uuuxuxhg" -> "dddxdxhg",
  "uuxccxhg" -> "uuxddxhg",
  "ccxddxhg" -> "uuxddxhg",
  "uuxhggg"  -> "ddxhggg"
};
