Commit e827df9b authored by Reyesrivera, Jose's avatar Reyesrivera, Jose
Browse files

add Hg frame

parent eeaa2906
(* Content-type: application/vnd.wolfram.mathematica *)
(*** Wolfram Notebook File ***)
(* http://www.wolfram.com/nb *)
(* CreatedBy='Mathematica 12.0' *)
(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[ 158, 7]
NotebookDataPosition[ 0, 0]
NotebookDataLength[ 35278, 1004]
NotebookOptionsPosition[ 33159, 967]
NotebookOutlinePosition[ 33500, 982]
CellTagsIndexPosition[ 33457, 979]
NotebookOptionsPosition[ 33001, 960]
NotebookOutlinePosition[ 33342, 975]
CellTagsIndexPosition[ 33299, 972]
WindowFrame->Normal*)
(* Beginning of Notebook Content *)
......@@ -982,30 +975,30 @@ CellTagsIndex->{}
(*NotebookFileOutline
Notebook[{
Cell[CellGroupData[{
Cell[580, 22, 381, 9, 94, "Input",ExpressionUUID->"b04d40e6-4afa-4e93-9c95-4b3cd75fd2b1"],
Cell[964, 33, 250, 4, 24, "Print",ExpressionUUID->"3c515af7-3093-4cbc-9e81-bae9d30a1fff"]
Cell[422, 15, 381, 9, 94, "Input",ExpressionUUID->"b04d40e6-4afa-4e93-9c95-4b3cd75fd2b1"],
Cell[806, 26, 250, 4, 24, "Print",ExpressionUUID->"3c515af7-3093-4cbc-9e81-bae9d30a1fff"]
}, Open ]],
Cell[1229, 40, 1453, 34, 94, "Input",ExpressionUUID->"e841b2dc-8a54-439b-9b13-01d2a06ac36d"],
Cell[2685, 76, 15323, 405, 1223, "Input",ExpressionUUID->"ff42d314-7c2b-4761-aa7e-599c36015f76"],
Cell[1071, 33, 1453, 34, 52, "Input",ExpressionUUID->"e841b2dc-8a54-439b-9b13-01d2a06ac36d"],
Cell[2527, 69, 15323, 405, 1102, "Input",ExpressionUUID->"ff42d314-7c2b-4761-aa7e-599c36015f76"],
Cell[CellGroupData[{
Cell[18033, 485, 1566, 45, 115, "Input",ExpressionUUID->"24c3be5d-915f-467b-a5ed-10fdfc2498c1"],
Cell[19602, 532, 855, 24, 71, "Output",ExpressionUUID->"ac91933d-1868-48fe-adb0-b683004fa38e"],
Cell[20460, 558, 1159, 34, 71, "Output",ExpressionUUID->"ef0276d8-9268-413c-8e6f-92b760328915"],
Cell[21622, 594, 1580, 48, 83, "Output",ExpressionUUID->"af9bb574-4935-430c-aa04-b07707d1bbb2"],
Cell[23205, 644, 1492, 45, 82, "Output",ExpressionUUID->"3ecb1a4a-8ea9-47b7-8d6a-35702d5ec0d4"],
Cell[24700, 691, 1257, 37, 71, "Output",ExpressionUUID->"a69717ad-15df-443f-be7e-64454fd927d8"]
Cell[17875, 478, 1566, 45, 115, "Input",ExpressionUUID->"24c3be5d-915f-467b-a5ed-10fdfc2498c1"],
Cell[19444, 525, 855, 24, 71, "Output",ExpressionUUID->"ac91933d-1868-48fe-adb0-b683004fa38e"],
Cell[20302, 551, 1159, 34, 71, "Output",ExpressionUUID->"ef0276d8-9268-413c-8e6f-92b760328915"],
Cell[21464, 587, 1580, 48, 83, "Output",ExpressionUUID->"af9bb574-4935-430c-aa04-b07707d1bbb2"],
Cell[23047, 637, 1492, 45, 82, "Output",ExpressionUUID->"3ecb1a4a-8ea9-47b7-8d6a-35702d5ec0d4"],
Cell[24542, 684, 1257, 37, 71, "Output",ExpressionUUID->"a69717ad-15df-443f-be7e-64454fd927d8"]
}, Open ]],
Cell[CellGroupData[{
Cell[25994, 733, 1377, 38, 157, "Input",ExpressionUUID->"5147ec4a-4274-4d72-970d-2267574ae770"],
Cell[27374, 773, 349, 9, 34, "Output",ExpressionUUID->"024dbb48-ad76-4857-864e-64ff96e7bfad"],
Cell[27726, 784, 675, 21, 83, "Output",ExpressionUUID->"e7f621bc-0a5e-4677-a7eb-57b0e69a64a2"],
Cell[28404, 807, 1081, 33, 83, "Output",ExpressionUUID->"f055ca76-164d-4617-bf25-ba5ab0582a93"],
Cell[29488, 842, 658, 22, 82, "Output",ExpressionUUID->"e13de4f4-5548-45da-9f95-d2ccee233029"],
Cell[30149, 866, 1140, 35, 91, "Output",ExpressionUUID->"ebabe41b-6bc2-4c9e-9d8e-2002b733a6aa"],
Cell[31292, 903, 571, 19, 79, "Output",ExpressionUUID->"e0ec2326-2ab2-4859-aca5-bc43c919a065"],
Cell[31866, 924, 1118, 35, 83, "Output",ExpressionUUID->"7dc37d65-2816-4074-9698-360aceb06a3a"]
Cell[25836, 726, 1377, 38, 157, "Input",ExpressionUUID->"5147ec4a-4274-4d72-970d-2267574ae770"],
Cell[27216, 766, 349, 9, 34, "Output",ExpressionUUID->"024dbb48-ad76-4857-864e-64ff96e7bfad"],
Cell[27568, 777, 675, 21, 83, "Output",ExpressionUUID->"e7f621bc-0a5e-4677-a7eb-57b0e69a64a2"],
Cell[28246, 800, 1081, 33, 83, "Output",ExpressionUUID->"f055ca76-164d-4617-bf25-ba5ab0582a93"],
Cell[29330, 835, 658, 22, 82, "Output",ExpressionUUID->"e13de4f4-5548-45da-9f95-d2ccee233029"],
Cell[29991, 859, 1140, 35, 91, "Output",ExpressionUUID->"ebabe41b-6bc2-4c9e-9d8e-2002b733a6aa"],
Cell[31134, 896, 571, 19, 79, "Output",ExpressionUUID->"e0ec2326-2ab2-4859-aca5-bc43c919a065"],
Cell[31708, 917, 1118, 35, 83, "Output",ExpressionUUID->"7dc37d65-2816-4074-9698-360aceb06a3a"]
}, Open ]],
Cell[32999, 962, 156, 3, 30, "Input",ExpressionUUID->"a220789d-793f-4ddf-a524-a1c63f94372b"]
Cell[32841, 955, 156, 3, 30, "Input",ExpressionUUID->"a220789d-793f-4ddf-a524-a1c63f94372b"]
}
]
*)
......
(* Content-type: application/vnd.wolfram.mathematica *)
(*** Wolfram Notebook File ***)
(* http://www.wolfram.com/nb *)
(* CreatedBy='Mathematica 12.0' *)
(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[ 158, 7]
NotebookDataLength[ 35278, 1004]
NotebookOptionsPosition[ 33159, 967]
NotebookOutlinePosition[ 33500, 982]
CellTagsIndexPosition[ 33457, 979]
WindowFrame->Normal*)
(* Beginning of Notebook Content *)
Notebook[{
Cell[CellGroupData[{
Cell[BoxData[{
RowBox[{"Needs", "[", "\"\<X`\>\"", "]"}], "\n",
RowBox[{"<<", "helicityvec`"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"SetDirectory", "[",
RowBox[{"NotebookDirectory", "[", "]"}], "]"}],
";"}], "\[IndentingNewLine]",
RowBox[{"NotebookSave", "[", "]"}]}], "Input",
CellLabel->
"In[1278]:=",ExpressionUUID->"b04d40e6-4afa-4e93-9c95-4b3cd75fd2b1"],
Cell[BoxData["\<\"helicityvec by J.G. Reyes v1. August 2019\"\>"], "Print",
CellChangeTimes->{{3.7822305769420233`*^9, 3.782230591171094*^9}},
CellLabel->
"During evaluation of \
In[1278]:=",ExpressionUUID->"3c515af7-3093-4cbc-9e81-bae9d30a1fff"]
}, Open ]],
Cell[BoxData[
RowBox[{
RowBox[{"$Assumptions", "=",
RowBox[{"{",
RowBox[{
RowBox[{"Element", "[",
RowBox[{
RowBox[{"{",
RowBox[{
"ki", ",", "E1", ",", "E2", ",", "k1", ",", "k2", ",", "k3", ",",
"k4", ",", " ", "\[Theta]1", ",", "\[Theta]2", ",", "\[Theta]3", ",",
"\[Theta]4", ",", " ", "\[Phi]1", ",", "\[Phi]2", ",", "\[Phi]3",
",", "\[Phi]4", ",", "p", ",", "MT2", ",", "GS", ",", "EL", ",",
"Alfas", ",", "a1", ",", "S", ",", "T", ",", "U", ",", "MT", ",",
"MH", ",", "MH2", ",", "\[Eta]4", ",", "\[Eta]3", ",", "\[Beta]"}],
"}"}], ",", "Reals"}], "]"}], ",",
RowBox[{"ki", ">", "0"}], " ", ",",
RowBox[{"MT2", ">", "0"}], ",",
RowBox[{"MT", ">", "0"}], ",",
RowBox[{"MH2", ">", "0"}], ",", " ",
RowBox[{"MH", ">", "0"}], ",",
RowBox[{"k4", ">", "0"}], ",",
RowBox[{"p", ">", "0"}], ",",
RowBox[{"k3", ">", "0"}], ",",
RowBox[{"E1", ">", "0"}], ",",
RowBox[{"E2", ">", "0"}], ",",
RowBox[{"k1", ">", "0"}], ",",
RowBox[{"k2", ">", "0"}], ",",
RowBox[{"rS", ">", "0"}], ",",
RowBox[{"kT4", ">", "0"}], ",",
RowBox[{"kT3", ">", "0"}], ",",
RowBox[{"\[Beta]", ">", "0"}]}], "}"}]}], ";"}]], "Input",
CellChangeTimes->{{3.782235612584735*^9, 3.78223563710957*^9}},
CellLabel->
"In[1530]:=",ExpressionUUID->"e841b2dc-8a54-439b-9b13-01d2a06ac36d"],
Cell[BoxData[
RowBox[{
RowBox[{"(*", " ",
RowBox[{
"Defines", " ", "the", " ", "center", " ", "of", " ", "mass", " ",
RowBox[{"frame", "."}]}], " ", "*)"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{
RowBox[{"CMSFrame", "=",
RowBox[{"{", "\[IndentingNewLine]", "\t\t",
RowBox[{
RowBox[{"k1x", "\[Rule]", " ", "0"}], ",",
RowBox[{"k1y", "\[Rule]", " ", "0"}], ",",
RowBox[{"k1z", "\[Rule]", " ",
RowBox[{"rS", "/", "2"}]}], ",", "\[IndentingNewLine]", "\t\t",
RowBox[{"k2x", "\[Rule]", " ", "0"}], ",",
RowBox[{"k2y", "\[Rule]", " ", "0"}], ",",
RowBox[{"k2z", "\[Rule]", " ",
RowBox[{
RowBox[{"-", "rS"}], "/", "2"}]}], ",", "\[IndentingNewLine]",
"\t\t",
RowBox[{"k5x", "\[Rule]", " ",
RowBox[{"-",
RowBox[{"(",
RowBox[{"k4x", "+", "k3x"}], ")"}]}]}], ",",
RowBox[{"k5y", "\[Rule]", " ",
RowBox[{"-",
RowBox[{"(",
RowBox[{"k4y", "+", "k3y"}], ")"}]}]}], ",",
RowBox[{"k5z", "\[Rule]", " ",
RowBox[{"-",
RowBox[{"(",
RowBox[{"k4z", " ", "+", "k3z"}], ")"}]}]}], ",",
"\[IndentingNewLine]", "\t\t",
RowBox[{"P", "->",
RowBox[{
RowBox[{"(",
RowBox[{"1", "/",
RowBox[{"(",
RowBox[{"Sqrt", "[", "S34", "]"}], ")"}]}], ")"}], "*",
RowBox[{"Sqrt", "[",
RowBox[{
RowBox[{"(",
RowBox[{"S34", "-",
RowBox[{"MH", "^", "2"}]}], ")"}], "*",
RowBox[{"(",
RowBox[{"S34", "-",
RowBox[{
RowBox[{"(",
RowBox[{"-", "MH"}], ")"}], "^", "2"}]}], ")"}]}], "]"}]}]}]}],
"}"}]}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]",
RowBox[{
RowBox[{"HiggsGluonFrame", "=",
RowBox[{"{", "\[IndentingNewLine]", "\t\t",
RowBox[{
RowBox[{"k1x", "\[Rule]", "0"}], ",",
RowBox[{"k1y", "\[Rule]", "0"}], ",",
RowBox[{"k1z", "\[Rule]", "E1"}], ",", "\[IndentingNewLine]", "\t\t",
RowBox[{"k2x", "\[Rule]",
RowBox[{"E2", "*",
RowBox[{"Sin", "[", "\[Theta]2", "]"}]}]}], ",",
RowBox[{"k2y", "->", "0"}], ",",
RowBox[{"k2z", "\[Rule]",
RowBox[{"E2", "*",
RowBox[{"Cos", "[", "\[Theta]2", "]"}]}]}], ",",
"\[IndentingNewLine]", "\t\t",
RowBox[{"k3x", "\[Rule]",
RowBox[{
RowBox[{"-", "P"}], "*",
RowBox[{"Sin", "[", "\[Theta]", "]"}], "*",
RowBox[{"Cos", "[", "\[Phi]", "]"}]}]}], ",",
RowBox[{"k3y", "\[Rule]",
RowBox[{
RowBox[{"-", "P"}], "*",
RowBox[{"Sin", "[", "\[Theta]", "]"}], "*",
RowBox[{"Sin", "[", "\[Phi]", "]"}]}]}], ",",
RowBox[{"k3z", "\[Rule]",
RowBox[{
RowBox[{"-", "P"}], "*",
RowBox[{"Cos", "[", "\[Theta]", "]"}]}]}], ",",
"\[IndentingNewLine]", "\t\t",
RowBox[{"k4x", "\[Rule]",
RowBox[{"P", "*",
RowBox[{"Sin", "[", "\[Theta]", "]"}], "*",
RowBox[{"Cos", "[", "\[Phi]", "]"}]}]}], ",",
RowBox[{"k4y", "\[Rule]",
RowBox[{"P", "*",
RowBox[{"Sin", "[", "\[Theta]", "]"}], "*",
RowBox[{"Sin", "[", "\[Phi]", "]"}]}]}], ",",
RowBox[{"k4z", "\[Rule]",
RowBox[{"P", "*",
RowBox[{"Cos", "[", "\[Theta]", "]"}]}]}], ",",
"\[IndentingNewLine]", "\t\t",
RowBox[{"k5x", "\[Rule]",
RowBox[{"E5", "*",
RowBox[{"Sin", "[", "\[Theta]5", "]"}]}]}], ",",
RowBox[{"k5y", "\[Rule]", "0"}], ",",
RowBox[{"k5z", "\[Rule]",
RowBox[{"E5", "*",
RowBox[{"Cos", "[", "\[Theta]5", "]"}]}]}], ",",
"\[IndentingNewLine]", "\t\t",
RowBox[{"P", "->",
RowBox[{
RowBox[{"(",
RowBox[{"1", "/",
RowBox[{"(",
RowBox[{"Sqrt", "[", "S34", "]"}], ")"}]}], ")"}], "*",
RowBox[{"Sqrt", "[",
RowBox[{
RowBox[{"(",
RowBox[{"S34", "-",
RowBox[{"MH", "^", "2"}]}], ")"}], "*",
RowBox[{"(",
RowBox[{"S34", "-",
RowBox[{
RowBox[{"(",
RowBox[{"-", "MH"}], ")"}], "^", "2"}]}], ")"}]}], "]"}]}]}]}],
"}"}]}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]",
RowBox[{"(*", " ",
RowBox[{
RowBox[{"Energies", " ", "in", " ", "terms", " ", "of", " ", "S"}], ",",
"Shat", ",",
RowBox[{"W", "^", "2"}]}], " ", "*)"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"SubEnergies", "=",
RowBox[{"{", "\[IndentingNewLine]", "\t\t",
RowBox[{
RowBox[{"E1", "\[Rule]",
RowBox[{
RowBox[{"(",
RowBox[{"1", "/",
RowBox[{"(",
RowBox[{"2", "*",
RowBox[{"Sqrt", "[", "S34", "]"}]}], ")"}]}], ")"}], "*",
RowBox[{"(",
RowBox[{"S", "-",
RowBox[{"W", "^", "2"}]}], ")"}]}]}], ",", "\[IndentingNewLine]",
"\t\t",
RowBox[{"E2", "\[Rule]", " ",
RowBox[{
RowBox[{"(",
RowBox[{"1", "/",
RowBox[{"(",
RowBox[{"2", "*",
RowBox[{"Sqrt", "[", "S34", "]"}]}], ")"}]}], ")"}], "*",
RowBox[{"(",
RowBox[{"S", "+",
RowBox[{"W", "^", "2"}]}], ")"}]}]}], ",", "\[IndentingNewLine]",
"\t\t",
RowBox[{"E5", "\[Rule]",
RowBox[{
RowBox[{"(",
RowBox[{"1", "/",
RowBox[{"(",
RowBox[{"2", "*",
RowBox[{"Sqrt", "[", "S34", "]"}]}], ")"}]}], ")"}], "*", "S",
"*",
RowBox[{"(",
RowBox[{"1", "-", "x"}], ")"}]}]}]}], "}"}]}], ";"}],
"\[IndentingNewLine]", "\[IndentingNewLine]",
RowBox[{"(*", " ",
RowBox[{"Returns", " ", "vector", " ",
RowBox[{"(",
RowBox[{"E", ",", "px", ",", "py", ",", "pz"}], ")"}], " ", "in", " ",
"Euclidean", " ", "coordinates"}], " ", "*)"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{
RowBox[{"vec4xyz", "[",
RowBox[{"m_", ",", "kx_", ",", "ky_", ",", "kz_"}], "]"}], ":=",
RowBox[{"Module", "[",
RowBox[{
RowBox[{"{",
RowBox[{"k", "=",
RowBox[{"Sqrt", "[",
RowBox[{
RowBox[{"kx", "^", "2"}], " ", "+",
RowBox[{"ky", "^", "2"}], " ", "+",
RowBox[{"kz", "^", "2"}]}], "]"}]}], "}"}], ",",
RowBox[{"{",
RowBox[{
RowBox[{"Sqrt", "[",
RowBox[{
RowBox[{"k", "^", "2"}], " ", "+", " ",
RowBox[{"m", "^", "2"}]}], "]"}], ",", "kx", ",", "ky", ",",
"kz"}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]",
"\[IndentingNewLine]",
RowBox[{"(*", " ",
RowBox[{"Returns", " ", "vector", " ",
RowBox[{"(",
RowBox[{"E", ",", "px", ",", "py", ",", "pz"}], ")"}], " ", "in", " ",
"Collider", " ", "coordinates"}], " ", "*)"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{
RowBox[{"vec4Tep", "[",
RowBox[{"m_", ",", "kT_", ",", "\[Eta]_", ",", "\[Phi]_"}], "]"}], ":=",
RowBox[{"Block", "[",
RowBox[{
RowBox[{"{",
RowBox[{
RowBox[{"k", "=",
RowBox[{"kT", " ",
RowBox[{"Cosh", "[", "\[Eta]", "]"}]}]}], ",",
RowBox[{"kx", "=", " ",
RowBox[{"kT", " ",
RowBox[{"Cos", "[", "\[Phi]", "]"}]}]}], ",",
RowBox[{"ky", "=", " ",
RowBox[{"kT", " ",
RowBox[{"Sin", "[", "\[Phi]", "]"}]}]}], ",",
RowBox[{"kz", " ", "=", " ",
RowBox[{"kT", " ",
RowBox[{"Sinh", "[", "\[Eta]", "]"}]}]}]}], "}"}], ",",
RowBox[{"{",
RowBox[{
RowBox[{"Sqrt", "[",
RowBox[{
RowBox[{"k", "^", "2"}], " ", "+", " ",
RowBox[{"m", "^", "2"}]}], "]"}], ",", "kx", ",", "ky", ",",
"kz"}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]",
"\[IndentingNewLine]",
RowBox[{"(*", " ",
RowBox[{"Returns", " ", "vector", " ",
RowBox[{"(",
RowBox[{"E", ",", "px", ",", "py", ",", "pz"}], ")"}], " ", "in", " ",
"Spherical", " ", "coordinates"}], " ", "*)"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{
RowBox[{"vec4ptp", "[",
RowBox[{"m_", ",", "p_", ",", "\[Theta]_", ",", "\[Phi]_"}], "]"}], ":=",
"\[IndentingNewLine]",
RowBox[{"Block", "[",
RowBox[{
RowBox[{"{",
RowBox[{
RowBox[{"k", "=", "p"}], ",",
RowBox[{"kx", "=", " ",
RowBox[{"p", "*",
RowBox[{"Sin", "[", "\[Theta]", "]"}], "*",
RowBox[{"Cos", "[", "\[Phi]", "]"}]}]}], ",",
RowBox[{"ky", "=", " ",
RowBox[{"p", "*",
RowBox[{"Sin", "[", "\[Phi]", "]"}], "*",
RowBox[{"Sin", "[", "\[Phi]", "]"}]}]}], ",",
RowBox[{"kz", " ", "=", " ",
RowBox[{"p", "*",
RowBox[{"Cos", "[", "\[Theta]", "]"}]}]}]}], "}"}], ",",
RowBox[{"{",
RowBox[{
RowBox[{"Sqrt", "[",
RowBox[{
RowBox[{"k", "^", "2"}], " ", "+", " ",
RowBox[{"m", "^", "2"}]}], "]"}], ",", "kx", ",", "ky", ",",
"kz"}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]",
"\[IndentingNewLine]",
RowBox[{"(*", " ",
RowBox[{
"Definition", " ", "of", " ", "the", " ", "four", " ", "vector", " ",
"of", " ", "particle", " ", "3", " ",
RowBox[{"(", "Higgs", ")"}]}], " ", "*)"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"vecetaphi3", "=",
RowBox[{
RowBox[{"vec4Tep", "[",
RowBox[{"MH", ",", "kT3", ",", "\[Eta]3", ",", "\[Phi]3"}], "]"}], "//",
"Simplify"}]}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]",
RowBox[{"(*", " ",
RowBox[{
"Definition", " ", "of", " ", "the", " ", "four", " ", "vector", " ",
"of", " ", "particle", " ", "4", " ",
RowBox[{"(",
RowBox[{"outgoing", " ", "gluon"}], ")"}]}], " ", "*)"}],
"\[IndentingNewLine]",
RowBox[{
RowBox[{"vecetaphi4", "=",
RowBox[{"(",
RowBox[{
RowBox[{"vec4Tep", "[",
RowBox[{"0", ",", "kT4", ",", "\[Eta]4", ",", "\[Phi]4"}], "]"}], "//",
"Simplify"}], ")"}]}], ";"}], "\[IndentingNewLine]",
"\[IndentingNewLine]",
RowBox[{"(*", " ",
RowBox[{
RowBox[{
"Transformation", " ", "of", " ", "momentum", " ", "variables", " ",
"from", " ", "Euclidean", " ", "to", " ", "Collider", " ",
"coordinates", " ", "for", " ", "particles", " ", "3"}], ",", "4."}],
" ", "*)"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"vecsubetaphi", "=",
RowBox[{"{",
RowBox[{
RowBox[{"k4x", "\[Rule]", " ",
RowBox[{"vecetaphi4", "[",
RowBox[{"[", "2", "]"}], "]"}]}], ",",
RowBox[{"k4y", "\[Rule]", " ",
RowBox[{"vecetaphi4", "[",
RowBox[{"[", "3", "]"}], "]"}]}], ",",
RowBox[{"k4z", "\[Rule]", " ",
RowBox[{"vecetaphi4", "[",
RowBox[{"[", "4", "]"}], "]"}]}], ",",
RowBox[{"k3x", "\[Rule]", " ",
RowBox[{"vecetaphi3", "[",
RowBox[{"[", "2", "]"}], "]"}]}], ",",
RowBox[{"k3y", "\[Rule]", " ",
RowBox[{"vecetaphi3", "[",
RowBox[{"[", "3", "]"}], "]"}]}], ",",
RowBox[{"k3z", "\[Rule]", " ",
RowBox[{"vecetaphi3", "[",
RowBox[{"[", "4", "]"}], "]"}]}]}], "}"}]}], ";"}],
"\[IndentingNewLine]", "\[IndentingNewLine]",
RowBox[{"(*", " ",
RowBox[{
"Four", " ", "momentum", " ", "vectors", " ", "for", " ", "all", " ",
RowBox[{"particles", "."}]}], " ", "*)"}], "\[IndentingNewLine]",
RowBox[{"(*", " ",
RowBox[{
RowBox[{
RowBox[{"g",
RowBox[{"(", "k1", ")"}], " ", "g",
RowBox[{"(", "k2", ")"}]}], " ", "\[Rule]", " ",
RowBox[{"h",
RowBox[{"(", "k3", ")"}]}]}], ",",
RowBox[{"g",
RowBox[{"(", "k4", ")"}]}], ",",
RowBox[{"g",
RowBox[{"(", "k5", ")"}]}]}], " ", "*)"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"k1vec", "=",
RowBox[{
RowBox[{"vec4xyz", "[",
RowBox[{"0", ",", "k1x", ",", "k1y", ",", "k1z"}], "]"}], "//",
"Simplify"}]}], ";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"k2vec", "=",
RowBox[{
RowBox[{"vec4xyz", "[",
RowBox[{"0", ",", "k2x", ",", "k2y", ",", "k2z"}], "]"}], "//",
"Simplify"}]}], ";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"k3vec", "=",
RowBox[{
RowBox[{"vec4xyz", "[",
RowBox[{"MH", ",", "k3x", ",", "k3y", ",", "k3z"}], "]"}], "//",
"Simplify"}]}], ";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"k4vec", "=",
RowBox[{
RowBox[{"vec4xyz", "[",
RowBox[{"0", ",", "k4x", ",", "k4y", ",", "k4z"}], "]"}], "//",
"Simplify"}]}], ";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"k5vec", "=",
RowBox[{
RowBox[{"vec4xyz", "[",
RowBox[{"0", ",", "k5x", ",", "k5y", ",", "k5z"}], "]"}], "//",
"Simplify"}]}], ";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"kvec1", "=",
RowBox[{
RowBox[{
RowBox[{"k1vec", "/.", " ", "vecsubetaphi"}], "/.",
RowBox[{"rS", "\[Rule]",
RowBox[{"Sqrt", "[", "S", "]"}]}]}], "//", "Simplify"}]}], ";"}],
"\[IndentingNewLine]",
RowBox[{
RowBox[{"kvec2", "=",
RowBox[{
RowBox[{
RowBox[{"k2vec", "/.", " ", "vecsubetaphi"}], "/.",
RowBox[{"rS", "\[Rule]",
RowBox[{"Sqrt", "[", "S", "]"}]}]}], "//", "Simplify"}]}], ";"}],
"\[IndentingNewLine]",
RowBox[{
RowBox[{"kvec3", "=",
RowBox[{
RowBox[{
RowBox[{"k3vec", "/.", " ", "vecsubetaphi"}], "/.",
RowBox[{"rS", "\[Rule]",
RowBox[{"Sqrt", "[", "S", "]"}]}]}], "//", "Simplify"}]}], ";"}],
"\[IndentingNewLine]",
RowBox[{
RowBox[{"kvec4", "=",
RowBox[{
RowBox[{
RowBox[{"k4vec", "/.", " ", "vecsubetaphi"}], "/.",
RowBox[{"rS", "\[Rule]",
RowBox[{"Sqrt", "[", "S", "]"}]}]}], "//", "Simplify"}]}], ";"}],
"\[IndentingNewLine]",
RowBox[{
RowBox[{"kvec5", "=",
RowBox[{
RowBox[{
RowBox[{"k5vec", "/.", " ", "vecsubetaphi"}], "/.",
RowBox[{"rS", "\[Rule]",
RowBox[{"Sqrt", "[", "S", "]"}]}]}], "//", "Simplify"}]}],
";"}]}]}]], "Input",
CellChangeTimes->{{3.78223063866676*^9, 3.782230646358025*^9}, {
3.7822306976524982`*^9, 3.782230735632852*^9}, {3.782230816029155*^9,
3.78223083197003*^9}, {3.78223086668169*^9, 3.782231008566526*^9}, {
3.78223109826705*^9, 3.7822311853511333`*^9}, {3.782231249051858*^9,
3.7822312498812933`*^9}, {3.782231532373032*^9, 3.782231875765934*^9}, {
3.78223197509763*^9, 3.782232118677737*^9}, {3.7822330868501797`*^9,
3.7822330905511*^9}, {3.782235895252698*^9, 3.7822359017283*^9}, {
3.78240263025075*^9, 3.782402789080329*^9}, 3.7824029254166937`*^9, {
3.782403181768487*^9, 3.78240318846362*^9}},
CellLabel->
"In[1623]:=",ExpressionUUID->"ff42d314-7c2b-4761-aa7e-599c36015f76"],
Cell[CellGroupData[{
Cell[BoxData[{
RowBox[{"k1vec", "=",
RowBox[{
RowBox[{
RowBox[{
RowBox[{"vec4xyz", "[",
RowBox[{"0", ",", "k1x", ",", "k1y", ",", "k1z"}], "]"}], "//.",
"HiggsGluonFrame"}], "//.", "SubEnergies"}], "//",
"Simplify"}]}], "\[IndentingNewLine]",
RowBox[{"k2vec", "=",
RowBox[{
RowBox[{
RowBox[{
RowBox[{"vec4xyz", "[",
RowBox[{"0", ",", "k2x", ",", "k2y", ",", "k2z"}], "]"}], "//.",
"HiggsGluonFrame"}], "//.", "SubEnergies"}], "//",
"Simplify"}]}], "\[IndentingNewLine]",
RowBox[{"k3vec", "=",
RowBox[{
RowBox[{
RowBox[{
RowBox[{"vec4xyz", "[",
RowBox[{"MH", ",", "k3x", ",", "k3y", ",", "k3z"}], "]"}], "//.",
"HiggsGluonFrame"}], "//.", "SubEnergies"}], "//",
"Simplify"}]}], "\[IndentingNewLine]",
RowBox[{"k4vec", "=",
RowBox[{
RowBox[{
RowBox[{
RowBox[{"vec4xyz", "[",
RowBox[{"0", ",", "k4x", ",", "k4y", ",", "k4z"}], "]"}], "//.",
"HiggsGluonFrame"}], "//.", "SubEnergies"}], "//",
"Simplify"}]}], "\[IndentingNewLine]",
RowBox[{"k5vec", "=",
RowBox[{
RowBox[{
RowBox[{
RowBox[{"vec4xyz", "[",
RowBox[{"0", ",", "k5x", ",", "k5y", ",", "k5z"}], "]"}], "//.",
"HiggsGluonFrame"}], "//.", "SubEnergies"}], "//",
"Simplify"}]}]}], "Input",
CellChangeTimes->{{3.782230653325596*^9, 3.7822306553843822`*^9}, {
3.7822319188263807`*^9, 3.782231929610517*^9}, {3.7824029339221983`*^9,
3.782402948629674*^9}},
CellLabel->
"In[1642]:=",ExpressionUUID->"24c3be5d-915f-467b-a5ed-10fdfc2498c1"],
Cell[BoxData[
RowBox[{"{",
RowBox[{
RowBox[{
FractionBox["1", "2"], " ",
SqrtBox[
FractionBox[
SuperscriptBox[
RowBox[{"(",
RowBox[{"S", "-",
SuperscriptBox["W", "2"]}], ")"}], "2"], "S34"]]}], ",", "0", ",",
"0", ",",
FractionBox[
RowBox[{"S", "-",
SuperscriptBox["W", "2"]}],
RowBox[{"2", " ",
SqrtBox["S34"]}]]}], "}"}]], "Output",
CellChangeTimes->{
3.7822306561860027`*^9, 3.782231012372971*^9, 3.782231191940976*^9,
3.782231930205492*^9, 3.782232021116662*^9, {3.782232072247366*^9,
3.782232120997113*^9}, 3.78223309583575*^9, 3.78223464180264*^9,
3.782235909183299*^9, 3.782235972840762*^9, {3.7824029295974216`*^9,
3.7824029494431334`*^9}, 3.782403192579413*^9},
CellLabel->
"Out[1642]=",ExpressionUUID->"ac91933d-1868-48fe-adb0-b683004fa38e"],
Cell[BoxData[
RowBox[{"{",
RowBox[{
RowBox[{
FractionBox["1", "2"], " ",
SqrtBox[
FractionBox[
SuperscriptBox[
RowBox[{"(",
RowBox[{"S", "+",
SuperscriptBox["W", "2"]}], ")"}], "2"], "S34"]]}], ",",
FractionBox[
RowBox[{
RowBox[{"(",
RowBox[{"S", "+",
SuperscriptBox["W", "2"]}], ")"}], " ",
RowBox[{"Sin", "[", "\[Theta]2", "]"}]}],
RowBox[{"2", " ",
SqrtBox["S34"]}]], ",", "0", ",",
FractionBox[
RowBox[{
RowBox[{"(",