FeynArts 3.9; Last update: 2 Dec 14;

Added a few functions and a manual section (7.4) for model-file debugging.

Author: t-hahn

FeynArts.m

@@ -2,7 +2,7 @@
 
 This is FeynArts, Version 3.9
 Copyright by Sepp Kueblbeck, Hagen Eck, and Thomas Hahn 1991-2014
-last modified 24 Oct 14 by Thomas Hahn
+last modified 2 Dec 14 by Thomas Hahn
 
 Release notes:
 
@@ -46,7 +46,7 @@
 Print[""];
 Print["FeynArts 3.9"];
 Print["by Hagen Eck, Sepp Kueblbeck, and Thomas Hahn"];
-Print["last revised 24 Oct 14"]
+Print["last revised 2 Dec 14"]
 
 
 BeginPackage["FeynArts`"]
@@ -544,6 +544,16 @@ symbols which must (or must not, if preceded by \"!\") appear in the
 "Reinitialize is an option of InitializeModel.  InitializeModel will
 reinitialize the current model only if Reinitialize is set to True."
 
+TagCouplings::usage =
+"TagCouplings is an option of InitializeModel.  If True, all couplings
+will be tagged by Coupling[i], where i is the coupling's position in
+M$CouplingMatrices, such that parts of the amplitude can be traced back
+to specific couplings later."
+
+Coupling::usage =
+"Coupling[i...] is the identifier used to tag the coupling at
+M$CouplingMatrices[[i]]."
+
 GenericModelEdit::usage =
 "GenericModelEdit is an option of InitializeModel.  It specifies code
 that will be executed directly after loading the generic model, i.e.
@@ -1054,8 +1064,12 @@ as arguments the quantum numbers of the involved fields (times -1 for
 to add graph information to the amplitude."
 
 VertexDebug::usage =
-"VertexDebug[debuginfo] is a function invoked whenever a vertex cannot
-be resolved.  It is used for debugging FeynArts."
+"VertexDebug[info] is a function invoked whenever a vertex cannot be
+resolved.  It is used for debugging FeynArts."
+
+VertexMonitor::usage =
+"VertexMonitor[info] is a function invoked whenever a vertex is resolved. 
+It is used for debugging FeynArts."
 
 VertexFunction::usage =
 "VertexFunction[o][f1, f2, ...] represents the 1PI vertex function of

FeynArts/Analytic.m

@@ -2,7 +2,7 @@
 	Analytic.m
 		Translation of InsertFields output into
 		analytic expressions
-		last modified 10 Apr 14 th
+		last modified 2 Dec 14 th
 *)
 
 Begin["`Analytic`"]
@@ -570,6 +570,9 @@ kinematical expression (for a G[-]).  If neither method resolves TheC,
 
   cv = SignResolve[sym, cvr = TheC[kin]@@ vert];
 
+  VertexMonitor[{"vert" -> vert, "kin" -> kin, "cv" -> cv,
+    "fi" -> {fi}, "cto" -> cto}];
+
   If[ !FreeQ[cv, TheC],
     Message[CreateFeynAmp::nocoupl, vert, kin];
     VertexDebug[{"vert" -> vert, "kin" -> kin, "cv" -> cvr,

FeynArts/Initialize.m

@@ -1,7 +1,7 @@
 (*
 	Initialize.m
 		Functions for the initialization of models
-		last modified 19 Jul 13 th
+		last modified 2 Dec 14 th
 *)
 
 Begin["`Initialize`"]
@@ -208,6 +208,7 @@
 Options[ InitializeModel ] = {
   GenericModel -> "Lorentz",
   Reinitialize -> True,
+  TagCouplings -> False,
   GenericModelEdit :> Null,
   ModelEdit :> Null
 }
@@ -260,9 +261,10 @@ assigned values somewhere in your model file (most often \"i\" or \"j\").  \
 InitializeModel[ opt:P$Options ] := InitializeModel[{}, opt]
 
 InitializeModel[ model_, options:P$Options ] :=
-Block[ {reini, genmod, mod,
+Block[ {reini, tagcoup, genmod, mod,
 opt = ActualOptions[InitializeModel, options]},
   reini = TrueQ[Reinitialize /. opt];
+  tagcoup = TrueQ[TagCouplings /. opt];
 
   genmod = ToModelName[GenericModel /. opt];
   If[ reini || genmod =!= $GenericModel,
@@ -525,6 +527,15 @@ assigned values somewhere in your model file (most often \"i\" or \"j\").  \
 Assign[ fi__, a_ -> b_ ] := a[fi] = b
 
 
+TagCoupling[ x_, {_, 1, __} ] := x
+
+TagCoupling[ x_, {n__} ] := x Coupling[n]
+
+
+TagCoupling[ coup_ == expr_, {n_} ] := coup == 
+  MapIndexed[#1 Coupling[n, #2[[1]] - 1]&, expr, {2}]
+
+
 (* Initialization of a classes model: *)
 
 Attributes[ InitModel ] = {HoldAll}
@@ -584,7 +595,10 @@ assigned values somewhere in your model file (most often \"i\" or \"j\").  \
     InitCoupling@@@ Block[ {Coup},
       _Coup = 0;
       TheCoeff[ fi_ ] := Throw[Message[InitializeModel::unknown, fi]];
-      Scan[SetCoupling, M$CouplingMatrices];
+      Scan[SetCoupling, If[ tagcoup, 
+        MapIndexed[TagCoupling, M$CouplingMatrices, {4}],
+      (* else *)
+        M$CouplingMatrices ]];
       TheCoeff[ fi_ ] =.;
       _Coup =.;
       DownValues[Coup]

FeynArts39.m

@@ -2,7 +2,7 @@
 
 This is FeynArts, Version 3.9
 Copyright by Sepp Kueblbeck, Hagen Eck, and Thomas Hahn 1991-2014
-last modified 24 Oct 14 by Thomas Hahn
+last modified 2 Dec 14 by Thomas Hahn
 
 Release notes:
 
@@ -46,7 +46,7 @@
 Print[""];
 Print["FeynArts 3.9"];
 Print["by Hagen Eck, Sepp Kueblbeck, and Thomas Hahn"];
-Print["last revised 24 Oct 14"]
+Print["last revised 2 Dec 14"]
 
 
 BeginPackage["FeynArts`"]
@@ -544,6 +544,16 @@ symbols which must (or must not, if preceded by \"!\") appear in the
 "Reinitialize is an option of InitializeModel.  InitializeModel will
 reinitialize the current model only if Reinitialize is set to True."
 
+TagCouplings::usage =
+"TagCouplings is an option of InitializeModel.  If True, all couplings
+will be tagged by Coupling[i], where i is the coupling's position in
+M$CouplingMatrices, such that parts of the amplitude can be traced back
+to specific couplings later."
+
+Coupling::usage =
+"Coupling[i...] is the identifier used to tag the coupling at
+M$CouplingMatrices[[i]]."
+
 GenericModelEdit::usage =
 "GenericModelEdit is an option of InitializeModel.  It specifies code
 that will be executed directly after loading the generic model, i.e.
@@ -1054,8 +1064,12 @@ as arguments the quantum numbers of the involved fields (times -1 for
 to add graph information to the amplitude."
 
 VertexDebug::usage =
-"VertexDebug[debuginfo] is a function invoked whenever a vertex cannot
-be resolved.  It is used for debugging FeynArts."
+"VertexDebug[info] is a function invoked whenever a vertex cannot be
+resolved.  It is used for debugging FeynArts."
+
+VertexMonitor::usage =
+"VertexMonitor[info] is a function invoked whenever a vertex is resolved. 
+It is used for debugging FeynArts."
 
 VertexFunction::usage =
 "VertexFunction[o][f1, f2, ...] represents the 1PI vertex function of

Models/ModelMaker.m

@@ -2,7 +2,7 @@
 	ModelMaker.m
 		derives the Feynman rules from the Lagrangian
 		and splices them into a template model file
-		last modified 28 Sep 07 th
+		last modified 2 Dec 14 th
 *)
 
 
@@ -87,7 +87,7 @@
 FuncD[expr_, f_List] := ReduceDeltas[Fold[FuncD, expr, f]]
 
 FuncD[expr_, f_Field] :=
-Block[ {NCSign, s, times, texpr},
+Block[ {NCSign, s, times, texpr},
   If[ !FreeQ[f, P$NonCommuting],
     s = -1;
     NCSign[P$NonCommuting] := s = -s ];