Added check if edmParticle is a fermion. Otherwise we don't know how …

…to calculate the diagrams.
FlexibleSUSY · Jan 26, 2017 · 1f25afc · 1f25afc
1 parent 9b0f34e
commit 1f25afc
Showing 1 changed file with 3 additions and 5 deletions.
diff --git a/meta/EDM.m b/meta/EDM.m
@@ -221,6 +221,9 @@ If you add new kinds of vertices (e.g for new diagram types):
                             {SARAH`FieldToInsert[1], SARAH`FieldToInsert[2]}}],
             subtypedDiagrams, uniqueDiagrams},
 
+           If[TreeMasses`IsFermion[edmField] =!= True,
+              Return[{edmField,{}}]];
+
            subtypedDiagrams = (Module[{photonEmitter = #[[2,1]],
                                        exchangeField = #[[2,2]],
                                        subType},
@@ -466,11 +469,6 @@ If you add new kinds of vertices (e.g for new diagram types):
            Return[cachedContributingDiagrams];
           ];
 
-(* Returns a list of all concrete diagram evaluators
- format: {{edmField1, {"DiagramEvaluator<OneLoopDiagram<1>, Fe, VP>", "...", ... }},
-          {edmField2, {"...", ... }},
-          ...}
- that need to be invoked in our calculation *)
 ConcreteDiagramEvaluators[] :=
      ({#[[1]],
          (("DiagramEvaluator<" <> SymbolName @ Head @ #[[1]] <> "<" <>