FeynArts 3.3

Author: t-hahn

FeynArts.m

@@ -2,7 +2,7 @@
 	Analytic.m
 		Translation of InsertFields output into
 		analytic expressions
-		last modified 22 Mar 06 th
+		last modified 10 Jan 08 th
 *)
 
 Begin["`Analytic`"]
@@ -507,7 +507,7 @@ need another (-1)^(Length[ext]/2).
   ins = ReplacePart[gm, sgen/s, -1] /. {ru} /.
     anti -> AntiParticle /.
     app[ x_. (fi:P$Generic)[n__], k__ ] :> x fi[n, k];
-  deltas = DeleteCases[ Union@@ Diagonal/@
+  deltas = DeleteCases[ Union@@ CouplingDeltas/@
     Union[ Cases[ins, G[_][cto_][fi__][__] :> FieldPoint[cto][fi]] ],
     _Integer ];
   ins = ins /. G -> GtoC /. Mass -> TheMass /. gaugeru /.
@@ -708,6 +708,7 @@ kinematical expression (for a G[-]).  If neither method resolves TheC,
       Range@@ Sort[Floor[{b, a}]] ]] ]
 
 
+(*
 DiagramSelect[ tops:TopologyList[info__][__], crit_ ] :=
 Block[ {lev, Rule},
   lev = ResolveLevel[InsertionLevel /. {info}][[-1]];
@@ -717,10 +718,33 @@ kinematical expression (for a G[-]).  If neither method resolves TheC,
     (Graph[__][__] -> _[]) :> Seq[] /.
     (Topology[__][__] -> _[]) :> Seq[]
 ]
+*)
+
+DiagramSelect[ tops:TopologyList[info__][__], crit_ ] :=
+Block[ {lev, Rule},
+  lev = ResolveLevel[InsertionLevel /. {info}][[-1]];
+  Rule[_] := Sequence[];
+  Apply[ #1 -> (#2 /. Graph[__, lev == _][fi__] :> Seq[] /;
+      crit[{fi}, #1] =!= True)&, tops, 1 ] /.
+    (Graph[__][__] -> _[]) :> Seq[] /.
+    (Topology[__][__] -> _[]) :> Seq[]
+]
 
 DiagramSelect[ amp_, crit_ ] := Select[amp, crit]
 
 
+Attributes[prop] = {Orderless}
+
+Attributes[merge] = {Flat, Orderless}
+
+merge[ prop[i_, j_], prop[j_, k_] ] := prop[i, k]
+
+FermionRouting[ fields_:{}, top:P$Topology ] := Level[
+  merge@@ Apply[ prop[ #1[[1]], #2[[1]] ]&,
+    Select[AddFieldNo[top] /. fields, !FreeQ[#, P$NonCommuting]&], 1 ],
+  {-1} ]
+
+
 DiagramComplement[ tops:TopologyList[info__][___],
   more:TopologyList[__][___].. ] :=
 Block[ {lev = ResolveLevel[InsertionLevel /. {info}][[-1]]},
@@ -769,7 +793,7 @@ kinematical expression (for a G[-]).  If neither method resolves TheC,
   Global`PolarizationVector[ _, mom_, li_ ] =
     Global`PolarizationVector[mom, li];
 
-  Global`DiracSpinor[ mom_, mass_, ___ ] := FeynArts`Spinor[mom, mass];
+  (* Global`DiracSpinor[ mom_, mass_, ___ ] := FeynArts`Spinor[mom, mass]; *)
 
   Index[ Global`Lorentz, n_ ] := Index[Global`Lorentz, n] =
     ToExpression["li" <> ToString[n]];

FeynArts/Analytic.m

@@ -1,7 +1,7 @@
 (*
 	Graphics.m
 		Graphics routines for FeynArts
-		last modified 10 Jan 05 th
+		last modified 4 Sep 07 th
 *)
 
 Begin["`Graphics`"]
@@ -73,7 +73,9 @@
   SheetHeader -> Automatic,
   Numbering -> Full,
   FieldNumbers -> False,
-  DisplayFunction :> $DisplayFunction
+  If[ $VersionNumber < 6,
+    DisplayFunction :> $DisplayFunction,
+    DisplayFunction :> (Print/@ Render[#]&) ]
 }
 
 Paint::nolevel =
@@ -227,7 +229,8 @@
   opt___?OptionQ ] :=
 Block[ {rg},
   rg = Render[g, InferFormat[chan, format],
-    ImageSize -> (ImageSize /. {opt} /. Options[Display])];
+    ImageSize -> (ImageSize /. {opt} /. Options[Display] /.
+                   ImageSize -> Automatic)];
   MapThread[Display[##, format, opt]&,
     {FilePerSheet[chan, Length[rg]], rg}]
 ]
@@ -238,7 +241,8 @@
   opt___?OptionQ ] :=
 Block[ {rg},
   rg = Render[g, InferFormat[chan, format],
-    ImageSize -> (ImageSize /. {opt} /. Options[Export])];
+    ImageSize -> (ImageSize /. {opt} /. Options[Export] /.
+                   ImageSize -> Automatic)];
   MapThread[Export[##, format, opt]&,
     {FilePerSheet[chan, Length[rg]], rg}]
 ]
@@ -388,11 +392,16 @@
     AspectRatio -> rows/cols]
 ]
 
+
 DiagramBox[ Null, _ ] = {}
 
-DiagramBox[ g_, {yoff_, xoff_} ] :=
-  Rectangle[ {xoff - 1, rows - yoff} DiagramSize,
-             {xoff, rows - yoff + 1} DiagramSize, MmaRender[g] ]
+DiagramBox[ g_, {yoff_, xoff_} ] := Inset[ MmaRender[g],
+  {xoff - 1, rows - yoff} DiagramSize, {0, 0}, {1, 1} DiagramSize ]
+
+If[ $VersionNumber < 6,
+  Inset[obj_, pos_, _, size_] := Rectangle[pos, pos + size, obj]
+]
+
 
 Title[ ] = {}