Only keep vertices involving a particle and its antiparticle

FlexibleSUSY · Jan 25, 2016 · a0d3542 · a0d3542
1 parent 432b08f
commit a0d3542
Showing 1 changed file with 2 additions and 1 deletion.
diff --git a/meta/EffectiveCouplings.m b/meta/EffectiveCouplings.m
@@ -168,8 +168,9 @@ previous results (e.g. define along the lines of f[p] := f[p] = ...)    *)
                 neededCoups = Select[neededTwoBodyDecays[[2]],
                                     (MemberQ[neededVectorBosonInteractions[[2]], #[[1]]] ||
                                      MemberQ[neededVectorBosonInteractions[[2]], #[[2]]])&];
+                (* only keep vertices of the form pD -> p AntiParticle[p] *)
+                neededCoups = Cases[neededCoups, {p1_, p2_, _} /; p1 === AntiParticle[p2]];
                 (* filter out massless states and Goldstones *)
-                Print["neededCoupgs = ", neededCoups];
                 neededCoups = Select[neededCoups, (!IsMasslessOrGoldstone[#[[1]]] && !IsMasslessOrGoldstone[#[[2]]])&];
                 result = Append[result, {couplings[[i]], #[[3]]& /@ neededCoups}];
               ];