(* Content-type: application/vnd.wolfram.mathematica *) (*** Wolfram Notebook File ***) (* http://www.wolfram.com/nb *) (* CreatedBy='Mathematica 9.0' *) (*CacheID: 234*) (* Internal cache information: NotebookFileLineBreakTest NotebookFileLineBreakTest NotebookDataPosition[ 157, 7] NotebookDataLength[ 14072, 424] NotebookOptionsPosition[ 13240, 392] NotebookOutlinePosition[ 13599, 408] CellTagsIndexPosition[ 13556, 405] WindowFrame->Normal*) (* Beginning of Notebook Content *) Notebook[{ Cell[BoxData[{ RowBox[{ RowBox[{ "FIESTAPath", " ", "=", " ", "\"\\""}], ";"}], "\n", RowBox[{ RowBox[{"Get", "[", RowBox[{"FIESTAPath", " ", "<>", " ", "\"\\""}], "]"}], ";"}], "\n", RowBox[{ RowBox[{"NumberOfSubkernels", " ", "=", " ", "8"}], ";"}], "\n", RowBox[{ RowBox[{"NumberOfLinks", " ", "=", " ", "16"}], ";"}], "\n", RowBox[{ RowBox[{"CurrentIntegrator", " ", "=", " ", "\"\ \""}], ";"}], "\n", RowBox[{ RowBox[{"CurrentIntegratorSettings", " ", "=", " ", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"\"\ \"", ",", " ", "\"\<1.000000E-06\>\""}], "}"}], ",", " ", RowBox[{"{", RowBox[{"\"\ \"", ",", " ", "\"\<75000000\>\""}], "}"}]}], "}"}]}], ";"}]}], "Input", CellChangeTimes->{{3.810250834416864*^9, 3.8102508571677*^9}, { 3.8102509089504223`*^9, 3.810250956142045*^9}, {3.810260472223688*^9, 3.8102604740595417`*^9}, {3.830819714369892*^9, 3.8308197205196733`*^9}, 3.830998728241206*^9}], Cell[BoxData[ RowBox[{ RowBox[{ RowBox[{"ComplexMode", " ", "=", " ", "False"}], ";"}], RowBox[{"(*", RowBox[{ RowBox[{"In", " ", "physical", " ", "region"}], ",", " ", RowBox[{ "should", " ", "be", " ", "set", " ", "to", " ", "\"\ \""}]}], "*)"}]}]], "Input", CellChangeTimes->{{3.830819737820456*^9, 3.830819778382249*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"(*", RowBox[{"define", " ", "propagators"}], "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"den1", " ", "=", " ", RowBox[{ RowBox[{"-", RowBox[{ RowBox[{"(", RowBox[{"k1", " ", "+", " ", "p1"}], ")"}], "^", "2"}]}], " ", "+", " ", RowBox[{"m", "^", "2"}]}]}], ";", " ", RowBox[{"den2", " ", "=", " ", RowBox[{ RowBox[{"-", RowBox[{ RowBox[{"(", RowBox[{"k2", " ", "+", " ", "p1"}], ")"}], "^", "2"}]}], " ", "+", " ", RowBox[{"m", "^", "2"}]}]}], ";", " ", RowBox[{"den3", " ", "=", " ", RowBox[{ RowBox[{"-", RowBox[{ RowBox[{"(", RowBox[{"k3", " ", "+", " ", "p1"}], ")"}], "^", "2"}]}], " ", "+", " ", RowBox[{"m", "^", "2"}]}]}], ";", " ", RowBox[{"den4", " ", "=", " ", RowBox[{"-", RowBox[{ RowBox[{"(", RowBox[{"k3", " ", "+", " ", "p2"}], ")"}], "^", "2"}]}]}], ";", " ", RowBox[{"den5", " ", "=", " ", RowBox[{"-", RowBox[{ RowBox[{"(", RowBox[{"k2", " ", "+", " ", "p2"}], ")"}], "^", "2"}]}]}], ";", " ", RowBox[{"den6", " ", "=", " ", RowBox[{"-", RowBox[{ RowBox[{"(", RowBox[{"k1", " ", "+", " ", "p2"}], ")"}], "^", "2"}]}]}], ";", " ", RowBox[{"den7", " ", "=", " ", RowBox[{"-", RowBox[{ RowBox[{"(", "k1", ")"}], "^", "2"}]}]}], ";", " ", RowBox[{"den8", " ", "=", " ", RowBox[{"-", RowBox[{ RowBox[{"(", RowBox[{"k1", " ", "-", " ", "k2"}], ")"}], "^", "2"}]}]}], ";", " ", RowBox[{"den9", " ", "=", " ", RowBox[{"-", RowBox[{ RowBox[{"(", RowBox[{"k2", " ", "-", " ", "\n", "k3"}], ")"}], "^", "2"}]}]}], ";", " ", RowBox[{"den10", " ", "=", " ", RowBox[{"-", RowBox[{ RowBox[{"(", RowBox[{"k1", " ", "-", " ", "k3"}], ")"}], "^", "2"}]}]}], ";", " ", RowBox[{"den11", " ", "=", " ", RowBox[{"-", RowBox[{ RowBox[{"(", "k2", ")"}], "^", "2"}]}]}], ";", " ", RowBox[{"den12", " ", "=", " ", RowBox[{"-", RowBox[{ RowBox[{"(", "k3", ")"}], "^", "2"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"prop", "=", RowBox[{"{", RowBox[{ "den1", ",", "den2", ",", "den3", ",", "den4", ",", "den5", ",", "den6", ",", "den7", ",", "den8", ",", "den9", ",", "den10", ",", "den11", ",", "den12"}], "}"}]}], ";"}]}]}]], "Input", CellChangeTimes->{{3.830819811826806*^9, 3.83081984022077*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"MIs", "=", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ "1", ",", "1", ",", "0", ",", "1", ",", "1", ",", "0", ",", "1", ",", "1", ",", "1", ",", "1", ",", RowBox[{"-", "1"}], ",", "1"}], "}"}], ",", RowBox[{"{", RowBox[{"1", ",", "1", ",", "0", ",", "1", ",", "1", ",", RowBox[{"-", "1"}], ",", "1", ",", "1", ",", "1", ",", "1", ",", RowBox[{"-", "1"}], ",", "1"}], "}"}]}], "}"}]}], ";"}]], "Input", CellChangeTimes->{{3.830819852663513*^9, 3.830819858065852*^9}, 3.830827535250352*^9}], Cell[BoxData[ RowBox[{ RowBox[{"rep", "=", RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"p1", "^", "2"}], "\[Rule]", RowBox[{"m", "^", "2"}]}], ",", RowBox[{ RowBox[{"p2", "^", "2"}], "\[Rule]", "0"}], ",", RowBox[{ RowBox[{"p1", "*", "p2"}], "\[Rule]", RowBox[{"-", RowBox[{"(", RowBox[{ RowBox[{"(", RowBox[{"s", "-", RowBox[{"m", "^", "2"}]}], ")"}], "/", "2"}], ")"}]}]}], ",", RowBox[{"s", "\[Rule]", RowBox[{"-", "1.3"}]}], ",", RowBox[{"m", "\[Rule]", "1.0"}]}], "}"}]}], ";"}]], "Input", CellChangeTimes->{{3.830819847582056*^9, 3.8308198694263763`*^9}, { 3.830828089480948*^9, 3.8308281036225243`*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"SDEvaluate", "[", RowBox[{ RowBox[{"UF", "[", RowBox[{ RowBox[{"{", RowBox[{"k1", ",", "k2", ",", "k3"}], "}"}], ",", "prop", ",", "rep"}], "]"}], ",", RowBox[{"MIs", "[", RowBox[{"[", "1", "]"}], "]"}], ",", "0"}], "]"}]], "Input", CellChangeTimes->{{3.830819877183343*^9, 3.830819931691865*^9}}], Cell["\<\ FIESTA 4.1 Current integrator: vegasCuba CurrentIntegratorSettings: \ {{\"epsrel\",\"1.000000E-06\"},{\"epsabs\",\"1.000000E-12\"},{\"mineval\",\"0\ \"},{\"maxeval\",\"75000000\"},{\"nstart\",\"1000\"},{\"nincrease\",\"500\"},{\ \"seed\",\"0\"},{\"rng\",\"0\"}} Integration test passed Starting 8 subkernels UsingC: True NumberOfLinks: 16 UsingQLink: True Strategy: STRATEGY_S Integration has to be performed up to order 0 KLink created (2013 version)! You can read information on QOpen, QRead, \ QRemoveDatabase, QClose, QList, QSize, QPut, QGet, QSafeGet, QCheck and QRemove Sector decomposition - 9 sectors Primary sector 1 resulted in 30 sectors. Primary sector 2 resulted in 30 sectors. Primary sector 3 resulted in 150 sectors. Primary sector 4 resulted in 184 sectors. Primary sector 6 resulted in 152 sectors. Primary sector 8 resulted in 222 sectors. Primary sector 7 resulted in 238 sectors. Primary sector 9 resulted in 940 sectors. Primary sector 5 resulted in 2461 sectors. Totally: 45.0665 seconds; 18 sectors. Preparing database: 12.1874 seconds. Variable substitution..........110.7695 seconds; 8814 terms. Pole resolution..........43.503 seconds; 53338 terms. Expression preparation..........99.2294 seconds; 53338 terms. Epsilon expansion..........172.8283 seconds; 118367 terms. Preparing integration strings..........384.9944 seconds; 102576 terms. Database ready for integration. Terms of order -6: 62, max vars: 2 Integrating..........13.0251 seconds. Returned answer: 0.000788 + pm* 0 (0.000788)*ep^(-6) Terms of order -5: 366, max vars: 3 Integrating..........126.2503 seconds. Returned answer: -0.006299 + pm* 0 (0.000788)*ep^(-6)+(-0.003936)*ep^(-5) Terms of order -4: 1115, max vars: 4 Integrating..........654.2228 seconds. Returned answer: 0.034069 + pm* 1.*^-6 (0.000788)*ep^(-6)+(-0.003936)*ep^(-5)+(0.021002 + 1.*^-6*pm6)*ep^(-4) Terms of order -3: 2417, max vars: 5 Integrating..........2529.641 seconds. Returned answer: -0.090317 + pm* 0.000028 (0.000788)*ep^(-6)+(-0.003936)*ep^(-5)+(0.021002 + \ 1.*^-6*pm9)*ep^(-4)+(-0.025767 + 0.000028*pm10)*ep^(-3) Terms of order -2: 4267, max vars: 6 Integrating..........8053.6642 seconds. Returned answer: 0.42222 + pm* 0.000611 (0.000788)*ep^(-6)+(-0.003936)*ep^(-5)+(0.021002 + \ 1.*^-6*pm13)*ep^(-4)+(-0.025767 + 0.000028*pm14)*ep^(-3)+(0.344999 + \ 0.000617*pm15)*ep^(-2) Terms of order -1: 6397, max vars: 7 Integrating..........23106.8338 seconds. Returned answer: 0.404303 + pm* 0.012763 (0.000788)*ep^(-6)+(-0.003936)*ep^(-5)+(0.021002 + \ 1.*^-6*pm18)*ep^(-4)+(-0.025767 + 0.000028*pm19)*ep^(-3)+(0.344999 + \ 0.000617*pm20)*ep^(-2)+(1.297968 + 0.012895*pm21)*ep^(-1) Terms of order 0: 8143, max vars: 8 Integrating..........64314.2292 seconds. Returned answer: 1.848032 + pm* 0.256157 (0.000788)*ep^(-6)+(-0.003936)*ep^(-5)+(0.021002 + \ 1.*^-6*pm24)*ep^(-4)+(-0.025767 + 0.000028*pm25)*ep^(-3)+(0.344999 + \ 0.000617*pm26)*ep^(-2)+(1.297968 + 0.012895*pm27)*ep^(-1)+(5.638255 + \ 0.259042*pm28)*1 Total integration time: 98799.9828 Total time used: 99683.2 seconds.\ \>", "Print", CellChangeTimes->{{3.830819920360607*^9, 3.830819946792556*^9}, 3.8308275524038467`*^9}], Cell[BoxData[ RowBox[{"5.638255`", "\[VeryThinSpace]", "+", FractionBox["0.000788`", SuperscriptBox["ep", "6"]], "-", FractionBox["0.003936`", SuperscriptBox["ep", "5"]], "+", FractionBox[ RowBox[{"0.021002`", "\[VeryThinSpace]", "+", RowBox[{"1.`*^-6", " ", "pm31"}]}], SuperscriptBox["ep", "4"]], "+", FractionBox[ RowBox[{ RowBox[{"-", "0.025767`"}], "+", RowBox[{"0.000028`", " ", "pm32"}]}], SuperscriptBox["ep", "3"]], "+", FractionBox[ RowBox[{"0.344999`", "\[VeryThinSpace]", "+", RowBox[{"0.000617`", " ", "pm33"}]}], SuperscriptBox["ep", "2"]], "+", FractionBox[ RowBox[{"1.297968`", "\[VeryThinSpace]", "+", RowBox[{"0.012895`", " ", "pm34"}]}], "ep"], "+", RowBox[{"0.259042`", " ", "pm35"}]}]], "Output", CellChangeTimes->{3.83092724393983*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"SDEvaluate", "[", RowBox[{ RowBox[{"UF", "[", RowBox[{ RowBox[{"{", RowBox[{"k1", ",", "k2", ",", "k3"}], "}"}], ",", "prop", ",", "rep"}], "]"}], ",", RowBox[{"MIs", "[", RowBox[{"[", "2", "]"}], "]"}], ",", "0"}], "]"}]], "Input", CellChangeTimes->{{3.8308198786216497`*^9, 3.830819878622467*^9}, { 3.830819914178754*^9, 3.83081993374749*^9}}], Cell["\<\ FIESTA 4.1 Current integrator: vegasCuba CurrentIntegratorSettings: \ {{\"epsrel\",\"1.000000E-06\"},{\"epsabs\",\"1.000000E-12\"},{\"mineval\",\"0\ \"},{\"maxeval\",\"75000000\"},{\"nstart\",\"1000\"},{\"nincrease\",\"500\"},{\ \"seed\",\"0\"},{\"rng\",\"0\"}} Integration test passed Starting 8 subkernels UsingC: True NumberOfLinks: 16 UsingQLink: True Strategy: STRATEGY_S Integration has to be performed up to order 0 Sector decomposition - 9 sectors Primary sector 1 resulted in 30 sectors. Primary sector 2 resulted in 30 sectors. Primary sector 3 resulted in 150 sectors. Primary sector 4 resulted in 184 sectors. Primary sector 6 resulted in 152 sectors. Primary sector 8 resulted in 222 sectors. Primary sector 7 resulted in 238 sectors. Primary sector 9 resulted in 940 sectors. Primary sector 5 resulted in 2461 sectors. Totally: 46.2186 seconds; 18 sectors. Preparing database: 126.8654 seconds. Variable substitution..........1473.58 seconds; 26442 terms. Pole resolution..........142.3332 seconds; 31932 terms. Expression preparation..........142.5814 seconds; 31932 terms. Epsilon expansion..........127.6014 seconds; 37422 terms. Preparing integration strings..........224.0041 seconds; 37422 terms. Database ready for integration. Terms of order -1: 5490, max vars: 7 Integrating..........5154.5312 seconds. Returned answer: -6.694279 + pm* 0.000015 (-6.694279 + 0.000015*pm36)*ep^(-1) Terms of order 0: 26442, max vars: 8 Integrating..........59522.8412 seconds. Returned answer: -63.121114 + pm* 0.000286 (-6.694279 + 0.000015*pm37)*ep^(-1)+(-63.121114 + 0.000286*pm38)*1 Total integration time: 64678.4034 Total time used: 66967.3 seconds.\ \>", "Print", CellChangeTimes->{3.830927244778428*^9}], Cell[BoxData[ RowBox[{ RowBox[{"-", "63.121114`"}], "+", FractionBox[ RowBox[{ RowBox[{"-", "6.694279`"}], "+", RowBox[{"0.000015`", " ", "pm39"}]}], "ep"], "+", RowBox[{"0.000286`", " ", "pm40"}]}]], "Output", CellChangeTimes->{3.8309942206903477`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"66967.3", "+", "99683.2"}], ")"}], "/", "24"}], "/", "3600"}]], "Input", CellChangeTimes->{{3.830998683117448*^9, 3.830998699582486*^9}}], Cell[BoxData["1.9288252314814813`"], "Output", CellChangeTimes->{3.830998700563266*^9}] }, Open ]] }, WindowSize->{2560, 1337}, WindowMargins->{{-1, Automatic}, {Automatic, -1}}, ShowSelection->True, FrontEndVersion->"9.0 for Linux x86 (64-bit) (February 7, 2013)", StyleDefinitions->"Default.nb" ] (* End of Notebook Content *) (* Internal cache information *) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[557, 20, 1084, 29, 143, "Input"], Cell[1644, 51, 362, 10, 32, "Input"], Cell[2009, 63, 2684, 85, 99, "Input"], Cell[4696, 150, 600, 15, 32, "Input"], Cell[5299, 167, 737, 22, 32, "Input"], Cell[CellGroupData[{ Cell[6061, 193, 359, 10, 32, "Input"], Cell[6423, 205, 3173, 74, 951, "Print"], Cell[9599, 281, 836, 23, 56, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[10472, 309, 409, 11, 32, "Input"], Cell[10884, 322, 1728, 43, 615, "Print"], Cell[12615, 367, 275, 8, 51, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[12927, 380, 206, 6, 32, "Input"], Cell[13136, 388, 88, 1, 65, "Output"] }, Open ]] } ] *) (* End of internal cache information *)

map