(* Content-type: application/mathematica *)

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

(* CreatedBy='Mathematica 6.0' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       145,          7]
NotebookDataLength[     50813,       1545]
NotebookOptionsPosition[     45605,       1361]
NotebookOutlinePosition[     45963,       1377]
CellTagsIndexPosition[     45920,       1374]
WindowFrame->Normal
ContainsDynamic->False*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell["\<\
Tree-level decays of the Higgs boson\
\>", "Title",
 CellChangeTimes->{{3.42892672627883*^9, 3.428926824323577*^9}, {
  3.428937703319875*^9, 3.428937712432496*^9}}],

Cell[CellGroupData[{

Cell["Input FeynCalc", "Subsection"],

Cell[BoxData[
 RowBox[{
  RowBox[{"<<", "HighEnergyPhysics`fc`"}], ";"}]], "Input"]
}, Open  ]],

Cell[CellGroupData[{

Cell["H \[RightArrow] ff", "Section",
 CellChangeTimes->{{3.4289377207239313`*^9, 3.428937731049849*^9}}],

Cell[CellGroupData[{

Cell["Kinematics 1->2", "Subsection",
 CellChangeTimes->{{3.428926835790217*^9, 3.428926837886043*^9}}],

Cell[BoxData[
 RowBox[{"\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{
     RowBox[{"ScalarProduct", "[", 
      RowBox[{"P", ",", "P"}], "]"}], "=", 
     RowBox[{"mh", "^", "2"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{
     RowBox[{"ScalarProduct", "[", 
      RowBox[{"p1", ",", "p1"}], "]"}], "=", 
     RowBox[{"mf", "^", "2"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{
     RowBox[{"ScalarProduct", "[", 
      RowBox[{"p2", ",", "p2"}], "]"}], "=", 
     RowBox[{"mf", "^", "2"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{
     RowBox[{"ScalarProduct", "[", 
      RowBox[{"P", ",", "p1"}], "]"}], "=", 
     RowBox[{
      RowBox[{"mh", "^", "2"}], "/", "2"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{
     RowBox[{"ScalarProduct", "[", 
      RowBox[{"P", ",", "p2"}], "]"}], "=", 
     RowBox[{
      RowBox[{"mh", "^", "2"}], "/", "2"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{
     RowBox[{"ScalarProduct", "[", 
      RowBox[{"p1", ",", "p2"}], "]"}], "=", 
     RowBox[{
      RowBox[{
       RowBox[{"mh", "^", "2"}], "/", "2"}], "-", 
      RowBox[{"mf", "^", "2"}]}]}], ";"}], "\[IndentingNewLine]"}]}]], "Input",\

 CellChangeTimes->{{3.4289268403011503`*^9, 3.428926946965074*^9}, {
  3.428937810100576*^9, 3.428937873020487*^9}, {3.428938071207912*^9, 
  3.428938078071608*^9}, {3.428938213984881*^9, 3.428938224130307*^9}, {
  3.4289384819556913`*^9, 3.428938484248788*^9}}]
}, Closed]],

Cell[CellGroupData[{

Cell["Amplitude (1 diagram)", "Subsection",
 CellChangeTimes->{{3.4289292421393127`*^9, 3.428929244675111*^9}}],

Cell["The summed and averaged square matrix element", "Text",
 CellChangeTimes->{{3.4289271420095243`*^9, 3.4289271454722757`*^9}, {
  3.4289273354903717`*^9, 3.428927351049911*^9}, {3.4289275105378*^9, 
  3.42892751343362*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"M2", "=", 
  RowBox[{
   RowBox[{
    RowBox[{
     RowBox[{
      RowBox[{"(", 
       RowBox[{"mf", "/", "v"}], ")"}], "^", "2"}], 
     RowBox[{"Tr", "[", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{
         RowBox[{"GSD", "[", "p1", "]"}], "+", "mf"}], ")"}], ".", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"GSD", "[", "p2", "]"}], "-", "mf"}], ")"}]}], "]"}]}], "//",
     "Expand"}], "//", "Simplify"}]}]], "Input",
 CellChangeTimes->{{3.428927033246361*^9, 3.4289271401967373`*^9}, {
   3.42892727787171*^9, 3.428927413135103*^9}, {3.4289274432481937`*^9, 
   3.4289274930463*^9}, {3.428929121395401*^9, 3.428929122088388*^9}, 
   3.428929375657975*^9, {3.428929541763723*^9, 3.428929566521463*^9}, {
   3.428937860468915*^9, 3.428937864548534*^9}, {3.4289380364317427`*^9, 
   3.428938062823503*^9}, {3.42893809477054*^9, 3.4289381246774197`*^9}, {
   3.4289382285018063`*^9, 3.428938230286139*^9}}],

Cell[BoxData[
 FormBox[
  FractionBox[
   RowBox[{
    RowBox[{"2", " ", 
     SuperscriptBox["mf", "2"], " ", 
     SuperscriptBox["mh", "2"]}], "-", 
    RowBox[{"8", " ", 
     SuperscriptBox["mf", "4"]}]}], 
   SuperscriptBox["v", "2"]], TraditionalForm]], "Output",
 CellChangeTimes->{
  3.4289272789512367`*^9, 3.428927318445057*^9, {3.4289273534290743`*^9, 
   3.428927404797317*^9}, {3.4289274516566753`*^9, 3.42892749447148*^9}, 
   3.4289291227183*^9, {3.428929370895793*^9, 3.428929376751039*^9}, 
   3.428929567605669*^9, {3.428938112386113*^9, 3.428938125511932*^9}, 
   3.4289382307430983`*^9, {3.428938486376433*^9, 3.4289384975076303`*^9}, 
   3.428939015665769*^9}]
}, Open  ]],

Cell["We insert he definition of the Fermi constant", "Text",
 CellChangeTimes->{{3.428927552171399*^9, 3.428927561923086*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"M2", "=", 
  RowBox[{"M2", "/.", 
   RowBox[{"v", "\[Rule]", 
    RowBox[{"1", "/", 
     RowBox[{"Sqrt", "[", 
      RowBox[{"Gf", " ", 
       RowBox[{"Sqrt", "[", "2", "]"}]}], "]"}]}]}]}]}]], "Input",
 CellChangeTimes->{{3.428927568404044*^9, 3.42892758705121*^9}, {
  3.4289381595073957`*^9, 3.4289381716760798`*^9}}],

Cell[BoxData[
 FormBox[
  RowBox[{
   SqrtBox["2"], " ", "Gf", " ", 
   RowBox[{"(", 
    RowBox[{
     RowBox[{"2", " ", 
      SuperscriptBox["mf", "2"], " ", 
      SuperscriptBox["mh", "2"]}], "-", 
     RowBox[{"8", " ", 
      SuperscriptBox["mf", "4"]}]}], ")"}]}], TraditionalForm]], "Output",
 CellChangeTimes->{3.428927588058874*^9, 3.4289291256255703`*^9, 
  3.428929571106011*^9, 3.428938172246612*^9, 3.428938233053894*^9, 
  3.4289384978556747`*^9, 3.428939015694807*^9}]
}, Open  ]]
}, Closed]],

Cell[CellGroupData[{

Cell["Phase Space", "Subsection",
 CellChangeTimes->{{3.4289292202493353`*^9, 3.428929221917066*^9}}],

Cell["The two particle phase space is", "Text",
 CellChangeTimes->{{3.4289286541859837`*^9, 3.428928662992796*^9}, {
  3.428938253498197*^9, 3.428938253866067*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"d\[CapitalPhi]2", "=", 
  RowBox[{
   RowBox[{
    RowBox[{"(", 
     RowBox[{"2", "Pi"}], ")"}], "^", "4"}], 
   RowBox[{"\[Delta]", "[", 
    RowBox[{"P", "-", "p1", "-", "p2"}], "]"}], 
   RowBox[{
    RowBox[{
     RowBox[{"d3p1", "/", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{"2", "Pi"}], ")"}], "^", "3"}]}], "/", "2"}], "/", "E1"}], 
   " ", 
   RowBox[{
    RowBox[{
     RowBox[{"d3p2", "/", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{"2", "Pi"}], ")"}], "^", "3"}]}], "/", "2"}], "/", 
    "E2"}]}]}]], "Input",
 CellChangeTimes->{{3.428928668808182*^9, 3.428928748243292*^9}, {
  3.4289288594542103`*^9, 3.428928869884198*^9}, {3.42893818360579*^9, 
  3.428938207781085*^9}, {3.4289382412163973`*^9, 3.4289382417010736`*^9}}],

Cell[BoxData[
 FormBox[
  FractionBox[
   RowBox[{"d3p1", " ", "d3p2", " ", 
    RowBox[{"\[Delta]", "(", 
     RowBox[{"P", "-", "p1", "-", "p2"}], ")"}]}], 
   RowBox[{"16", " ", "E1", " ", "E2", " ", 
    SuperscriptBox["\[Pi]", "2"]}]], TraditionalForm]], "Output",
 CellChangeTimes->{{3.428928862716354*^9, 3.428928870072042*^9}, 
   3.428929126861431*^9, 3.428929574710586*^9, {3.428938234752118*^9, 
   3.428938241930442*^9}, 3.428938501869619*^9, 3.4289390157281837`*^9}]
}, Open  ]],

Cell["\<\
Performing the integration over d3p2 in the rest frame of the H gives\
\>", "Text",
 CellChangeTimes->{{3.428928749821361*^9, 3.428928784904477*^9}, {
  3.4289288233373938`*^9, 3.428928824472897*^9}, {3.428938246362335*^9, 
  3.428938248774928*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"PS", "=", 
  RowBox[{
   RowBox[{
    RowBox[{"1", "/", "4"}], "/", 
    RowBox[{
     RowBox[{"(", 
      RowBox[{"2", "Pi"}], ")"}], "^", "2"}]}], " ", 
   RowBox[{"\[Delta]", "[", 
    RowBox[{"mh", "-", 
     RowBox[{"2", "E1"}]}], "]"}], 
   RowBox[{"d3p1", "/", 
    RowBox[{"E1", "^", "2"}]}]}]}]], "Input",
 CellChangeTimes->{{3.4289287922794847`*^9, 3.4289288819309*^9}, {
  3.42893825962851*^9, 3.42893826484441*^9}}],

Cell[BoxData[
 FormBox[
  FractionBox[
   RowBox[{"d3p1", " ", 
    RowBox[{"\[Delta]", "(", 
     RowBox[{"mh", "-", 
      RowBox[{"2", " ", "E1"}]}], ")"}]}], 
   RowBox[{"16", " ", 
    SuperscriptBox["E1", "2"], " ", 
    SuperscriptBox["\[Pi]", "2"]}]], TraditionalForm]], "Output",
 CellChangeTimes->{{3.428928849656254*^9, 3.428928882523725*^9}, 
   3.428929068231016*^9, 3.428929127845537*^9, 3.428929574829381*^9, 
   3.42893826621164*^9, 3.428938504167341*^9, 3.4289386058116293`*^9, 
   3.4289386914875193`*^9, 3.4289390157737207`*^9}]
}, Open  ]],

Cell["\<\
Going to spherical coordinates, and usging p1 dp1 = E1 dE1, we get\
\>", "Text",
 CellChangeTimes->{{3.428928883567082*^9, 3.4289289264749002`*^9}, {
  3.4289382922025423`*^9, 3.428938296218473*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"PS", "=", 
  RowBox[{
   RowBox[{"PS", "/.", 
    RowBox[{
     RowBox[{"\[Delta]", "[", 
      RowBox[{"mh", "-", 
       RowBox[{"2", "E1"}]}], "]"}], "\[Rule]", 
     RowBox[{
      RowBox[{"1", "/", "2"}], 
      RowBox[{"\[Delta]", "[", 
       RowBox[{"E1", "-", 
        RowBox[{"mh", "/", "2"}]}], "]"}]}]}]}], "/.", " ", 
   RowBox[{"d3p1", "\[Rule]", 
    RowBox[{"4", "Pi", " ", 
     RowBox[{"Sqrt", "[", 
      RowBox[{
       RowBox[{"E1", "^", "2"}], "-", 
       RowBox[{"mf", "^", "2"}]}], "]"}], "E1", " ", "dE1"}]}]}]}]], "Input",
 CellChangeTimes->{{3.4289289441444197`*^9, 3.428929022421103*^9}, 
   3.428929070151713*^9, {3.4289382799967327`*^9, 3.428938290196628*^9}, {
   3.428938584824209*^9, 3.428938599062233*^9}}],

Cell[BoxData[
 FormBox[
  FractionBox[
   RowBox[{"dE1", " ", 
    SqrtBox[
     RowBox[{
      SuperscriptBox["E1", "2"], "-", 
      SuperscriptBox["mf", "2"]}]], " ", 
    RowBox[{"\[Delta]", "(", 
     RowBox[{"E1", "-", 
      FractionBox["mh", "2"]}], ")"}]}], 
   RowBox[{"8", " ", "E1", " ", "\[Pi]"}]], TraditionalForm]], "Output",
 CellChangeTimes->{
  3.428929025150517*^9, 3.428929070325218*^9, 3.4289291311060667`*^9, 
   3.428929574875102*^9, 3.428938297415737*^9, 3.428938505065439*^9, {
   3.428938600499318*^9, 3.428938606660725*^9}, {3.4289386861864023`*^9, 
   3.428938692218741*^9}, 3.4289390158068933`*^9}]
}, Open  ]],

Cell["Performing the integration", "Text",
 CellChangeTimes->{{3.428929031938621*^9, 3.428929037282918*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"PS", "=", 
  RowBox[{
   RowBox[{"PS", "/.", 
    RowBox[{
     RowBox[{"dE1", " ", 
      RowBox[{"\[Delta]", "[", 
       RowBox[{"E1", "-", 
        RowBox[{"mh", "/", "2"}]}], "]"}]}], "\[Rule]", "1"}]}], "/.", 
   RowBox[{"E1", "\[Rule]", 
    RowBox[{"mh", "/", "2"}]}]}]}]], "Input",
 CellChangeTimes->{{3.428929044586605*^9, 3.428929056603554*^9}, {
  3.428938304472186*^9, 3.428938307564096*^9}, {3.428938621479186*^9, 
  3.428938627852805*^9}, {3.4289386823375387`*^9, 3.428938682419998*^9}}],

Cell[BoxData[
 FormBox[
  FractionBox[
   SqrtBox[
    RowBox[{
     FractionBox[
      SuperscriptBox["mh", "2"], "4"], "-", 
     SuperscriptBox["mf", "2"]}]], 
   RowBox[{"4", " ", "mh", " ", "\[Pi]"}]], TraditionalForm]], "Output",
 CellChangeTimes->{{3.428929056861126*^9, 3.428929071114458*^9}, 
   3.4289291321148567`*^9, 3.428929574912538*^9, 3.428938307856097*^9, 
   3.428938506325878*^9, 3.4289386297502537`*^9, {3.4289386829460487`*^9, 
   3.428938694746279*^9}, 3.42893901583956*^9}]
}, Open  ]]
}, Closed]],

Cell[CellGroupData[{

Cell["Result", "Subsection",
 CellChangeTimes->{{3.4289292627868023`*^9, 3.4289292633475027`*^9}}],

Cell["The decay rate of the H then becomes", "Text",
 CellChangeTimes->{{3.428929073874954*^9, 3.428929081291203*^9}, 
   3.4289383237898083`*^9}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"\[CapitalGamma]", "=", 
  RowBox[{
   RowBox[{
    RowBox[{
     RowBox[{"1", "/", "2"}], "/", "mh"}], " ", "PS", " ", "M2"}], "//", 
   "Simplify"}]}]], "Input",
 CellChangeTimes->{{3.428929087038114*^9, 3.428929107627933*^9}, 
   3.428929346599481*^9, {3.428938314727579*^9, 3.428938334043777*^9}, 
   3.42893864433534*^9, 3.4289389462936773`*^9}],

Cell[BoxData[
 FormBox[
  FractionBox[
   RowBox[{"Gf", " ", 
    SqrtBox[
     RowBox[{
      SuperscriptBox["mh", "2"], "-", 
      RowBox[{"4", " ", 
       SuperscriptBox["mf", "2"]}]}]], " ", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{
       SuperscriptBox["mf", "2"], " ", 
       SuperscriptBox["mh", "2"]}], "-", 
      RowBox[{"4", " ", 
       SuperscriptBox["mf", "4"]}]}], ")"}]}], 
   RowBox[{"4", " ", 
    SqrtBox["2"], " ", 
    SuperscriptBox["mh", "2"], " ", "\[Pi]"}]], TraditionalForm]], "Output",
 CellChangeTimes->{{3.428929108658009*^9, 3.428929133003119*^9}, 
   3.4289293476527367`*^9, 3.4289295771229353`*^9, {3.428938319056734*^9, 
   3.4289383343919973`*^9}, 3.4289385085419197`*^9, {3.4289386319695272`*^9, 
   3.428938644554867*^9}, 3.428938696988208*^9, 3.428938947199583*^9, 
   3.42893901589124*^9}]
}, Open  ]]
}, Closed]]
}, Open  ]],

Cell[CellGroupData[{

Cell["H \[RightArrow] WW", "Section",
 CellChangeTimes->{{3.4289377207239313`*^9, 3.428937731049849*^9}, {
  3.428939028678343*^9, 3.428939029987813*^9}}],

Cell[CellGroupData[{

Cell["Kinematics 1->2", "Subsection",
 CellChangeTimes->{{3.428926835790217*^9, 3.428926837886043*^9}}],

Cell[BoxData[
 RowBox[{"\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{
     RowBox[{"ScalarProduct", "[", 
      RowBox[{"P", ",", "P"}], "]"}], "=", 
     RowBox[{"mh", "^", "2"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{
     RowBox[{"ScalarProduct", "[", 
      RowBox[{"pw1", ",", "pw1"}], "]"}], "=", 
     RowBox[{"mw", "^", "2"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{
     RowBox[{"ScalarProduct", "[", 
      RowBox[{"pw2", ",", "pw2"}], "]"}], "=", 
     RowBox[{"mw", "^", "2"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{
     RowBox[{"ScalarProduct", "[", 
      RowBox[{"P", ",", "pw1"}], "]"}], "=", 
     RowBox[{
      RowBox[{"mh", "^", "2"}], "/", "2"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{
     RowBox[{"ScalarProduct", "[", 
      RowBox[{"P", ",", "pw2"}], "]"}], "=", 
     RowBox[{
      RowBox[{"mh", "^", "2"}], "/", "2"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{
     RowBox[{"ScalarProduct", "[", 
      RowBox[{"pw1", ",", "pw2"}], "]"}], "=", 
     RowBox[{
      RowBox[{
       RowBox[{"mh", "^", "2"}], "/", "2"}], "-", 
      RowBox[{"mw", "^", "2"}]}]}], ";"}], "\[IndentingNewLine]"}]}]], "Input",\

 CellChangeTimes->{{3.4289268403011503`*^9, 3.428926946965074*^9}, {
  3.428937810100576*^9, 3.428937873020487*^9}, {3.428938071207912*^9, 
  3.428938078071608*^9}, {3.428938213984881*^9, 3.428938224130307*^9}, {
  3.4289384819556913`*^9, 3.428938484248788*^9}, {3.428939073598165*^9, 
  3.428939082348546*^9}, {3.42893937775745*^9, 3.4289393838581676`*^9}}]
}, Closed]],

Cell[CellGroupData[{

Cell["Amplitude (1 diagram)", "Subsection",
 CellChangeTimes->{{3.4289292421393127`*^9, 3.428929244675111*^9}}],

Cell["The summed and averaged square matrix element", "Text",
 CellChangeTimes->{{3.4289271420095243`*^9, 3.4289271454722757`*^9}, {
  3.4289273354903717`*^9, 3.428927351049911*^9}, {3.4289275105378*^9, 
  3.42892751343362*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"M2", "=", 
  RowBox[{
   RowBox[{
    RowBox[{"(", 
     RowBox[{"v", " ", 
      RowBox[{
       RowBox[{"gw", "^", "2"}], "/", "2"}]}], ")"}], "^", "2"}], 
   RowBox[{"(", 
    RowBox[{
     RowBox[{"-", 
      RowBox[{"MetricTensor", "[", 
       RowBox[{"mu", ",", "nu"}], "]"}]}], "+", 
     RowBox[{
      RowBox[{"FV", "[", 
       RowBox[{"pw1", ",", "mu"}], "]"}], 
      RowBox[{
       RowBox[{"FV", "[", 
        RowBox[{"pw1", ",", "nu"}], "]"}], "/", 
       RowBox[{"mw", "^", "2"}]}]}]}], ")"}], 
   RowBox[{"(", 
    RowBox[{
     RowBox[{"-", 
      RowBox[{"MetricTensor", "[", 
       RowBox[{"mu", ",", "nu"}], "]"}]}], "+", 
     RowBox[{
      RowBox[{"FV", "[", 
       RowBox[{"pw2", ",", "mu"}], "]"}], 
      RowBox[{
       RowBox[{"FV", "[", 
        RowBox[{"pw2", ",", "nu"}], "]"}], "/", 
       RowBox[{"mw", "^", "2"}]}]}]}], ")"}]}]}]], "Input",
 CellChangeTimes->{{3.428927033246361*^9, 3.4289271401967373`*^9}, {
   3.42892727787171*^9, 3.428927413135103*^9}, {3.4289274432481937`*^9, 
   3.4289274930463*^9}, {3.428929121395401*^9, 3.428929122088388*^9}, 
   3.428929375657975*^9, {3.428929541763723*^9, 3.428929566521463*^9}, {
   3.428937860468915*^9, 3.428937864548534*^9}, {3.4289380364317427`*^9, 
   3.428938062823503*^9}, {3.42893809477054*^9, 3.4289381246774197`*^9}, {
   3.4289382285018063`*^9, 3.428938230286139*^9}, {3.4289390981304398`*^9, 
   3.428939148746867*^9}, {3.428939388292745*^9, 3.428939394515564*^9}, 
   3.428939614471209*^9}],

Cell[BoxData[
 FormBox[
  RowBox[{
   FractionBox["1", "4"], " ", 
   SuperscriptBox["gw", "4"], " ", 
   SuperscriptBox["v", "2"], " ", 
   RowBox[{"(", 
    RowBox[{
     FractionBox[
      RowBox[{
       SuperscriptBox["pw1", "mu"], " ", 
       SuperscriptBox["pw1", "nu"]}], 
      SuperscriptBox["mw", "2"]], "-", 
     SuperscriptBox["g", 
      RowBox[{"mu", "nu"}]]}], ")"}], " ", 
   RowBox[{"(", 
    RowBox[{
     FractionBox[
      RowBox[{
       SuperscriptBox["pw2", "mu"], " ", 
       SuperscriptBox["pw2", "nu"]}], 
      SuperscriptBox["mw", "2"]], "-", 
     SuperscriptBox["g", 
      RowBox[{"mu", "nu"}]]}], ")"}]}], TraditionalForm]], "Output",
 CellChangeTimes->{
  3.4289272789512367`*^9, 3.428927318445057*^9, {3.4289273534290743`*^9, 
   3.428927404797317*^9}, {3.4289274516566753`*^9, 3.42892749447148*^9}, 
   3.4289291227183*^9, {3.428929370895793*^9, 3.428929376751039*^9}, 
   3.428929567605669*^9, {3.428938112386113*^9, 3.428938125511932*^9}, 
   3.4289382307430983`*^9, {3.428938486376433*^9, 3.4289384975076303`*^9}, 
   3.428939015665769*^9, 3.428939149066834*^9, 3.4289393600066032`*^9, 
   3.428939394733592*^9, {3.428939598302187*^9, 3.42893961837188*^9}}]
}, Open  ]],

Cell["We constrat the Lorentz indices", "Text",
 CellChangeTimes->{{3.428939152271905*^9, 3.428939164364278*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"M2", "=", 
  RowBox[{
   RowBox[{"Contract", "[", "M2", "]"}], "//", "Simplify"}]}]], "Input",
 CellChangeTimes->{{3.428939165250428*^9, 3.428939180998008*^9}}],

Cell[BoxData[
 FormBox[
  FractionBox[
   RowBox[{
    SuperscriptBox["gw", "4"], " ", 
    RowBox[{"(", 
     RowBox[{
      SuperscriptBox["mh", "4"], "-", 
      RowBox[{"4", " ", 
       SuperscriptBox["mw", "2"], " ", 
       SuperscriptBox["mh", "2"]}], "+", 
      RowBox[{"12", " ", 
       SuperscriptBox["mw", "4"]}]}], ")"}], " ", 
    SuperscriptBox["v", "2"]}], 
   RowBox[{"16", " ", 
    SuperscriptBox["mw", "4"]}]], TraditionalForm]], "Output",
 CellChangeTimes->{{3.428939173219203*^9, 3.428939181353017*^9}, 
   3.428939361929927*^9, 3.428939396412657*^9, 3.428939619873785*^9}]
}, Open  ]],

Cell["We insert he definition of the Fermi constant", "Text",
 CellChangeTimes->{{3.428927552171399*^9, 3.428927561923086*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"M2", "=", 
  RowBox[{
   RowBox[{"M2", "/.", 
    RowBox[{"v", "\[Rule]", 
     RowBox[{"1", "/", 
      RowBox[{"Sqrt", "[", 
       RowBox[{"Gf", " ", 
        RowBox[{"Sqrt", "[", "2", "]"}]}], "]"}]}]}]}], "/.", 
   RowBox[{"gw", "\[Rule]", 
    RowBox[{"Sqrt", "[", 
     RowBox[{"Gf", " ", "8", 
      RowBox[{
       RowBox[{"mw", "^", "2"}], "/", 
       RowBox[{"Sqrt", "[", "2", "]"}]}]}], "]"}]}]}]}]], "Input",
 CellChangeTimes->{{3.428927568404044*^9, 3.42892758705121*^9}, {
  3.4289381595073957`*^9, 3.4289381716760798`*^9}, {3.428939195889434*^9, 
  3.428939216735014*^9}}],

Cell[BoxData[
 FormBox[
  RowBox[{
   SqrtBox["2"], " ", "Gf", " ", 
   RowBox[{"(", 
    RowBox[{
     SuperscriptBox["mh", "4"], "-", 
     RowBox[{"4", " ", 
      SuperscriptBox["mw", "2"], " ", 
      SuperscriptBox["mh", "2"]}], "+", 
     RowBox[{"12", " ", 
      SuperscriptBox["mw", "4"]}]}], ")"}]}], TraditionalForm]], "Output",
 CellChangeTimes->{3.428927588058874*^9, 3.4289291256255703`*^9, 
  3.428929571106011*^9, 3.428938172246612*^9, 3.428938233053894*^9, 
  3.4289384978556747`*^9, 3.428939015694807*^9, 3.428939217535921*^9, 
  3.428939398074815*^9, 3.428939621035542*^9}]
}, Open  ]]
}, Closed]],

Cell[CellGroupData[{

Cell["Phase Space", "Subsection",
 CellChangeTimes->{{3.4289292202493353`*^9, 3.428929221917066*^9}}],

Cell["The two particle phase space is", "Text",
 CellChangeTimes->{{3.4289286541859837`*^9, 3.428928662992796*^9}, {
  3.428938253498197*^9, 3.428938253866067*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"d\[CapitalPhi]2", "=", 
  RowBox[{
   RowBox[{
    RowBox[{"(", 
     RowBox[{"2", "Pi"}], ")"}], "^", "4"}], 
   RowBox[{"\[Delta]", "[", 
    RowBox[{"P", "-", "pw1", "-", "pw2"}], "]"}], 
   RowBox[{
    RowBox[{
     RowBox[{"d3pw1", "/", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{"2", "Pi"}], ")"}], "^", "3"}]}], "/", "2"}], "/", "E1"}], 
   " ", 
   RowBox[{
    RowBox[{
     RowBox[{"d3pw2", "/", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{"2", "Pi"}], ")"}], "^", "3"}]}], "/", "2"}], "/", 
    "E2"}]}]}]], "Input",
 CellChangeTimes->{{3.428928668808182*^9, 3.428928748243292*^9}, {
  3.4289288594542103`*^9, 3.428928869884198*^9}, {3.42893818360579*^9, 
  3.428938207781085*^9}, {3.4289382412163973`*^9, 3.4289382417010736`*^9}, {
  3.428939402360241*^9, 3.428939407583867*^9}}],

Cell[BoxData[
 FormBox[
  FractionBox[
   RowBox[{"d3pw1", " ", "d3pw2", " ", 
    RowBox[{"\[Delta]", "(", 
     RowBox[{"P", "-", "pw1", "-", "pw2"}], ")"}]}], 
   RowBox[{"16", " ", "E1", " ", "E2", " ", 
    SuperscriptBox["\[Pi]", "2"]}]], TraditionalForm]], "Output",
 CellChangeTimes->{{3.428928862716354*^9, 3.428928870072042*^9}, 
   3.428929126861431*^9, 3.428929574710586*^9, {3.428938234752118*^9, 
   3.428938241930442*^9}, 3.428938501869619*^9, 3.4289390157281837`*^9, 
   3.428939224535726*^9, 3.4289394090399523`*^9, 3.428939622836937*^9}]
}, Open  ]],

Cell["\<\
Performing the integration over d3pw2 in the rest frame of the H gives\
\>", "Text",
 CellChangeTimes->{{3.428928749821361*^9, 3.428928784904477*^9}, {
   3.4289288233373938`*^9, 3.428928824472897*^9}, {3.428938246362335*^9, 
   3.428938248774928*^9}, 3.428939418781519*^9}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"PS", "=", 
  RowBox[{
   RowBox[{
    RowBox[{"1", "/", "4"}], "/", 
    RowBox[{
     RowBox[{"(", 
      RowBox[{"2", "Pi"}], ")"}], "^", "2"}]}], " ", 
   RowBox[{"\[Delta]", "[", 
    RowBox[{"mh", "-", 
     RowBox[{"2", "E1"}]}], "]"}], 
   RowBox[{"d3pw1", "/", 
    RowBox[{"E1", "^", "2"}]}]}]}]], "Input",
 CellChangeTimes->{{3.4289287922794847`*^9, 3.4289288819309*^9}, {
   3.42893825962851*^9, 3.42893826484441*^9}, 3.428939414455236*^9}],

Cell[BoxData[
 FormBox[
  FractionBox[
   RowBox[{"d3pw1", " ", 
    RowBox[{"\[Delta]", "(", 
     RowBox[{"mh", "-", 
      RowBox[{"2", " ", "E1"}]}], ")"}]}], 
   RowBox[{"16", " ", 
    SuperscriptBox["E1", "2"], " ", 
    SuperscriptBox["\[Pi]", "2"]}]], TraditionalForm]], "Output",
 CellChangeTimes->{{3.428928849656254*^9, 3.428928882523725*^9}, 
   3.428929068231016*^9, 3.428929127845537*^9, 3.428929574829381*^9, 
   3.42893826621164*^9, 3.428938504167341*^9, 3.4289386058116293`*^9, 
   3.4289386914875193`*^9, 3.4289390157737207`*^9, {3.428939227658648*^9, 
   3.428939241509007*^9}, 3.428939424632287*^9, 3.428939623873477*^9}]
}, Open  ]],

Cell["\<\
Going to spherical coordinates, and usging pw1 dpw1 = E1 dE1, we get\
\>", "Text",
 CellChangeTimes->{{3.428928883567082*^9, 3.4289289264749002`*^9}, {
  3.4289382922025423`*^9, 3.428938296218473*^9}, {3.428939420508946*^9, 
  3.4289394222210503`*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"PS", "=", 
  RowBox[{
   RowBox[{"PS", "/.", 
    RowBox[{
     RowBox[{"\[Delta]", "[", 
      RowBox[{"mh", "-", 
       RowBox[{"2", "E1"}]}], "]"}], "\[Rule]", 
     RowBox[{
      RowBox[{"1", "/", "2"}], 
      RowBox[{"\[Delta]", "[", 
       RowBox[{"E1", "-", 
        RowBox[{"mh", "/", "2"}]}], "]"}]}]}]}], "/.", " ", 
   RowBox[{"d3pw1", "\[Rule]", 
    RowBox[{"4", "Pi", " ", 
     RowBox[{"Sqrt", "[", 
      RowBox[{
       RowBox[{"E1", "^", "2"}], "-", 
       RowBox[{"mw", "^", "2"}]}], "]"}], "E1", " ", "dE1"}]}]}]}]], "Input",
 CellChangeTimes->{{3.4289289441444197`*^9, 3.428929022421103*^9}, 
   3.428929070151713*^9, {3.4289382799967327`*^9, 3.428938290196628*^9}, {
   3.428938584824209*^9, 3.428938599062233*^9}, {3.428939238657551*^9, 
   3.4289392398953457`*^9}, 3.428939430416103*^9}],

Cell[BoxData[
 FormBox[
  FractionBox[
   RowBox[{"dE1", " ", 
    SqrtBox[
     RowBox[{
      SuperscriptBox["E1", "2"], "-", 
      SuperscriptBox["mw", "2"]}]], " ", 
    RowBox[{"\[Delta]", "(", 
     RowBox[{"E1", "-", 
      FractionBox["mh", "2"]}], ")"}]}], 
   RowBox[{"8", " ", "E1", " ", "\[Pi]"}]], TraditionalForm]], "Output",
 CellChangeTimes->{
  3.428929025150517*^9, 3.428929070325218*^9, 3.4289291311060667`*^9, 
   3.428929574875102*^9, 3.428938297415737*^9, 3.428938505065439*^9, {
   3.428938600499318*^9, 3.428938606660725*^9}, {3.4289386861864023`*^9, 
   3.428938692218741*^9}, 3.4289390158068933`*^9, {3.428939229150958*^9, 
   3.428939242434461*^9}, 3.4289394312423964`*^9, 3.428939624808075*^9}]
}, Open  ]],

Cell["Performing the integration", "Text",
 CellChangeTimes->{{3.428929031938621*^9, 3.428929037282918*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"PS", "=", 
  RowBox[{
   RowBox[{"PS", "/.", 
    RowBox[{
     RowBox[{"dE1", " ", 
      RowBox[{"\[Delta]", "[", 
       RowBox[{"E1", "-", 
        RowBox[{"mh", "/", "2"}]}], "]"}]}], "\[Rule]", "1"}]}], "/.", 
   RowBox[{"E1", "\[Rule]", 
    RowBox[{"mh", "/", "2"}]}]}]}]], "Input",
 CellChangeTimes->{{3.428929044586605*^9, 3.428929056603554*^9}, {
  3.428938304472186*^9, 3.428938307564096*^9}, {3.428938621479186*^9, 
  3.428938627852805*^9}, {3.4289386823375387`*^9, 3.428938682419998*^9}}],

Cell[BoxData[
 FormBox[
  FractionBox[
   SqrtBox[
    RowBox[{
     FractionBox[
      SuperscriptBox["mh", "2"], "4"], "-", 
     SuperscriptBox["mw", "2"]}]], 
   RowBox[{"4", " ", "mh", " ", "\[Pi]"}]], TraditionalForm]], "Output",
 CellChangeTimes->{{3.428929056861126*^9, 3.428929071114458*^9}, 
   3.4289291321148567`*^9, 3.428929574912538*^9, 3.428938307856097*^9, 
   3.428938506325878*^9, 3.4289386297502537`*^9, {3.4289386829460487`*^9, 
   3.428938694746279*^9}, 3.42893901583956*^9, 3.428939246350231*^9, 
   3.4289394342927837`*^9, 3.428939625835762*^9}]
}, Open  ]]
}, Closed]],

Cell[CellGroupData[{

Cell["Result", "Subsection",
 CellChangeTimes->{{3.4289292627868023`*^9, 3.4289292633475027`*^9}}],

Cell["The decay rate of the H then becomes", "Text",
 CellChangeTimes->{{3.428929073874954*^9, 3.428929081291203*^9}, 
   3.4289383237898083`*^9}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"\[CapitalGamma]", "=", 
  RowBox[{
   RowBox[{
    RowBox[{
     RowBox[{
      RowBox[{
       RowBox[{"1", "/", "2"}], "/", "mh"}], " ", "PS", " ", "M2"}], "/.", 
     RowBox[{"mw", "\[Rule]", 
      RowBox[{"x", " ", "mh"}]}]}], "//", "Simplify"}], "//", 
   "PowerExpand"}]}]], "Input",
 CellChangeTimes->{{3.428929087038114*^9, 3.428929107627933*^9}, 
   3.428929346599481*^9, {3.428938314727579*^9, 3.428938334043777*^9}, 
   3.42893864433534*^9, 3.4289389462936773`*^9, {3.428939734911998*^9, 
   3.428939761347455*^9}}],

Cell[BoxData[
 FormBox[
  FractionBox[
   RowBox[{"Gf", " ", 
    SuperscriptBox["mh", "3"], " ", 
    SqrtBox[
     RowBox[{"1", "-", 
      RowBox[{"4", " ", 
       SuperscriptBox["x", "2"]}]}]], " ", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"12", " ", 
       SuperscriptBox["x", "4"]}], "-", 
      RowBox[{"4", " ", 
       SuperscriptBox["x", "2"]}], "+", "1"}], ")"}]}], 
   RowBox[{"8", " ", 
    SqrtBox["2"], " ", "\[Pi]"}]], TraditionalForm]], "Output",
 CellChangeTimes->{{3.428929108658009*^9, 3.428929133003119*^9}, 
   3.4289293476527367`*^9, 3.4289295771229353`*^9, {3.428938319056734*^9, 
   3.4289383343919973`*^9}, 3.4289385085419197`*^9, {3.4289386319695272`*^9, 
   3.428938644554867*^9}, 3.428938696988208*^9, 3.428938947199583*^9, 
   3.42893901589124*^9, 3.428939251053282*^9, 3.4289394364600153`*^9, 
   3.4289396272755547`*^9, {3.42893974381048*^9, 3.428939762040904*^9}}]
}, Open  ]]
}, Closed]]
}, Open  ]],

Cell[CellGroupData[{

Cell["H \[RightArrow] ZZ", "Section",
 CellChangeTimes->{{3.4289377207239313`*^9, 3.428937731049849*^9}, {
  3.428939028678343*^9, 3.428939029987813*^9}, {3.428939649559754*^9, 
  3.4289396513572483`*^9}}],

Cell[CellGroupData[{

Cell["Kinematics 1->2", "Subsection",
 CellChangeTimes->{{3.428926835790217*^9, 3.428926837886043*^9}}],

Cell[BoxData[
 RowBox[{"\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{
     RowBox[{"ScalarProduct", "[", 
      RowBox[{"P", ",", "P"}], "]"}], "=", 
     RowBox[{"mh", "^", "2"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{
     RowBox[{"ScalarProduct", "[", 
      RowBox[{"pz1", ",", "pz1"}], "]"}], "=", 
     RowBox[{"mz", "^", "2"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{
     RowBox[{"ScalarProduct", "[", 
      RowBox[{"pz2", ",", "pz2"}], "]"}], "=", 
     RowBox[{"mz", "^", "2"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{
     RowBox[{"ScalarProduct", "[", 
      RowBox[{"P", ",", "pz1"}], "]"}], "=", 
     RowBox[{
      RowBox[{"mh", "^", "2"}], "/", "2"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{
     RowBox[{"ScalarProduct", "[", 
      RowBox[{"P", ",", "pz2"}], "]"}], "=", 
     RowBox[{
      RowBox[{"mh", "^", "2"}], "/", "2"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{
     RowBox[{"ScalarProduct", "[", 
      RowBox[{"pz1", ",", "pz2"}], "]"}], "=", 
     RowBox[{
      RowBox[{
       RowBox[{"mh", "^", "2"}], "/", "2"}], "-", 
      RowBox[{"mz", "^", "2"}]}]}], ";"}], "\[IndentingNewLine]"}]}]], "Input",\

 CellChangeTimes->{{3.4289268403011503`*^9, 3.428926946965074*^9}, {
  3.428937810100576*^9, 3.428937873020487*^9}, {3.428938071207912*^9, 
  3.428938078071608*^9}, {3.428938213984881*^9, 3.428938224130307*^9}, {
  3.4289384819556913`*^9, 3.428938484248788*^9}, {3.428939073598165*^9, 
  3.428939082348546*^9}, {3.42893937775745*^9, 3.4289393838581676`*^9}, {
  3.428939657707238*^9, 3.428939671283299*^9}}]
}, Closed]],

Cell[CellGroupData[{

Cell["Amplitude (1 diagram)", "Subsection",
 CellChangeTimes->{{3.4289292421393127`*^9, 3.428929244675111*^9}}],

Cell["The summed and averaged square matrix element", "Text",
 CellChangeTimes->{{3.4289271420095243`*^9, 3.4289271454722757`*^9}, {
  3.4289273354903717`*^9, 3.428927351049911*^9}, {3.4289275105378*^9, 
  3.42892751343362*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"M2", "=", 
  RowBox[{
   RowBox[{
    RowBox[{"(", 
     RowBox[{"v", " ", 
      RowBox[{
       RowBox[{
        RowBox[{"gw", "^", "2"}], "/", "2"}], "/", 
       RowBox[{"cw", "^", "2"}]}]}], ")"}], "^", "2"}], 
   RowBox[{"(", 
    RowBox[{
     RowBox[{"-", 
      RowBox[{"MetricTensor", "[", 
       RowBox[{"mu", ",", "nu"}], "]"}]}], "+", 
     RowBox[{
      RowBox[{"FV", "[", 
       RowBox[{"pz1", ",", "mu"}], "]"}], 
      RowBox[{
       RowBox[{"FV", "[", 
        RowBox[{"pz1", ",", "nu"}], "]"}], "/", 
       RowBox[{"mz", "^", "2"}]}]}]}], ")"}], 
   RowBox[{"(", 
    RowBox[{
     RowBox[{"-", 
      RowBox[{"MetricTensor", "[", 
       RowBox[{"mu", ",", "nu"}], "]"}]}], "+", 
     RowBox[{
      RowBox[{"FV", "[", 
       RowBox[{"pz2", ",", "mu"}], "]"}], 
      RowBox[{
       RowBox[{"FV", "[", 
        RowBox[{"pz2", ",", "nu"}], "]"}], "/", 
       RowBox[{"mz", "^", "2"}]}]}]}], ")"}]}]}]], "Input",
 CellChangeTimes->{{3.428927033246361*^9, 3.4289271401967373`*^9}, {
   3.42892727787171*^9, 3.428927413135103*^9}, {3.4289274432481937`*^9, 
   3.4289274930463*^9}, {3.428929121395401*^9, 3.428929122088388*^9}, 
   3.428929375657975*^9, {3.428929541763723*^9, 3.428929566521463*^9}, {
   3.428937860468915*^9, 3.428937864548534*^9}, {3.4289380364317427`*^9, 
   3.428938062823503*^9}, {3.42893809477054*^9, 3.4289381246774197`*^9}, {
   3.4289382285018063`*^9, 3.428938230286139*^9}, {3.4289390981304398`*^9, 
   3.428939148746867*^9}, {3.428939388292745*^9, 3.428939394515564*^9}, 
   3.428939614471209*^9, {3.4289396974209337`*^9, 3.428939729732409*^9}, {
   3.42893977829564*^9, 3.4289397787246437`*^9}}],

Cell[BoxData[
 FormBox[
  FractionBox[
   RowBox[{
    SuperscriptBox["gw", "4"], " ", 
    SuperscriptBox["v", "2"], " ", 
    RowBox[{"(", 
     RowBox[{
      FractionBox[
       RowBox[{
        SuperscriptBox["pz1", "mu"], " ", 
        SuperscriptBox["pz1", "nu"]}], 
       SuperscriptBox["mz", "2"]], "-", 
      SuperscriptBox["g", 
       RowBox[{"mu", "nu"}]]}], ")"}], " ", 
    RowBox[{"(", 
     RowBox[{
      FractionBox[
       RowBox[{
        SuperscriptBox["pz2", "mu"], " ", 
        SuperscriptBox["pz2", "nu"]}], 
       SuperscriptBox["mz", "2"]], "-", 
      SuperscriptBox["g", 
       RowBox[{"mu", "nu"}]]}], ")"}]}], 
   RowBox[{"4", " ", 
    SuperscriptBox["cw", "4"]}]], TraditionalForm]], "Output",
 CellChangeTimes->{
  3.4289272789512367`*^9, 3.428927318445057*^9, {3.4289273534290743`*^9, 
   3.428927404797317*^9}, {3.4289274516566753`*^9, 3.42892749447148*^9}, 
   3.4289291227183*^9, {3.428929370895793*^9, 3.428929376751039*^9}, 
   3.428929567605669*^9, {3.428938112386113*^9, 3.428938125511932*^9}, 
   3.4289382307430983`*^9, {3.428938486376433*^9, 3.4289384975076303`*^9}, 
   3.428939015665769*^9, 3.428939149066834*^9, 3.4289393600066032`*^9, 
   3.428939394733592*^9, {3.428939598302187*^9, 3.42893961837188*^9}, 
   3.428939778926634*^9, 3.4289398622949247`*^9}]
}, Open  ]],

Cell["We constrat the Lorentz indices", "Text",
 CellChangeTimes->{{3.428939152271905*^9, 3.428939164364278*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"M2", "=", 
  RowBox[{
   RowBox[{"Contract", "[", "M2", "]"}], "//", "Simplify"}]}]], "Input",
 CellChangeTimes->{{3.428939165250428*^9, 3.428939180998008*^9}}],

Cell[BoxData[
 FormBox[
  FractionBox[
   RowBox[{
    SuperscriptBox["gw", "4"], " ", 
    RowBox[{"(", 
     RowBox[{
      SuperscriptBox["mh", "4"], "-", 
      RowBox[{"4", " ", 
       SuperscriptBox["mz", "2"], " ", 
       SuperscriptBox["mh", "2"]}], "+", 
      RowBox[{"12", " ", 
       SuperscriptBox["mz", "4"]}]}], ")"}], " ", 
    SuperscriptBox["v", "2"]}], 
   RowBox[{"16", " ", 
    SuperscriptBox["cw", "4"], " ", 
    SuperscriptBox["mz", "4"]}]], TraditionalForm]], "Output",
 CellChangeTimes->{{3.428939173219203*^9, 3.428939181353017*^9}, 
   3.428939361929927*^9, 3.428939396412657*^9, 3.428939619873785*^9, 
   3.428939782429967*^9, 3.428939869387093*^9}]
}, Open  ]],

Cell["We insert he definition of the Fermi constant", "Text",
 CellChangeTimes->{{3.428927552171399*^9, 3.428927561923086*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"M2", "=", 
  RowBox[{
   RowBox[{"M2", "/.", 
    RowBox[{"v", "\[Rule]", 
     RowBox[{"1", "/", 
      RowBox[{"Sqrt", "[", 
       RowBox[{"Gf", " ", 
        RowBox[{"Sqrt", "[", "2", "]"}]}], "]"}]}]}]}], "/.", 
   RowBox[{"gw", "\[Rule]", 
    RowBox[{"Sqrt", "[", 
     RowBox[{"Gf", " ", "8", 
      RowBox[{
       RowBox[{"mw", "^", "2"}], "/", 
       RowBox[{"Sqrt", "[", "2", "]"}]}]}], "]"}]}]}]}]], "Input",
 CellChangeTimes->{{3.428927568404044*^9, 3.42892758705121*^9}, {
  3.4289381595073957`*^9, 3.4289381716760798`*^9}, {3.428939195889434*^9, 
  3.428939216735014*^9}}],

Cell[BoxData[
 FormBox[
  FractionBox[
   RowBox[{
    SqrtBox["2"], " ", "Gf", " ", 
    SuperscriptBox["mw", "4"], " ", 
    RowBox[{"(", 
     RowBox[{
      SuperscriptBox["mh", "4"], "-", 
      RowBox[{"4", " ", 
       SuperscriptBox["mz", "2"], " ", 
       SuperscriptBox["mh", "2"]}], "+", 
      RowBox[{"12", " ", 
       SuperscriptBox["mz", "4"]}]}], ")"}]}], 
   RowBox[{
    SuperscriptBox["cw", "4"], " ", 
    SuperscriptBox["mz", "4"]}]], TraditionalForm]], "Output",
 CellChangeTimes->{3.428927588058874*^9, 3.4289291256255703`*^9, 
  3.428929571106011*^9, 3.428938172246612*^9, 3.428938233053894*^9, 
  3.4289384978556747`*^9, 3.428939015694807*^9, 3.428939217535921*^9, 
  3.428939398074815*^9, 3.428939621035542*^9, 3.428939786171034*^9, 
  3.4289398708976583`*^9}]
}, Open  ]]
}, Closed]],

Cell[CellGroupData[{

Cell["Phase Space", "Subsection",
 CellChangeTimes->{{3.4289292202493353`*^9, 3.428929221917066*^9}}],

Cell["The two particle phase space is", "Text",
 CellChangeTimes->{{3.4289286541859837`*^9, 3.428928662992796*^9}, {
  3.428938253498197*^9, 3.428938253866067*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"d\[CapitalPhi]2", "=", 
  RowBox[{
   RowBox[{"1", "/", "2"}], 
   RowBox[{
    RowBox[{"(", 
     RowBox[{"2", "Pi"}], ")"}], "^", "4"}], 
   RowBox[{"\[Delta]", "[", 
    RowBox[{"P", "-", "pz1", "-", "pz2"}], "]"}], 
   RowBox[{
    RowBox[{
     RowBox[{"d3pz1", "/", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{"2", "Pi"}], ")"}], "^", "3"}]}], "/", "2"}], "/", "E1"}], 
   " ", 
   RowBox[{
    RowBox[{
     RowBox[{"d3pz2", "/", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{"2", "Pi"}], ")"}], "^", "3"}]}], "/", "2"}], "/", 
    "E2"}]}]}]], "Input",
 CellChangeTimes->{{3.428928668808182*^9, 3.428928748243292*^9}, {
  3.4289288594542103`*^9, 3.428928869884198*^9}, {3.42893818360579*^9, 
  3.428938207781085*^9}, {3.4289382412163973`*^9, 3.4289382417010736`*^9}, {
  3.428939402360241*^9, 3.428939407583867*^9}, {3.4289397947484217`*^9, 
  3.4289398003366203`*^9}, {3.42893988381181*^9, 3.428939884041622*^9}}],

Cell[BoxData[
 FormBox[
  FractionBox[
   RowBox[{"d3pz1", " ", "d3pz2", " ", 
    RowBox[{"\[Delta]", "(", 
     RowBox[{"P", "-", "pz1", "-", "pz2"}], ")"}]}], 
   RowBox[{"16", " ", "E1", " ", "E2", " ", 
    SuperscriptBox["\[Pi]", "2"]}]], TraditionalForm]], "Output",
 CellChangeTimes->{{3.428928862716354*^9, 3.428928870072042*^9}, 
   3.428929126861431*^9, 3.428929574710586*^9, {3.428938234752118*^9, 
   3.428938241930442*^9}, 3.428938501869619*^9, 3.4289390157281837`*^9, 
   3.428939224535726*^9, 3.4289394090399523`*^9, 3.428939622836937*^9, 
   3.428939802399288*^9, 3.4289398723729563`*^9}]
}, Open  ]],

Cell["\<\
Notice the additional factor 1/2 in the phase spac, coming form the fact that \
we have identical particles in the final state.\
\>", "Text",
 CellChangeTimes->{{3.428939885985323*^9, 3.428939907485496*^9}}],

Cell["\<\
Performing the integration over d3pz2 in the rest frame of the H gives\
\>", "Text",
 CellChangeTimes->{{3.428928749821361*^9, 3.428928784904477*^9}, {
   3.4289288233373938`*^9, 3.428928824472897*^9}, {3.428938246362335*^9, 
   3.428938248774928*^9}, 3.428939418781519*^9, {3.428939816157881*^9, 
   3.4289398163019114`*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"PS", "=", 
  RowBox[{
   RowBox[{"1", "/", "2"}], " ", 
   RowBox[{
    RowBox[{"1", "/", "4"}], "/", 
    RowBox[{
     RowBox[{"(", 
      RowBox[{"2", "Pi"}], ")"}], "^", "2"}]}], " ", 
   RowBox[{"\[Delta]", "[", 
    RowBox[{"mh", "-", 
     RowBox[{"2", "E1"}]}], "]"}], 
   RowBox[{"d3pz1", "/", 
    RowBox[{"E1", "^", "2"}]}]}]}]], "Input",
 CellChangeTimes->{{3.4289287922794847`*^9, 3.4289288819309*^9}, {
   3.42893825962851*^9, 3.42893826484441*^9}, 3.428939414455236*^9, {
   3.4289398057924757`*^9, 3.42893980595168*^9}, {3.428939917120672*^9, 
   3.428939918020076*^9}}],

Cell[BoxData[
 FormBox[
  FractionBox[
   RowBox[{"d3pz1", " ", 
    RowBox[{"\[Delta]", "(", 
     RowBox[{"mh", "-", 
      RowBox[{"2", " ", "E1"}]}], ")"}]}], 
   RowBox[{"32", " ", 
    SuperscriptBox["E1", "2"], " ", 
    SuperscriptBox["\[Pi]", "2"]}]], TraditionalForm]], "Output",
 CellChangeTimes->{{3.428928849656254*^9, 3.428928882523725*^9}, 
   3.428929068231016*^9, 3.428929127845537*^9, 3.428929574829381*^9, 
   3.42893826621164*^9, 3.428938504167341*^9, 3.4289386058116293`*^9, 
   3.4289386914875193`*^9, 3.4289390157737207`*^9, {3.428939227658648*^9, 
   3.428939241509007*^9}, 3.428939424632287*^9, 3.428939623873477*^9, 
   3.4289398062227793`*^9, 3.428939873429842*^9, 3.428939918388637*^9}]
}, Open  ]],

Cell["\<\
Going to spherical coordinates, and usging pz1 dpz1 = E1 dE1, we get\
\>", "Text",
 CellChangeTimes->{{3.428928883567082*^9, 3.4289289264749002`*^9}, {
  3.4289382922025423`*^9, 3.428938296218473*^9}, {3.428939420508946*^9, 
  3.4289394222210503`*^9}, {3.428939811789967*^9, 3.4289398135420017`*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"PS", "=", 
  RowBox[{
   RowBox[{"PS", "/.", 
    RowBox[{
     RowBox[{"\[Delta]", "[", 
      RowBox[{"mh", "-", 
       RowBox[{"2", "E1"}]}], "]"}], "\[Rule]", 
     RowBox[{
      RowBox[{"1", "/", "2"}], 
      RowBox[{"\[Delta]", "[", 
       RowBox[{"E1", "-", 
        RowBox[{"mh", "/", "2"}]}], "]"}]}]}]}], "/.", " ", 
   RowBox[{"d3pz1", "\[Rule]", 
    RowBox[{"4", "Pi", " ", 
     RowBox[{"Sqrt", "[", 
      RowBox[{
       RowBox[{"E1", "^", "2"}], "-", 
       RowBox[{"mz", "^", "2"}]}], "]"}], "E1", " ", "dE1"}]}]}]}]], "Input",
 CellChangeTimes->{{3.4289289441444197`*^9, 3.428929022421103*^9}, 
   3.428929070151713*^9, {3.4289382799967327`*^9, 3.428938290196628*^9}, {
   3.428938584824209*^9, 3.428938599062233*^9}, {3.428939238657551*^9, 
   3.4289392398953457`*^9}, 3.428939430416103*^9, {3.428939810323925*^9, 
   3.428939819169014*^9}}],

Cell[BoxData[
 FormBox[
  FractionBox[
   RowBox[{"dE1", " ", 
    SqrtBox[
     RowBox[{
      SuperscriptBox["E1", "2"], "-", 
      SuperscriptBox["mz", "2"]}]], " ", 
    RowBox[{"\[Delta]", "(", 
     RowBox[{"E1", "-", 
      FractionBox["mh", "2"]}], ")"}]}], 
   RowBox[{"16", " ", "E1", " ", "\[Pi]"}]], TraditionalForm]], "Output",
 CellChangeTimes->{
  3.428929025150517*^9, 3.428929070325218*^9, 3.4289291311060667`*^9, 
   3.428929574875102*^9, 3.428938297415737*^9, 3.428938505065439*^9, {
   3.428938600499318*^9, 3.428938606660725*^9}, {3.4289386861864023`*^9, 
   3.428938692218741*^9}, 3.4289390158068933`*^9, {3.428939229150958*^9, 
   3.428939242434461*^9}, 3.4289394312423964`*^9, 3.428939624808075*^9, 
   3.4289398195382347`*^9, 3.428939874626782*^9, 3.428939921689764*^9}]
}, Open  ]],

Cell["Performing the integration", "Text",
 CellChangeTimes->{{3.428929031938621*^9, 3.428929037282918*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"PS", "=", 
  RowBox[{
   RowBox[{"PS", "/.", 
    RowBox[{
     RowBox[{"dE1", " ", 
      RowBox[{"\[Delta]", "[", 
       RowBox[{"E1", "-", 
        RowBox[{"mh", "/", "2"}]}], "]"}]}], "\[Rule]", "1"}]}], "/.", 
   RowBox[{"E1", "\[Rule]", 
    RowBox[{"mh", "/", "2"}]}]}]}]], "Input",
 CellChangeTimes->{{3.428929044586605*^9, 3.428929056603554*^9}, {
  3.428938304472186*^9, 3.428938307564096*^9}, {3.428938621479186*^9, 
  3.428938627852805*^9}, {3.4289386823375387`*^9, 3.428938682419998*^9}}],

Cell[BoxData[
 FormBox[
  FractionBox[
   SqrtBox[
    RowBox[{
     FractionBox[
      SuperscriptBox["mh", "2"], "4"], "-", 
     SuperscriptBox["mz", "2"]}]], 
   RowBox[{"8", " ", "mh", " ", "\[Pi]"}]], TraditionalForm]], "Output",
 CellChangeTimes->{{3.428929056861126*^9, 3.428929071114458*^9}, 
   3.4289291321148567`*^9, 3.428929574912538*^9, 3.428938307856097*^9, 
   3.428938506325878*^9, 3.4289386297502537`*^9, {3.4289386829460487`*^9, 
   3.428938694746279*^9}, 3.42893901583956*^9, 3.428939246350231*^9, 
   3.4289394342927837`*^9, 3.428939625835762*^9, 3.428939823598246*^9, 
   3.428939875651293*^9, 3.4289399228618383`*^9}]
}, Open  ]]
}, Closed]],

Cell[CellGroupData[{

Cell["Result", "Subsection",
 CellChangeTimes->{{3.4289292627868023`*^9, 3.4289292633475027`*^9}}],

Cell["The decay rate of the H then becomes", "Text",
 CellChangeTimes->{{3.428929073874954*^9, 3.428929081291203*^9}, 
   3.4289383237898083`*^9}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"\[CapitalGamma]", "=", 
  RowBox[{
   RowBox[{
    RowBox[{
     RowBox[{
      RowBox[{
       RowBox[{
        RowBox[{"1", "/", "2"}], "/", "mh"}], " ", "PS", " ", "M2"}], "/.", 
      RowBox[{"mw", "\[Rule]", 
       RowBox[{"mz", " ", "cw"}]}]}], "/.", 
     RowBox[{"mz", "\[Rule]", 
      RowBox[{"x", " ", "mh"}]}]}], "//", "Simplify"}], "//", 
   "PowerExpand"}]}]], "Input",
 CellChangeTimes->{{3.428929087038114*^9, 3.428929107627933*^9}, 
   3.428929346599481*^9, {3.428938314727579*^9, 3.428938334043777*^9}, 
   3.42893864433534*^9, 3.4289389462936773`*^9, {3.42893982964193*^9, 
   3.4289398502038717`*^9}}],

Cell[BoxData[
 FormBox[
  FractionBox[
   RowBox[{"Gf", " ", 
    SuperscriptBox["mh", "3"], " ", 
    SqrtBox[
     RowBox[{"1", "-", 
      RowBox[{"4", " ", 
       SuperscriptBox["x", "2"]}]}]], " ", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"12", " ", 
       SuperscriptBox["x", "4"]}], "-", 
      RowBox[{"4", " ", 
       SuperscriptBox["x", "2"]}], "+", "1"}], ")"}]}], 
   RowBox[{"16", " ", 
    SqrtBox["2"], " ", "\[Pi]"}]], TraditionalForm]], "Output",
 CellChangeTimes->{{3.428929108658009*^9, 3.428929133003119*^9}, 
   3.4289293476527367`*^9, 3.4289295771229353`*^9, {3.428938319056734*^9, 
   3.4289383343919973`*^9}, 3.4289385085419197`*^9, {3.4289386319695272`*^9, 
   3.428938644554867*^9}, 3.428938696988208*^9, 3.428938947199583*^9, 
   3.42893901589124*^9, 3.428939251053282*^9, 3.4289394364600153`*^9, 
   3.4289396272755547`*^9, {3.4289398454713497`*^9, 3.4289398770178833`*^9}, 
   3.428939924241026*^9}]
}, Open  ]]
}, Closed]]
}, Open  ]]
}, Open  ]]
},
WindowSize->{710, 706},
WindowMargins->{{Automatic, 230}, {Automatic, 4}},
ShowSelection->True,
FrontEndVersion->"6.0 for Mac OS X x86 (32-bit) (April 20, 2007)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[CellGroupData[{
Cell[590, 23, 175, 4, 76, "Title"],
Cell[CellGroupData[{
Cell[790, 31, 36, 0, 34, "Subsection"],
Cell[829, 33, 83, 2, 27, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[949, 40, 105, 1, 67, "Section"],
Cell[CellGroupData[{
Cell[1079, 45, 103, 1, 34, "Subsection"],
Cell[1185, 48, 1494, 42, 133, "Input"]
}, Closed]],
Cell[CellGroupData[{
Cell[2716, 95, 111, 1, 26, "Subsection"],
Cell[2830, 98, 228, 3, 26, "Text"],
Cell[CellGroupData[{
Cell[3083, 105, 955, 23, 27, "Input"],
Cell[4041, 130, 682, 16, 70, "Output"]
}, Open  ]],
Cell[4738, 149, 127, 1, 26, "Text"],
Cell[CellGroupData[{
Cell[4890, 154, 346, 9, 27, "Input"],
Cell[5239, 165, 485, 13, 70, "Output"]
}, Open  ]]
}, Closed]],
Cell[CellGroupData[{
Cell[5773, 184, 101, 1, 26, "Subsection"],
Cell[5877, 187, 164, 2, 26, "Text"],
Cell[CellGroupData[{
Cell[6066, 193, 782, 24, 27, "Input"],
Cell[6851, 219, 479, 10, 70, "Output"]
}, Open  ]],
Cell[7345, 232, 259, 5, 26, "Text"],
Cell[CellGroupData[{
Cell[7629, 241, 451, 14, 27, "Input"],
Cell[8083, 257, 547, 13, 70, "Output"]
}, Open  ]],
Cell[8645, 273, 209, 4, 26, "Text"],
Cell[CellGroupData[{
Cell[8879, 281, 765, 21, 27, "Input"],
Cell[9647, 304, 627, 16, 70, "Output"]
}, Open  ]],
Cell[10289, 323, 108, 1, 26, "Text"],
Cell[CellGroupData[{
Cell[10422, 328, 526, 13, 27, "Input"],
Cell[10951, 343, 496, 12, 70, "Output"]
}, Open  ]]
}, Closed]],
Cell[CellGroupData[{
Cell[11496, 361, 98, 1, 26, "Subsection"],
Cell[11597, 364, 146, 2, 26, "Text"],
Cell[CellGroupData[{
Cell[11768, 370, 373, 9, 27, "Input"],
Cell[12144, 381, 837, 23, 60, "Output"]
}, Open  ]]
}, Closed]]
}, Open  ]],
Cell[CellGroupData[{
Cell[13042, 411, 154, 2, 67, "Section"],
Cell[CellGroupData[{
Cell[13221, 417, 103, 1, 34, "Subsection"],
Cell[13327, 420, 1598, 43, 133, "Input"]
}, Closed]],
Cell[CellGroupData[{
Cell[14962, 468, 111, 1, 26, "Subsection"],
Cell[15076, 471, 228, 3, 26, "Text"],
Cell[CellGroupData[{
Cell[15329, 478, 1514, 40, 43, "Input"],
Cell[16846, 520, 1199, 31, 48, "Output"]
}, Open  ]],
Cell[18060, 554, 113, 1, 26, "Text"],
Cell[CellGroupData[{
Cell[18198, 559, 184, 4, 27, "Input"],
Cell[18385, 565, 597, 17, 51, "Output"]
}, Open  ]],
Cell[18997, 585, 127, 1, 26, "Text"],
Cell[CellGroupData[{
Cell[19149, 590, 613, 17, 27, "Input"],
Cell[19765, 609, 593, 15, 39, "Output"]
}, Open  ]]
}, Closed]],
Cell[CellGroupData[{
Cell[20407, 630, 101, 1, 26, "Subsection"],
Cell[20511, 633, 164, 2, 26, "Text"],
Cell[CellGroupData[{
Cell[20700, 639, 835, 25, 27, "Input"],
Cell[21538, 666, 555, 11, 48, "Output"]
}, Open  ]],
Cell[22108, 680, 284, 5, 26, "Text"],
Cell[CellGroupData[{
Cell[22417, 689, 475, 14, 27, "Input"],
Cell[22895, 705, 642, 14, 48, "Output"]
}, Open  ]],
Cell[23552, 722, 262, 5, 26, "Text"],
Cell[CellGroupData[{
Cell[23839, 731, 840, 22, 27, "Input"],
Cell[24682, 755, 723, 17, 59, "Output"]
}, Open  ]],
Cell[25420, 775, 108, 1, 26, "Text"],
Cell[CellGroupData[{
Cell[25553, 780, 526, 13, 27, "Input"],
Cell[26082, 795, 568, 13, 71, "Output"]
}, Open  ]]
}, Closed]],
Cell[CellGroupData[{
Cell[26699, 814, 98, 1, 26, "Subsection"],
Cell[26800, 817, 146, 2, 26, "Text"],
Cell[CellGroupData[{
Cell[26971, 823, 551, 14, 27, "Input"],
Cell[27525, 839, 905, 22, 60, "Output"]
}, Open  ]]
}, Closed]]
}, Open  ]],
Cell[CellGroupData[{
Cell[28491, 868, 205, 3, 67, "Section"],
Cell[CellGroupData[{
Cell[28721, 875, 103, 1, 34, "Subsection"],
Cell[28827, 878, 1647, 44, 133, "Input"]
}, Closed]],
Cell[CellGroupData[{
Cell[30511, 927, 111, 1, 26, "Subsection"],
Cell[30625, 930, 228, 3, 26, "Text"],
Cell[CellGroupData[{
Cell[30878, 937, 1671, 43, 43, "Input"],
Cell[32552, 982, 1310, 34, 59, "Output"]
}, Open  ]],
Cell[33877, 1019, 113, 1, 26, "Text"],
Cell[CellGroupData[{
Cell[34015, 1024, 184, 4, 27, "Input"],
Cell[34202, 1030, 682, 19, 51, "Output"]
}, Open  ]],
Cell[34899, 1052, 127, 1, 26, "Text"],
Cell[CellGroupData[{
Cell[35051, 1057, 613, 17, 27, "Input"],
Cell[35667, 1076, 788, 21, 53, "Output"]
}, Open  ]]
}, Closed]],
Cell[CellGroupData[{
Cell[36504, 1103, 101, 1, 26, "Subsection"],
Cell[36608, 1106, 164, 2, 26, "Text"],
Cell[CellGroupData[{
Cell[36797, 1112, 962, 27, 27, "Input"],
Cell[37762, 1141, 605, 12, 48, "Output"]
}, Open  ]],
Cell[38382, 1156, 217, 4, 26, "Text"],
Cell[38602, 1162, 336, 6, 26, "Text"],
Cell[CellGroupData[{
Cell[38963, 1172, 610, 17, 27, "Input"],
Cell[39576, 1191, 714, 15, 48, "Output"]
}, Open  ]],
Cell[40305, 1209, 310, 5, 26, "Text"],
Cell[CellGroupData[{
Cell[40640, 1218, 890, 23, 27, "Input"],
Cell[41533, 1243, 796, 18, 59, "Output"]
}, Open  ]],
Cell[42344, 1264, 108, 1, 26, "Text"],
Cell[CellGroupData[{
Cell[42477, 1269, 526, 13, 27, "Input"],
Cell[43006, 1284, 640, 14, 71, "Output"]
}, Open  ]]
}, Closed]],
Cell[CellGroupData[{
Cell[43695, 1304, 98, 1, 26, "Subsection"],
Cell[43796, 1307, 146, 2, 26, "Text"],
Cell[CellGroupData[{
Cell[43967, 1313, 646, 17, 27, "Input"],
Cell[44616, 1332, 937, 23, 60, "Output"]
}, Open  ]]
}, Closed]]
}, Open  ]]
}, Open  ]]
}
]
*)

(* End of internal cache information *)
