
(* description = "pp -> VV+jets|pp -> Z/A(->ll)A+jets|gg -> llAj|Loop-squared matrix elements for off-shell Z/A boson with leptonic decay (l+l- and nn) plus on-shell photon and jet production." *)

SelectInterference = {eQED -> 6};
CreateTopologiesOptions = {ExcludeTopologies -> {Snails, WFCorrectionCTs, TadpoleCTs, Loops[6]}};
SetParameters = JoinOptions[{ME -> 0, YE -> 0}];

spliqq = {
  SelectTreeDiagrams -> (False &),
  SelectLoopDiagrams -> ContainsFermionLoop,
  SelectCTDiagrams -> NFieldPropagators[V[5],1]
  (* Approximation where the vector boson only couples to the loop: *)
  (* SelectLoopDiagrams -> (ContainsFermionLoop[##] && NFieldPropagators[F[3|4],0,"Type"->Internal][##] &) *)
};

(* eexaggg   *) AddProcess[FeynArtsProcess -> {V[5], V[5]} -> {F[2,{1}], -F[2,{1}], V[1], V[5]}];
(* nenexaggg *) AddProcess[FeynArtsProcess -> {V[5], V[5]} -> {F[1,{1}], -F[1,{1}], V[1], V[5]}];

(* eexuuxag   *) AddProcess[FeynArtsProcess -> {F[3,{1}], -F[3,{1}]} -> {F[2,{1}], -F[2,{1}], V[1], V[5]}, Sequence @@ spliqq];
(* nenexuuxag *) AddProcess[FeynArtsProcess -> {F[3,{1}], -F[3,{1}]} -> {F[1,{1}], -F[1,{1}], V[1], V[5]}, Sequence @@ spliqq];
(* eexddxag   *) AddProcess[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {F[2,{1}], -F[2,{1}], V[1], V[5]}, Sequence @@ spliqq];
(* nenexddxag *) AddProcess[FeynArtsProcess -> {F[4,{1}], -F[4,{1}]} -> {F[1,{1}], -F[1,{1}], V[1], V[5]}, Sequence @@ spliqq];
(* eexbbxag   *) AddProcess[FeynArtsProcess -> {F[4,{3}], -F[4,{3}]} -> {F[2,{1}], -F[2,{1}], V[1], V[5]}, Sequence @@ spliqq];
(* nenexbbxag *) AddProcess[FeynArtsProcess -> {F[4,{3}], -F[4,{3}]} -> {F[1,{1}], -F[1,{1}], V[1], V[5]}, Sequence @@ spliqq];

ChannelMap = {
  "eexbbxag"   -> "eexddxag MB=0 YB=0",
  "nenexbbxag" -> "nenexddxag MB=0= YB=0"
};
