FeynArts 3.9; Last update: 4 May 15;

Added `Charge' quantum number to S[5] (charged Higgs) in MSSM.mod (MSSMCT.mod not affected).

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 2 Dec 14 by Thomas Hahn
+last modified 4 Apr 15 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 2 Dec 14"]
+Print["last revised 4 Apr 15"]
 
 
 BeginPackage["FeynArts`"]
@@ -544,16 +544,6 @@ 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.
@@ -1064,12 +1054,8 @@ as arguments the quantum numbers of the involved fields (times -1 for
 to add graph information to the amplitude."
 
 VertexDebug::usage =
-"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."
+"VertexDebug[debuginfo] is a function invoked whenever a vertex cannot
+be resolved.  It is used for debugging FeynArts."
 
 VertexFunction::usage =
 "VertexFunction[o][f1, f2, ...] represents the 1PI vertex function of
@@ -1362,7 +1348,9 @@ topologies where it specifies whether the field numbers (the n in
 "P$Generic is the pattern for generic fields."
 
 P$NonCommuting::usage =
-"P$NonCommuting is the pattern for the non-commuting generic fields."
+"P$NonCommuting is the pattern for the non-commuting generic fields.
+A list {patt1, patt2, ...} (e.g. {F, U}) means that chains are built
+first particles matching patt1, then for patt2, etc."
 
 P$InsertionObjects::usage =
 "P$InsertionObjects matches the objects in the generic amplitude that
@@ -1373,7 +1361,7 @@ topologies where it specifies whether the field numbers (the n in
 
 P$Generic = F | S | V | T | U | _Mix | _Rev
 
-P$NonCommuting = F | U
+P$NonCommuting = {F, U}
 
 P$InsertionObjects = G[_][_][__][__] | _Mass | _GaugeXi |
   VertexFunction[_][__]

FeynArts/Analytic.m

@@ -2,7 +2,7 @@
 	Analytic.m
 		Translation of InsertFields output into
 		analytic expressions
-		last modified 2 Dec 14 th
+		last modified 4 Apr 15 th
 *)
 
 Begin["`Analytic`"]
@@ -56,7 +56,8 @@
 
 CreateFeynAmp[ tops:TopologyList[info___][___], options___Rule ] :=
 Block[ {alevel, pref, next, gaugeru, truncru, momcons, graphinfo, toplist,
-amps, head, topnr = 1, opt = ActualOptions[CreateFeynAmp, options]},
+amps, head, topnr = 1, p$nc = Level[P$NonCommuting, {-1}, Alternatives],
+opt = ActualOptions[CreateFeynAmp, options]},
 
   If[ (alevel = ResolveLevel[AmplitudeLevel /. opt /. {info} /.
         Options[InsertFields]]) === $Failed,
@@ -329,20 +330,31 @@ explicit spinor index with which it is possible (outside of FeynArts)
 	   e.g. when computing counter terms from the self energies *)
 Fixgmc[ c1:_[__, _?SelfConjugate], r___, c2:_[__, _?SelfConjugate] ] :=
   dot[c1, r, c2] /;
-  OrderedQ[{c2[[0, 1]][ c2[[2, 1]] ], c1[[0, 1]][ c1[[1, 1]] ]}]
+  OrderedQ[{c2[[0, 1]] @ c2[[2, 1]], c1[[0, 1]] @ c1[[1, 1]]}]
 
 Fixgmc[ c__ ] := Reverse[ ReverseProp/@ dot[c] ]
 
 
-MakeFermionChains[ top_ ] := top /; FreeQ[top, P$NonCommuting]
+Chkgmc[ c_ ] := Chkgmc[Level[c, {2}][[{1, -2}]], c]
+
+Chkgmc[ {v_, v_}, c_tr ] := c
+
+Chkgmc[ {Vertex[1][_], Vertex[1][_]}, c_dot ] := c
+
+Chkgmc[ _, c_ ] := List@@ c
+
+
+NCSelect[ top_, nc_ ] := (
+  ch = {ch, BuildChain@@ gmc/@ Select[top, !FreeQ[#[[3]], nc]&]};
+  Select[top, FreeQ[#[[3]], nc]&] )
+
+
+MakeFermionChains[ top_ ] := top /; FreeQ[top, p$nc]
 
 MakeFermionChains[ top_ ] :=
-Block[ {res, ext},
-  res = Append[
-    Select[top, FreeQ[#[[3]], P$NonCommuting]&],
-    BuildChain@@ gmc/@ Select[top, !FreeQ[#[[3]], P$NonCommuting]&] /.
-      gmc -> Fixgmc
-  ] /. BuildChain -> Sequence;
+Block[ {ch = {}, res, ext},
+  res = Fold[NCSelect, top, Flatten[{P$NonCommuting}]];
+  res = Flatten[{res, Cases[ch, gmc[c__] :> Chkgmc[Fixgmc[c]], Infinity]}];
 
 	(* Since fermion chains are always traversed opposite to the
 	   fermion flow, we need the sign of the permutation that gets
@@ -403,7 +415,7 @@ need another (-1)^(Length[ext]/2).
   If[ perm === $Failed, Return[{}] ];
   v = v[[perm]];
   If[ cto < 0,
-    I PV[ If[FreeQ[v, P$NonCommuting], Identity, NonCommutative][
+    I PV[ If[FreeQ[v, p$nc], Identity, NonCommutative][
             VertexFunction[-cto]@@ v ] ],
   (* else *)
     AnalyticalCoupling[cto]@@ v ]
@@ -802,8 +814,7 @@ kinematical expression (for a G[-]).  If neither method resolves TheC,
 
 FermionRouting[ gr_:{}, top:P$Topology, ___ ] := Level[
   merge@@ Apply[ prop[ #1[[1]], #2[[1]] ]&,
-    Select[AddFieldNo[top] /. List@@ gr, !FreeQ[#, P$NonCommuting]&],
-    1 ],
+    Select[AddFieldNo[top] /. List@@ gr, !FreeQ[#, p$nc]&], 1 ],
   {-1} ]