(************** Content-type: application/mathematica **************

                    Mathematica-Compatible Notebook

This notebook can be used with any Mathematica-compatible
application, such as Mathematica, MathReader or Publicon. The data
for the notebook starts with the line containing stars above.

To get the notebook into a Mathematica-compatible application, do
one of the following:

* Save the data starting with the line of stars above into a file
  with a name ending in .nb, then open the file inside the
  application;

* Copy the data starting with the line of stars above to the
  clipboard, then use the Paste menu command inside the application.

Data for notebooks contains only printable 7-bit ASCII and can be
sent directly in email or through ftp in text mode.  Newlines can be
CR, LF or CRLF (Unix, Macintosh or MS-DOS style).

NOTE: If you modify the data for this notebook not in a Mathematica-
compatible application, you must delete the line below containing
the word CacheID, otherwise Mathematica-compatible applications may
try to use invalid cache data.

For more information on notebooks and Mathematica-compatible 
applications, contact Wolfram Research:
  web: http://www.wolfram.com
  email: info@wolfram.com
  phone: +1-217-398-0700 (U.S.)

Notebook reader applications are available free of charge from 
Wolfram Research.
*******************************************************************)

(*CacheID: 232*)


(*NotebookFileLineBreakTest
NotebookFileLineBreakTest*)
(*NotebookOptionsPosition[     32589,        988]*)
(*NotebookOutlinePosition[     33288,       1012]*)
(*  CellTagsIndexPosition[     33244,       1008]*)
(*WindowFrame->Normal*)



Notebook[{

Cell[CellGroupData[{
Cell["gg \[Rule] H at NLO in the EFT", "Title",
  PageWidth->PaperWidth],

Cell[CellGroupData[{

Cell["Input FeynCalc", "Subsection",
  PageWidth->PaperWidth],

Cell[BoxData[
    \(<< HighEnergyPhysics`fc`\)], "Input",
  PageWidth->PaperWidth],

Cell[TextData[{
  StyleBox["FeynCalc",
    FontWeight->"Bold"],
  " ",
  "4.1.0.3b",
  "    ",
  " Evaluate ?FeynCalc for help or visit ",
  ButtonBox["www.feyncalc.org",
    ButtonData:>{
      URL[ "http://www.feyncalc.org"], None},
    ButtonStyle->"Hyperlink",
    ButtonNote->"http://www.feyncalc.org"]
}], "Text",
  GeneratedCell->True,
  CellAutoOverwrite->True],

Cell[BoxData[
    \(\($LimitTo4 = False;\)\)], "Input",
  PageWidth->PaperWidth]
}, Closed]],

Cell[CellGroupData[{

Cell["Virtual diagrams :preliminaries", "Section",
  PageWidth->PaperWidth],

Cell[CellGroupData[{

Cell["Kinematics", "Subsection",
  PageWidth->PaperWidth],

Cell["\<\
-----I take all momenta outgoing  

     p1 + p2 + p3  = 0     
     
     p1^2=0
     p2^2=0
     p3^3=Q^2
     

        \
\>", "Text",
  PageWidth->PaperWidth],

Cell[BoxData[{
    \(\(ScalarProduct[p1, p1] = 0;\)\), "\[IndentingNewLine]", 
    \(\(ScalarProduct[p2, p2] = 0;\)\), "\[IndentingNewLine]", 
    \(\(ScalarProduct[p3, p3] = Q2;\)\), "\[IndentingNewLine]", 
    \(\(ScalarProduct[p1, p3] = \(-Q2\)/2;\)\), "\[IndentingNewLine]", 
    \(\(ScalarProduct[p1, p2] = Q2/2;\)\), "\[IndentingNewLine]", 
    \(\(ScalarProduct[p2, p3] = \(-\ Q2\)/2;\)\), "\[IndentingNewLine]", 
    \(\(ScalarProduct[p1, e1] = 0;\)\), "\[IndentingNewLine]", 
    \(\(ScalarProduct[p2, e2] = 0;\)\), "\[IndentingNewLine]", 
    \(\(ScalarProduct[p1, e2] = 0;\)\), "\[IndentingNewLine]", 
    \(\(\(ScalarProduct[p2, e1] = 0;\)\(\[IndentingNewLine]\)
    \)\), "\[IndentingNewLine]", 
    \(\(ScalarProduct[p1, p1, Dimension \[Rule] D] = 
        0;\)\), "\[IndentingNewLine]", 
    \(\(ScalarProduct[p2, p2, Dimension \[Rule] D] = 
        0;\)\), "\[IndentingNewLine]", 
    \(\(ScalarProduct[p3, p3, Dimension \[Rule] D] = 
        Q2;\)\), "\[IndentingNewLine]", 
    \(\(ScalarProduct[p1, p3, Dimension \[Rule] D] = \(-Q2\)/
          2;\)\), "\[IndentingNewLine]", 
    \(\(ScalarProduct[p1, p2, Dimension \[Rule] D] = 
        Q2/2;\)\), "\[IndentingNewLine]", 
    \(\(ScalarProduct[p2, p3, Dimension \[Rule] D] = \(-\ Q2\)/
          2;\)\), "\[IndentingNewLine]", 
    \(\(ScalarProduct[p1, e1, Dimension \[Rule] D] = 
        0;\)\), "\[IndentingNewLine]", 
    \(\(ScalarProduct[p2, e2, Dimension \[Rule] D] = 
        0;\)\), "\[IndentingNewLine]", 
    \(\(ScalarProduct[p1, e2, Dimension \[Rule] D] = 
        0;\)\), "\[IndentingNewLine]", 
    \(\(ScalarProduct[p2, e1, Dimension \[Rule] D] = 
        0;\)\), "\[IndentingNewLine]", 
    \(\)}], "Input",
  PageWidth->PaperWidth]
}, Closed]],

Cell[CellGroupData[{

Cell["Verteces and Propagators", "Subsection",
  PageWidth->PaperWidth],

Cell[TextData[{
  StyleBox["GGG is the kinematic part of the three-gluon vtx (momenta \
outgoing, clockwise ordering):\nVTX(ggg) = (-\[ImaginaryI] ",
    FontSize->14],
  Cell[BoxData[
      \(TraditionalForm\`g\_s\)],
    FontSize->14],
  StyleBox[" ) (\[ImaginaryI] ",
    FontSize->14],
  Cell[BoxData[
      \(TraditionalForm\`f\^abc\)],
    FontSize->14],
  StyleBox[") GGG\nVTX(qqg) = ( -\[ImaginaryI] ",
    FontSize->14],
  Cell[BoxData[
      \(TraditionalForm\`g\_s\)],
    FontSize->14],
  ")",
  StyleBox[" ",
    FontSize->14],
  Cell[BoxData[
      \(TraditionalForm\`\((T\^a)\)\_ij\)],
    FontSize->14],
  StyleBox[" ",
    FontSize->14],
  Cell[BoxData[
      \(TraditionalForm\`\[Gamma]\^\[Mu]\)],
    FontSize->14],
  StyleBox["\nGluon Propagator= ",
    FontSize->14],
  Cell[BoxData[
      FormBox[
        FractionBox[
          StyleBox[\(\(-\[ImaginaryI]\)\ g\^\[Mu]\[Nu]\),
            FontSize->16], \(p\^2\)], TraditionalForm]],
    FontSize->14],
  "\nQuark ",
  StyleBox["Propagator= ",
    FontSize->14],
  Cell[BoxData[
      FormBox[
        FractionBox[
          StyleBox[\(\(\[ImaginaryI]\)\(\ \)\),
            FontSize->16], \(p\&^\)], TraditionalForm]],
    FontSize->14]
}], "Text",
  PageWidth->PaperWidth],

Cell[BoxData[{
    \(\(GGG[p1_, p2_, p3_, m1_, m2_, m3_] := 
        FV[p1 - p2, m3]\ MT[m1, m2] + FV[p2 - p3, m1]\ MT[m2, m3] + 
          FV[p3 - p1, m2]\ MT[m1, m3];\)\), "\[IndentingNewLine]", 
    \(\(GGGD[p1_, p2_, p3_, m1_, m2_, m3_] := 
        FVD[p1 - p2, m3]\ MTD[m1, m2] + FVD[p2 - p3, m1]\ MTD[m2, m3] + 
          FVD[p3 - p1, m2]\ MTD[m1, m3];\)\), "\[IndentingNewLine]", 
    \(\(GGGG[m1_, m2_, m3_, m4_] := 
        2\ MT[m1, m2]\ MT[m3, m4] - MT[m1, m3]\ MT[m2, m4] - 
          MT[m1, m4]\ MT[m2, m3];\)\), "\[IndentingNewLine]", 
    \(\(GGGGD[m1_, m2_, m3_, m4_] := 
        2\ MTD[m1, m2]\ MTD[m3, m4] - MTD[m1, m3]\ MTD[m2, m4] - 
          MTD[m1, m4]\ MTD[m2, m3];\)\), "\[IndentingNewLine]", 
    \(\(PropQuark = I;\)\), "\[IndentingNewLine]", 
    \(\(PropGluon\  = \(-I\);\)\), "\[IndentingNewLine]", 
    \(\(vtx = \(-I\)\ gs;\)\)}], "Input",
  PageWidth->PaperWidth]
}, Closed]],

Cell[CellGroupData[{

Cell["Born Matrix element ", "Subsection",
  PageWidth->PaperWidth],

Cell[CellGroupData[{

Cell[BoxData[{
    \(factborn = gs2\), "\[IndentingNewLine]", 
    \(\(Born = 
        factborn \((\ 
              MTD[mu, nu]\ SPD[p1, p2] - FVD[p1, nu]\ FVD[p2, mu])\)\ FVD[e1, 
              mu]\ FVD[e2, nu] // Contract;\)\), "\[IndentingNewLine]", 
    \(\(Born4 = Born /. \ D \[Rule] 4;\)\), "\[IndentingNewLine]", 
    \(Born = Born\)}], "Input",
  PageWidth->PaperWidth],

Cell[BoxData[
    \(TraditionalForm\`gs2\)], "Output"],

Cell[BoxData[
    FormBox[
      RowBox[{\(1\/2\), " ", "gs2", " ", "Q2", " ", 
        RowBox[{
          FormBox["e1",
            "TraditionalForm"], "\[NoBreak]", "\[CenterDot]", "\[NoBreak]", 
          FormBox["e2",
            "TraditionalForm"]}]}], TraditionalForm]], "Output"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
    \(BornSq = 
      Normal[Series[\((\[IndentingNewLine]gs2^2\ /\((D - 2)\)^2*\((\ 
                      MTD[mu, nu]\ SPD[p1, p2] - 
                        FVD[p1, nu]\ FVD[p2, mu])\) \((\ 
                      MTD[mu, nu]\ SPD[p1, p2] - 
                        FVD[p1, nu]\ FVD[p2, mu])\)\  // Contract)\)\  /. \ 
              D \[Rule] 4 - 2  e, {e, 0, 2}]] // 
        Simplify\), "\[IndentingNewLine]", 
    \(\(BornSq4 = BornSq /. \ e \[Rule] 0;\)\)}], "Input"],

Cell[BoxData[
    \(TraditionalForm\`1\/8\ \((e\^2 + e + 1)\)\ gs2\^2\ Q2\^2\)], "Output"]
}, Open  ]]
}, Closed]],

Cell[CellGroupData[{

Cell["My Scalar Integrals", "Subsection",
  PageWidth->PaperWidth],

Cell[BoxData[{
    \(\(subInt[expr_] := 
        expr /. \ {\[IndentingNewLine]B0[x__] \[Rule] \ 
              DUPI/\((I\ Pi^2)\)\ MyB0[x], \[IndentingNewLine]C0[
                x__] \[Rule] \ 
              DUPI/\((I\ Pi^2)\)\ MyC0[x], \[IndentingNewLine]D0[
                x__] \[Rule] \ 
              DUPI/\((I\ Pi^2)\)\ MyD0[x]};\)\), "\[IndentingNewLine]", 
    \(\(MyC0[0, 0, Q2, 0, 0, 
          0] = \((c\[CapitalGamma]*\((2/e^2 - Pi^2 - 2/e*Log[Q2/mu2] + 
                  Log[Q2/mu2]^2)\))\)/\((2*
              Q2)\);\)\), "\[IndentingNewLine]", 
    \(\(MyB0[0, 0, 0] = 0;\)\), "\[IndentingNewLine]", 
    \(\(MyB0[Q2, 0, 0] = 
        c\[CapitalGamma]\ \((1/e + 2 - Log[Q2/mu2])\);\)\), "\n", 
    \(\)}], "Input",
  PageWidth->PaperWidth],

Cell[BoxData[""], "Input"],

Cell[BoxData[""], "Input"]
}, Closed]]
}, Closed]],

Cell[CellGroupData[{

Cell["CDR", "Section"],

Cell[CellGroupData[{

Cell["Vertex Diagram ", "Subsection",
  PageWidth->PaperWidth],

Cell[CellGroupData[{

Cell[BoxData[{
    \(\(\(colorvtx = SUNF[b, x, y] SUNF[y, a, x]\  // SUNSimplify\)\(\n\)
    \)\), "\[IndentingNewLine]", 
    \(\(num = 
        gs2^2*\((\(-gs2\))\)*
            I*\[IndentingNewLine]\((\ 
              MTD[al, be]\ SPD[\(-l\) + p1, p2 + l] - \ \ FVD[\(-l\) + p1, 
                    be]\ FVD[p2 + l, al])\)\ *\ 
            GGGD[\(-l\), p1, l - p1, ro, mu, al]*\[IndentingNewLine]GGGD[p2, 
              l, \(-l\) - p2, nu, ro, be]*\((\ 
              MTD[mu, nu]\ SPD[p1, p2] - FVD[p1, nu]\ FVD[p2, mu])\)\  // 
          Contract;\)\), "\n", 
    \(\(dens = FAD[l, p1 - l, l + p2];\)\), "\[IndentingNewLine]", 
    \(\(amp = CA\ dens*\ num/\((D - 2)\)^2;\)\), "\n", 
    \(inte1 = \((\(\((\(OneLoop[l, amp/DUPI] // PaVeReduce\) // Simplify)\) // 
                Factor\) // Simplify)\) /. \ B0[0, 0, 0] \[Rule] 0\  // 
        Simplify\)}], "Input",
  PageWidth->PaperWidth],

Cell[BoxData[
    FormBox[
      RowBox[{\(C\_A\), " ", 
        SubscriptBox["\[Delta]", 
          RowBox[{
            FormBox[
              FormBox["a",
                "TraditionalForm"],
              "TraditionalForm"], "\[NoBreak]", 
            FormBox[
              FormBox["b",
                "TraditionalForm"],
              "TraditionalForm"]}]]}], TraditionalForm]], "Output"],

Cell[BoxData[
    FormBox[
      RowBox[{"-", 
        RowBox[{
          RowBox[{"(", 
            
            RowBox[{\(C\_A\), " ", \(gs2\^3\), " ", \(\[Pi]\^2\), 
              " ", \(Q2\^2\), " ", 
              RowBox[{"(", 
                
                RowBox[{\(\((\(-20\)\ D\^2 + 73\ D - 52)\)\ \(\(B\_0\)(Q2, 0, 
                      0)\)\), "+", 
                  RowBox[{"8", " ", \((D\^2 - 3\ D + 2)\), " ", "Q2", " ", 
                    RowBox[{
                      FormBox[\("C"\_"0"\),
                        "TraditionalForm"], "\[NoBreak]", "(", "\[NoBreak]", 
                      FormBox["0",
                        "TraditionalForm"], "\[NoBreak]", ",", "\[NoBreak]", 
                      FormBox["0",
                        "TraditionalForm"], "\[NoBreak]", ",", "\[NoBreak]", 
                      "Q2", "\[NoBreak]", ",", "\[NoBreak]", 
                      FormBox["0",
                        "TraditionalForm"], "\[NoBreak]", ",", "\[NoBreak]", 
                      FormBox["0",
                        "TraditionalForm"], "\[NoBreak]", ",", "\[NoBreak]", 
                      FormBox["0",
                        "TraditionalForm"], "\[NoBreak]", ")"}]}]}], ")"}]}], 
            ")"}], "/", \((16\ \((D - 2)\)\^2\ \((D - 1)\)\ DUPI)\)}]}], 
      TraditionalForm]], "Output"]
}, Open  ]],

Cell[BoxData[""], "Input"]
}, Open  ]],

Cell[CellGroupData[{

Cell["4-gluon Diagram ", "Subsection",
  PageWidth->PaperWidth],

Cell[CellGroupData[{

Cell[BoxData[{
    \(\(color = SUNF[b, c, y] SUNF[y, a, c]\  // SUNSimplify;\)\), "\n", 
    \(\(num = 
        gs2^2*I*\ 
            gs2*\((\ MTD[ro, si]\ SPD[l, \(-l\) + p1 + p2] - \ \ FVD[l, 
                    si]\ FVD[\(-l\) + p1 + p2, ro])\)\ *\ 
            GGGGD[nu, mu, ro, si]\ *\((\ 
              MTD[mu, nu]\ SPD[p1, p2] - FVD[p1, nu]\ FVD[p2, mu])\)\  // 
          Contract;\)\), "\n", 
    \(\(dens = FAD[l, l - p1 - p2];\)\), "\n", 
    \(\(amp = CA/2\ \ dens*num/\((D - 2)\)^2;\)\), "\[IndentingNewLine]", 
    \(inte2 = \(\((\(OneLoop[l, amp/DUPI] // PaVeReduce\) // Simplify)\) // 
          Factor\) // Simplify\)}], "Input",
  PageWidth->PaperWidth],

Cell[BoxData[
    \(TraditionalForm\`\(-\(\(C\_A\ \((2\ D - 3)\)\ \((2\ D\^2 - 5\ D + 
                4)\)\ gs2\^3\ \[Pi]\^2\ Q2\^2\ \(\(B\_0\)(Q2, 0, 
              0)\)\)\/\(16\ \((D - 2)\)\^2\ \((D - 
                1)\)\ DUPI\)\)\)\)], "Output"]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Sum the virtuals", "Subsection"],

Cell[CellGroupData[{

Cell[BoxData[
    \(\(\(\[IndentingNewLine]\)\(res = \ \(\(\(-1\)/
                2*\[IndentingNewLine]Normal[
                Series[\(\(\((inte1 + inte2)\)/BornSq\)/
                          I\)/\((mu2/Q2)\)^\((e)\)\  /. \ 
                      D \[Rule] 4 - 2  e // \ subInt, {e, 0, 0}]] // 
            Simplify\) // PowerExpand\) // Expand\)\)\)], "Input",
  PageWidth->PaperWidth],

Cell[BoxData[
    \(TraditionalForm\`1\/2\ C\_A\ c\[CapitalGamma]\ gs2\ \[Pi]\^2 - \(C\_A\ \
c\[CapitalGamma]\ gs2\)\/e\^2\)], "Output"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
    \(\(add = 
        1/2 \((Normal[
                Series[\((\((1 + \ \(19/
                                4\)/\[Pi]\ gs2/\((4  \[Pi])\)\ )\)/\((1 + 
                          2\ \(gs2/\((4  \[Pi])\)\)/\[Pi])\))\), {gs2, 0, 
                    1}]] - 1)\) c\[CapitalGamma]\ CA/3\ 16 
           Pi^2\ ;\)\), "\[IndentingNewLine]", 
    \(\(virt = add + res;\)\[IndentingNewLine]\), "\[IndentingNewLine]", 
    \(add\)}], "Input",
  PageWidth->PaperWidth],

Cell[BoxData[
    \(TraditionalForm\`11\/6\ C\_A\ c\[CapitalGamma]\ gs2\)], "Output"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
    \(virt = \(\((\((virt*\(2/gs2\)/c\[CapitalGamma])\) // 
            Expand)\)\(*\)\(as\/\(2  \[Pi]\)\) \(c\[CapitalGamma]\)\(\ \
\)\)\)], "Input"],

Cell[BoxData[
    \(TraditionalForm\`\(as\ c\[CapitalGamma]\ \((\(-\(\(2\ C\_A\)\/e\^2\)\) \
+ \[Pi]\^2\ C\_A + \(11\ C\_A\)\/3)\)\)\/\(2\ \[Pi]\)\)], "Output"]
}, Open  ]],

Cell[BoxData[
    \(\[IndentingNewLine]\)], "Input"],

Cell[CellGroupData[{

Cell[BoxData[
    \(\(\(32/256\)/3\)/\((1 - e)\)^2\)], "Input"],

Cell[BoxData[
    \(TraditionalForm\`1\/\(24\ \((1 - e)\)\^2\)\)], "Output"]
}, Open  ]]
}, Open  ]]
}, Closed]],

Cell[CellGroupData[{

Cell["The real contributions", "Section",
  PageWidth->PaperWidth],

Cell[CellGroupData[{

Cell["Kinematics", "Subsection",
  PageWidth->PaperWidth],

Cell["\<\
-----I take all momenta outgoing  

     p1 + p2 + p3 + p4 = 0     
     
     p1^2=0
     p2^2=0
     p3^3=0
     p4^2=mh^2
     
-----invariants

     s  = (p1 + p2)^2 = (p3 + p4)^2=2 p1.p2=mh^2+2 p3.p4
     t  = (p1 + p3)^2 = (p2+  p4)^2 =2 p1.p3=mh^2+2 p2.p4
     u  = (p2 + p3)^2 = (p1 + p4)^2=2 p1.p4+mh^2=+2 p2.p3
    
     s + t + u = mh^2
    
     \[Sigma]+\[Tau]+\[Upsilon]=1
        
 -----scalar products    
     
     p1.p2=s/2
     p1.p3=(t)/2
     p1.p4= (u-mh^2)/2
     p2.p3= (u)/2
     p2.p4=(t-mh^2)/2
     p3.p4=(s-mh^2)/2
    
-----physical region for production   p1+p2=-p3-p4

      s>(mh)^2 ; 
      t<0; 
      u<0;
         
\
\>", "Text",
  PageWidth->PaperWidth],

Cell[BoxData[{
    \(\(ScalarProduct[p1, p1] = 0;\)\), "\[IndentingNewLine]", 
    \(\(ScalarProduct[p2, p2] = 0;\)\), "\[IndentingNewLine]", 
    \(\(ScalarProduct[p3, p3] = 0;\)\), "\[IndentingNewLine]", 
    \(\(ScalarProduct[q, p3] = u/2;\)\  (*q = 
        p2 + p3*) \), "\[IndentingNewLine]", 
    \(\(ScalarProduct[q, p2] = u/2;\)\  (*q = 
        p2 + p3*) \), "\[IndentingNewLine]", 
    \(\(ScalarProduct[q, q] = u;\)\  (*q = 
        p2 + p3*) \), "\[IndentingNewLine]", 
    \(\(ScalarProduct[p4, p4] = mh2;\)\), "\[IndentingNewLine]", 
    \(\(ScalarProduct[p1, p2] = s/2;\)\), "\[IndentingNewLine]", 
    \(\(ScalarProduct[p1, p3] = t/2;\)\), "\[IndentingNewLine]", 
    \(\(ScalarProduct[p1, p4] = \((u - mh2)\)/2;\)\), "\[IndentingNewLine]", 
    \(\(ScalarProduct[p2, p3] = u/2;\)\), "\[IndentingNewLine]", 
    \(\(ScalarProduct[p2, p4] = \ \((t - mh2)\)/
          2;\)\), "\[IndentingNewLine]", 
    \(\(ScalarProduct[p3, p4] = \((s - mh2)\)/
          2\ \ ;\)\), "\[IndentingNewLine]", 
    \(\(ScalarProduct[p3, e3] = 0\ ;\)\), "\[IndentingNewLine]", 
    \(\(ScalarProduct[p1, e1] = 0\ ;\)\), "\[IndentingNewLine]", 
    \(\(ScalarProduct[p2, e2] = 0\ ;\)\), "\[IndentingNewLine]", 
    \(\(s13 = t;\)\), "\[IndentingNewLine]", 
    \(\(s23 = u;\)\)}], "Input",
  PageWidth->PaperWidth]
}, Open  ]],

Cell[CellGroupData[{

Cell["Verteces and Propagators", "Subsection",
  PageWidth->PaperWidth],

Cell[TextData[{
  StyleBox["GGG is the kinematic part of the three-gluon vtx (momenta \
outgoing, clockwise ordering):\nVTX(ggg) = (-\[ImaginaryI] ",
    FontSize->14],
  Cell[BoxData[
      \(TraditionalForm\`g\_s\)],
    FontSize->14],
  StyleBox[" ) (\[ImaginaryI] ",
    FontSize->14],
  Cell[BoxData[
      \(TraditionalForm\`f\^abc\)],
    FontSize->14],
  StyleBox[") GGG\nVTX(qqg) = ( -\[ImaginaryI] ",
    FontSize->14],
  Cell[BoxData[
      \(TraditionalForm\`g\_s\)],
    FontSize->14],
  ")",
  StyleBox[" ",
    FontSize->14],
  Cell[BoxData[
      \(TraditionalForm\`\((T\^a)\)\_ij\)],
    FontSize->14],
  StyleBox[" ",
    FontSize->14],
  Cell[BoxData[
      \(TraditionalForm\`\[Gamma]\^\[Mu]\)],
    FontSize->14],
  StyleBox["\nGluon Propagator= ",
    FontSize->14],
  Cell[BoxData[
      FormBox[
        FractionBox[
          StyleBox[\(\(-\[ImaginaryI]\)\ g\^\[Mu]\[Nu]\),
            FontSize->16], \(p\^2\)], TraditionalForm]],
    FontSize->14],
  "\nQuark ",
  StyleBox["Propagator= ",
    FontSize->14],
  Cell[BoxData[
      FormBox[
        FractionBox[
          StyleBox[\(\(\[ImaginaryI]\)\(\ \)\),
            FontSize->16], \(p\&^\)], TraditionalForm]],
    FontSize->14]
}], "Text",
  PageWidth->PaperWidth],

Cell[BoxData[{
    \(\(GGGD[p1_, p2_, p3_, m1_, m2_, m3_] := 
        FourVector[p1 - p2, m3, \ Dimension\  \[Rule] D]\ MTD[m1, m2] + 
          FourVector[p2 - p3, m1, \ Dimension\  \[Rule] \ D]\ MTD[m2, m3] + 
          FourVector[p3 - p1, m2, \ Dimension\  \[Rule] \ D]\ MTD[m1, 
              m3];\)\), "\[IndentingNewLine]", 
    \(\(H[p1_, p2_, m1_, m2_] := 
        MTD[m1, m2]\ SPD[p1, p2]\  - \ 
          FourVector[p1, m2, \ Dimension \[Rule] D]*\[IndentingNewLine]\ 
            FourVector[p2, m1, \ 
              Dimension\  \[Rule] D];\)\), "\[IndentingNewLine]", 
    \(\(PropQuark = I;\)\), "\[IndentingNewLine]", 
    \(\(PropGluon\  = \(-I\);\)\), "\[IndentingNewLine]", 
    \(\(vtx = \(-I\)\ gs;\)\)}], "Input",
  PageWidth->PaperWidth]
}, Open  ]],

Cell[CellGroupData[{

Cell["Sum over the  four Feynman diagrams", "Subsection"],

Cell["\<\

uno= - gs (-f123) GGGD[p1,p2,-p1-p2,m1,m2,mu]  (-I MTD[mu,nu]/s) (I A) \
H[p3,p1+p2,m3,nu]//Contract;
tre= - gs (-f123) GGGD[p3,p1,-p3-p1,m3,m1,mu] (-I MTD[mu,nu]/t) (I A)  \
H[p2,p1+p3,m2,nu]//Contract;
qua= - gs (-f123) GGGD[p2,p3,-p3-p2,m2,m3,mu] (-I MTD[mu,nu]/u) (I A)  \
H[p1,p2+p3,m1,nu]//Contract;
due= -  A gs (-f123) GGGD[p1,p2,p3,m1,m2,m3];
res=uno+due+tre+qua//ExpandScalarProduct;\
\>", "Input"]
}, Open  ]],

Cell[CellGroupData[{

Cell["\<\
Now I have to square the amplitude.  The sum is performed over the physical \
polarizations  of the gluons:\
\>", "Subsection"],

Cell["\<\
res1=res;
res2=res /. m1->m1p /. m2->m2p /. m3-> m3p;
tot=res1*(-MTD[m1,m1p]+(FVD[p1,m1] FVD[p2,m1p]+FVD[p1,m1p] \
FVD[p2,m1])/SPD[p1,p2])//Contract;
tot=tot*(-MTD[m2,m2p]+(FVD[p2,m2] FVD[p3,m2p]+FVD[p2,m2p] \
FVD[p3,m2])/SPD[p3,p2])//Contract;
tot=tot*(-MTD[m3,m3p]+(FVD[p1,m3] FVD[p3,m3p]+FVD[p1,m3p] \
FVD[p3,m3])/SPD[p3,p1])//Contract;
tot=tot*res2//Contract//Expand;\
\>", "Input"],

Cell[CellGroupData[{

Cell["\<\
Amp2=((tot//Contract//Factor)/. A-> as/3/Pi/v /. gs^2-> 4 Pi as/. f123^2-> 24 \
//Simplify)//Factor\
\>", "Input"],

Cell[BoxData[
    \(TraditionalForm\`\(\(1\/\(3\ \[Pi]\ s\ t\ u\ v\^2\)\)\((32\ as\^3\ \((D\
\ s\^4 - 2\ s\^4 + 2\ D\ t\ s\^3 - 4\ t\ s\^3 + 2\ D\ u\ s\^3 - 4\ u\ s\^3 + 
            3\ D\ t\^2\ s\^2 - 6\ t\^2\ s\^2 + 3\ D\ u\^2\ s\^2 - 
            6\ u\^2\ s\^2 + 8\ D\ t\ u\ s\^2 - 20\ t\ u\ s\^2 + 
            2\ D\ t\^3\ s - 4\ t\^3\ s + 2\ D\ u\^3\ s - 4\ u\^3\ s + 
            8\ D\ t\ u\^2\ s - 20\ t\ u\^2\ s + 8\ D\ t\^2\ u\ s - 
            20\ t\^2\ u\ s + D\ t\^4 - 2\ t\^4 + D\ u\^4 - 2\ u\^4 + 
            2\ D\ t\ u\^3 - 4\ t\ u\^3 + 3\ D\ t\^2\ u\^2 - 6\ t\^2\ u\^2 + 
            2\ D\ t\^3\ u - 4\ t\^3\ u)\))\)\)\)], "Output"]
}, Open  ]],

Cell["\<\

Amp2=Amp2 /. D->4-2e//Simplify;\
\>", "Input"]
}, Open  ]],

Cell[CellGroupData[{

Cell["Real Amplitude Squared in D dimensions", "Subsection",
  PageWidth->PaperWidth],

Cell[BoxData[
    \(\(\(\[IndentingNewLine]\)\(\(Emme = 
        1/s \(\(\((\((mh2^4 + s^4 + t^4 + u^4)\)\ \((1 - 2\ e)\) + \ 
                    e/2\ \((mh2^2 + s^2 + t^2 + u^2)\)^2)\)/s\)/t\)/
            u;\)\[IndentingNewLine]
    \(RealD = Emme/\((1 - e)\)^2;\)\[IndentingNewLine]
    \(Real4 = RealD /. \ e \[Rule] 0;\)\)\)\)], "Input",
  PageWidth->PaperWidth]
}, Open  ]],

Cell[CellGroupData[{

Cell["Phase space in D dimensions", "Subsection",
  PageWidth->PaperWidth],

Cell[BoxData[{
    \(\[IndentingNewLine]\(PS = \(1\/\(8  \[Pi]\)\) \(\((\(\(4\)\(\ \)\(\[Pi]\
\)\(\ \)\)\/mh2)\)\^e\) 
          1\/Gamma[1 - e]\ \((mh2\/s)\)\^e\ \ \((1 - mh2\/s)\)\^\(1 - 2  e\)\ \
v\^\(-e\)\ \((omv)\)\^\(-e\);\)\), "\[IndentingNewLine]", 
    \(\(substu = {t \[Rule] \ \(-s\)\ \((1 - mh2\/s)\) \((omv)\), \ 
          u \[Rule] \ \(-s\)\ \((1 - mh2\/s)\) v\ , 
          s \[Rule] \ mh2/z};\)\), "\[IndentingNewLine]", 
    \(\(cGamma = \((1/16)\)/Pi^2\ mh2^\((\(-e\))\) \((4\ Pi)\)^e\ Gamma[
            1 + e]\ Gamma[1 - e]^2/Gamma[1 - 2\ e];\)\), "\n", 
    \(\(pgg = 
        2 \((z\ PlusDistribution[1/\((1 - z)\)] + \((1 - z)\)/z + 
              z \((1 - z)\) + 
              11/12\ DeltaFunction[1 - z])\);\)\), "\[IndentingNewLine]", 
    \(\(s0 = z;\)\)}], "Input",
  PageWidth->PaperWidth]
}, Closed]],

Cell[CellGroupData[{

Cell["CDR", "Subsection",
  PageWidth->PaperWidth],

Cell[CellGroupData[{

Cell[BoxData[{
    \(intando = 
      FullSimplify[
          FullSimplify[\(RealD\ PS\ c\[CapitalGamma]\)\/cGamma] \
//. \[InvisibleSpace]substu] /. \[InvisibleSpace]omv \[Rule] 
          1 - v; \), "\n", 
    \(intando = intando\/\((1 - z)\)\^\(\(-2\)\ e - 1\)\)}], "Input"],

Cell[BoxData[
    \(TraditionalForm\`\(-\(\((c\[CapitalGamma]\ \((1\/mh2)\)\^e\ mh2\^e\ \
\[Pi]\ \((1 - v)\)\^\(\(-e\) - 1\)\ v\^\(\(-e\) - 1\)\ z\^e\ \((e\ \((3\ \((1 \
- v)\)\^4\ \((z - 1)\)\^4 + 3\ v\^4\ \((z - 1)\)\^4 - 
                      2\ v\^2\ \((z\^2 + 1)\)\ \((z - 1)\)\^2 - 
                      2\ \((1 - v)\)\^2\ \((v\^2\ \((z - 1)\)\^2 + z\^2 + 
                            1)\)\ \((z - 1)\)\^2 + 3\ z\^4 - 2\ z\^2 + 3)\) - 
                2\ \((\((1 - v)\)\^4\ \((z - 1)\)\^4 + v\^4\ \((z - 1)\)\^4 + 
                      z\^4 + 1)\))\)\ \(\[CapitalGamma](
              1 - 2\ e)\))\)/\((\((e - 1)\)\^2\ \(\[CapitalGamma](1 - e)\)\^3\
\ \(\[CapitalGamma](e + 1)\))\)\)\)\)], "Output"]
}, Open  ]],

Cell[BoxData[{
    \(\(\[Sigma]r = 
        Integrate[intando, {v, 0, 1}, \ GenerateConditions \[Rule] False] // 
          PowerExpand;\)\), "\n", 
    \(\(\[Sigma]r = \(Normal[Series[\[Sigma]r, {e, 0, 1}]] // Factor\) // 
          FullSimplify;\)\), "\n", 
    \(\(\[Sigma]r\  = \[Sigma]r\ \ *\((\(\(-1\)\/\(2  e\)\) 
                DeltaFunction[1 - z] + PlusDistribution[1/\((1 - z)\)] - 
              2\ e\ PlusDistribution[Log[1 - z]/\((1 - z)\)])\);\)\), "\n", 
    \(\(\[Sigma]r = 
        Normal[Series[\[Sigma]r/\((1 + e + e^2)\), {e, 0, 0}]] // 
          Expand;\)\)}], "Input",
  PageWidth->PaperWidth]
}, Closed]]
}, Closed]],

Cell[CellGroupData[{

Cell["Results", "Section"],

Cell[CellGroupData[{

Cell["Virtual+Real+Counterterms  in CDR", "Subsection"],

Cell[BoxData[{
    \(\(\[Sigma]rn = \(\(\[Sigma]r/c\[CapitalGamma]\)/\[Pi]\)/2\ CA - 
          1/e\ 2\ Pgg\ z\ CA + 1/e\ 2\ pgg\ z\ CA;\)\), "\n", 
    \(\(\[Sigma]rn = \((\((\((Coefficient[\[Sigma]rn, DeltaFunction[1 - z]] /. 
                      z \[Rule] 1 // FullSimplify)\)*
                DeltaFunction[
                  1 - z])\) + \((\((\((Coefficient[\[Sigma]rn, 
                        PlusDistribution[1/\((1 - z)\)]] // 
                      FullSimplify)\))\) 
                PlusDistribution[
                  1/\((1 - z)\)])\) + \((\((\((Coefficient[\[Sigma]rn, 
                        PlusDistribution[Log[1 - z]/\((1 - z)\)]] // 
                      FullSimplify)\))\) 
                PlusDistribution[
                  Log[1 - z]/\((1 - z)\)])\) + \((\(\(\[Sigma]rn /. 
                    DeltaFunction[x_] \[Rule] 0\) /. 
                  PlusDistribution[Log[1 - z]/\((1 - z)\)] \[Rule] 0\) /. 
                PlusDistribution[1/\((1 - z)\)] \[Rule] 0)\))\);\)\), "\n", 
    \(\(\[Sigma]rn = \((\((\((Coefficient[\[Sigma]rn, 
                          1/e\ PlusDistribution[1/\((1 - z)\)]]/\((1 - z)\) // 
                      FullSimplify)\))\) 1/e)\) + \((\[Sigma]rn - 
                Coefficient[\[Sigma]rn, 1/e\ PlusDistribution[1/\((1 - z)\)]] 
                  1/e\ PlusDistribution[1/\((1 - z)\)])\) // 
          Expand;\)\)}], "Input",
  PageWidth->PaperWidth],

Cell[CellGroupData[{

Cell[BoxData[
    \(divD = 
      Coefficient[\[Sigma]rn, 
          DeltaFunction[1 - z]]*\(as/2\)/\[Pi]\ c\[CapitalGamma]\)], "Input"],

Cell[BoxData[
    \(TraditionalForm\`\(as\ c\[CapitalGamma]\ \((\(11\ C\_A\)\/\(3\ e\) + \
\(2\ C\_A\)\/e\^2 - \(\[Pi]\^2\ C\_A\)\/3)\)\)\/\(2\ \[Pi]\)\)], "Output"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
    \(UVC = \ \(\(1/
          e\)\(\ \)\(4\)\(*\)\(\(-11\)\/6\) \(CA\)\(\ \)\(as\/\(4  \[Pi]\)\) \
\(c\[CapitalGamma]\)\(\ \)\)\)], "Input"],

Cell[BoxData[
    \(TraditionalForm\`\(-\(\(11\ as\ C\_A\ c\[CapitalGamma]\)\/\(6\ e\ \[Pi]\
\)\)\)\)], "Output"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
    \(\((\(\((divD + UVC)\)/as\)/c\[CapitalGamma]\ 2\ \[Pi] // 
          Expand)\)\ as\ \(c\[CapitalGamma]/2\)/\[Pi]\)], "Input"],

Cell[BoxData[
    \(TraditionalForm\`\(as\ c\[CapitalGamma]\ \((\(2\ C\_A\)\/e\^2 - \(\[Pi]\
\^2\ C\_A\)\/3 + C\_A\/3)\)\)\/\(2\ \[Pi]\)\)], "Output"]
}, Open  ]],

Cell[BoxData[
    \(\(sally = 
        as/\[Pi]\ \((\(-11\)/2\ \((1 - z)\)^3 + 
              6\ \((1 + z^4 + \((1 - z)\)^4)\)\ PlusDistribution[
                  Log[1 - z]/\((1 - z)\)] - 
              6\ \((z^2\ PlusDistribution[1/\((1 - z)\)] + \((1 - z)\) + 
                    z^2 \((1 - z)\))\) Log[z])\);\)\)], "Input"],

Cell[BoxData[
    \(\((\(\(\((\[Sigma]rn*\(as/2\)/\[Pi] /. DeltaFunction[x_] \[Rule] 0)\) - 
                    sally /. CA \[Rule] 3\) /. Pgg \[Rule] 0 // Expand\) // 
            Simplify)\) /. 
        PlusDistribution[1/\((1 - z)\)] \[Rule] 1/\((1 - z)\) // 
      Factor\)], "Input"]
}, Open  ]]
}, Closed]],

Cell[CellGroupData[{

Cell["Gluon initiated process", "Section",
  PageWidth->PaperWidth],

Cell[CellGroupData[{

Cell["Kinematics for the Born", "Subsection",
  PageWidth->PaperWidth],

Cell["\<\
-----I take all momenta outgoing  

     p1 + p2 + p3  = 0     
     
     p1^2=0
     p2^2=0
     p3^3=Q^2
     

        \
\>", "Text",
  PageWidth->PaperWidth],

Cell[BoxData[{
    \(\(ScalarProduct[p1, p1] = 0;\)\), "\[IndentingNewLine]", 
    \(\(ScalarProduct[p2, p2] = 0;\)\), "\[IndentingNewLine]", 
    \(\(ScalarProduct[p3, p3] = Q2;\)\), "\[IndentingNewLine]", 
    \(\(ScalarProduct[p1, p3] = \(-Q2\)/2;\)\), "\[IndentingNewLine]", 
    \(\(ScalarProduct[p1, p2] = Q2/2;\)\), "\[IndentingNewLine]", 
    \(\(ScalarProduct[p2, p3] = \(-\ Q2\)/2;\)\), "\[IndentingNewLine]", 
    \(\(ScalarProduct[p1, e1] = 0;\)\), "\[IndentingNewLine]", 
    \(\(ScalarProduct[p2, e2] = 0;\)\), "\[IndentingNewLine]", 
    \(\(ScalarProduct[p1, e2] = 0;\)\), "\[IndentingNewLine]", 
    \(\(\(ScalarProduct[p2, e1] = 0;\)\(\[IndentingNewLine]\)
    \)\), "\[IndentingNewLine]", 
    \(\(ScalarProduct[p1, p1, Dimension \[Rule] D] = 
        0;\)\), "\[IndentingNewLine]", 
    \(\(ScalarProduct[p2, p2, Dimension \[Rule] D] = 
        0;\)\), "\[IndentingNewLine]", 
    \(\(ScalarProduct[p3, p3, Dimension \[Rule] D] = 
        Q2;\)\), "\[IndentingNewLine]", 
    \(\(ScalarProduct[p1, p3, Dimension \[Rule] D] = \(-Q2\)/
          2;\)\), "\[IndentingNewLine]", 
    \(\(ScalarProduct[p1, p2, Dimension \[Rule] D] = 
        Q2/2;\)\), "\[IndentingNewLine]", 
    \(\(ScalarProduct[p2, p3, Dimension \[Rule] D] = \(-\ Q2\)/
          2;\)\), "\[IndentingNewLine]", 
    \(\(ScalarProduct[p1, e1, Dimension \[Rule] D] = 
        0;\)\), "\[IndentingNewLine]", 
    \(\(ScalarProduct[p2, e2, Dimension \[Rule] D] = 
        0;\)\), "\[IndentingNewLine]", 
    \(\(ScalarProduct[p1, e2, Dimension \[Rule] D] = 
        0;\)\), "\[IndentingNewLine]", 
    \(\(ScalarProduct[p2, e1, Dimension \[Rule] D] = 
        0;\)\), "\[IndentingNewLine]", 
    \(\)}], "Input",
  PageWidth->PaperWidth]
}, Open  ]],

Cell[CellGroupData[{

Cell["Real Amplitude", "Subsection",
  PageWidth->PaperWidth],

Cell[BoxData[{
    \(\(Emme = 
        CF\ 1/s \((\((s^2 + u^2)\)\  - e\ \((s + u)\)^2)\)/
            t;\)\), "\[IndentingNewLine]", 
    \(\(RealD = Emme/\((1 - e)\);\)\), "\[IndentingNewLine]", 
    \(\(Real4 = RealD /. \ e \[Rule] 0;\)\)}], "Input"]
}, Open  ]],

Cell[CellGroupData[{

Cell["Integration over the phase space in D dimensions", "Subsection",
  PageWidth->PaperWidth],

Cell[BoxData[{
    \(\(PS = \(1\/\(8  \[Pi]\)\) \(\((\(\(4\)\(\ \)\(\[Pi]\)\(\ \
\)\)\/mh2)\)\^e\) 
          1\/Gamma[1 - e]\ \((mh2\/s)\)\^e\ \ \((1 - mh2\/s)\)\^\(1 - 2  e\)\ \
v\^\(-e\)\ \((1 - v)\)\^\(-e\);\)\), "\[IndentingNewLine]", 
    \(\(substu = {t \[Rule] \ \(-s\)\ \((1 - mh2\/s)\) \((1 - v)\), \ 
          u \[Rule] \ \(-s\)\ \((1 - mh2\/s)\) v\ , 
          s \[Rule] \ mh2/z};\)\), "\[IndentingNewLine]", 
    \(\(cGamma = \((1/16)\)/Pi^2\ mh2^\((\(-e\))\) \((4\ Pi)\)^e\ Gamma[
            1 + e]\ Gamma[1 - e]^2/Gamma[1 - 2\ e];\)\), "\n", 
    \(\(pgq = 
        CF \((\((1 - z)\)\^2 + 1)\)/z\  // 
          Factor;\)\), "\[IndentingNewLine]", 
    \(\)}], "Input",
  PageWidth->PaperWidth],

Cell[CellGroupData[{

Cell["dim-reg", "Subsubsection",
  PageWidth->PaperWidth],

Cell[CellGroupData[{

Cell[BoxData[{
    \(\(\[Sigma]r = 
        Integrate[\((RealD\ PS\ c\[CapitalGamma]/cGamma\  // 
                    Simplify)\) //. \ substu // Expand, {v, 0, 1}, \ 
            GenerateConditions \[Rule] False] // Simplify;\)\), "\n", 
    \(\(\[Sigma]r = 
        Normal[Series[\[Sigma]r/\((1 + \ e)\), {e, 0, 0}]] /. \ 
            gs^2\  \[Rule] \ as\ 4\ \[Pi]\  // FullSimplify;\)\), "\n", 
    \(\(\[Sigma]r = \(\(\(-\[Sigma]r\)/c\[CapitalGamma]\)/\[Pi]\)/2 - 
          1/e\ z\ \ \ Pgq\  + 1/e\ \ z\ \ pgq;\)\), "\n", 
    \(\[Sigma]r = \(\[Sigma]r // PowerExpand\) // Expand\)}], "Input",
  PageWidth->PaperWidth],

Cell[BoxData[
    \(TraditionalForm\`\(-\(1\/2\)\)\ C\_F\ z\^2 + 
      2\ C\_F\ \(log(1 - z)\)\ z\^2 - C\_F\ \(log(z)\)\ z\^2 + 
      3\ C\_F\ z - \(Pgq\ z\)\/e - 4\ C\_F\ \(log(1 - z)\)\ z + 
      2\ C\_F\ \(log(z)\)\ z - \(3\ C\_F\)\/2 + 4\ C\_F\ \(log(1 - z)\) - 
      2\ C\_F\ \(log(z)\)\)], "Output"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
    \(sally = 
      2*\((\(-\ \ \((1 - z)\)\)\ \((7 - 3  z)\)/3\  + 
            1/2\ z\ pgq\ \((1 + 2  Log[1 - z] - 
                  Log[z])\))\)\), "\[IndentingNewLine]", 
    \(\((\(\((\[Sigma]r - \ sally)\) /. \ CF -> 4/3\) /. \ Pgq \[Rule] 0)\) // 
      Simplify\)}], "Input"],

Cell[BoxData[
    \(TraditionalForm\`2\ \((1\/3\ \((7 - 3\ z)\)\ \((z - 1)\) + 
          1\/2\ C\_F\ \((z\^2 - 2\ z + 2)\)\ \((2\ \(log(1 - z)\) - log(z) + 
                1)\))\)\)], "Output"],

Cell[BoxData[
    \(TraditionalForm\`0\)], "Output"]
}, Open  ]],

Cell[BoxData[""], "Input"]
}, Open  ]]
}, Open  ]]
}, Closed]]
}, Open  ]]
},
FrontEndVersion->"4.1 for Macintosh",
ScreenRectangle->{{0, 1280}, {0, 832}},
WindowSize->{807, 655},
WindowMargins->{{Automatic, 126}, {Automatic, 0}},
PrintingCopies->1,
PrintingPageRange->{Automatic, Automatic}
]

(*******************************************************************
Cached data follows.  If you edit this Notebook file directly, not
using Mathematica, you must remove the line containing CacheID at
the top of  the file.  The cache data will then be recreated when
you save this file from within Mathematica.
*******************************************************************)

(*CellTagsOutline
CellTagsIndex->{}
*)

(*CellTagsIndex
CellTagsIndex->{}
*)

(*NotebookFileOutline
Notebook[{

Cell[CellGroupData[{
Cell[1727, 52, 72, 1, 108, "Title"],

Cell[CellGroupData[{
Cell[1824, 57, 61, 1, 46, "Subsection"],
Cell[1888, 60, 82, 2, 27, "Input"],
Cell[1973, 64, 369, 14, 32, "Text"],
Cell[2345, 80, 80, 2, 27, "Input"]
}, Closed]],

Cell[CellGroupData[{
Cell[2462, 87, 75, 1, 36, "Section"],

Cell[CellGroupData[{
Cell[2562, 92, 57, 1, 46, "Subsection"],
Cell[2622, 95, 172, 12, 204, "Text"],
Cell[2797, 109, 1720, 33, 450, "Input"]
}, Closed]],

Cell[CellGroupData[{
Cell[4554, 147, 71, 1, 30, "Subsection"],
Cell[4628, 150, 1246, 46, 152, "Text"],
Cell[5877, 198, 896, 16, 310, "Input"]
}, Closed]],

Cell[CellGroupData[{
Cell[6810, 219, 67, 1, 30, "Subsection"],

Cell[CellGroupData[{
Cell[6902, 224, 378, 8, 130, "Input"],
Cell[7283, 234, 54, 1, 70, "Output"],
Cell[7340, 237, 286, 7, 70, "Output"]
}, Open  ]],

Cell[CellGroupData[{
Cell[7663, 249, 488, 9, 110, "Input"],
Cell[8154, 260, 90, 1, 70, "Output"]
}, Open  ]]
}, Closed]],

Cell[CellGroupData[{
Cell[8293, 267, 66, 1, 30, "Subsection"],
Cell[8362, 270, 757, 16, 210, "Input"],
Cell[9122, 288, 26, 0, 30, "Input"],
Cell[9151, 290, 26, 0, 30, "Input"]
}, Closed]]
}, Closed]],

Cell[CellGroupData[{
Cell[9226, 296, 22, 0, 36, "Section"],

Cell[CellGroupData[{
Cell[9273, 300, 62, 1, 46, "Subsection"],

Cell[CellGroupData[{
Cell[9360, 305, 897, 17, 235, "Input"],
Cell[10260, 324, 394, 12, 26, "Output"],
Cell[10657, 338, 1327, 28, 29, "Output"]
}, Open  ]],
Cell[11999, 369, 26, 0, 27, "Input"]
}, Open  ]],

Cell[CellGroupData[{
Cell[12062, 374, 63, 1, 46, "Subsection"],

Cell[CellGroupData[{
Cell[12150, 379, 673, 13, 219, "Input"],
Cell[12826, 394, 251, 4, 48, "Output"]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{
Cell[13126, 404, 38, 0, 46, "Subsection"],

Cell[CellGroupData[{
Cell[13189, 408, 389, 7, 107, "Input"],
Cell[13581, 417, 136, 2, 43, "Output"]
}, Open  ]],

Cell[CellGroupData[{
Cell[13754, 424, 474, 10, 155, "Input"],
Cell[14231, 436, 85, 1, 43, "Output"]
}, Open  ]],

Cell[CellGroupData[{
Cell[14353, 442, 164, 3, 40, "Input"],
Cell[14520, 447, 160, 2, 52, "Output"]
}, Open  ]],
Cell[14695, 452, 52, 1, 43, "Input"],

Cell[CellGroupData[{
Cell[14772, 457, 63, 1, 27, "Input"],
Cell[14838, 460, 76, 1, 45, "Output"]
}, Open  ]]
}, Open  ]]
}, Closed]],

Cell[CellGroupData[{
Cell[14975, 468, 66, 1, 36, "Section"],

Cell[CellGroupData[{
Cell[15066, 473, 57, 1, 46, "Subsection"],
Cell[15126, 476, 702, 37, 574, "Text"],
Cell[15831, 515, 1310, 24, 299, "Input"]
}, Open  ]],

Cell[CellGroupData[{
Cell[17178, 544, 71, 1, 46, "Subsection"],
Cell[17252, 547, 1246, 46, 152, "Text"],
Cell[18501, 595, 757, 14, 187, "Input"]
}, Open  ]],

Cell[CellGroupData[{
Cell[19295, 614, 57, 0, 46, "Subsection"],
Cell[19355, 616, 419, 10, 147, "Input"]
}, Open  ]],

Cell[CellGroupData[{
Cell[19811, 631, 137, 3, 46, "Subsection"],
Cell[19951, 636, 396, 10, 147, "Input"],

Cell[CellGroupData[{
Cell[20372, 650, 124, 3, 42, "Input"],
Cell[20499, 655, 649, 9, 98, "Output"]
}, Open  ]],
Cell[21163, 667, 57, 3, 42, "Input"]
}, Open  ]],

Cell[CellGroupData[{
Cell[21257, 675, 85, 1, 46, "Subsection"],
Cell[21345, 678, 368, 7, 123, "Input"]
}, Open  ]],

Cell[CellGroupData[{
Cell[21750, 690, 74, 1, 46, "Subsection"],
Cell[21827, 693, 820, 15, 231, "Input"]
}, Closed]],

Cell[CellGroupData[{
Cell[22684, 713, 50, 1, 30, "Subsection"],

Cell[CellGroupData[{
Cell[22759, 718, 277, 6, 98, "Input"],
Cell[23039, 726, 708, 10, 121, "Output"]
}, Open  ]],
Cell[23762, 739, 618, 12, 171, "Input"]
}, Closed]]
}, Closed]],

Cell[CellGroupData[{
Cell[24429, 757, 26, 0, 36, "Section"],

Cell[CellGroupData[{
Cell[24480, 761, 55, 0, 46, "Subsection"],
Cell[24538, 763, 1409, 24, 299, "Input"],

Cell[CellGroupData[{
Cell[25972, 791, 136, 3, 27, "Input"],
Cell[26111, 796, 165, 2, 55, "Output"]
}, Open  ]],

Cell[CellGroupData[{
Cell[26313, 803, 155, 3, 42, "Input"],
Cell[26471, 808, 113, 2, 43, "Output"]
}, Open  ]],

Cell[CellGroupData[{
Cell[26621, 815, 144, 2, 27, "Input"],
Cell[26768, 819, 150, 2, 55, "Output"]
}, Open  ]],
Cell[26933, 824, 329, 6, 59, "Input"],
Cell[27265, 832, 289, 5, 43, "Input"]
}, Open  ]]
}, Closed]],

Cell[CellGroupData[{
Cell[27603, 843, 67, 1, 36, "Section"],

Cell[CellGroupData[{
Cell[27695, 848, 70, 1, 46, "Subsection"],
Cell[27768, 851, 172, 12, 174, "Text"],
Cell[27943, 865, 1720, 33, 363, "Input"]
}, Open  ]],

Cell[CellGroupData[{
Cell[29700, 903, 61, 1, 46, "Subsection"],
Cell[29764, 906, 253, 5, 59, "Input"]
}, Open  ]],

Cell[CellGroupData[{
Cell[30054, 916, 95, 1, 46, "Subsection"],
Cell[30152, 919, 712, 14, 186, "Input"],

Cell[CellGroupData[{
Cell[30889, 937, 57, 1, 42, "Subsubsection"],

Cell[CellGroupData[{
Cell[30971, 942, 623, 11, 139, "Input"],
Cell[31597, 955, 309, 5, 43, "Output"]
}, Open  ]],

Cell[CellGroupData[{
Cell[31943, 965, 300, 6, 43, "Input"],
Cell[32246, 973, 195, 3, 43, "Output"],
Cell[32444, 978, 52, 1, 26, "Output"]
}, Open  ]],
Cell[32511, 982, 26, 0, 27, "Input"]
}, Open  ]]
}, Open  ]]
}, Closed]]
}, Open  ]]
}
]
*)



(*******************************************************************
End of Mathematica Notebook file.
*******************************************************************)

