
(* description = "pp -> V+jets|pp -> A+jets|pp -> Ajj|Loop-squared photon plus two jets production." *)

SelectInterference = {eQED -> 2};

loop2approx = {
  SelectLoopDiagrams -> (ContainsFermionLoop[##] && ParticlesOnLoop[V[1],1][##] && NLegsOnLoop[4][##] &),
  SelectTreeDiagrams -> (False &),
  SelectCTDiagrams -> NFieldPropagators[V[5],1]
};

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

(* uuxagg *) AddProcess[FeynArtsProcess -> {F[3,{1}], -F[3,{1}]} -> {V[1], V[5], V[5]}, Sequence @@ loop2approx];
(* bbxagg *) AddProcess[FeynArtsProcess -> {F[4,{3}], -F[4,{3}]} -> {V[1], V[5], V[5]}, Sequence @@ loop2approx];

ChannelMap = {
  "ddxagg" -> "uuxagg",
  "bbxagg" -> "uuxagg MB=0"
};
