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

OpenLoopsModel = "HEFT";
SelectInterference = {eQED -> 2, DOI -> 0};
OnTheFlyMode = 0;

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

(* hggg  *) AddProcess[FeynArtsProcess -> {V[5]    ,  V[5]    } -> {S[1], V[5]}];
(* ddxhg *) AddProcess[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {S[1], V[5]}];

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

ChannelMap = {
  "uuxhg" -> "ddxhg",
  "bbxhg" -> "ddxhg MB=0",
  (* reals *)
  "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"
};
