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

add logs for 2 jets triangles

parent a89065a8
No preview for this file type
(* Created with the Wolfram Language for Students - Personal Use Only : www.wolfram.com *)
{}
(* Created with the Wolfram Language for Students - Personal Use Only : www.wolfram.com *)
{Log[(-MH^2 + 2*MT^2 + MH*Sqrt[MH^2 - 4*MT^2])/(2*MT^2)], Log[MT],
Log[Mu^2/MT^2], Log[Mu^2]}
(* Created with the Wolfram Language for Students - Personal Use Only : www.wolfram.com *)
{}
(* Created with the Wolfram Language for Students - Personal Use Only : www.wolfram.com *)
{}
(* Created with the Wolfram Language for Students - Personal Use Only : www.wolfram.com *)
{}
(* Created with the Wolfram Language for Students - Personal Use Only : www.wolfram.com *)
{Log[(-MH^2 + 2*MT^2 + MH*Sqrt[MH^2 - 4*MT^2])/(2*MT^2)],
Log[(2*MT^2 - S + Sqrt[S*(-4*MT^2 + S)])/(2*MT^2)],
Log[(2*MT^2 - T + Sqrt[T*(-4*MT^2 + T)])/(2*MT^2)],
Log[(2*MT^2 - U + Sqrt[U*(-4*MT^2 + U)])/(2*MT^2)], Log[Mu^2/MT^2]}
(* Created with the Wolfram Language for Students - Personal Use Only : www.wolfram.com *)
{}
(* Created with the Wolfram Language for Students - Personal Use Only : www.wolfram.com *)
{ScalarC0[0, 0, S, MT, MT, MT], ScalarC0[0, 0, T, MT, MT, MT],
ScalarC0[0, 0, U, MT, MT, MT], ScalarC0[0, MH^2, S, MT, MT, MT],
ScalarC0[0, MH^2, T, MT, MT, MT], ScalarC0[0, MH^2, U, MT, MT, MT],
ScalarD0[0, 0, 0, MH^2, S, T, MT, MT, MT, MT],
ScalarD0[0, 0, 0, MH^2, S, U, MT, MT, MT, MT],
ScalarD0[0, 0, 0, MH^2, T, U, MT, MT, MT, MT]}
(* Created with the Wolfram Language for Students - Personal Use Only : www.wolfram.com *)
{}
(* Created with the Wolfram Language for Students - Personal Use Only : www.wolfram.com *)
{Log[(-2*MH^2 + 2*MT^2 + S34 + T + U +
Sqrt[(-2*MH^2 + S34 + T + U)*(-2*MH^2 + 4*MT^2 + S34 + T + U)])/
(2*MT^2)], Log[(-MH^2 + 2*MT^2 + MH*Sqrt[MH^2 - 4*MT^2])/(2*MT^2)],
Log[(2*MT^2 - S34 + Sqrt[S34*(-4*MT^2 + S34)])/(2*MT^2)],
Log[(2*MT^2 - T + Sqrt[T*(-4*MT^2 + T)])/(2*MT^2)],
Log[(2*MT^2 - U + Sqrt[U*(-4*MT^2 + U)])/(2*MT^2)], Log[Mu^2/MT^2]}
Log[(-2*MH^2+2*MT^2+S34+T+U+Sqrt[(-2*MH^2+S34+T+U)*(-2*MH^2+4*MT^2+S34+T+U)])/(2*MT^2)]
Log[(-MH^2+2*MT^2+MH*Sqrt[MH^2-4*MT^2])/(2*MT^2)]
Log[(2*MT^2-S34+Sqrt[S34*(-4*MT^2+S34)])/(2*MT^2)]
Log[(2*MT^2-T+Sqrt[T*(-4*MT^2+T)])/(2*MT^2)]
Log[(2*MT^2-U+Sqrt[U*(-4*MT^2+U)])/(2*MT^2)]
Log[Mu^2/MT^2]
(* Created with the Wolfram Language for Students - Personal Use Only : www.wolfram.com *)
{}
(* Created with the Wolfram Language for Students - Personal Use Only : www.wolfram.com *)
{ScalarC0[MH^2, MH^2 - S - T - T14, T24, MT, MT, MT],
ScalarC0[MH^2, S, MH^2 - S34 - T14 - T24, MT, MT, MT],
ScalarC0[MH^2, T14, MH^2 - S - T24 - U, MT, MT, MT]}
ScalarC0[MH^2,MH^2-S-T-T14,T24,MT,MT,MT]
ScalarC0[MH^2,S,MH^2-S34-T14-T24,MT,MT,MT]
ScalarC0[MH^2,T14,MH^2-S-T24-U,MT,MT,MT]
No preview for this file type
(* Content-type: application/vnd.wolfram.mathematica *)
(*** Wolfram Notebook File ***)
(* http://www.wolfram.com/nb *)
(* CreatedBy='Mathematica 11.3' *)
(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[ 158, 7]
NotebookDataLength[ 13868, 362]
NotebookOptionsPosition[ 13055, 341]
NotebookOutlinePosition[ 13411, 357]
CellTagsIndexPosition[ 13368, 354]
WindowFrame->Normal*)
(* Beginning of Notebook Content *)
Notebook[{
Cell[BoxData[
RowBox[{"(*", " ",
RowBox[{
RowBox[{
RowBox[{
RowBox[{"LogsAnalysis", ".", "nb"}], " ", "\[IndentingNewLine]", " ",
RowBox[{"Process", ":", " ",
RowBox[{"g", " ", "+", " ", "g"}]}]}], " ", "\[Rule]", " ",
RowBox[{"H", " ", "+", " ", "njets"}]}], ",", " ",
RowBox[{"n", " ", "\[LessEqual]", " ",
RowBox[{"2", "\[IndentingNewLine]", " ",
RowBox[{"Model", ":", " ", "SMQCD"}]}]}], ",", " ",
RowBox[{
"Last", " ", "Modified", " ", "October", " ", "2019",
"\[IndentingNewLine]", "Created", " ",
RowBox[{"by", ":", " ",
RowBox[{
RowBox[{"J", ".", "G", ".", "Reyes"}], " ", "Rivera"}]}]}]}], " ",
"*)"}]], "Input",
CellChangeTimes->{{3.7483460640690107`*^9, 3.748346128736135*^9}, {
3.7505322064466867`*^9, 3.7505322096193037`*^9}, {3.750685818525031*^9,
3.7506858187985888`*^9}, {3.750685862467576*^9, 3.750685862527199*^9}, {
3.7513081531132383`*^9, 3.751308157253491*^9}, {3.761295781705222*^9,
3.761295787067747*^9}, {3.7640060179701233`*^9, 3.7640060182815123`*^9}, {
3.769416168255188*^9, 3.769416171052967*^9}, {3.7791271781333838`*^9,
3.779127187089452*^9}, {3.779550031652034*^9, 3.779550065586117*^9}, {
3.779811778477392*^9,
3.779811787712831*^9}},ExpressionUUID->"a9eb43f9-c437-4ff0-ae8b-\
eb31f4d06b38"],
Cell[CellGroupData[{
Cell[BoxData[{
RowBox[{"Needs", "[", "\"\<X`\>\"", "]"}], "\n",
RowBox[{
RowBox[{"SetDirectory", "[",
RowBox[{"NotebookDirectory", "[", "]"}], "]"}],
";"}], "\[IndentingNewLine]",
RowBox[{"NotebookSave", "[", "]"}]}], "Input",
CellChangeTimes->{{3.7694157567375174`*^9, 3.7694157651420507`*^9}, {
3.7694161379556293`*^9, 3.769416138416831*^9}, {3.7694161826886263`*^9,
3.769416192145475*^9}, {3.769416241287668*^9, 3.769416242139756*^9}, {
3.769421884524377*^9, 3.769421889789727*^9}, {3.7798117921235847`*^9,
3.779811792883956*^9}},
CellLabel->"In[1]:=",ExpressionUUID->"aa6e76d2-ec10-49c0-b6e0-1dd4d91b1a18"],
Cell[BoxData["\<\"\\!\\(\\*TemplateBox[List[\\\"\\\\\\\"Package-X v2.1.1, by \
Hiren H. Patel\\\\\\\\nFor more information, see the \\\\\\\"\\\", \
TemplateBox[List[\\\"\\\\\\\"guide\\\\\\\"\\\", \\\"paclet:X/guide/PackageX\\\
\"], \\\"HyperlinkPaclet\\\"]], \\\"RowDefault\\\"]\\)\"\>"], "Print",
CellChangeTimes->{3.779812578591907*^9},
CellLabel->
"During evaluation of \
In[1]:=",ExpressionUUID->"f8f2db75-7118-45a9-adb1-7fa65a7be2f2"]
}, Open ]],
Cell[BoxData[
RowBox[{
RowBox[{"$Assumptions", "=",
RowBox[{"{",
RowBox[{
RowBox[{"Element", "[",
RowBox[{
RowBox[{"{",
RowBox[{
"ki", ",", "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[{"k1", ">", "0"}], ",",
RowBox[{"k2", ">", "0"}], ",",
RowBox[{"rS", ">", "0"}], ",",
RowBox[{"kT4", ">", "0"}], ",",
RowBox[{"kT3", ">", "0"}], ",",
RowBox[{"\[Beta]", ">", "0"}]}], "}"}]}], ";"}]], "Input",
CellChangeTimes->{
3.7767696725954723`*^9, {3.7795525243916817`*^9, 3.7795525271701508`*^9}},
CellLabel->"In[4]:=",ExpressionUUID->"1c12ba0e-fd6c-4b7e-8de6-2b784fa72d92"],
Cell[BoxData[
RowBox[{"\[IndentingNewLine]",
RowBox[{"(*", " ",
RowBox[{
RowBox[{"Here", " ", "the", " ", "program", " ",
RowBox[{"logs", ".", "py"}]}], ",", " ",
RowBox[{"dilogs", ".", "py"}], ",", " ",
RowBox[{"polylogs", ".", "py"}], ",", " ",
RowBox[{
RowBox[{"scalar", ".", "py"}], " ", "is", " ", "used", " ", "to", " ",
"produce", " ", "the", " ", "lists", " ", "of", " ", "all", " ", "log",
" ", "containing", " ", "terms"}]}], " ", "*)"}], "\[IndentingNewLine]",
RowBox[{"(*", " ", "\[IndentingNewLine]", " ",
RowBox[{"Do", ":", "\[IndentingNewLine]", " ",
RowBox[{
RowBox[{"python", " ",
RowBox[{"logs", ".", "py"}], " ", "\"\<../feynamps/ggH_PV.m\>\""}], " ",
"|", " ",
RowBox[{
RowBox[{
RowBox[{
RowBox[{
RowBox[{
RowBox[{"sort", " ", "-", "u"}], " ", ">"}], " ", ".."}], "/",
"logs"}], "/", "0"}],
RowBox[{"jets", "/",
RowBox[{"logs", ".", "txt"}]}], "\n", " ", "python", " ",
RowBox[{"logs", ".", "py"}], " ", "\"\<../feynamps/ggHg_PV.m\>\""}],
" ", "|", " ",
RowBox[{
RowBox[{
RowBox[{
RowBox[{
RowBox[{
RowBox[{"sort", " ", "-", "u"}], " ", ">"}], " ", ".."}], "/",
"logs"}], "/", "1"}],
RowBox[{"jets", "/",
RowBox[{"logs", ".", "txt"}]}], "\[IndentingNewLine]", " ",
"python", " ",
RowBox[{"logs", ".", "py"}], " ", "\"\<../feynamps/ggHgg_T_PV.m\>\""}],
" ", "|", " ",
RowBox[{
RowBox[{
RowBox[{
RowBox[{
RowBox[{
RowBox[{"sort", " ", "-", "u"}], " ", ">"}], " ", ".."}], "/",
"logs"}], "/", "2"}],
RowBox[{"jets", "/",
RowBox[{"logs", ".", "txt"}]}]}]}]}], "\[IndentingNewLine]",
"*)"}]}]], "Input",
CellChangeTimes->{{3.769951408543744*^9, 3.769951516031106*^9}, {
3.770009459536234*^9, 3.770009522564651*^9}, {3.775826780721582*^9,
3.775826781197331*^9}, {3.778510819803602*^9, 3.778510820329376*^9}, {
3.7791280023926783`*^9, 3.779128054274996*^9}, {3.779202869762745*^9,
3.779202874883494*^9}, {3.7798118395001507`*^9, 3.77981193775534*^9}, {
3.779812260719655*^9, 3.779812263113318*^9}, {3.779812529369989*^9,
3.779812563752565*^9}},ExpressionUUID->"476a37c6-482c-4f61-9ced-\
5d80a625711e"],
Cell[BoxData[
RowBox[{
RowBox[{"(*", " ",
RowBox[{"H", " ", "+", " ",
RowBox[{"0", " ", "jets"}]}], " ", "*)"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{
RowBox[{"SetDirectory", "[",
RowBox[{"NotebookDirectory", "[", "]"}], "]"}], ";"}],
"\[IndentingNewLine]",
RowBox[{
RowBox[{"SetDirectory", "[",
RowBox[{"StringJoin", "[",
RowBox[{
RowBox[{"ParentDirectory", "[", "]"}], ",", "\"\</logs/0jets\>\""}],
"]"}], "]"}], ";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"logs0", "=",
RowBox[{"ToExpression", "[",
RowBox[{"Import", "[",
RowBox[{"\"\<logs.txt\>\"", ",", "\"\<Lines\>\""}], "]"}], "]"}]}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"dilogs0", "=",
RowBox[{"ToExpression", "[",
RowBox[{"Import", "[",
RowBox[{"\"\<dilogs.txt\>\"", ",", "\"\<Lines\>\""}], "]"}], "]"}]}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"polylogs0", "=",
RowBox[{"ToExpression", "[",
RowBox[{"Import", "[",
RowBox[{"\"\<polylogs.txt\>\"", ",", "\"\<Lines\>\""}], "]"}], "]"}]}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"scalar0", "=",
RowBox[{"ToExpression", "[",
RowBox[{"Import", "[",
RowBox[{"\"\<scalar.txt\>\"", ",", "\"\<Lines\>\""}], "]"}], "]"}]}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"Export", "[",
RowBox[{"\"\<logs.m\>\"", ",", "logs0"}], "]"}], ";"}],
"\[IndentingNewLine]",
RowBox[{
RowBox[{"Export", "[",
RowBox[{"\"\<dilogs.m\>\"", ",", "dilogs0"}], "]"}], ";"}],
"\[IndentingNewLine]",
RowBox[{
RowBox[{"Export", "[",
RowBox[{"\"\<polylogs.m\>\"", ",", "polylogs0"}], "]"}], ";"}],
"\[IndentingNewLine]",
RowBox[{
RowBox[{"Export", "[",
RowBox[{"\"\<scalar.m\>\"", ",", "scalar0"}], "]"}], ";"}],
"\[IndentingNewLine]", "\[IndentingNewLine]",
RowBox[{"(*", " ",
RowBox[{"H", " ", "+", " ",
RowBox[{"1", " ", "jets"}]}], " ", "*)"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"SetDirectory", "[",
RowBox[{"NotebookDirectory", "[", "]"}], "]"}], ";"}],
"\[IndentingNewLine]",
RowBox[{
RowBox[{"SetDirectory", "[",
RowBox[{"StringJoin", "[",
RowBox[{
RowBox[{"ParentDirectory", "[", "]"}], ",", "\"\</logs/1jets\>\""}],
"]"}], "]"}], ";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"logs1", "=",
RowBox[{"ToExpression", "[",
RowBox[{"Import", "[",
RowBox[{"\"\<logs.txt\>\"", ",", "\"\<Lines\>\""}], "]"}], "]"}]}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"dilogs1", "=",
RowBox[{"ToExpression", "[",
RowBox[{"Import", "[",
RowBox[{"\"\<dilogs.txt\>\"", ",", "\"\<Lines\>\""}], "]"}], "]"}]}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"polylogs1", "=",
RowBox[{"ToExpression", "[",
RowBox[{"Import", "[",
RowBox[{"\"\<polylogs.txt\>\"", ",", "\"\<Lines\>\""}], "]"}], "]"}]}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"scalar1", "=",
RowBox[{"ToExpression", "[",
RowBox[{"Import", "[",
RowBox[{"\"\<scalar.txt\>\"", ",", "\"\<Lines\>\""}], "]"}], "]"}]}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"Export", "[",
RowBox[{"\"\<logs.m\>\"", ",", "logs1"}], "]"}], ";"}],
"\[IndentingNewLine]",
RowBox[{
RowBox[{"Export", "[",
RowBox[{"\"\<dilogs.m\>\"", ",", "dilogs1"}], "]"}], ";"}],
"\[IndentingNewLine]",
RowBox[{
RowBox[{"Export", "[",
RowBox[{"\"\<polylogs.m\>\"", ",", "polylogs1"}], "]"}], ";"}],
"\[IndentingNewLine]",
RowBox[{
RowBox[{"Export", "[",
RowBox[{"\"\<scalar.m\>\"", ",", "scalar1"}], "]"}], ";"}],
"\[IndentingNewLine]", "\[IndentingNewLine]",
RowBox[{"(*", " ",
RowBox[{"H", " ", "+", " ",
RowBox[{"2", " ", "jets"}]}], " ", "*)"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"SetDirectory", "[",
RowBox[{"NotebookDirectory", "[", "]"}], "]"}], ";"}],
"\[IndentingNewLine]",
RowBox[{
RowBox[{"SetDirectory", "[",
RowBox[{"StringJoin", "[",
RowBox[{
RowBox[{"ParentDirectory", "[", "]"}], ",", "\"\</logs/2jets\>\""}],
"]"}], "]"}], ";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"logs2", "=",
RowBox[{"ToExpression", "[",
RowBox[{"Import", "[",
RowBox[{"\"\<logs.txt\>\"", ",", "\"\<Lines\>\""}], "]"}], "]"}]}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"dilogs2", "=",
RowBox[{"ToExpression", "[",
RowBox[{"Import", "[",
RowBox[{"\"\<dilogs.txt\>\"", ",", "\"\<Lines\>\""}], "]"}], "]"}]}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"polylogs2", "=",
RowBox[{"ToExpression", "[",
RowBox[{"Import", "[",
RowBox[{"\"\<polylogs.txt\>\"", ",", "\"\<Lines\>\""}], "]"}], "]"}]}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"scalar2", "=",
RowBox[{"ToExpression", "[",
RowBox[{"Import", "[",
RowBox[{"\"\<scalar.txt\>\"", ",", "\"\<Lines\>\""}], "]"}], "]"}]}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"Export", "[",
RowBox[{"\"\<logs.m\>\"", ",", "logs2"}], "]"}], ";"}],
"\[IndentingNewLine]",
RowBox[{
RowBox[{"Export", "[",
RowBox[{"\"\<dilogs.m\>\"", ",", "dilogs2"}], "]"}], ";"}],
"\[IndentingNewLine]",
RowBox[{
RowBox[{"Export", "[",
RowBox[{"\"\<polylogs.m\>\"", ",", "polylogs2"}], "]"}], ";"}],
"\[IndentingNewLine]",
RowBox[{
RowBox[{"Export", "[",
RowBox[{"\"\<scalar.m\>\"", ",", "scalar2"}], "]"}], ";"}],
"\[IndentingNewLine]"}]}]], "Input",
CellChangeTimes->{{3.779550546825959*^9, 3.779550604971881*^9}, {
3.7796243610532722`*^9, 3.7796243844412127`*^9}, {3.779624469152783*^9,
3.779624472292886*^9}, {3.779624508680273*^9, 3.779624508897086*^9}, {
3.779812156719249*^9, 3.779812164746841*^9}, {3.779812231784854*^9,
3.779812238364973*^9}, 3.7798122716966133`*^9, {3.779812466296836*^9,
3.779812510809091*^9}, {3.779812589963299*^9, 3.7798125922402687`*^9}},
CellLabel->"In[35]:=",ExpressionUUID->"dfd6f63d-512e-4dca-8e27-b6dc7e095594"]
},
WindowSize->{900, 1034},
WindowMargins->{{Automatic, -900}, {Automatic, 0}},
FrontEndVersion->"11.3 for Mac OS X x86 (32-bit, 64-bit Kernel) (March 5, \
2018)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)
(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[558, 20, 1325, 28, 94, "Input",ExpressionUUID->"a9eb43f9-c437-4ff0-ae8b-eb31f4d06b38"],
Cell[CellGroupData[{
Cell[1908, 52, 636, 12, 73, "Input",ExpressionUUID->"aa6e76d2-ec10-49c0-b6e0-1dd4d91b1a18"],
Cell[2547, 66, 443, 7, 46, "Print",ExpressionUUID->"f8f2db75-7118-45a9-adb1-7fa65a7be2f2"]
}, Open ]],
Cell[3005, 76, 1381, 32, 94, "Input",ExpressionUUID->"1c12ba0e-fd6c-4b7e-8de6-2b784fa72d92"],
Cell[4389, 110, 2376, 59, 199, "Input",ExpressionUUID->"476a37c6-482c-4f61-9ced-5d80a625711e"],
Cell[6768, 171, 6283, 168, 766, "Input",ExpressionUUID->"dfd6f63d-512e-4dca-8e27-b6dc7e095594"]
}
]
*)
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment