(* description = "pp -> jets|pp -> jets|pp -> jjj|All gluon gg ->gg(g)." *)

SelectInterference = {eQED -> 0};

OpenLoopsModel = "SM_no4g";
SetParameters = JoinOptions[{MB -> 0, nf -> 5}];

InsertFieldsOptions = {Restrictions -> {ExcludeParticles -> {F[3,{3}]}, NoQuarkMixing}}



(*
(* gggg   *) AddProcess[FeynArtsProcess -> {V[5]    ,  V[5]    } -> {V[5]    ,  V[5]    }, Sequence @@ nf0, OnTheFlyMode -> 1];
(* ggggg   *) AddProcess[FeynArtsProcess -> {V[5]    ,  V[5]    } -> {V[5]    ,  V[5],  V[5]   }, Sequence @@ nf0, OnTheFlyMode -> 1];
*)


(* gggg   *) AddProcess[FeynArtsProcess -> {V[5]    ,  V[5]    } -> {V[5]    ,  V[5]    }, OnTheFlyMode -> 0];
(* ggggg   *) AddProcess[FeynArtsProcess -> {V[5]    ,  V[5]    } -> {V[5]    ,  V[5],  V[5]   }, OnTheFlyMode -> 0];


(* qqxgg   *) AddProcess[FeynArtsProcess -> {V[5]    ,  V[5]    } -> {F[3,{1}], -F[3,{1}]}, OnTheFlyMode -> 0];
(* qqxggg   *) AddProcess[FeynArtsProcess -> {V[5]    ,  V[5]    } -> {F[3,{1}], -F[3,{1}],  V[5]   }, OnTheFlyMode -> 0];
(* qqqxqxg   *) AddProcess[FeynArtsProcess -> {F[3,{1}], -F[3,{1}]   } -> {F[3,{1}], -F[3,{1}],  V[5]   }, OnTheFlyMode -> 0];