(* Content-type: application/vnd.wolfram.mathematica *)

(*** Wolfram Notebook File ***)
(* https://www.wolfram.com/nb *)

(* CreatedBy='Mathematica 8.0' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       157,          7]
NotebookDataLength[     41811,       1273]
NotebookOptionsPosition[     40238,       1215]
NotebookOutlinePosition[     40582,       1230]
CellTagsIndexPosition[     40539,       1227]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{
Cell[BoxData[
 RowBox[{"<<", "HighEnergyPhysics`FeynCalc`"}]], "Input"],

Cell["\<\
Loading FeynCalc from C:\\Archivos de programa\\Wolfram Research\\Mathematica\
\\8.0\\AddOns\\Applications\\HighEnergyPhysics\
\>", "Text",
 CellMargins->{{Inherited, Inherited}, {1, 0}},
 GeneratedCell->True,
 CellAutoOverwrite->True,
 CellChangeTimes->{3.537705500515625*^9}],

Cell[TextData[{
 StyleBox["FeynCalc",
  FontWeight->"Bold"],
 " ",
 "8.1.0",
 " For help, type ?FeynCalc, open ",
 ButtonBox["FeynCalcRef8.nb",
  BaseStyle->"AddOnsLink",
  ButtonFunction:>NotebookOpen[
   "C:\\Archivos de programa\\Wolfram \
Research\\Mathematica\\8.0\\AddOns\\Applications\\HighEnergyPhysics\\\
Documentation\\English\\FeynCalcRef8.nb"],
  ButtonData:>{"Short Overview", "intro"},
  ButtonNote->"Open FeynCalcRef8.nb"],
 " or visit ",
 ButtonBox["www.feyncalc.org",
  BaseStyle->"Hyperlink",
  ButtonData:>{
    URL["https://www.feyncalc.org/"], None},
  ButtonNote->"https://www.feyncalc.org/"]
}], "Text",
 CellMargins->{{Inherited, Inherited}, {1, 0}},
 GeneratedCell->True,
 CellAutoOverwrite->True,
 CellChangeTimes->{3.537705506671875*^9}],

Cell[TextData[{
 "Loading FeynArts, see ",
 ButtonBox["www.feynarts.de",
  BaseStyle->"Hyperlink",
  ButtonData:>{
    URL["https://www.feynarts.de/"], None},
  ButtonNote->"https://www.feynarts.de/"],
 " for documentation"
}], "Text",
 CellMargins->{{Inherited, Inherited}, {1, 0}},
 GeneratedCell->True,
 CellAutoOverwrite->True,
 CellChangeTimes->{3.5377055170625*^9}],

Cell["FeynArts 3.4 patched for use with FeynCalc", "Text",
 CellMargins->{{Inherited, Inherited}, {1, 0}},
 GeneratedCell->True,
 CellAutoOverwrite->True,
 CellChangeTimes->{3.53770551746875*^9}],

Cell[BoxData[
 RowBox[{
  RowBox[{"SetMandelstam", "[",
   RowBox[{"s", ",", "t", ",", "u", ",", "p1", ",", "p2", ",",
    RowBox[{"-", "p3"}], ",",
    RowBox[{"-", "p4"}], ",", "m", ",", "m", ",", "0", ",", "0"}], "]"}],
  ";"}]], "Input",
 CellChangeTimes->{{3.537289952578125*^9, 3.53728999478125*^9}}],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"M1", "=",
   RowBox[{
    RowBox[{"g", "^", "2"}], "*",
    RowBox[{
     RowBox[{"Spinor", "[",
      RowBox[{
       RowBox[{"-", "p2"}], ",", "m"}], "]"}], ".",
     RowBox[{"GA", "[", "\[Nu]", "]"}], ".",
     RowBox[{"(",
      RowBox[{
       RowBox[{"GS", "[",
        RowBox[{"p1", "-", "p3"}], "]"}], "+", "m"}], ")"}], ".",
     RowBox[{"GA", "[", "\[Mu]", "]"}], ".",
     RowBox[{"Spinor", "[",
      RowBox[{"p1", ",", "m"}], "]"}]}], "*",
    RowBox[{"Conjugate", "[",
     RowBox[{"PolarizationVector", "[",
      RowBox[{"p3", ",", "\[Mu]"}], "]"}], "]"}], "*",
    RowBox[{
     RowBox[{"Conjugate", "[",
      RowBox[{"PolarizationVector", "[",
       RowBox[{"p4", ",", "\[Nu]"}], "]"}], "]"}], "/",
     RowBox[{"(",
      RowBox[{
       RowBox[{"SP", "[",
        RowBox[{
         RowBox[{"p1", "-", "p3"}], ",",
         RowBox[{"p1", "-", "p3"}]}], "]"}], "-",
       RowBox[{"m", "^", "2"}]}], ")"}]}]}]}], ";"}], "\n",
 RowBox[{
  RowBox[{"M1A", "=",
   RowBox[{
    RowBox[{
     RowBox[{"ComplexConjugate", "[", "M1", "]"}], " ", "/.", " ",
     RowBox[{"\[Mu]", "\[Rule]", " ", "\[Mu]2"}]}], " ", "/.", "  ",
    RowBox[{"\[Nu]", "\[Rule]", " ", "\[Nu]2"}]}]}], ";"}], "\n",
 RowBox[{
  RowBox[{"F1", "=",
   RowBox[{
    RowBox[{
     RowBox[{"1", "/", "9"}], "*",
     RowBox[{"SUNTrace", "[",
      RowBox[{"SUNT", "[",
       RowBox[{"a3", ",", "a4", ",", "a4", ",", "a3"}], "]"}], "]"}]}], " ", "/.",
     " ",
    RowBox[{"SUNN", "\[Rule]", " ", "3"}]}]}], ";"}], "\n",
 RowBox[{
  RowBox[{"T11", "=",
   RowBox[{
    RowBox[{"(",
     RowBox[{
      RowBox[{
       RowBox[{
        RowBox[{"1", "/",
         RowBox[{"(",
          RowBox[{"16", "Pi", " ",
           RowBox[{"s", "^", "2"}]}], ")"}]}], "*",
        RowBox[{"F1", "/", "4"}], "*",
        RowBox[{"(",
         RowBox[{
          RowBox[{
           RowBox[{
            RowBox[{"M1A", "  ", "M1"}], " ", "//", " ",
            "DoPolarizationSums"}], " ", "//", " ", "Contract"}], " ", "//",
          " ", "FermionSpinSum"}], ")"}]}], " ", "/.", " ",
       RowBox[{"DiracTrace", "\[Rule]", "TR"}]}], " ", "//", " ",
      "Simplify"}], " ", ")"}], "/.", " ",
    RowBox[{"g", "\[Rule]",
     RowBox[{"Sqrt", "[",
      RowBox[{"4", " ", "Pi", " ", "\[Alpha]"}], "]"}]}]}]}], ";"}], "\n",
 RowBox[{
  RowBox[{"T11exp", "=",
   RowBox[{
    RowBox[{
     RowBox[{"1", "/",
      RowBox[{"(",
       RowBox[{"16", "Pi", " ",
        RowBox[{"s", "^", "2"}]}], ")"}]}], "*",
     RowBox[{"g", "^", "4"}], "*",
     RowBox[{"(",
      RowBox[{
       RowBox[{"2", "^", "5"}], "/",
       RowBox[{"3", "^", "3"}]}], ")"}], "*",
     RowBox[{
      RowBox[{"(",
       RowBox[{
        RowBox[{"t", "*", "u"}], "+",
        RowBox[{
         RowBox[{"m", "^", "2"}], "*",
         RowBox[{"(",
          RowBox[{"s", "-",
           RowBox[{"2", "t"}], "-",
           RowBox[{"3",
            RowBox[{"m", "^", "2"}]}]}], ")"}]}]}], ")"}], "/",
      RowBox[{
       RowBox[{"(",
        RowBox[{"t", "-",
         RowBox[{"m", "^", "2"}]}], ")"}], "^", "2"}]}]}], "  ", "/.", " ",
    RowBox[{"g", "\[Rule]",
     RowBox[{"Sqrt", "[",
      RowBox[{"4", " ", "Pi", " ", "\[Alpha]"}], "]"}]}]}]}], ";"}], "\n",
 RowBox[{"TrickMandelstam", "[",
  RowBox[{
   RowBox[{"T11", "-", "T11exp"}], ",",
   RowBox[{"{",
    RowBox[{"s", ",", "t", ",", "u", ",",
     RowBox[{"2",
      RowBox[{"m", "^", "2"}]}]}], "}"}]}], "]"}]}], "Input",
 CellChangeTimes->{{3.53729006521875*^9, 3.53729007975*^9}, {
   3.537291055203125*^9, 3.5372910589375*^9}, {3.537294459890625*^9,
   3.53729446984375*^9}, 3.53729468453125*^9, {3.53729659815625*^9,
   3.537296623359375*^9}, {3.53729667146875*^9, 3.537296758109375*^9}}],

Cell[BoxData[
 FormBox["0", TraditionalForm]], "Output",
 CellChangeTimes->{{3.53729669040625*^9, 3.537296759265625*^9},
   3.537302445125*^9, 3.537532348421875*^9, 3.53761063946875*^9,
   3.53761105490625*^9, 3.53761147303125*^9, 3.537705535625*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"M2", "=",
   RowBox[{
    RowBox[{"g", "^", "2"}], "*",
    RowBox[{
     RowBox[{"Spinor", "[",
      RowBox[{
       RowBox[{"-", "p2"}], ",", "m"}], "]"}], ".",
     RowBox[{"GA", "[", "\[Nu]", "]"}], ".",
     RowBox[{"(",
      RowBox[{
       RowBox[{"GS", "[",
        RowBox[{"p1", "-", "p4"}], "]"}], "+", "m"}], ")"}], ".",
     RowBox[{"GA", "[", "\[Mu]", "]"}], ".",
     RowBox[{"Spinor", "[",
      RowBox[{"p1", ",", "m"}], "]"}]}], "*",
    RowBox[{"Conjugate", "[",
     RowBox[{"PolarizationVector", "[",
      RowBox[{"p4", ",", "\[Mu]"}], "]"}], "]"}], "*",
    RowBox[{
     RowBox[{"Conjugate", "[",
      RowBox[{"PolarizationVector", "[",
       RowBox[{"p3", ",", "\[Nu]"}], "]"}], "]"}], "/",
     RowBox[{"(",
      RowBox[{
       RowBox[{"SP", "[",
        RowBox[{
         RowBox[{"p1", "-", "p4"}], ",",
         RowBox[{"p1", "-", "p4"}]}], "]"}], "-",
       RowBox[{"m", "^", "2"}]}], ")"}]}]}]}], ";"}], "\n",
 RowBox[{
  RowBox[{"M2A", "=",
   RowBox[{
    RowBox[{
     RowBox[{"ComplexConjugate", "[", "M2", "]"}], " ", "/.", " ",
     RowBox[{"\[Mu]", "\[Rule]", " ", "\[Mu]2"}]}], " ", "/.", "  ",
    RowBox[{"\[Nu]", "\[Rule]", " ", "\[Nu]2"}]}]}], ";"}], "\n",
 RowBox[{
  RowBox[{"F2", "=",
   RowBox[{
    RowBox[{
     RowBox[{"1", "/", "9"}], "*",
     RowBox[{"SUNTrace", "[",
      RowBox[{"SUNT", "[",
       RowBox[{"a4", ",", "a3", ",", "a3", ",", "a4"}], "]"}], "]"}]}], " ", "/.",
     " ",
    RowBox[{"SUNN", "\[Rule]", " ", "3"}]}]}], ";"}], "\n",
 RowBox[{
  RowBox[{"T22", "=",
   RowBox[{
    RowBox[{"(",
     RowBox[{
      RowBox[{
       RowBox[{
        RowBox[{"1", "/",
         RowBox[{"(",
          RowBox[{"16", "Pi", " ",
           RowBox[{"s", "^", "2"}]}], ")"}]}], "*",
        RowBox[{"F2", "/", "4"}], "*",
        RowBox[{"(",
         RowBox[{
          RowBox[{
           RowBox[{
            RowBox[{"M2A", "  ", "M2"}], " ", "//", " ",
            "DoPolarizationSums"}], " ", "//", " ", "Contract"}], " ", "//",
          " ", "FermionSpinSum"}], ")"}]}], " ", "/.", " ",
       RowBox[{"DiracTrace", "\[Rule]", "TR"}]}], " ", "//", " ",
      "Simplify"}], " ", ")"}], "/.", " ",
    RowBox[{"g", "\[Rule]",
     RowBox[{"Sqrt", "[",
      RowBox[{"4", " ", "Pi", " ", "\[Alpha]"}], "]"}]}]}]}], ";"}], "\n",
 RowBox[{
  RowBox[{"T22exp", "=",
   RowBox[{
    RowBox[{
     RowBox[{"1", "/",
      RowBox[{"(",
       RowBox[{"16", "Pi", " ",
        RowBox[{"s", "^", "2"}]}], ")"}]}], "*",
     RowBox[{"g", "^", "4"}], "*",
     RowBox[{"(",
      RowBox[{
       RowBox[{"2", "^", "5"}], "/",
       RowBox[{"3", "^", "3"}]}], ")"}], "*",
     RowBox[{
      RowBox[{"(",
       RowBox[{
        RowBox[{"t", "*", "u"}], "+",
        RowBox[{
         RowBox[{"m", "^", "2"}], "*",
         RowBox[{"(",
          RowBox[{"s", "-",
           RowBox[{"2", "u"}], "-",
           RowBox[{"3",
            RowBox[{"m", "^", "2"}]}]}], ")"}]}]}], ")"}], "/",
      RowBox[{
       RowBox[{"(",
        RowBox[{"u", "-",
         RowBox[{"m", "^", "2"}]}], ")"}], "^", "2"}]}]}], "  ", "/.", " ",
    RowBox[{"g", "\[Rule]",
     RowBox[{"Sqrt", "[",
      RowBox[{"4", " ", "Pi", " ", "\[Alpha]"}], "]"}]}]}]}], ";"}], "\n",
 RowBox[{"TrickMandelstam", "[",
  RowBox[{
   RowBox[{"T22", "-", "T22exp"}], ",",
   RowBox[{"{",
    RowBox[{"s", ",", "t", ",", "u", ",",
     RowBox[{"2",
      RowBox[{"m", "^", "2"}]}]}], "}"}]}], "]"}]}], "Input",
 CellChangeTimes->{{3.5372953049375*^9, 3.5372953164375*^9}, {
   3.537296693953125*^9, 3.53729669671875*^9}, 3.537296769984375*^9, {
   3.53729681096875*^9, 3.5372968891875*^9}}],

Cell[BoxData[
 FormBox["0", TraditionalForm]], "Output",
 CellChangeTimes->{{3.5372968764375*^9, 3.53729689046875*^9},
   3.5373000729375*^9, 3.537302447296875*^9, 3.53753235028125*^9,
   3.537610642109375*^9, 3.53761105753125*^9, 3.537611474375*^9,
   3.537705536671875*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"F12", "=",
   RowBox[{
    RowBox[{
     RowBox[{"1", "/", "9"}], "*",
     RowBox[{"SUNTrace", "[",
      RowBox[{"SUNT", "[",
       RowBox[{"a3", ",", "a4", ",", "a3", ",", "a4"}], "]"}], "]"}]}], " ", "/.",
     " ",
    RowBox[{"SUNN", "\[Rule]", " ", "3"}]}]}], ";"}], "\n",
 RowBox[{
  RowBox[{"T12", "=",
   RowBox[{
    RowBox[{"(",
     RowBox[{
      RowBox[{
       RowBox[{
        RowBox[{"1", "/",
         RowBox[{"(",
          RowBox[{"16", "Pi", " ",
           RowBox[{"s", "^", "2"}]}], ")"}]}], "*",
        RowBox[{"F12", "/", "4"}], "*",
        RowBox[{"(",
         RowBox[{
          RowBox[{
           RowBox[{
            RowBox[{"M1A", "  ", "M2"}], " ", "//", " ",
            "DoPolarizationSums"}], " ", "//", " ", "Contract"}], " ", "//",
          " ", "FermionSpinSum"}], ")"}]}], " ", "/.", " ",
       RowBox[{"DiracTrace", "\[Rule]", "TR"}]}], " ", "//", " ",
      "Simplify"}], " ", ")"}], "/.", " ",
    RowBox[{"g", "\[Rule]",
     RowBox[{"Sqrt", "[",
      RowBox[{"4", " ", "Pi", " ", "\[Alpha]"}], "]"}]}]}]}], ";"}], "\n",
 RowBox[{
  RowBox[{"T12exp", "=",
   RowBox[{
    RowBox[{
     RowBox[{
      RowBox[{
       RowBox[{"-", "4"}], "/", "27"}], "/",
      RowBox[{"(",
       RowBox[{"16", "Pi", " ",
        RowBox[{"s", "^", "2"}]}], ")"}]}], "*",
     RowBox[{"g", "^", "4"}], "*",
     RowBox[{"m", "^", "2"}], "*",
     RowBox[{
      RowBox[{
       RowBox[{"(",
        RowBox[{"s", "-",
         RowBox[{"4",
          RowBox[{"m", "^", "2"}]}]}], ")"}], "/",
       RowBox[{"(",
        RowBox[{"t", "-",
         RowBox[{"m", "^", "2"}]}], ")"}]}], "/",
      RowBox[{"(",
       RowBox[{"u", "-",
        RowBox[{"m", "^", "2"}]}], ")"}]}]}], " ", "/.", " ",
    RowBox[{"g", "\[Rule]",
     RowBox[{"Sqrt", "[",
      RowBox[{"4", " ", "Pi", " ", "\[Alpha]"}], "]"}]}]}]}], ";"}], "\n",
 RowBox[{"TrickMandelstam", "[",
  RowBox[{
   RowBox[{"T12", "-", "T12exp"}], ",",
   RowBox[{"{",
    RowBox[{"s", ",", "t", ",", "u", ",",
     RowBox[{"2",
      RowBox[{"m", "^", "2"}]}]}], "}"}]}], "]"}]}], "Input",
 CellChangeTimes->{{3.537296905125*^9, 3.537296949359375*^9}, {
  3.537296995984375*^9, 3.537297013796875*^9}, {3.537297054953125*^9,
  3.537297100453125*^9}, {3.5372971809375*^9, 3.537297276703125*^9}, {
  3.537297359734375*^9, 3.53729738140625*^9}}],

Cell[BoxData[
 FormBox["0", TraditionalForm]], "Output",
 CellChangeTimes->{{3.537297006546875*^9, 3.537297084875*^9}, {
   3.5372972189375*^9, 3.537297277578125*^9}, {3.537297364828125*^9,
   3.53729738209375*^9}, 3.537302449875*^9, 3.537532353796875*^9,
   3.5376106443125*^9, 3.53761105984375*^9, 3.537611476375*^9,
   3.53770553903125*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"M3", "=",
   RowBox[{
    RowBox[{
     RowBox[{
      RowBox[{"I", " ",
       RowBox[{"g", "^", "2"}], "*",
       RowBox[{
        RowBox[{"Spinor", "[",
         RowBox[{
          RowBox[{"-", "p2"}], ",", "m"}], "]"}], ".",
        RowBox[{"GA", "[", "\[Mu]", "]"}], ".",
        RowBox[{"Spinor", "[",
         RowBox[{"p1", ",", "m"}], "]"}]}], "*",
       RowBox[{"Conjugate", "[",
        RowBox[{"PolarizationVector", "[",
         RowBox[{"p3", ",", "\[Nu]"}], "]"}], "]"}], "*",
       RowBox[{"Conjugate", "[",
        RowBox[{"PolarizationVector", "[",
         RowBox[{"p4", ",", "\[Sigma]"}], "]"}], "]"}], "*",
       RowBox[{
        RowBox[{"(",
         RowBox[{"(",
          RowBox[{
           RowBox[{
            RowBox[{"MT", "[",
             RowBox[{"\[Mu]", ",", "\[Nu]"}], "]"}], " ",
            RowBox[{"(",
             RowBox[{
              RowBox[{"FV", "[",
               RowBox[{"p", ",", "\[Sigma]"}], "]"}], "-",
              RowBox[{"FV", "[",
               RowBox[{"q", ",", "\[Sigma]"}], "]"}]}], ")"}]}], "+",
           RowBox[{
            RowBox[{"MT", "[",
             RowBox[{"\[Nu]", ",", "\[Sigma]"}], "]"}], " ",
            RowBox[{"(",
             RowBox[{
              RowBox[{"FV", "[",
               RowBox[{"q", ",", "\[Mu]"}], "]"}], "-",
              RowBox[{"FV", "[",
               RowBox[{"r", ",", "\[Mu]"}], "]"}]}], ")"}]}], "+",
           RowBox[{
            RowBox[{"MT", "[",
             RowBox[{"\[Mu]", ",", "\[Sigma]"}], "]"}], " ",
            RowBox[{"(",
             RowBox[{
              RowBox[{"FV", "[",
               RowBox[{"r", ",", "\[Nu]"}], "]"}], "-",
              RowBox[{"FV", "[",
               RowBox[{"p", ",", "\[Nu]"}], "]"}]}], ")"}]}]}], ")"}], ")"}],
        "/",
        RowBox[{"SP", "[",
         RowBox[{
          RowBox[{"p3", "+", "p4"}], ",",
          RowBox[{"p3", "+", "p4"}]}], "]"}]}]}], "/.",
      RowBox[{
       RowBox[{"FV", "[",
        RowBox[{"q", ",", "i_"}], "]"}], "\[Rule]",
       RowBox[{"-",
        RowBox[{"FV", "[",
         RowBox[{"p3", ",", "i"}], "]"}]}]}]}], "/.",
     RowBox[{
      RowBox[{"FV", "[",
       RowBox[{"r", ",", "i_"}], "]"}], "\[Rule]",
      RowBox[{"-",
       RowBox[{"FV", "[",
        RowBox[{"p4", ",", "i"}], "]"}]}]}]}], "/.",
    RowBox[{
     RowBox[{"FV", "[",
      RowBox[{"p", ",", "i_"}], "]"}], "\[Rule]",
     RowBox[{
      RowBox[{"FV", "[",
       RowBox[{"p3", ",", "i"}], "]"}], "+",
      RowBox[{"FV", "[",
       RowBox[{"p4", ",", "i"}], "]"}]}]}]}]}], ";"}], "\n",
 RowBox[{
  RowBox[{"M3A", "=",
   RowBox[{
    RowBox[{
     RowBox[{
      RowBox[{"ComplexConjugate", "[", "M3", "]"}], " ", "/.", " ",
      RowBox[{"\[Mu]", "\[Rule]", " ", "\[Mu]2"}]}], " ", "/.", "  ",
     RowBox[{"\[Nu]", "\[Rule]", " ", "\[Nu]2"}]}], " ", "/.",
    RowBox[{"\[Sigma]", "\[Rule]", "\[Sigma]2"}]}]}], ";"}], "\n",
 RowBox[{
  RowBox[{"F3", "=",
   RowBox[{
    RowBox[{"SUNSimplify", "[",
     RowBox[{
      RowBox[{
       RowBox[{"1", "/", "9"}], "*",
       RowBox[{"SUNF", "[",
        RowBox[{"b2", ",", "a3", ",", "a4"}], "]"}], "*",
       RowBox[{"SUNF", "[",
        RowBox[{"b", ",", "a3", ",", "a4"}], "]"}], "*",
       RowBox[{"SUNTrace", "[",
        RowBox[{"SUNT", "[",
         RowBox[{"b", ",", "b2"}], "]"}], "]"}]}], ",",
      RowBox[{"SUNNToCACF", "\[Rule]", "False"}]}], "]"}], " ", "/.", " ",
    RowBox[{"SUNN", "\[Rule]", " ", "3"}]}]}], ";"}], "\n",
 RowBox[{
  RowBox[{"T33", "=",
   RowBox[{
    RowBox[{"(",
     RowBox[{
      RowBox[{
       RowBox[{
        RowBox[{"1", "/",
         RowBox[{"(",
          RowBox[{"16", "Pi", " ",
           RowBox[{"s", "^", "2"}]}], ")"}]}], "*",
        RowBox[{"F3", "/", "4"}], "*",
        RowBox[{"(",
         RowBox[{
          RowBox[{
           RowBox[{
            RowBox[{"M3A", "  ", "M3"}], " ", "//", " ",
            "DoPolarizationSums"}], " ", "//", " ", "Contract"}], " ", "//",
          " ", "FermionSpinSum"}], ")"}]}], " ", "/.", " ",
       RowBox[{"DiracTrace", "\[Rule]", "TR"}]}], " ", "//", " ",
      "Simplify"}], " ", ")"}], "/.", " ",
    RowBox[{"g", "\[Rule]",
     RowBox[{"Sqrt", "[",
      RowBox[{"4", " ", "Pi", " ", "\[Alpha]"}], "]"}]}]}]}], ";"}], "\n",
 RowBox[{
  RowBox[{"T33exp", "=",
   RowBox[{
    RowBox[{
     RowBox[{"1", "/",
      RowBox[{"(",
       RowBox[{"16", "Pi", " ",
        RowBox[{"s", "^", "2"}]}], ")"}]}], "*",
     RowBox[{"g", "^", "4"}], "*",
     RowBox[{"(",
      RowBox[{
       RowBox[{"-", "2"}], "/", "3"}], ")"}], "*",
     RowBox[{
      RowBox[{"(",
       RowBox[{
        RowBox[{"s", "^", "2"}], "+",
        RowBox[{"7", "*",
         RowBox[{"(",
          RowBox[{
           RowBox[{"t", "^", "2"}], "+",
           RowBox[{"u", "^", "2"}]}], ")"}]}], "+",
        RowBox[{"4", "t", "*", "u"}], "+",
        RowBox[{"34",
         RowBox[{"m", "^", "2"}], " ", "s"}], "-",
        RowBox[{"18",
         RowBox[{"m", "^", "4"}]}]}], ")"}], "/",
      RowBox[{"s", "^", "2"}]}]}], "  ", "/.", " ",
    RowBox[{"g", "\[Rule]",
     RowBox[{"Sqrt", "[",
      RowBox[{"4", " ", "Pi", " ", "\[Alpha]"}], "]"}]}]}]}], ";"}], "\n",
 RowBox[{
  RowBox[{"TrickMandelstam", "[",
   RowBox[{
    RowBox[{"T33", "-", "T33exp"}], ",",
    RowBox[{"{",
     RowBox[{"s", ",", "t", ",", "u", ",",
      RowBox[{"2",
       RowBox[{"m", "^", "2"}]}]}], "}"}]}], "]"}], " "}]}], "Input",
 CellChangeTimes->CompressedData["
1:eJxTTMoPSmViYGAQA2IQzdCgz3Un6pUjQ8MWHjDN8EAAQt8Wg4h7SkPoL8oQ
WswHIm8bAOELh0P456IhfKsUCD0nFULXZUPo9mII/b4KQj9og9B3eyH63fog
9PWpEPrwNIg863IIf+kKCJ20FiKevgNC5+6FiDcegrrnIoT/G0rXdd0F07sh
dEP2Tr9oEM21D0wz7D0E4U85A+HH3YTQLz9A6H3fIPLHf0P4LMz+YNocSi/m
hNBmvBD6khCEzhcG0QBnpJjl
  "]],

Cell[BoxData[
 FormBox["0", TraditionalForm]], "Output",
 CellChangeTimes->{
  3.53753289628125*^9, 3.53753295765625*^9, {3.537533011015625*^9,
   3.53753303509375*^9}, 3.53761064803125*^9, 3.537611063140625*^9,
   3.53761148015625*^9, 3.537705542109375*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"F13", "=",
   RowBox[{
    RowBox[{"SUNSimplify", "[",
     RowBox[{
      RowBox[{"CalcColorFactor", "[",
       RowBox[{
        RowBox[{"1", "/", "9"}], "*",
        RowBox[{"SUNF", "[",
         RowBox[{"b", ",", "a3", ",", "a4"}], "]"}], "*",
        RowBox[{"SUNTrace", "[",
         RowBox[{"SUNT", "[",
          RowBox[{"a3", ",", "a4", ",", "b"}], "]"}], "]"}]}], "]"}], ",",
      RowBox[{"SUNNToCACF", "\[Rule]", "False"}]}], "]"}], " ", "/.", " ",
    RowBox[{"SUNN", "\[Rule]", " ", "3"}]}]}], ";"}], "\n",
 RowBox[{
  RowBox[{"T13", "=",
   RowBox[{
    RowBox[{"(",
     RowBox[{
      RowBox[{
       RowBox[{
        RowBox[{"1", "/",
         RowBox[{"(",
          RowBox[{"16", "Pi", " ",
           RowBox[{"s", "^", "2"}]}], ")"}]}], "*",
        RowBox[{"F13", "/", "4"}], "*",
        RowBox[{"(",
         RowBox[{
          RowBox[{
           RowBox[{
            RowBox[{"M1A", "  ", "M3"}], " ", "//", " ",
            "DoPolarizationSums"}], " ", "//", " ", "Contract"}], " ", "//",
          " ", "FermionSpinSum"}], ")"}]}], " ", "/.", " ",
       RowBox[{"DiracTrace", "\[Rule]", "TR"}]}], " ", "//", " ",
      "Simplify"}], " ", ")"}], "/.", " ",
    RowBox[{"g", "\[Rule]",
     RowBox[{"Sqrt", "[",
      RowBox[{"4", " ", "Pi", " ", "\[Alpha]"}], "]"}]}]}]}], ";"}], "\n",
 RowBox[{
  RowBox[{"T13exp", "=",
   RowBox[{
    RowBox[{
     RowBox[{
      RowBox[{
       RowBox[{"-", "4"}], "/", "3"}], "/",
      RowBox[{"(",
       RowBox[{"16", "Pi", " ",
        RowBox[{"s", "^", "2"}]}], ")"}]}], "*",
     RowBox[{"g", "^", "4"}], "*",
     RowBox[{
      RowBox[{
       RowBox[{"(",
        RowBox[{
         RowBox[{"t", "^", "2"}], "+",
         RowBox[{
          RowBox[{"m", "^", "2"}], "*",
          RowBox[{"(",
           RowBox[{
            RowBox[{"m", "^", "2"}], "+", "s", "-",
            RowBox[{"2", "t"}]}], ")"}]}]}], ")"}], "/",
       RowBox[{"(",
        RowBox[{"t", "-",
         RowBox[{"m", "^", "2"}]}], ")"}]}], "/", "s"}]}], " ", "/.", " ",
    RowBox[{"g", "\[Rule]",
     RowBox[{"Sqrt", "[",
      RowBox[{"4", " ", "Pi", " ", "\[Alpha]"}], "]"}]}]}]}], ";"}], "\n",
 RowBox[{"TrickMandelstam", "[",
  RowBox[{
   RowBox[{"T13", "-", "T13exp"}], ",",
   RowBox[{"{",
    RowBox[{"s", ",", "t", ",", "u", ",",
     RowBox[{"2",
      RowBox[{"m", "^", "2"}]}]}], "}"}]}], "]"}]}], "Input",
 CellChangeTimes->{{3.53729901403125*^9, 3.537299032125*^9}, {
   3.537299126046875*^9, 3.537299184453125*^9}, {3.537299296140625*^9,
   3.537299308703125*^9}, {3.53729933875*^9, 3.53729946165625*^9}, {
   3.53729951040625*^9, 3.5372995545625*^9}, {3.537299746*^9,
   3.537299755921875*^9}, {3.537300120171875*^9, 3.537300121703125*^9}, {
   3.53730245803125*^9, 3.53730246690625*^9}, {3.53760792071875*^9,
   3.53760792271875*^9}, 3.537610359109375*^9}],

Cell[BoxData[
 FormBox["0", TraditionalForm]], "Output",
 CellChangeTimes->{{3.537299372984375*^9, 3.5372994625*^9}, {
   3.53729951128125*^9, 3.53729955565625*^9}, {3.5372997484375*^9,
   3.537299760078125*^9}, 3.537300123359375*^9, {3.537302454171875*^9,
   3.537302468703125*^9}, 3.53753256621875*^9, 3.53753296225*^9,
   3.53753303821875*^9, 3.53760792459375*^9, 3.53761036159375*^9,
   3.53761065484375*^9, 3.5376110658125*^9, 3.5376114834375*^9,
   3.537705545515625*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"F23", "=",
   RowBox[{
    RowBox[{"SUNSimplify", "[",
     RowBox[{
      RowBox[{"CalcColorFactor", "[",
       RowBox[{
        RowBox[{"1", "/", "9"}], "*",
        RowBox[{"SUNF", "[",
         RowBox[{"b", ",", "a3", ",", "a4"}], "]"}], "*",
        RowBox[{"SUNTrace", "[",
         RowBox[{"SUNT", "[",
          RowBox[{"a4", ",", "a3", ",", "b"}], "]"}], "]"}]}], "]"}], ",",
      RowBox[{"SUNNToCACF", "\[Rule]", "False"}]}], "]"}], " ", "/.", " ",
    RowBox[{"SUNN", "\[Rule]", " ", "3"}]}]}], ";"}], "\n",
 RowBox[{
  RowBox[{"T23", "=",
   RowBox[{
    RowBox[{"(",
     RowBox[{
      RowBox[{
       RowBox[{
        RowBox[{"1", "/",
         RowBox[{"(",
          RowBox[{"16", "Pi", " ",
           RowBox[{"s", "^", "2"}]}], ")"}]}], "*",
        RowBox[{"F23", "/", "4"}], "*",
        RowBox[{"(",
         RowBox[{
          RowBox[{
           RowBox[{
            RowBox[{"M2A", "  ", "M3"}], " ", "//", " ",
            "DoPolarizationSums"}], " ", "//", " ", "Contract"}], " ", "//",
          " ", "FermionSpinSum"}], ")"}]}], " ", "/.", " ",
       RowBox[{"DiracTrace", "\[Rule]", "TR"}]}], " ", "//", " ",
      "Simplify"}], " ", ")"}], "/.", " ",
    RowBox[{"g", "\[Rule]",
     RowBox[{"Sqrt", "[",
      RowBox[{"4", " ", "Pi", " ", "\[Alpha]"}], "]"}]}]}]}], ";"}], "\n",
 RowBox[{
  RowBox[{"T23exp", "=",
   RowBox[{
    RowBox[{
     RowBox[{
      RowBox[{
       RowBox[{"-", "4"}], "/", "3"}], "/",
      RowBox[{"(",
       RowBox[{"16", "Pi", " ",
        RowBox[{"s", "^", "2"}]}], ")"}]}], "*",
     RowBox[{"g", "^", "4"}], "*",
     RowBox[{
      RowBox[{
       RowBox[{"(",
        RowBox[{
         RowBox[{"u", "^", "2"}], "+",
         RowBox[{
          RowBox[{"m", "^", "2"}], "*",
          RowBox[{"(",
           RowBox[{
            RowBox[{"m", "^", "2"}], "+", "s", "-",
            RowBox[{"2", "u"}]}], ")"}]}]}], ")"}], "/",
       RowBox[{"(",
        RowBox[{"u", "-",
         RowBox[{"m", "^", "2"}]}], ")"}]}], "/", "s"}]}], " ", "/.", " ",
    RowBox[{"g", "\[Rule]",
     RowBox[{"Sqrt", "[",
      RowBox[{"4", " ", "Pi", " ", "\[Alpha]"}], "]"}]}]}]}], ";"}], "\n",
 RowBox[{"TrickMandelstam", "[",
  RowBox[{
   RowBox[{"T23", "-", "T23exp"}], ",",
   RowBox[{"{",
    RowBox[{"s", ",", "t", ",", "u", ",",
     RowBox[{"2",
      RowBox[{"m", "^", "2"}]}]}], "}"}]}], "]"}]}], "Input",
 CellChangeTimes->{{3.53729981834375*^9, 3.53729985040625*^9}, {
  3.537299916375*^9, 3.537299931984375*^9}, {3.53729999546875*^9,
  3.537300016765625*^9}, {3.53730010834375*^9, 3.537300129*^9}, {
  3.537300162734375*^9, 3.53730020834375*^9}, {3.537300296890625*^9,
  3.537300303578125*^9}, {3.5373003368125*^9, 3.537300379046875*^9}}],

Cell[BoxData[
 FormBox["0", TraditionalForm]], "Output",
 CellChangeTimes->{
  3.53729997190625*^9, {3.537300004921875*^9, 3.537300029140625*^9},
   3.537300078015625*^9, 3.537300131359375*^9, {3.537300178953125*^9,
   3.537300211546875*^9}, 3.5373003051875*^9, {3.53730034071875*^9,
   3.537300380609375*^9}, 3.537302473578125*^9, 3.537532567875*^9,
   3.53753296503125*^9, 3.5375330399375*^9, 3.537610378109375*^9,
   3.53761065703125*^9, 3.537611068890625*^9, 3.537611485875*^9,
   3.537705549*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"G", "=",
   RowBox[{
    RowBox[{"-", "I"}], " ",
    RowBox[{"g", "^", "2"}], "*",
    RowBox[{
     RowBox[{
      RowBox[{"Spinor", "[",
       RowBox[{
        RowBox[{"-", "p2"}], ",", "m"}], "]"}], ".",
      RowBox[{"GS", "[", "p3", "]"}], ".",
      RowBox[{"Spinor", "[",
       RowBox[{"p1", ",", "m"}], "]"}]}], "/", "s"}]}]}], ";"}], "\n",
 RowBox[{
  RowBox[{"Ga", "=",
   RowBox[{"ComplexConjugate", "[", "G", "]"}]}], ";"}], "\n",
 RowBox[{
  RowBox[{"FG", "=",
   RowBox[{
    RowBox[{"SUNSimplify", "[",
     RowBox[{
      RowBox[{
       RowBox[{"1", "/", "9"}], "*",
       RowBox[{"SUNF", "[",
        RowBox[{"b2", ",", "a3", ",", "a4"}], "]"}], "*",
       RowBox[{"SUNF", "[",
        RowBox[{"b", ",", "a3", ",", "a4"}], "]"}], "*",
       RowBox[{"SUNTrace", "[",
        RowBox[{"SUNT", "[",
         RowBox[{"b2", ",", "b"}], "]"}], "]"}]}], ",",
      RowBox[{"SUNNToCACF", "\[Rule]", "False"}]}], "]"}], " ", "/.", " ",
    RowBox[{"SUNN", "\[Rule]", " ", "3"}]}]}], ";"}], "\n",
 RowBox[{
  RowBox[{"TG", "=",
   RowBox[{
    RowBox[{"-",
     RowBox[{"(",
      RowBox[{
       RowBox[{
        RowBox[{
         RowBox[{"1", "/",
          RowBox[{"(",
           RowBox[{"16", "Pi", " ",
            RowBox[{"s", "^", "2"}]}], ")"}]}], "*",
         RowBox[{"FG", "/", "4"}], "*",
         RowBox[{"(",
          RowBox[{
           RowBox[{
            RowBox[{"Ga", " ", "G"}], " ", "//", " ", "Contract"}], " ", "//",
            " ", "FermionSpinSum"}], ")"}]}], " ", "/.", " ",
        RowBox[{"DiracTrace", "\[Rule]", "TR"}]}], " ", "//", " ",
       "Simplify"}], " ", ")"}]}], "/.", " ",
    RowBox[{"g", "\[Rule]",
     RowBox[{"Sqrt", "[",
      RowBox[{"4", " ", "Pi", " ", "\[Alpha]"}], "]"}]}]}]}], ";"}], "\n",
 RowBox[{
  RowBox[{"TGexp", "=",
   RowBox[{
    RowBox[{
     RowBox[{
      RowBox[{"-", "1"}], "/",
      RowBox[{"(",
       RowBox[{"16", "Pi", " ",
        RowBox[{"s", "^", "2"}]}], ")"}]}], "*",
     RowBox[{"g", "^", "4"}], "*",
     RowBox[{"(",
      RowBox[{"2", "/", "3"}], ")"}], "*",
     RowBox[{"(",
      RowBox[{
       RowBox[{"m", "^", "2"}], "-", "t"}], ")"}], "*",
     RowBox[{
      RowBox[{"(",
       RowBox[{
        RowBox[{"m", "^", "2"}], "-", "u"}], ")"}], "/",
      RowBox[{"s", "^", "2"}]}]}], "  ", "/.", " ",
    RowBox[{"g", "\[Rule]",
     RowBox[{"Sqrt", "[",
      RowBox[{"4", " ", "Pi", " ", "\[Alpha]"}], "]"}]}]}]}], ";"}], "\n",
 RowBox[{"TrickMandelstam", "[",
  RowBox[{
   RowBox[{"TG", "-", "TGexp"}], ",",
   RowBox[{"{",
    RowBox[{"s", ",", "t", ",", "u", ",",
     RowBox[{"2",
      RowBox[{"m", "^", "2"}]}]}], "}"}]}], "]"}]}], "Input",
 CellChangeTimes->{{3.53730071371875*^9, 3.5373009011875*^9}, {
   3.5373009396875*^9, 3.537300956578125*^9}, {3.537301056140625*^9,
   3.537301197453125*^9}, {3.53730123828125*^9, 3.537301308703125*^9}, {
   3.53730134125*^9, 3.537301345484375*^9}, {3.537301427375*^9,
   3.5373014384375*^9}, {3.537301473984375*^9, 3.53730147496875*^9}, {
   3.537301618265625*^9, 3.5373016850625*^9}, {3.537302883390625*^9,
   3.537302894015625*^9}, {3.537302937703125*^9, 3.537303006640625*^9},
   3.537303211828125*^9, 3.53730524946875*^9, {3.537532579078125*^9,
   3.537532581390625*^9}}],

Cell[BoxData[
 FormBox["0", TraditionalForm]], "Output",
 CellChangeTimes->{{3.537301159*^9, 3.537301198421875*^9},
   3.53730124228125*^9, {3.53730130653125*^9, 3.53730130928125*^9}, {
   3.537301346234375*^9, 3.537301366140625*^9}, {3.53730142840625*^9,
   3.537301439203125*^9}, 3.537301475671875*^9, 3.537301545515625*^9, {
   3.537301644734375*^9, 3.53730168590625*^9}, 3.537302495890625*^9, {
   3.537302886734375*^9, 3.537303007125*^9}, 3.537303212359375*^9,
   3.537305251171875*^9, 3.53753258371875*^9, 3.537532970125*^9,
   3.537533046109375*^9, 3.537610406546875*^9, 3.53761065865625*^9,
   3.5376110705625*^9, 3.53761148709375*^9, 3.537705556890625*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"fin", "=",
  RowBox[{
   RowBox[{
    RowBox[{"1", "/",
     RowBox[{"(",
      RowBox[{"16", "Pi", " ",
       RowBox[{"s", "^", "2"}]}], ")"}]}], "*",
    RowBox[{"2", "^", "5"}], "*",
    RowBox[{
     RowBox[{"g", "^", "4"}], "/",
     RowBox[{"3", "^", "3"}]}], "*",
    RowBox[{"(",
     RowBox[{
      RowBox[{
       RowBox[{"(",
        RowBox[{
         RowBox[{"t", "*", "u"}], "+",
         RowBox[{
          RowBox[{"m", "^", "2"}], "*",
          RowBox[{"(",
           RowBox[{"s", "-",
            RowBox[{"2", "t"}], "-",
            RowBox[{"3",
             RowBox[{"m", "^", "2"}]}]}], ")"}]}]}], ")"}], "/",
       RowBox[{
        RowBox[{"(",
         RowBox[{"t", "-",
          RowBox[{"m", "^", "2"}]}], ")"}], "^", "2"}]}], "+",
      RowBox[{
       RowBox[{"(",
        RowBox[{
         RowBox[{"t", "*", "u"}], "+",
         RowBox[{
          RowBox[{"m", "^", "2"}], "*",
          RowBox[{"(",
           RowBox[{"s", "-",
            RowBox[{"2", "u"}], "-",
            RowBox[{"3",
             RowBox[{"m", "^", "2"}]}]}], ")"}]}]}], ")"}], "/",
       RowBox[{
        RowBox[{"(",
         RowBox[{"u", "-",
          RowBox[{"m", "^", "2"}]}], ")"}], "^", "2"}]}], "+",
      RowBox[{
       RowBox[{"m", "^", "2"}], "*",
       RowBox[{
        RowBox[{
         RowBox[{"(",
          RowBox[{
           RowBox[{"m", "^", "2"}], "+",
           RowBox[{"2", "*", "s"}]}], ")"}], "/",
         RowBox[{"(",
          RowBox[{"t", "-",
           RowBox[{"m", "^", "2"}]}], ")"}]}], "/",
        RowBox[{"(",
         RowBox[{"u", "-",
          RowBox[{"m", "^", "2"}]}], ")"}]}]}], "-",
      RowBox[{
       RowBox[{"9", "/", "4"}], "*",
       RowBox[{
        RowBox[{"(",
         RowBox[{
          RowBox[{"t", "^", "2"}], "+",
          RowBox[{"u", "^", "2"}], "+",
          RowBox[{"6",
           RowBox[{"m", "^", "2"}], "s"}], "-",
          RowBox[{"2",
           RowBox[{"m", "^", "4"}]}]}], ")"}], "/",
        RowBox[{"s", "^", "2"}]}]}]}], ")"}]}], "/.", " ",
   RowBox[{"g", "\[Rule]",
    RowBox[{"Sqrt", "[",
     RowBox[{"4", " ", "Pi", " ", "\[Alpha]"}], "]"}]}]}]}]], "Input",
 CellChangeTimes->{{3.5377119749375*^9, 3.537712015140625*^9}, {
  3.53771211590625*^9, 3.537712186171875*^9}, {3.53771223128125*^9,
  3.5377122335625*^9}, {3.537712301921875*^9, 3.537712314890625*^9}, {
  3.537712371375*^9, 3.53771237915625*^9}}],

Cell[BoxData[
 FormBox[
  FractionBox[
   RowBox[{"32", " ", "\[Pi]", " ",
    SuperscriptBox["\[Alpha]", "2"], " ",
    RowBox[{"(",
     RowBox[{
      FractionBox[
       RowBox[{
        SuperscriptBox["m", "2"], " ",
        RowBox[{"(",
         RowBox[{
          SuperscriptBox["m", "2"], "+",
          RowBox[{"2", " ", "s"}]}], ")"}]}],
       RowBox[{
        RowBox[{"(",
         RowBox[{"t", "-",
          SuperscriptBox["m", "2"]}], ")"}], " ",
        RowBox[{"(",
         RowBox[{"u", "-",
          SuperscriptBox["m", "2"]}], ")"}]}]], "+",
      FractionBox[
       RowBox[{
        RowBox[{
         SuperscriptBox["m", "2"], " ",
         RowBox[{"(",
          RowBox[{
           RowBox[{"-",
            RowBox[{"3", " ",
             SuperscriptBox["m", "2"]}]}], "+", "s", "-",
           RowBox[{"2", " ", "t"}]}], ")"}]}], "+",
        RowBox[{"t", " ", "u"}]}],
       SuperscriptBox[
        RowBox[{"(",
         RowBox[{"t", "-",
          SuperscriptBox["m", "2"]}], ")"}], "2"]], "+",
      FractionBox[
       RowBox[{
        RowBox[{
         SuperscriptBox["m", "2"], " ",
         RowBox[{"(",
          RowBox[{
           RowBox[{"-",
            RowBox[{"3", " ",
             SuperscriptBox["m", "2"]}]}], "+", "s", "-",
           RowBox[{"2", " ", "u"}]}], ")"}]}], "+",
        RowBox[{"t", " ", "u"}]}],
       SuperscriptBox[
        RowBox[{"(",
         RowBox[{"u", "-",
          SuperscriptBox["m", "2"]}], ")"}], "2"]], "-",
      FractionBox[
       RowBox[{"9", " ",
        RowBox[{"(",
         RowBox[{
          RowBox[{"-",
           RowBox[{"2", " ",
            SuperscriptBox["m", "4"]}]}], "+",
          RowBox[{"6", " ",
           SuperscriptBox["m", "2"], " ", "s"}], "+",
          SuperscriptBox["t", "2"], "+",
          SuperscriptBox["u", "2"]}], ")"}]}],
       RowBox[{"4", " ",
        SuperscriptBox["s", "2"]}]]}], ")"}]}],
   RowBox[{"27", " ",
    SuperscriptBox["s", "2"]}]], TraditionalForm]], "Output",
 CellChangeTimes->{{3.5377121816875*^9, 3.537712186953125*^9},
   3.537712242140625*^9, 3.53771231709375*^9, 3.537712388015625*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"finb", "=",
  RowBox[{"T11", "+", "T22", "+",
   RowBox[{"2", "*", "T12exp"}], "+",
   RowBox[{"2", "*", "T13exp"}], "+",
   RowBox[{"2", "*", "T23exp"}], "+", "T33exp", "+",
   RowBox[{"2", "*", "TGexp"}]}]}]], "Input",
 CellChangeTimes->{{3.537712192109375*^9, 3.537712240078125*^9}, {
  3.537712365359375*^9, 3.53771238665625*^9}}],

Cell[BoxData[
 FormBox[
  RowBox[{
   RowBox[{"-",
    FractionBox[
     RowBox[{"4", " ", "\[Pi]", " ",
      SuperscriptBox["\[Alpha]", "2"], " ",
      RowBox[{"(",
       RowBox[{
        SuperscriptBox["m", "2"], "-", "t"}], ")"}], " ",
      RowBox[{"(",
       RowBox[{
        SuperscriptBox["m", "2"], "-", "u"}], ")"}]}],
     RowBox[{"3", " ",
      SuperscriptBox["s", "4"]}]]}], "-",
   FractionBox[
    RowBox[{"8", " ", "\[Pi]", " ",
     SuperscriptBox["\[Alpha]", "2"], " ",
     RowBox[{"(",
      RowBox[{
       RowBox[{
        SuperscriptBox["m", "2"], " ",
        RowBox[{"(",
         RowBox[{
          SuperscriptBox["m", "2"], "+", "s", "-",
          RowBox[{"2", " ", "t"}]}], ")"}]}], "+",
       SuperscriptBox["t", "2"]}], ")"}]}],
    RowBox[{"3", " ",
     SuperscriptBox["s", "3"], " ",
     RowBox[{"(",
      RowBox[{"t", "-",
       SuperscriptBox["m", "2"]}], ")"}]}]], "-",
   FractionBox[
    RowBox[{"8", " ", "\[Pi]", " ",
     SuperscriptBox["\[Alpha]", "2"], " ",
     RowBox[{"(",
      RowBox[{
       RowBox[{
        SuperscriptBox["m", "2"], " ",
        RowBox[{"(",
         RowBox[{
          SuperscriptBox["m", "2"], "+", "s", "-",
          RowBox[{"2", " ", "u"}]}], ")"}]}], "+",
       SuperscriptBox["u", "2"]}], ")"}]}],
    RowBox[{"3", " ",
     SuperscriptBox["s", "3"], " ",
     RowBox[{"(",
      RowBox[{"u", "-",
       SuperscriptBox["m", "2"]}], ")"}]}]], "-",
   FractionBox[
    RowBox[{"8", " ", "\[Pi]", " ",
     SuperscriptBox["\[Alpha]", "2"], " ",
     SuperscriptBox["m", "2"], " ",
     RowBox[{"(",
      RowBox[{"s", "-",
       RowBox[{"4", " ",
        SuperscriptBox["m", "2"]}]}], ")"}]}],
    RowBox[{"27", " ",
     SuperscriptBox["s", "2"], " ",
     RowBox[{"(",
      RowBox[{"t", "-",
       SuperscriptBox["m", "2"]}], ")"}], " ",
     RowBox[{"(",
      RowBox[{"u", "-",
       SuperscriptBox["m", "2"]}], ")"}]}]], "+",
   FractionBox[
    RowBox[{"32", " ", "\[Pi]", " ",
     SuperscriptBox["\[Alpha]", "2"], " ",
     RowBox[{"(",
      RowBox[{
       RowBox[{"3", " ",
        SuperscriptBox["m", "4"]}], "-",
       RowBox[{
        SuperscriptBox["m", "2"], " ",
        RowBox[{"(",
         RowBox[{
          RowBox[{"2", " ", "s"}], "+",
          RowBox[{"5", " ", "t"}], "+",
          RowBox[{"3", " ", "u"}]}], ")"}]}], "+",
       RowBox[{"t", " ", "u"}]}], ")"}]}],
    RowBox[{"27", " ",
     SuperscriptBox["s", "2"], " ",
     SuperscriptBox[
      RowBox[{"(",
       RowBox[{
        SuperscriptBox["m", "2"], "-", "t"}], ")"}], "2"]}]], "+",
   FractionBox[
    RowBox[{"32", " ", "\[Pi]", " ",
     SuperscriptBox["\[Alpha]", "2"], " ",
     RowBox[{"(",
      RowBox[{
       RowBox[{"3", " ",
        SuperscriptBox["m", "4"]}], "-",
       RowBox[{
        SuperscriptBox["m", "2"], " ",
        RowBox[{"(",
         RowBox[{
          RowBox[{"2", " ", "s"}], "+",
          RowBox[{"3", " ", "t"}], "+",
          RowBox[{"5", " ", "u"}]}], ")"}]}], "+",
       RowBox[{"t", " ", "u"}]}], ")"}]}],
    RowBox[{"27", " ",
     SuperscriptBox["s", "2"], " ",
     SuperscriptBox[
      RowBox[{"(",
       RowBox[{
        SuperscriptBox["m", "2"], "-", "u"}], ")"}], "2"]}]], "-",
   FractionBox[
    RowBox[{"2", " ", "\[Pi]", " ",
     SuperscriptBox["\[Alpha]", "2"], " ",
     RowBox[{"(",
      RowBox[{
       RowBox[{"-",
        RowBox[{"18", " ",
         SuperscriptBox["m", "4"]}]}], "+",
       RowBox[{"34", " ",
        SuperscriptBox["m", "2"], " ", "s"}], "+",
       SuperscriptBox["s", "2"], "+",
       RowBox[{"7", " ",
        RowBox[{"(",
         RowBox[{
          SuperscriptBox["t", "2"], "+",
          SuperscriptBox["u", "2"]}], ")"}]}], "+",
       RowBox[{"4", " ", "t", " ", "u"}]}], ")"}]}],
    RowBox[{"3", " ",
     SuperscriptBox["s", "4"]}]]}], TraditionalForm]], "Output",
 CellChangeTimes->{{3.53771222409375*^9, 3.5377122433125*^9}, {
   3.537712319671875*^9, 3.53771234421875*^9}, 3.5377123895*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"TrickMandelstam", "[",
  RowBox[{
   RowBox[{"fin", "-", "finb"}], ",",
   RowBox[{"{",
    RowBox[{"s", ",", "t", ",", "u", ",",
     RowBox[{"2",
      RowBox[{"m", "^", "2"}]}]}], "}"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.537712244875*^9, 3.5377122804375*^9}, {
  3.53771233996875*^9, 3.537712341359375*^9}}],

Cell[BoxData[
 FormBox["0", TraditionalForm]], "Output",
 CellChangeTimes->{{3.537712246515625*^9, 3.5377122810625*^9}, {
   3.537712320765625*^9, 3.537712345078125*^9}, 3.53771239078125*^9}]
}, Open  ]]
},
WindowSize->{1016, 668},
WindowMargins->{{-5, Automatic}, {Automatic, 2}},
FrontEndVersion->"8.0 for Microsoft Windows (32-bit) (November 7, 2010)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[557, 20, 71, 1, 31, "Input"],
Cell[631, 23, 287, 7, 16, "Text"],
Cell[921, 32, 762, 24, 16, "Text"],
Cell[1686, 58, 369, 12, 16, "Text"],
Cell[2058, 72, 195, 4, 16, "Text"],
Cell[2256, 78, 310, 7, 31, "Input"],
Cell[CellGroupData[{
Cell[2591, 89, 3819, 113, 172, "Input"],
Cell[6413, 204, 252, 4, 30, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[6702, 213, 3726, 112, 172, "Input"],
Cell[10431, 327, 278, 5, 30, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[10746, 337, 2412, 72, 112, "Input"],
Cell[13161, 411, 346, 6, 30, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[13544, 422, 5882, 170, 212, "Input"],
Cell[19429, 594, 260, 5, 30, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[19726, 604, 2910, 81, 112, "Input"],
Cell[22639, 687, 482, 8, 30, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[23158, 700, 2790, 79, 112, "Input"],
Cell[25951, 781, 507, 9, 30, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[26495, 795, 3331, 94, 132, "Input"],
Cell[29829, 891, 670, 10, 30, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[30536, 906, 2472, 76, 92, "Input"],
Cell[33011, 984, 2169, 67, 72, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[35217, 1056, 363, 8, 31, "Input"],
Cell[35583, 1066, 4062, 127, 105, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[39682, 1198, 346, 9, 31, "Input"],
Cell[40031, 1209, 191, 3, 30, "Output"]
}, Open  ]]
}
]
*)

(* End of internal cache information *)