
(* description = "pp -> VV+jets|pp -> AA+jets|pp -> AAjj|Loop-squared photon pair plus two jets production." *)

SelectInterference = {eQED -> 4};
ColourCorrelations = False;

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

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

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

ChannelMap = {
  "ddxaagg" -> "uuxaagg",
  "bbxaagg" -> "uuxaagg MB=0"
};
