Select Git revision
helper_functions.py
RDownScript.nb 130.78 KiB
(* 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[ 133756, 2738]
NotebookOptionsPosition[ 128157, 2637]
NotebookOutlinePosition[ 128497, 2652]
CellTagsIndexPosition[ 128454, 2649]
WindowFrame->Normal*)
(* Beginning of Notebook Content *)
Notebook[{
Cell[CellGroupData[{
Cell["RDown excitation time", "Title",
CellChangeTimes->{{3.7171355710697517`*^9,
3.7171356028895473`*^9}},ExpressionUUID->"b53761d4-1b24-4b86-b93e-\
9c44d4b29d65"],
Cell[BoxData[
RowBox[{"Quit", "[", "]"}]], "Input",
CellChangeTimes->{{3.772628159141247*^9, 3.772628161249263*^9}, {
3.785052589392028*^9, 3.785052592823535*^9}},
CellLabel->
"In[101]:=",ExpressionUUID->"980e79f6-9c07-4072-bd1e-300ac946e8c3"],
Cell[CellGroupData[{
Cell["Setup", "Section",
CellChangeTimes->{{3.717145376675954*^9, 3.7171453778419027`*^9}, {
3.725360092733674*^9,
3.725360093603505*^9}},ExpressionUUID->"fe9e23be-b684-46cc-9b54-\
957e00ada2e3"],
Cell[CellGroupData[{
Cell["Load package and set root paths", "Subsection",
CellChangeTimes->{{3.799119833730782*^9, 3.7991198361011353`*^9}, {
3.801906007839068*^9,
3.8019060131649723`*^9}},ExpressionUUID->"e952a5b9-bbfe-43a9-9668-\
8567ca76bf42"],
Cell[BoxData[
RowBox[{
RowBox[{"rootpath", "=", "\"\<~/git/rdstackingproject/\>\""}],
";"}]], "Code",
CellChangeTimes->{{3.801893836989571*^9, 3.8018938554497147`*^9}, {
3.8018991673772917`*^9, 3.801899171272743*^9}, {3.801899220625408*^9,
3.801899220979004*^9}, {3.801899255251523*^9, 3.8018992567883883`*^9}, {
3.801900861015358*^9, 3.801900867609252*^9}, {3.8019010062879868`*^9,
3.801901006499617*^9}, {3.801905973467177*^9, 3.8019060065780773`*^9}},
CellLabel->
"In[152]:=",ExpressionUUID->"ded0ca57-00ba-4008-906d-ee5041d06f91"],
Cell[BoxData[{
RowBox[{
RowBox[{"SetDirectory", "[",
RowBox[{"rootpath", "<>", "\"\<code\>\""}], "]"}], ";"}], "\n",
RowBox[{"<<", "RDown.m"}], "\n",
RowBox[{"<<", "mcmc.m"}]}], "Code",
CellChangeTimes->{{3.799119881021441*^9, 3.799119911381034*^9}, {
3.80077965736977*^9, 3.800779661525262*^9}, {3.801736519951005*^9, 3.8017365270768414`*^9}, {3.801750413991219*^9, 3.80175045124013*^9}, {
3.801754031026664*^9, 3.8017540356351223`*^9}, {3.801893846779317*^9,
3.8018938856249447`*^9}, {3.80189915075161*^9, 3.801899152360096*^9}, {
3.801899239104719*^9, 3.801899268047777*^9}, 3.801899313824655*^9},
CellLabel->
"In[153]:=",ExpressionUUID->"df18942a-15e3-4e8b-880e-923fca9e18d7"],
Cell[BoxData[{
RowBox[{
RowBox[{
"mysxscase", "=",
"\"\</Users/xisco/SXS/BBH_SKS_d14.3_q1.22_sA_0_0_0.330_sB_0_0_-0.440\>\""}]\
, ";"}], "\n",
RowBox[{
RowBox[{"tshift", "=", "0"}], ";"}], "\n",
RowBox[{
RowBox[{"npoints", "=", "10"}], ";"}], "\n",
RowBox[{
RowBox[{"ntones", "=", "1"}], ";"}], "\n",
RowBox[{
RowBox[{"\[Omega]fact", "=", "0.05"}], ";"}]}], "Code",
CellChangeTimes->{{3.8019064095845547`*^9, 3.801906411959391*^9}, {
3.801906985378249*^9, 3.801906993215358*^9}},
CellLabel->
"In[228]:=",ExpressionUUID->"d0e5d8fb-ddc3-4ce5-93a5-2a83bdf621d4"]
}, Closed]],
Cell[CellGroupData[{
Cell["Setup paths and other notebook folders", "Subsection",
CellChangeTimes->{{3.76257774139565*^9,
3.7625777648770647`*^9}},ExpressionUUID->"46366112-0a9f-4a64-97d3-\
72f08bceb964"],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{
RowBox[{"(*", " ",
RowBox[{"Change", " ", "them", " ", "accordingly"}], " ", "*)"}], "\n",
RowBox[{"(*",
RowBox[{"Notebook", " ", "directory"}], "*)"}], "\n",
RowBox[{
RowBox[{
RowBox[{"notdir", "=",
RowBox[{"rootpath", "<>", "\"\<code/\>\""}]}], ";"}], "\n",
RowBox[{"(*",
RowBox[{"Directory", " ", "of", " ", "the", " ",
RowBox[{"l", "'"}], "s", " ", "Ringdown", " ", "data"}], " ", "*)"}],
"\n",
RowBox[{
RowBox[{"modedir", "=", "rootpath"}], ";"}], "\n",
RowBox[{"(*", " ",
RowBox[{
"True", " ", "if", " ", "you", " ", "want", " ", "to", " ", "export", " ",
"the", " ", "plots"}], " ", "*)"}], "\n",
RowBox[{"ExportQ", "=", "False"}], "\n",
RowBox[{"ExportDir", "=", "notdir"}]}]}]], "Code",
CellChangeTimes->{
3.7625787067637377`*^9, {3.7625787621100388`*^9, 3.762578766498042*^9}, {
3.7625788039303837`*^9, 3.76257880843638*^9}, {3.762579176451227*^9,
3.762579205498322*^9}, 3.783842026281891*^9, {3.785048517032866*^9,
3.785048568634179*^9}, 3.785052064855216*^9, 3.785052098677169*^9, {
3.785052662782302*^9, 3.785052711482973*^9}, {3.799112511472287*^9,
3.799112532529355*^9}, 3.799119917854619*^9, {3.8007685960415173`*^9,
3.800768686902137*^9}, {3.8007688081677713`*^9, 3.80076883021767*^9}, {
3.8007688671018257`*^9, 3.800768867313781*^9}, {3.8017363103133783`*^9,
3.80173641269767*^9}, {3.801736444416366*^9, 3.801736507343718*^9}, {
3.801736668985515*^9, 3.801736672338735*^9}, {3.801887712719111*^9,
3.8018877185807247`*^9}, {3.8018939370001*^9, 3.801893986031884*^9}, 3.801899307537121*^9, {3.801899347317583*^9, 3.8018994021627417`*^9},
3.801899549788517*^9, {3.801900970284079*^9, 3.8019010086863937`*^9}},
CellLabel->
"In[121]:=",ExpressionUUID->"2c37e1b9-5cb9-4980-a108-163d386515de"],
Cell[BoxData["False"], "Output",
CellChangeTimes->{
3.8019008734562407`*^9, {3.801900996088092*^9, 3.801901009074788*^9}},
CellLabel->
"Out[123]=",ExpressionUUID->"04cad810-7ff2-4fa4-88b1-a3cb2e1eaf60"],
Cell[BoxData["\<\"~/git/rdstackingproject/code/\"\>"], "Output",
CellChangeTimes->{
3.8019008734562407`*^9, {3.801900996088092*^9, 3.801901009076292*^9}},
CellLabel->
"Out[124]=",ExpressionUUID->"3a282303-63af-4a85-87c0-523e68dad528"]
}, Open ]]
}, Closed]],
Cell[CellGroupData[{
Cell["Plot Formatting", "Subsection",
CellChangeTimes->{{3.801736555565851*^9,
3.80173655943112*^9}},ExpressionUUID->"2b160476-27ad-417e-ba6b-\
aba2f29dfd9c"],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{
RowBox[{"(*", " ", "PlotStyle", " ", "*)"}], "\n",
RowBox[{
RowBox[{"markers", "=",
RowBox[{"{",
RowBox[{
"\"\<\[FilledCircle]\>\"", ",", "\"\<\[FilledSquare]\>\"", ",",
"\"\<\[FilledDiamond]\>\"", ",", "\"\<\[FilledUpTriangle]\>\"", ",",
"\"\<\[FilledDownTriangle]\>\"", ",", "\"\<\[EmptyCircle]\>\"", ",",
"\"\<\[EmptySquare]\>\"", ",", "\"\<\[EmptyDiamond]\>\"", ",",
"\"\<\[EmptyUpTriangle]\>\"", ",", "\"\<\[EmptyDownTriangle]\>\""}],
"}"}]}], "\n",
RowBox[{"colors", "=",
RowBox[{"ColorData", "[",
RowBox[{"97", ",", "\"\<ColorList\>\""}], "]"}]}], "\n",
RowBox[{
RowBox[{"l1", "=", ".025"}], ";"}], "\n",
RowBox[{
RowBox[{"l2", "=", ".001"}], ";"}], "\n",
RowBox[{
RowBox[{"s", "=",
RowBox[{"l1", "/", "2"}]}], ";"}]}]}]], "Code",
CellChangeTimes->{{3.770258655384822*^9, 3.770258658711926*^9}, {
3.785048597172567*^9, 3.785048611784268*^9}, {3.7850518765860777`*^9,
3.785051879447013*^9}},
CellLabel->
"In[102]:=",ExpressionUUID->"91d0d186-3694-4095-b0ad-f5476cf3486b"],
Cell[BoxData[
RowBox[{"{",
RowBox[{"\<\"\[FilledCircle]\"\>", ",", "\<\"\[FilledSquare]\"\>",
",", "\<\"\[FilledDiamond]\"\>", ",", "\<\"\[FilledUpTriangle]\"\>",
",", "\<\"\[FilledDownTriangle]\"\>", ",", "\<\"\[EmptyCircle]\"\>",
",", "\<\"\[EmptySquare]\"\>", ",", "\<\"\[EmptyDiamond]\"\>",
",", "\<\"\[EmptyUpTriangle]\"\>", ",", "\<\"\[EmptyDownTriangle]\"\>"}],
"}"}]], "Output",
CellChangeTimes->{{3.785048598454445*^9, 3.7850486121743803`*^9},
3.7850492534864817`*^9, 3.7850518797807217`*^9, 3.7851593453111877`*^9,
3.788525055881117*^9, 3.788588576142179*^9, 3.789016362482636*^9,
3.7891011551402693`*^9, 3.7971326144003468`*^9, 3.799112667826545*^9,
3.799113935678566*^9, 3.799119920100224*^9, 3.79912099675138*^9, 3.799121678020488*^9, 3.79912196315336*^9, 3.7991220841956873`*^9, {
3.799122272245305*^9, 3.799122300399253*^9}, 3.8000746708421507`*^9,
3.8000885291890907`*^9, 3.800088987611136*^9, {3.800089240452045*^9,
3.800089269881946*^9}, 3.800089321591539*^9, 3.800089431077919*^9,
3.800089495187615*^9, 3.800089595616612*^9, 3.800148380233251*^9,
3.800539327345366*^9, 3.800610134971425*^9, 3.800705385664913*^9,
3.80070556987372*^9, 3.800768504913447*^9, 3.800798370679039*^9,
3.8008772883571987`*^9, 3.801281244944542*^9, 3.801285594509903*^9,
3.801567580742796*^9, {3.801723018647998*^9, 3.801723046894418*^9},
3.801736561561638*^9, 3.801738393523629*^9, 3.801744291978118*^9,
3.8017504625330067`*^9, 3.801754015153778*^9, 3.801754114168932*^9,
3.8018007150796423`*^9, 3.8018040233399773`*^9, 3.801882349914794*^9,
3.8018951501756773`*^9, {3.80190087812988*^9, 3.8019008850418577`*^9}},
CellLabel->
"Out[102]=",ExpressionUUID->"4661e259-0c1d-4e90-b4af-b76e6ee6e71c"],
Cell[BoxData[
RowBox[{"{",
RowBox[{
InterpretationBox[
ButtonBox[
TooltipBox[
GraphicsBox[{
{GrayLevel[0], RectangleBox[{0, 0}]},
{GrayLevel[0], RectangleBox[{1, -1}]},
{RGBColor[0.368417, 0.506779, 0.709798],
RectangleBox[{0, -1}, {2, 1}]}},
AspectRatio->1,
DefaultBaseStyle->"ColorSwatchGraphics",
Frame->True,
FrameStyle->RGBColor[
0.24561133333333335`, 0.3378526666666667, 0.4731986666666667],
FrameTicks->None,
ImageSize->
Dynamic[{
Automatic, 1.35 CurrentValue["FontCapHeight"]/AbsoluteCurrentValue[
Magnification]}],
PlotRangePadding->None],
StyleBox[
RowBox[{"RGBColor", "[",
RowBox[{"0.368417`", ",", "0.506779`", ",", "0.709798`"}], "]"}],
NumberMarks -> False]],
Appearance->None,
BaseStyle->{},
BaselinePosition->Baseline,
ButtonFunction:>With[{Typeset`box$ = EvaluationBox[]},
If[
Not[
AbsoluteCurrentValue["Deployed"]],
SelectionMove[Typeset`box$, All, Expression];
FrontEnd`Private`$ColorSelectorInitialAlpha = 1;
FrontEnd`Private`$ColorSelectorInitialColor =
RGBColor[0.368417, 0.506779, 0.709798];
FrontEnd`Private`$ColorSelectorUseMakeBoxes = True;
MathLink`CallFrontEnd[
FrontEnd`AttachCell[Typeset`box$,
FrontEndResource["RGBColorValueSelector"], {0, {Left, Bottom}}, {
Left, Top},
"ClosingActions" -> {
"SelectionDeparture", "ParentChanged", "EvaluatorQuit"}]]]],
DefaultBaseStyle->{},
Evaluator->Automatic,
Method->"Preemptive"],
RGBColor[0.368417, 0.506779, 0.709798],
Editable->False,
Selectable->False], ",",
InterpretationBox[
ButtonBox[
TooltipBox[
GraphicsBox[{ {GrayLevel[0], RectangleBox[{0, 0}]},
{GrayLevel[0], RectangleBox[{1, -1}]},
{RGBColor[0.880722, 0.611041, 0.142051],
RectangleBox[{0, -1}, {2, 1}]}},
AspectRatio->1,
DefaultBaseStyle->"ColorSwatchGraphics",
Frame->True,
FrameStyle->RGBColor[
0.587148, 0.40736066666666665`, 0.09470066666666668],
FrameTicks->None,
ImageSize->
Dynamic[{
Automatic, 1.35 CurrentValue["FontCapHeight"]/AbsoluteCurrentValue[
Magnification]}],
PlotRangePadding->None],
StyleBox[
RowBox[{"RGBColor", "[",
RowBox[{"0.880722`", ",", "0.611041`", ",", "0.142051`"}], "]"}],
NumberMarks -> False]],
Appearance->None,
BaseStyle->{},
BaselinePosition->Baseline,
ButtonFunction:>With[{Typeset`box$ = EvaluationBox[]},
If[
Not[
AbsoluteCurrentValue["Deployed"]],
SelectionMove[Typeset`box$, All, Expression];
FrontEnd`Private`$ColorSelectorInitialAlpha = 1;
FrontEnd`Private`$ColorSelectorInitialColor =
RGBColor[0.880722, 0.611041, 0.142051];
FrontEnd`Private`$ColorSelectorUseMakeBoxes = True;
MathLink`CallFrontEnd[
FrontEnd`AttachCell[Typeset`box$,
FrontEndResource["RGBColorValueSelector"], {0, {Left, Bottom}}, {
Left, Top},
"ClosingActions" -> {
"SelectionDeparture", "ParentChanged", "EvaluatorQuit"}]]]],
DefaultBaseStyle->{},
Evaluator->Automatic,
Method->"Preemptive"],
RGBColor[0.880722, 0.611041, 0.142051],
Editable->False,
Selectable->False], ",",
InterpretationBox[
ButtonBox[
TooltipBox[
GraphicsBox[{
{GrayLevel[0], RectangleBox[{0, 0}]},
{GrayLevel[0], RectangleBox[{1, -1}]},
{RGBColor[0.560181, 0.691569, 0.194885],
RectangleBox[{0, -1}, {2, 1}]}},
AspectRatio->1,
DefaultBaseStyle->"ColorSwatchGraphics",
Frame->True,
FrameStyle->RGBColor[
0.37345400000000006`, 0.461046, 0.12992333333333334`],
FrameTicks->None,
ImageSize->
Dynamic[{
Automatic, 1.35 CurrentValue["FontCapHeight"]/AbsoluteCurrentValue[
Magnification]}],
PlotRangePadding->None],
StyleBox[
RowBox[{"RGBColor", "[",
RowBox[{"0.560181`", ",", "0.691569`", ",", "0.194885`"}], "]"}],
NumberMarks -> False]],
Appearance->None,
BaseStyle->{},
BaselinePosition->Baseline,
ButtonFunction:>With[{Typeset`box$ = EvaluationBox[]}, If[
Not[
AbsoluteCurrentValue["Deployed"]],
SelectionMove[Typeset`box$, All, Expression];
FrontEnd`Private`$ColorSelectorInitialAlpha = 1;
FrontEnd`Private`$ColorSelectorInitialColor =
RGBColor[0.560181, 0.691569, 0.194885];
FrontEnd`Private`$ColorSelectorUseMakeBoxes = True;
MathLink`CallFrontEnd[
FrontEnd`AttachCell[Typeset`box$,
FrontEndResource["RGBColorValueSelector"], {0, {Left, Bottom}}, {
Left, Top},
"ClosingActions" -> {
"SelectionDeparture", "ParentChanged", "EvaluatorQuit"}]]]],
DefaultBaseStyle->{},
Evaluator->Automatic,
Method->"Preemptive"],
RGBColor[0.560181, 0.691569, 0.194885],
Editable->False,
Selectable->False], ",",
InterpretationBox[
ButtonBox[
TooltipBox[
GraphicsBox[{
{GrayLevel[0], RectangleBox[{0, 0}]},
{GrayLevel[0], RectangleBox[{1, -1}]},
{RGBColor[0.922526, 0.385626, 0.209179],
RectangleBox[{0, -1}, {2, 1}]}},
AspectRatio->1,
DefaultBaseStyle->"ColorSwatchGraphics",
Frame->True,
FrameStyle->RGBColor[
0.6150173333333333, 0.25708400000000003`, 0.13945266666666667`],
FrameTicks->None,
ImageSize->
Dynamic[{
Automatic, 1.35 CurrentValue["FontCapHeight"]/AbsoluteCurrentValue[
Magnification]}],
PlotRangePadding->None],
StyleBox[
RowBox[{"RGBColor", "[",
RowBox[{"0.922526`", ",", "0.385626`", ",", "0.209179`"}], "]"}],
NumberMarks -> False]],
Appearance->None,
BaseStyle->{},
BaselinePosition->Baseline,
ButtonFunction:>With[{Typeset`box$ = EvaluationBox[]},
If[
Not[
AbsoluteCurrentValue["Deployed"]],
SelectionMove[Typeset`box$, All, Expression];
FrontEnd`Private`$ColorSelectorInitialAlpha = 1;
FrontEnd`Private`$ColorSelectorInitialColor =
RGBColor[0.922526, 0.385626, 0.209179];
FrontEnd`Private`$ColorSelectorUseMakeBoxes = True;
MathLink`CallFrontEnd[
FrontEnd`AttachCell[Typeset`box$,
FrontEndResource["RGBColorValueSelector"], {0, {Left, Bottom}}, {
Left, Top},
"ClosingActions" -> {
"SelectionDeparture", "ParentChanged", "EvaluatorQuit"}]]]],
DefaultBaseStyle->{},
Evaluator->Automatic,
Method->"Preemptive"],
RGBColor[0.922526, 0.385626, 0.209179],
Editable->False,
Selectable->False], ",",
InterpretationBox[
ButtonBox[
TooltipBox[ GraphicsBox[{
{GrayLevel[0], RectangleBox[{0, 0}]},
{GrayLevel[0], RectangleBox[{1, -1}]},
{RGBColor[0.528488, 0.470624, 0.701351],
RectangleBox[{0, -1}, {2, 1}]}},
AspectRatio->1,
DefaultBaseStyle->"ColorSwatchGraphics",
Frame->True,
FrameStyle->RGBColor[
0.3523253333333333, 0.3137493333333333, 0.46756733333333333`],
FrameTicks->None,
ImageSize->
Dynamic[{
Automatic, 1.35 CurrentValue["FontCapHeight"]/AbsoluteCurrentValue[
Magnification]}],
PlotRangePadding->None],
StyleBox[
RowBox[{"RGBColor", "[",
RowBox[{"0.528488`", ",", "0.470624`", ",", "0.701351`"}], "]"}],
NumberMarks -> False]],
Appearance->None,
BaseStyle->{},
BaselinePosition->Baseline,
ButtonFunction:>With[{Typeset`box$ = EvaluationBox[]},
If[
Not[
AbsoluteCurrentValue["Deployed"]],
SelectionMove[Typeset`box$, All, Expression];
FrontEnd`Private`$ColorSelectorInitialAlpha = 1;
FrontEnd`Private`$ColorSelectorInitialColor =
RGBColor[0.528488, 0.470624, 0.701351];
FrontEnd`Private`$ColorSelectorUseMakeBoxes = True;
MathLink`CallFrontEnd[
FrontEnd`AttachCell[Typeset`box$,
FrontEndResource["RGBColorValueSelector"], {0, {Left, Bottom}}, {
Left, Top},
"ClosingActions" -> {
"SelectionDeparture", "ParentChanged", "EvaluatorQuit"}]]]],
DefaultBaseStyle->{},
Evaluator->Automatic,
Method->"Preemptive"],
RGBColor[0.528488, 0.470624, 0.701351],
Editable->False,
Selectable->False], ",",
InterpretationBox[
ButtonBox[
TooltipBox[
GraphicsBox[{
{GrayLevel[0], RectangleBox[{0, 0}]},
{GrayLevel[0], RectangleBox[{1, -1}]},
{RGBColor[0.772079, 0.431554, 0.102387],
RectangleBox[{0, -1}, {2, 1}]}},
AspectRatio->1,
DefaultBaseStyle->"ColorSwatchGraphics",
Frame->True,
FrameStyle->RGBColor[
0.5147193333333333, 0.28770266666666666`, 0.06825800000000001],
FrameTicks->None,
ImageSize->
Dynamic[{
Automatic, 1.35 CurrentValue["FontCapHeight"]/AbsoluteCurrentValue[
Magnification]}],
PlotRangePadding->None],
StyleBox[
RowBox[{"RGBColor", "[",
RowBox[{"0.772079`", ",", "0.431554`", ",", "0.102387`"}], "]"}],
NumberMarks -> False]],
Appearance->None,
BaseStyle->{},
BaselinePosition->Baseline, ButtonFunction:>With[{Typeset`box$ = EvaluationBox[]},
If[
Not[
AbsoluteCurrentValue["Deployed"]],
SelectionMove[Typeset`box$, All, Expression];
FrontEnd`Private`$ColorSelectorInitialAlpha = 1;
FrontEnd`Private`$ColorSelectorInitialColor =
RGBColor[0.772079, 0.431554, 0.102387];
FrontEnd`Private`$ColorSelectorUseMakeBoxes = True;
MathLink`CallFrontEnd[
FrontEnd`AttachCell[Typeset`box$,
FrontEndResource["RGBColorValueSelector"], {0, {Left, Bottom}}, {
Left, Top},
"ClosingActions" -> {
"SelectionDeparture", "ParentChanged", "EvaluatorQuit"}]]]],
DefaultBaseStyle->{},
Evaluator->Automatic,
Method->"Preemptive"],
RGBColor[0.772079, 0.431554, 0.102387],
Editable->False,
Selectable->False], ",",
InterpretationBox[
ButtonBox[
TooltipBox[
GraphicsBox[{
{GrayLevel[0], RectangleBox[{0, 0}]},
{GrayLevel[0], RectangleBox[{1, -1}]},
{RGBColor[0.363898, 0.618501, 0.782349],
RectangleBox[{0, -1}, {2, 1}]}},
AspectRatio->1,
DefaultBaseStyle->"ColorSwatchGraphics",
Frame->True,
FrameStyle->RGBColor[0.24259866666666668`, 0.412334, 0.521566],
FrameTicks->None,
ImageSize->
Dynamic[{
Automatic, 1.35 CurrentValue["FontCapHeight"]/AbsoluteCurrentValue[
Magnification]}],
PlotRangePadding->None],
StyleBox[
RowBox[{"RGBColor", "[",
RowBox[{"0.363898`", ",", "0.618501`", ",", "0.782349`"}], "]"}],
NumberMarks -> False]],
Appearance->None,
BaseStyle->{},
BaselinePosition->Baseline,
ButtonFunction:>With[{Typeset`box$ = EvaluationBox[]},
If[
Not[
AbsoluteCurrentValue["Deployed"]],
SelectionMove[Typeset`box$, All, Expression];
FrontEnd`Private`$ColorSelectorInitialAlpha = 1;
FrontEnd`Private`$ColorSelectorInitialColor =
RGBColor[0.363898, 0.618501, 0.782349];
FrontEnd`Private`$ColorSelectorUseMakeBoxes = True;
MathLink`CallFrontEnd[
FrontEnd`AttachCell[Typeset`box$,
FrontEndResource["RGBColorValueSelector"], {0, {Left, Bottom}}, {
Left, Top},
"ClosingActions" -> {
"SelectionDeparture", "ParentChanged", "EvaluatorQuit"}]]]],
DefaultBaseStyle->{},
Evaluator->Automatic,
Method->"Preemptive"],
RGBColor[0.363898, 0.618501, 0.782349],
Editable->False,
Selectable->False], ",",
InterpretationBox[
ButtonBox[
TooltipBox[ GraphicsBox[{
{GrayLevel[0], RectangleBox[{0, 0}]},
{GrayLevel[0], RectangleBox[{1, -1}]},
{RGBColor[1, 0.75, 0], RectangleBox[{0, -1}, {2, 1}]}},
AspectRatio->1,
DefaultBaseStyle->"ColorSwatchGraphics",
Frame->True,
FrameStyle->RGBColor[0.6666666666666666, 0.5, 0.],
FrameTicks->None,
ImageSize->
Dynamic[{
Automatic, 1.35 CurrentValue["FontCapHeight"]/AbsoluteCurrentValue[
Magnification]}],
PlotRangePadding->None],
StyleBox[
RowBox[{"RGBColor", "[",
RowBox[{"1", ",", "0.75`", ",", "0"}], "]"}], NumberMarks -> False]],
Appearance->None,
BaseStyle->{},
BaselinePosition->Baseline,
ButtonFunction:>With[{Typeset`box$ = EvaluationBox[]},
If[
Not[
AbsoluteCurrentValue["Deployed"]],
SelectionMove[Typeset`box$, All, Expression];
FrontEnd`Private`$ColorSelectorInitialAlpha = 1;
FrontEnd`Private`$ColorSelectorInitialColor = RGBColor[1, 0.75, 0];
FrontEnd`Private`$ColorSelectorUseMakeBoxes = True;
MathLink`CallFrontEnd[
FrontEnd`AttachCell[Typeset`box$,
FrontEndResource["RGBColorValueSelector"], {0, {Left, Bottom}}, {
Left, Top},
"ClosingActions" -> {
"SelectionDeparture", "ParentChanged", "EvaluatorQuit"}]]]],
DefaultBaseStyle->{},
Evaluator->Automatic,
Method->"Preemptive"],
RGBColor[1, 0.75, 0],
Editable->False,
Selectable->False], ",",
InterpretationBox[
ButtonBox[
TooltipBox[
GraphicsBox[{
{GrayLevel[0], RectangleBox[{0, 0}]},
{GrayLevel[0], RectangleBox[{1, -1}]},
{RGBColor[0.647624, 0.37816, 0.614037],
RectangleBox[{0, -1}, {2, 1}]}},
AspectRatio->1,
DefaultBaseStyle->"ColorSwatchGraphics",
Frame->True,
FrameStyle->RGBColor[
0.4317493333333333, 0.2521066666666667, 0.40935800000000006`],
FrameTicks->None,
ImageSize->
Dynamic[{
Automatic, 1.35 CurrentValue["FontCapHeight"]/AbsoluteCurrentValue[
Magnification]}],
PlotRangePadding->None],
StyleBox[
RowBox[{"RGBColor", "[",
RowBox[{"0.647624`", ",", "0.37816`", ",", "0.614037`"}], "]"}],
NumberMarks -> False]],
Appearance->None,
BaseStyle->{},
BaselinePosition->Baseline,
ButtonFunction:>With[{Typeset`box$ = EvaluationBox[]},
If[
Not[ AbsoluteCurrentValue["Deployed"]],
SelectionMove[Typeset`box$, All, Expression];
FrontEnd`Private`$ColorSelectorInitialAlpha = 1;
FrontEnd`Private`$ColorSelectorInitialColor =
RGBColor[0.647624, 0.37816, 0.614037];
FrontEnd`Private`$ColorSelectorUseMakeBoxes = True;
MathLink`CallFrontEnd[
FrontEnd`AttachCell[Typeset`box$,
FrontEndResource["RGBColorValueSelector"], {0, {Left, Bottom}}, {
Left, Top},
"ClosingActions" -> {
"SelectionDeparture", "ParentChanged", "EvaluatorQuit"}]]]],
DefaultBaseStyle->{},
Evaluator->Automatic,
Method->"Preemptive"],
RGBColor[0.647624, 0.37816, 0.614037],
Editable->False,
Selectable->False], ",",
InterpretationBox[
ButtonBox[
TooltipBox[
GraphicsBox[{
{GrayLevel[0], RectangleBox[{0, 0}]},
{GrayLevel[0], RectangleBox[{1, -1}]},
{RGBColor[0.571589, 0.586483, 0.], RectangleBox[{0, -1}, {2, 1}]}},
AspectRatio->1,
DefaultBaseStyle->"ColorSwatchGraphics",
Frame->True,
FrameStyle->RGBColor[0.38105933333333336`, 0.39098866666666665`, 0.],
FrameTicks->None,
ImageSize->
Dynamic[{
Automatic, 1.35 CurrentValue["FontCapHeight"]/AbsoluteCurrentValue[
Magnification]}],
PlotRangePadding->None],
StyleBox[
RowBox[{"RGBColor", "[",
RowBox[{"0.571589`", ",", "0.586483`", ",", "0.`"}], "]"}],
NumberMarks -> False]],
Appearance->None,
BaseStyle->{},
BaselinePosition->Baseline,
ButtonFunction:>With[{Typeset`box$ = EvaluationBox[]},
If[
Not[
AbsoluteCurrentValue["Deployed"]],
SelectionMove[Typeset`box$, All, Expression];
FrontEnd`Private`$ColorSelectorInitialAlpha = 1;
FrontEnd`Private`$ColorSelectorInitialColor =
RGBColor[0.571589, 0.586483, 0.];
FrontEnd`Private`$ColorSelectorUseMakeBoxes = True;
MathLink`CallFrontEnd[
FrontEnd`AttachCell[Typeset`box$,
FrontEndResource["RGBColorValueSelector"], {0, {Left, Bottom}}, {
Left, Top},
"ClosingActions" -> {
"SelectionDeparture", "ParentChanged", "EvaluatorQuit"}]]]],
DefaultBaseStyle->{},
Evaluator->Automatic,
Method->"Preemptive"],
RGBColor[0.571589, 0.586483, 0.],
Editable->False,
Selectable->False], ",",
InterpretationBox[
ButtonBox[
TooltipBox[
GraphicsBox[{
{GrayLevel[0], RectangleBox[{0, 0}]},
{GrayLevel[0], RectangleBox[{1, -1}]},
{RGBColor[0.915, 0.3325, 0.2125], RectangleBox[{0, -1}, {2, 1}]}}, AspectRatio->1,
DefaultBaseStyle->"ColorSwatchGraphics",
Frame->True,
FrameStyle->RGBColor[
0.6100000000000001, 0.22166666666666668`, 0.14166666666666666`],
FrameTicks->None,
ImageSize->
Dynamic[{
Automatic, 1.35 CurrentValue["FontCapHeight"]/AbsoluteCurrentValue[
Magnification]}],
PlotRangePadding->None],
StyleBox[
RowBox[{"RGBColor", "[",
RowBox[{"0.915`", ",", "0.3325`", ",", "0.2125`"}], "]"}],
NumberMarks -> False]],
Appearance->None,
BaseStyle->{},
BaselinePosition->Baseline,
ButtonFunction:>With[{Typeset`box$ = EvaluationBox[]},
If[
Not[
AbsoluteCurrentValue["Deployed"]],
SelectionMove[Typeset`box$, All, Expression];
FrontEnd`Private`$ColorSelectorInitialAlpha = 1;
FrontEnd`Private`$ColorSelectorInitialColor =
RGBColor[0.915, 0.3325, 0.2125];
FrontEnd`Private`$ColorSelectorUseMakeBoxes = True;
MathLink`CallFrontEnd[
FrontEnd`AttachCell[Typeset`box$,
FrontEndResource["RGBColorValueSelector"], {0, {Left, Bottom}}, {
Left, Top},
"ClosingActions" -> {
"SelectionDeparture", "ParentChanged", "EvaluatorQuit"}]]]],
DefaultBaseStyle->{},
Evaluator->Automatic,
Method->"Preemptive"],
RGBColor[0.915, 0.3325, 0.2125],
Editable->False,
Selectable->False], ",",
InterpretationBox[
ButtonBox[
TooltipBox[
GraphicsBox[{
{GrayLevel[0], RectangleBox[{0, 0}]},
{GrayLevel[0], RectangleBox[{1, -1}]},
{RGBColor[0.40082222609352647`, 0.5220066643438841, 0.85],
RectangleBox[{0, -1}, {2, 1}]}},
AspectRatio->1,
DefaultBaseStyle->"ColorSwatchGraphics",
Frame->True,
FrameStyle->RGBColor[
0.2672148173956843, 0.34800444289592275`, 0.5666666666666667],
FrameTicks->None,
ImageSize->
Dynamic[{
Automatic, 1.35 CurrentValue["FontCapHeight"]/AbsoluteCurrentValue[
Magnification]}],
PlotRangePadding->None],
StyleBox[
RowBox[{"RGBColor", "[",
RowBox[{
"0.40082222609352647`", ",", "0.5220066643438841`", ",", "0.85`"}],
"]"}], NumberMarks -> False]],
Appearance->None,
BaseStyle->{},
BaselinePosition->Baseline,
ButtonFunction:>With[{Typeset`box$ = EvaluationBox[]},
If[
Not[
AbsoluteCurrentValue["Deployed"]], SelectionMove[Typeset`box$, All, Expression];
FrontEnd`Private`$ColorSelectorInitialAlpha = 1;
FrontEnd`Private`$ColorSelectorInitialColor =
RGBColor[0.40082222609352647`, 0.5220066643438841, 0.85];
FrontEnd`Private`$ColorSelectorUseMakeBoxes = True;
MathLink`CallFrontEnd[
FrontEnd`AttachCell[Typeset`box$,
FrontEndResource["RGBColorValueSelector"], {0, {Left, Bottom}}, {
Left, Top},
"ClosingActions" -> {
"SelectionDeparture", "ParentChanged", "EvaluatorQuit"}]]]],
DefaultBaseStyle->{},
Evaluator->Automatic,
Method->"Preemptive"],
RGBColor[0.40082222609352647`, 0.5220066643438841, 0.85],
Editable->False,
Selectable->False], ",",
InterpretationBox[
ButtonBox[
TooltipBox[
GraphicsBox[{
{GrayLevel[0], RectangleBox[{0, 0}]},
{GrayLevel[0], RectangleBox[{1, -1}]},
{RGBColor[0.9728288904374106, 0.621644452187053, 0.07336199581899142],
RectangleBox[{0, -1}, {2, 1}]}},
AspectRatio->1,
DefaultBaseStyle->"ColorSwatchGraphics",
Frame->True,
FrameStyle->RGBColor[
0.6485525936249404, 0.4144296347913687, 0.048907997212660946`],
FrameTicks->None,
ImageSize->
Dynamic[{
Automatic, 1.35 CurrentValue["FontCapHeight"]/AbsoluteCurrentValue[
Magnification]}],
PlotRangePadding->None],
StyleBox[
RowBox[{"RGBColor", "[",
RowBox[{
"0.9728288904374106`", ",", "0.621644452187053`", ",",
"0.07336199581899142`"}], "]"}], NumberMarks -> False]],
Appearance->None,
BaseStyle->{},
BaselinePosition->Baseline,
ButtonFunction:>With[{Typeset`box$ = EvaluationBox[]},
If[
Not[
AbsoluteCurrentValue["Deployed"]],
SelectionMove[Typeset`box$, All, Expression];
FrontEnd`Private`$ColorSelectorInitialAlpha = 1;
FrontEnd`Private`$ColorSelectorInitialColor =
RGBColor[0.9728288904374106, 0.621644452187053, 0.07336199581899142];
FrontEnd`Private`$ColorSelectorUseMakeBoxes = True;
MathLink`CallFrontEnd[
FrontEnd`AttachCell[Typeset`box$,
FrontEndResource["RGBColorValueSelector"], {0, {Left, Bottom}}, {
Left, Top},
"ClosingActions" -> {
"SelectionDeparture", "ParentChanged", "EvaluatorQuit"}]]]],
DefaultBaseStyle->{},
Evaluator->Automatic,
Method->"Preemptive"],
RGBColor[0.9728288904374106, 0.621644452187053, 0.07336199581899142],
Editable->False,
Selectable->False], ",",
InterpretationBox[
ButtonBox[
TooltipBox[
GraphicsBox[{
{GrayLevel[0], RectangleBox[{0, 0}]}, {GrayLevel[0], RectangleBox[{1, -1}]},
{RGBColor[0.736782672705901, 0.358, 0.5030266573755369],
RectangleBox[{0, -1}, {2, 1}]}},
AspectRatio->1,
DefaultBaseStyle->"ColorSwatchGraphics",
Frame->True,
FrameStyle->RGBColor[
0.4911884484706007, 0.23866666666666667`, 0.3353511049170246],
FrameTicks->None,
ImageSize->
Dynamic[{
Automatic, 1.35 CurrentValue["FontCapHeight"]/AbsoluteCurrentValue[
Magnification]}],
PlotRangePadding->None],
StyleBox[
RowBox[{"RGBColor", "[",
RowBox[{
"0.736782672705901`", ",", "0.358`", ",", "0.5030266573755369`"}],
"]"}], NumberMarks -> False]],
Appearance->None,
BaseStyle->{},
BaselinePosition->Baseline,
ButtonFunction:>With[{Typeset`box$ = EvaluationBox[]},
If[
Not[
AbsoluteCurrentValue["Deployed"]],
SelectionMove[Typeset`box$, All, Expression];
FrontEnd`Private`$ColorSelectorInitialAlpha = 1;
FrontEnd`Private`$ColorSelectorInitialColor =
RGBColor[0.736782672705901, 0.358, 0.5030266573755369];
FrontEnd`Private`$ColorSelectorUseMakeBoxes = True;
MathLink`CallFrontEnd[
FrontEnd`AttachCell[Typeset`box$,
FrontEndResource["RGBColorValueSelector"], {0, {Left, Bottom}}, {
Left, Top},
"ClosingActions" -> {
"SelectionDeparture", "ParentChanged", "EvaluatorQuit"}]]]],
DefaultBaseStyle->{},
Evaluator->Automatic,
Method->"Preemptive"],
RGBColor[0.736782672705901, 0.358, 0.5030266573755369],
Editable->False,
Selectable->False], ",",
InterpretationBox[
ButtonBox[
TooltipBox[
GraphicsBox[{
{GrayLevel[0], RectangleBox[{0, 0}]},
{GrayLevel[0], RectangleBox[{1, -1}]},
{RGBColor[0.28026441037696703`, 0.715, 0.4292089322474965],
RectangleBox[{0, -1}, {2, 1}]}},
AspectRatio->1,
DefaultBaseStyle->"ColorSwatchGraphics",
Frame->True,
FrameStyle->RGBColor[
0.18684294025131137`, 0.4766666666666667, 0.2861392881649977],
FrameTicks->None,
ImageSize->
Dynamic[{
Automatic, 1.35 CurrentValue["FontCapHeight"]/AbsoluteCurrentValue[
Magnification]}],
PlotRangePadding->None],
StyleBox[
RowBox[{"RGBColor", "[",
RowBox[{
"0.28026441037696703`", ",", "0.715`", ",", "0.4292089322474965`"}],
"]"}], NumberMarks -> False]],
Appearance->None,
BaseStyle->{},
BaselinePosition->Baseline, ButtonFunction:>With[{Typeset`box$ = EvaluationBox[]},
If[
Not[
AbsoluteCurrentValue["Deployed"]],
SelectionMove[Typeset`box$, All, Expression];
FrontEnd`Private`$ColorSelectorInitialAlpha = 1;
FrontEnd`Private`$ColorSelectorInitialColor =
RGBColor[0.28026441037696703`, 0.715, 0.4292089322474965];
FrontEnd`Private`$ColorSelectorUseMakeBoxes = True;
MathLink`CallFrontEnd[
FrontEnd`AttachCell[Typeset`box$,
FrontEndResource["RGBColorValueSelector"], {0, {Left, Bottom}}, {
Left, Top},
"ClosingActions" -> {
"SelectionDeparture", "ParentChanged", "EvaluatorQuit"}]]]],
DefaultBaseStyle->{},
Evaluator->Automatic,
Method->"Preemptive"],
RGBColor[0.28026441037696703`, 0.715, 0.4292089322474965],
Editable->False,
Selectable->False]}], "}"}]], "Output",
CellChangeTimes->{{3.785048598454445*^9, 3.7850486121743803`*^9},
3.7850492534864817`*^9, 3.7850518797807217`*^9, 3.7851593453111877`*^9,
3.788525055881117*^9, 3.788588576142179*^9, 3.789016362482636*^9,
3.7891011551402693`*^9, 3.7971326144003468`*^9, 3.799112667826545*^9,
3.799113935678566*^9, 3.799119920100224*^9, 3.79912099675138*^9,
3.799121678020488*^9, 3.79912196315336*^9, 3.7991220841956873`*^9, {
3.799122272245305*^9, 3.799122300399253*^9}, 3.8000746708421507`*^9,
3.8000885291890907`*^9, 3.800088987611136*^9, {3.800089240452045*^9,
3.800089269881946*^9}, 3.800089321591539*^9, 3.800089431077919*^9,
3.800089495187615*^9, 3.800089595616612*^9, 3.800148380233251*^9,
3.800539327345366*^9, 3.800610134971425*^9, 3.800705385664913*^9,
3.80070556987372*^9, 3.800768504913447*^9, 3.800798370679039*^9,
3.8008772883571987`*^9, 3.801281244944542*^9, 3.801285594509903*^9,
3.801567580742796*^9, {3.801723018647998*^9, 3.801723046894418*^9},
3.801736561561638*^9, 3.801738393523629*^9, 3.801744291978118*^9,
3.8017504625330067`*^9, 3.801754015153778*^9, 3.801754114168932*^9,
3.8018007150796423`*^9, 3.8018040233399773`*^9, 3.801882349914794*^9,
3.8018951501756773`*^9, {3.80190087812988*^9, 3.801900885048917*^9}},
CellLabel->
"Out[103]=",ExpressionUUID->"d642568b-341d-4984-ae2d-30a8eeb30577"]
}, Open ]]
}, Closed]],
Cell[CellGroupData[{
Cell["Other functions", "Subsection",
CellChangeTimes->{{3.8017505701278543`*^9,
3.801750573727201*^9}},ExpressionUUID->"44e4b7cd-d026-4079-8fd8-\
3f63cf14cfa7"],
Cell[BoxData[{
RowBox[{
RowBox[{
RowBox[{"Options", "[", "pdf", "]"}], "=",
RowBox[{"{",
RowBox[{"\"\<SNR\>\"", "\[Rule]", "1"}], "}"}]}], ";"}], "\n",
RowBox[{
RowBox[{"pdf", "[",
RowBox[{"ansatz_", ",", "data_", ",",
RowBox[{"tend_", "?", "NumericQ"}], ",",
RowBox[{"x_", "?",
RowBox[{"(",
RowBox[{
RowBox[{"VectorQ", "[",
RowBox[{"#", ",", "NumericQ"}], "]"}], "&"}], ")"}]}], ",",
RowBox[{"OptionsPattern", "[", "]"}]}], "]"}], ":=",
RowBox[{"Block", "[",
RowBox[{
RowBox[{"{", RowBox[{
"ansl", ",", "cfit", ",", "cfitd", ",", "h1red", ",", "h2red", ",",
"norm1", ",", "norm2", ",", "myTable", ",", "snr", ",", "tmax"}], "}"}],
",", "\[IndentingNewLine]", "\n",
RowBox[{
RowBox[{"snr", "=",
RowBox[{"OptionValue", "[", "\"\<SNR\>\"", "]"}]}], ";", "\n",
"\[IndentingNewLine]",
RowBox[{"tmax", "=",
RowBox[{"data", "[",
RowBox[{"[",
RowBox[{"1", ",", "1"}], "]"}], "]"}]}], ";", "\n",
RowBox[{"ansl", "=", "ansatz"}], ";", "\[IndentingNewLine]",
RowBox[{"cfitd", "=",
RowBox[{"Transpose", "[",
RowBox[{"{",
RowBox[{
RowBox[{"data", "[",
RowBox[{"[",
RowBox[{"All", ",", "1"}], "]"}], "]"}], ",",
RowBox[{"ansl", "/.",
RowBox[{"t", "->",
RowBox[{"data", "[",
RowBox[{"[",
RowBox[{"All", ",", "1"}], "]"}], "]"}]}]}]}], "}"}], "]"}]}],
";", "\[IndentingNewLine]", "\[IndentingNewLine]",
RowBox[{"h1red", "=",
RowBox[{
RowBox[{"Select", "[",
RowBox[{"data", ",",
RowBox[{
RowBox[{"tmax", "<=",
RowBox[{"#", "[",
RowBox[{"[", "1", "]"}], "]"}], "<=",
RowBox[{"tmax", "+", "tend"}]}], "&"}]}], "]"}], "[",
RowBox[{"[",
RowBox[{"All", ",", "2"}], "]"}], "]"}]}], ";", "\[IndentingNewLine]",
RowBox[{"h2red", "=",
RowBox[{
RowBox[{"Select", "[",
RowBox[{"cfitd", ",",
RowBox[{
RowBox[{"tmax", "<=",
RowBox[{"#", "[",
RowBox[{"[", "1", "]"}], "]"}], "<=",
RowBox[{"tmax", "+", "tend"}]}], "&"}]}], "]"}], "[",
RowBox[{"[",
RowBox[{"All", ",", "2"}], "]"}], "]"}]}], ";", "\[IndentingNewLine]",
RowBox[{"norm1", "=",
RowBox[{"Total", "[",
RowBox[{
RowBox[{"(",
RowBox[{"Abs", "@", "h1red"}], ")"}], "^", "2"}], " ", "]"}]}], ";",
"\n",
RowBox[{"norm2", "=",
RowBox[{"Total", "[",
RowBox[{
RowBox[{"(",
RowBox[{"Abs", "@", "h2red"}], ")"}], "^", "2"}], " ", "]"}]}], ";",
"\n", "\n",
RowBox[{"myTable", "=",
RowBox[{"h1red", " ",
RowBox[{"Conjugate", "@", "h2red"}]}]}], ";", "\n",
RowBox[{
RowBox[{"-",
RowBox[{"snr", "^", "2"}]}],
RowBox[{"(",
RowBox[{"1", "-",
RowBox[{"(",
RowBox[{ RowBox[{"Re", "@",
RowBox[{"Total", "[", "myTable", "]"}]}], "/",
RowBox[{"Sqrt", "[",
RowBox[{"norm1", " ", "norm2"}], "]"}]}], ")"}]}], ")"}]}]}]}],
"]"}]}]}], "Code",
CellChangeTimes->{{3.777193638014988*^9, 3.7771937755781918`*^9},
3.777194208058053*^9, {3.7771943136764383`*^9, 3.777194320159019*^9}, {
3.777194706435961*^9, 3.77719471195916*^9}, 3.777194844196529*^9,
3.777194878365808*^9, {3.80175364495243*^9, 3.801753680033052*^9}, {
3.801753838001347*^9, 3.801753865156281*^9}, {3.801753993285178*^9,
3.801753994197847*^9}, {3.801754228905569*^9, 3.801754233024242*^9}, {
3.801754997079472*^9, 3.801755001571432*^9}, {3.801755848922812*^9,
3.8017558520508947`*^9}, {3.801800854608481*^9, 3.801800861689166*^9}, {
3.801803635098641*^9, 3.801803732376952*^9}, {3.80180382005339*^9,
3.8018038226174717`*^9}, {3.801803864490718*^9, 3.801803875455367*^9}, {
3.801803950318574*^9, 3.801804017270652*^9}},
CellLabel->
"In[107]:=",ExpressionUUID->"48958061-652d-46b3-9aca-2ec5e658ec2a"],
Cell[BoxData[{
RowBox[{
RowBox[{
RowBox[{"Options", "[", "pdfFit", "]"}], "=",
RowBox[{"{",
RowBox[{"\"\<SNR\>\"", "\[Rule]", "1"}], "}"}]}], ";"}], "\n",
RowBox[{
RowBox[{"pdfFit", "[",
RowBox[{"ansatz_", ",", "data_", ",",
RowBox[{"tend_", "?", "NumericQ"}], ",",
RowBox[{"x_", "?",
RowBox[{"(",
RowBox[{
RowBox[{"VectorQ", "[",
RowBox[{"#", ",", "NumericQ"}], "]"}], "&"}], ")"}]}], ",",
RowBox[{"OptionsPattern", "[", "]"}]}], "]"}], ":=",
RowBox[{"Block", "[",
RowBox[{
RowBox[{"{",
RowBox[{
"ansl", ",", "cfit", ",", "cfitd", ",", "h1red", ",", "h2red", ",",
"norm1", ",", "norm2", ",", "myTable", ",", "snr", ",", "tmax"}], "}"}],
",", "\[IndentingNewLine]", "\n",
RowBox[{
RowBox[{"snr", "=",
RowBox[{"OptionValue", "[", "\"\<SNR\>\"", "]"}]}], ";", "\n",
"\[IndentingNewLine]",
RowBox[{"tmax", "=",
RowBox[{"data", "[",
RowBox[{"[",
RowBox[{"1", ",", "1"}], "]"}], "]"}]}], ";", "\n",
RowBox[{"ansl", "=", "ansatz"}], ";", "\[IndentingNewLine]",
RowBox[{"cfitd", "=",
RowBox[{"Transpose", "[",
RowBox[{"{",
RowBox[{
RowBox[{"data", "[",
RowBox[{"[",
RowBox[{"All", ",", "1"}], "]"}], "]"}], ",",
RowBox[{"ansl", "/.",
RowBox[{"t", "->",
RowBox[{"data", "[",
RowBox[{"[",
RowBox[{"All", ",", "1"}], "]"}], "]"}]}]}]}], "}"}], "]"}]}],
";", "\[IndentingNewLine]",
RowBox[{"h1red", "=",
RowBox[{
RowBox[{"Select", "[",
RowBox[{"data", ",",
RowBox[{
RowBox[{"tmax", "<=", RowBox[{"#", "[",
RowBox[{"[", "1", "]"}], "]"}], "<=",
RowBox[{"tmax", "+", "tend"}]}], "&"}]}], "]"}], "[",
RowBox[{"[",
RowBox[{"All", ",", "2"}], "]"}], "]"}]}], ";", "\[IndentingNewLine]",
RowBox[{"h2red", "=",
RowBox[{
RowBox[{"Select", "[",
RowBox[{"cfitd", ",",
RowBox[{
RowBox[{"tmax", "<=",
RowBox[{"#", "[",
RowBox[{"[", "1", "]"}], "]"}], "<=",
RowBox[{"tmax", "+", "tend"}]}], "&"}]}], "]"}], "[",
RowBox[{"[",
RowBox[{"All", ",", "2"}], "]"}], "]"}]}], ";", "\[IndentingNewLine]",
"\n",
RowBox[{"Total", "[",
RowBox[{"(",
RowBox[{
RowBox[{
RowBox[{"Re", "[",
RowBox[{"h1red", "-", "h2red"}], "]"}], "^", "2"}], "+",
RowBox[{
RowBox[{"Im", "[",
RowBox[{"h1red", "-", "h2red"}], "]"}], "^", "2"}]}], ")"}],
"]"}]}]}], "\n", "]"}]}]}], "Code",
CellChangeTimes->{{3.777193638014988*^9, 3.7771937755781918`*^9},
3.777194208058053*^9, {3.7771943136764383`*^9, 3.777194320159019*^9}, {
3.777194706435961*^9, 3.77719471195916*^9}, 3.777194844196529*^9,
3.777194878365808*^9, {3.80175364495243*^9, 3.801753680033052*^9}, {
3.801753838001347*^9, 3.801753865156281*^9}, {3.801753993285178*^9,
3.801753994197847*^9}, {3.801754228905569*^9, 3.801754233024242*^9}, {
3.801754997079472*^9, 3.801755001571432*^9}, {3.801755848922812*^9,
3.8017558520508947`*^9}, {3.801800854608481*^9, 3.801800861689166*^9}, {
3.801803635098641*^9, 3.801803732376952*^9}, {3.80180382005339*^9,
3.8018038226174717`*^9}, {3.801803864490718*^9, 3.801803875455367*^9}, {
3.801803950318574*^9, 3.801804017270652*^9}, 3.8018087211731253`*^9, {
3.8018087834835567`*^9, 3.801808901759162*^9}},
CellLabel->
"In[109]:=",ExpressionUUID->"0204b879-5f51-44f5-8a4f-95cc490da72f"]
}, Closed]],
Cell[CellGroupData[{
Cell["Imports", "Subsection",
CellChangeTimes->{{3.7991159558445053`*^9,
3.799115957130929*^9}},ExpressionUUID->"8488942b-dc04-47e9-8ca6-\
119712f58bc1"],
Cell[CellGroupData[{
Cell["RDown data", "Subsubsection",
CellChangeTimes->{{3.801736631218856*^9,
3.8017366327690687`*^9}},ExpressionUUID->"f0a39a1f-5153-4be0-b03d-\
fd24806eca3d"],
Cell[BoxData[
RowBox[{"(*", " ",
RowBox[{
RowBox[{
RowBox[{
"The", " ", "QNM", " ", "spectrum", " ", "is", " ", "obtained", " ",
"from", " ",
RowBox[{"https", ":"}]}], "//",
RowBox[{
RowBox[{
RowBox[{
RowBox[{"pages", ".", "jh", ".", "edu"}], "/",
RowBox[{"~", "eberti2"}]}], "/", "ringdown"}], "/.", " ", RowBox[{"Up", " ", "to", " ", "now"}]}]}], ",", " ",
RowBox[{
RowBox[{
"we", " ", "can", " ", "load", " ", "any", " ", "lmn", " ", "combination",
" ", "up", " ", "to", " ", "l"}], "=",
RowBox[{
RowBox[{"4", " ", "and", " ", "n"}], "=",
RowBox[{
RowBox[{
"7.", " ", "To", " ", "extend", " ", "the", " ", "l", " ", "number",
" ", "you", " ", "need", " ", "to", " ", "download", " ", "the", " ",
"l"}], ">",
RowBox[{"4", " ",
RowBox[{"files", "."}]}]}]}]}]}], " ", "*)"}]], "Input",
CellChangeTimes->{{3.801737163175469*^9, 3.801737257507339*^9}},
CellLabel->
"In[159]:=",ExpressionUUID->"7694c13e-a765-44ce-9adc-09e59459c3d2"],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{
RowBox[{"(*", " ",
RowBox[{
"Import", " ", "all", " ", "the", " ", "modes", " ", "you", " ", "want",
" ", "to", " ", "use"}], " ", "*)"}], "\n",
RowBox[{"modelist", "=",
RowBox[{"{",
RowBox[{
RowBox[{"{",
RowBox[{"2", ",", "2", ",", "0"}], "}"}], ",",
RowBox[{"{",
RowBox[{"2", ",", "2", ",", "1"}], "}"}], ",",
RowBox[{"{",
RowBox[{"2", ",", "2", ",", "2"}], "}"}], ",",
RowBox[{"{",
RowBox[{"2", ",", "2", ",", "3"}], "}"}], ",",
RowBox[{"{",
RowBox[{"2", ",", "2", ",", "4"}], "}"}], ",",
RowBox[{"{",
RowBox[{"2", ",", "2", ",", "5"}], "}"}], ",",
RowBox[{"{",
RowBox[{"2", ",", "2", ",", "6"}], "}"}], ",",
RowBox[{"{",
RowBox[{"2", ",", "2", ",", "7"}], "}"}]}], "}"}]}]}]], "Code",
CellChangeTimes->{{3.799116007730418*^9, 3.7991160596093693`*^9}, {
3.80173672888538*^9, 3.801736729854124*^9}, {3.801737022977298*^9,
3.8017370304868717`*^9}, {3.801741821615449*^9, 3.801741825149316*^9}, {
3.8018891925929213`*^9, 3.801889199423003*^9}, {3.801906226368847*^9,
3.801906250821456*^9}},
CellLabel->
"In[171]:=",ExpressionUUID->"763f62ae-036b-47b4-86e4-d1b0e24fc92b"],
Cell[BoxData[
RowBox[{"{",
RowBox[{
RowBox[{"{",
RowBox[{"2", ",", "2", ",", "0"}], "}"}], ",",
RowBox[{"{",
RowBox[{"2", ",", "2", ",", "1"}], "}"}], ",",
RowBox[{"{",
RowBox[{"2", ",", "2", ",", "2"}], "}"}], ",",
RowBox[{"{",
RowBox[{"2", ",", "2", ",", "3"}], "}"}], ",",
RowBox[{"{",
RowBox[{"2", ",", "2", ",", "4"}], "}"}], ",",
RowBox[{"{",
RowBox[{"2", ",", "2", ",", "5"}], "}"}], ",",
RowBox[{"{",
RowBox[{"2", ",", "2", ",", "6"}], "}"}], ",", RowBox[{"{",
RowBox[{"2", ",", "2", ",", "7"}], "}"}]}], "}"}]], "Output",
CellChangeTimes->{
3.801737030858226*^9, 3.801738395377088*^9, {3.801741822713435*^9,
3.801741825490649*^9}, 3.801744292163055*^9, 3.801750462612952*^9,
3.801754015269165*^9, 3.80175411432535*^9, 3.8018007152800283`*^9,
3.801804024324876*^9, 3.801882350073197*^9, 3.801889169966961*^9,
3.801889200128565*^9, 3.8018951503870993`*^9, 3.801900897849856*^9,
3.801906099558012*^9, 3.801906251101836*^9},
CellLabel->
"Out[171]=",ExpressionUUID->"81efd3bd-daab-4062-8b73-d2a7435e97df"]
}, Open ]],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{
RowBox[{"(*", " ",
RowBox[{
RowBox[{
RowBox[{"Files", " ", "for", " ", "the", " ",
RowBox[{"modes", ".", " ", "Note"}], " ", "that", " ", "in", " ",
"this", " ", "file", " ", "notation", " ", "the", " ", "fundamental",
" ", "tone", " ", "corresponds", " ", "to", " ", "n"}], "=",
RowBox[{"1.", " ", "However"}]}], ",", " ",
RowBox[{
RowBox[{
"for", " ", "the", " ", "code", " ", "notation", " ", "we", " ", "will",
" ", "always", " ", "use", " ", "n"}], "=", "0."}]}], " ", "*)"}],
"\n",
RowBox[{"Modefiles", "=",
RowBox[{"Table", "[",
RowBox[{
RowBox[{"modedir", "<>", "\"\<l\>\"", "<>",
RowBox[{"ToString", "[",
RowBox[{"modelist", "[",
RowBox[{"[",
RowBox[{"j", ",", "1"}], "]"}], "]"}], "]"}], "<>", "\"\</n\>\"", "<>",
RowBox[{"ToString", "[",
RowBox[{
RowBox[{"modelist", "[",
RowBox[{"[",
RowBox[{"j", ",", "3"}], "]"}], "]"}], "+", "1"}], "]"}], "<>",
"\"\<l\>\"", "<>",
RowBox[{"ToString", "[",
RowBox[{"modelist", "[",
RowBox[{"[",
RowBox[{"j", ",", "1"}], "]"}], "]"}], "]"}], "<>", "\"\<m\>\"", "<>",
RowBox[{"ToString", "[",
RowBox[{"modelist", "[",
RowBox[{"[",
RowBox[{"j", ",", "2"}], "]"}], "]"}], "]"}], "<>", "\"\<.dat\>\""}],
",",
RowBox[{"{",
RowBox[{"j", ",",
RowBox[{"Length", "@", "modelist"}]}], "}"}]}], "]"}]}]}]], "Code",
CellChangeTimes->{{3.799116066718706*^9, 3.799116183969667*^9}, {
3.799116221329751*^9, 3.799116223632145*^9}, {3.7991162786627827`*^9,
3.799116285885654*^9}, {3.7991163534374104`*^9, 3.79911639784995*^9}, {
3.799116694401691*^9, 3.7991167077116737`*^9}, 3.8007687513426323`*^9, {
3.800768840870967*^9, 3.800768851624331*^9}, 3.801736580303645*^9, {
3.8017367250936527`*^9, 3.8017367258167553`*^9}, {3.801737036972227*^9,
3.8017371320562687`*^9}, {3.8017418292956*^9, 3.8017418338861313`*^9},
3.801741882215084*^9},
CellLabel->
"In[172]:=",ExpressionUUID->"37d6369f-02c3-4af4-91b6-0490e703d64a"],
Cell[BoxData[
RowBox[{"{",
RowBox[{"\<\"~/git/rdstackingproject/l2/n1l2m2.dat\"\>", ",", "\<\"~/git/rdstackingproject/l2/n2l2m2.dat\"\>",
",", "\<\"~/git/rdstackingproject/l2/n3l2m2.dat\"\>",
",", "\<\"~/git/rdstackingproject/l2/n4l2m2.dat\"\>",
",", "\<\"~/git/rdstackingproject/l2/n5l2m2.dat\"\>",
",", "\<\"~/git/rdstackingproject/l2/n6l2m2.dat\"\>",
",", "\<\"~/git/rdstackingproject/l2/n7l2m2.dat\"\>",
",", "\<\"~/git/rdstackingproject/l2/n8l2m2.dat\"\>"}], "}"}]], "Output",
CellChangeTimes->{
3.799116184826573*^9, 3.799116224061639*^9, 3.7991162873873873`*^9,
3.799116355402153*^9, 3.799116398123296*^9, {3.799116696222416*^9,
3.799116708064439*^9}, 3.799119929262824*^9, 3.799121680129674*^9,
3.7991219634626427`*^9, 3.799122084375304*^9, {3.799122272555999*^9,
3.799122300874164*^9}, 3.800074671099461*^9, 3.8000885294715147`*^9,
3.8000889878447647`*^9, {3.800089240559409*^9, 3.800089270099875*^9},
3.8000893217612467`*^9, 3.800089431239541*^9, 3.800089495352811*^9,
3.800089595851823*^9, 3.800148405206428*^9, 3.800539327493127*^9,
3.800610135108379*^9, 3.8007055994616632`*^9, 3.8007685050622683`*^9,
3.800768736245232*^9, {3.800768841340745*^9, 3.800768873828692*^9},
3.800798370835478*^9, 3.800877288561569*^9, 3.801281245216522*^9,
3.801285597638463*^9, 3.80156758216507*^9, {3.801723018844685*^9,
3.8017230487687283`*^9}, {3.801736575720043*^9, 3.801736580659892*^9},
3.80173673506961*^9, 3.801737031790604*^9, 3.8017371395428753`*^9,
3.801738396238654*^9, 3.801741884414989*^9, 3.801744292170224*^9,
3.801750462618342*^9, 3.801754015286237*^9, 3.801754114331324*^9,
3.801800715386989*^9, 3.8018040243921223`*^9, 3.8018823500782957`*^9, {
3.801889182274848*^9, 3.801889203976694*^9}, 3.801895150442212*^9,
3.801900897858267*^9, 3.801906099628516*^9, 3.8019062528095083`*^9},
CellLabel->
"Out[172]=",ExpressionUUID->"1fd520ec-5fee-43d1-967d-a3a48f981509"]
}, Open ]],
Cell[BoxData[
RowBox[{
RowBox[{"(*", " ",
RowBox[{"Import", " ", "modes", " ", "data"}], " ", "*)"}], "\n",
RowBox[{
RowBox[{"Modedata", "=",
RowBox[{"Import", "/@", "Modefiles"}]}], ";"}]}]], "Code",
CellChangeTimes->{{3.799117127671782*^9, 3.799117133283991*^9}, {
3.801736719816572*^9, 3.801736720908989*^9}},
CellLabel->
"In[173]:=",ExpressionUUID->"827b1be5-400a-4c10-97cf-df0471536b52"]
}, Closed]],
Cell[CellGroupData[{
Cell["NR data + EOBFits", "Subsubsection",
CellChangeTimes->{{3.762592741285844*^9,
3.762592747256867*^9}},ExpressionUUID->"e69b91c9-6c5f-47bb-9540-\
07f42335cb1a"],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"modes", "=",
RowBox[{"{",
RowBox[{"{",
RowBox[{"2", ",", "2"}], "}"}], "}"}]}]], "Code",
CellChangeTimes->{{3.71721646622243*^9, 3.717216472004093*^9}, {
3.7172238332578897`*^9, 3.7172238699527483`*^9}, {3.7172239989120007`*^9,
3.717223999908312*^9}, {3.717230827569207*^9, 3.7172308284476624`*^9}, {
3.717233966560224*^9, 3.717233966743574*^9}, 3.762057923479889*^9, {
3.762681745081925*^9, 3.7626817794296007`*^9}, {3.7991133898969717`*^9,
3.799113390253549*^9}, {3.799113545789241*^9, 3.799113608581436*^9},
3.8018995918649387`*^9, {3.8019009292101917`*^9, 3.8019009565374517`*^9}, {
3.80190121855807*^9, 3.8019012197457848`*^9}},
CellLabel->
"In[130]:=",ExpressionUUID->"b361f1b5-e01e-4264-bf04-6b665f472fd9"],
Cell[BoxData[ RowBox[{"{",
RowBox[{"{",
RowBox[{"2", ",", "2"}], "}"}], "}"}]], "Output",
CellChangeTimes->{{3.7172238612766*^9, 3.717223870939024*^9},
3.717224000107691*^9, 3.717230828736228*^9, 3.717233929169869*^9,
3.7172339670405207`*^9, 3.717234264957567*^9, 3.7173898088553953`*^9,
3.717397795928178*^9, 3.7174090874909067`*^9, 3.719042266653441*^9,
3.719725716094987*^9, 3.7204979816442013`*^9, 3.720601795601798*^9,
3.7620579265503798`*^9, 3.762171305194047*^9, 3.762185874430632*^9,
3.762230695496544*^9, 3.7625796880016623`*^9, 3.762602459472851*^9,
3.7626038253288*^9, {3.762681745570024*^9, 3.762681784730443*^9},
3.7643108351004267`*^9, 3.764392213490958*^9, 3.764406815613749*^9,
3.7644071500828*^9, 3.764416940095358*^9, 3.764471032983673*^9,
3.764502309669742*^9, 3.76451318465303*^9, 3.764567390742053*^9,
3.768568010349392*^9, 3.768572297815539*^9, 3.768645827003063*^9, {
3.768649825183477*^9, 3.768649837869774*^9}, 3.768701845001706*^9,
3.768706532972953*^9, 3.768882352986763*^9, 3.768909643770508*^9,
3.7691352413787107`*^9, 3.769137120092252*^9, 3.769138853253389*^9,
3.7691738608203077`*^9, 3.769228308058159*^9, 3.769488529403152*^9,
3.7694887741511297`*^9, 3.769754030511262*^9, 3.769849331321813*^9,
3.769927143176673*^9, 3.7699298467484093`*^9, {3.769934942968072*^9,
3.7699349540967703`*^9}, 3.770001458398054*^9, 3.770089588112008*^9,
3.7700911499652567`*^9, 3.770091211153775*^9, 3.770096723859185*^9,
3.770175902879304*^9, 3.770355811386429*^9, 3.770368573284638*^9,
3.770431012241467*^9, 3.770452812170587*^9, 3.7712189390391197`*^9,
3.7716732302402697`*^9, 3.7719070307003307`*^9, 3.771907438653419*^9,
3.771915621511821*^9, 3.772190839282982*^9, 3.7721965104057217`*^9,
3.772247855130703*^9, 3.772256628537322*^9, 3.772259081090678*^9,
3.772348454834467*^9, 3.772367800626193*^9, 3.7723708897053223`*^9,
3.77241999957555*^9, 3.772428158523457*^9, 3.7724306638823853`*^9,
3.7724404539315863`*^9, 3.772529270837591*^9, 3.77287591387982*^9,
3.772878049981765*^9, 3.772885185257471*^9, 3.772942819052163*^9,
3.772943402403473*^9, 3.772943541727195*^9, 3.772944053630115*^9,
3.772949568223918*^9, 3.772950462032238*^9, 3.7729505152950287`*^9,
3.7730492848622227`*^9, 3.7730496269622517`*^9, 3.773144075606292*^9,
3.773460623067893*^9, 3.773463998185182*^9, 3.773471311405634*^9,
3.7734788836585703`*^9, 3.77347963931531*^9, 3.773481967361264*^9,
3.773645277259672*^9, 3.77364680378931*^9, 3.773986244987804*^9,
3.776510085446184*^9, 3.7765790810668097`*^9, 3.776662520028639*^9, {
3.776668739189172*^9, 3.7766687414143867`*^9}, 3.7766699887032843`*^9,
3.77667547130352*^9, 3.776676365915367*^9, 3.776677554777287*^9,
3.7766781567455072`*^9, 3.776679727486134*^9, 3.7766797993092937`*^9,
3.776679922316041*^9, 3.7766814663166523`*^9, 3.776762885603821*^9,
3.777024814636004*^9, 3.777117695799756*^9, 3.777184967225831*^9,
3.777625642290798*^9, 3.7776907191949377`*^9, 3.777725375557233*^9,
3.777786888566283*^9, 3.777953025346758*^9, 3.777964327509592*^9,
3.778220763424073*^9, 3.778241132739333*^9, 3.778242215252248*^9,
3.7783898515967083`*^9, 3.7784896886312447`*^9, 3.77883439339159*^9,
3.77891303655623*^9, 3.778918648410945*^9, 3.7789996814424963`*^9,
3.779014946820965*^9, 3.779016358474268*^9, 3.779091498367682*^9,
3.7791090633146276`*^9, 3.7794380428020563`*^9, 3.7794390353822803`*^9,
3.779446129755233*^9, 3.779512776856132*^9, 3.7795447057568073`*^9,
3.7795459143084173`*^9, 3.7796202554136963`*^9, 3.779967327286769*^9,
3.78005038675193*^9, 3.780061151611998*^9, 3.780115034323454*^9,
3.780115441240199*^9, {3.780997648127914*^9, 3.78099766425393*^9},
3.7824639793355637`*^9, 3.78394716725804*^9, 3.785049266720189*^9,
3.785052848198894*^9, 3.785057075575471*^9, 3.785159356530836*^9,
3.788525065848336*^9, 3.788588588179237*^9, 3.789016374435004*^9,
3.7891011634419317`*^9, 3.797132624620028*^9, 3.797135025936468*^9,
3.7991133916199923`*^9, 3.799113509759405*^9, {3.7991135463122883`*^9,
3.799113551635894*^9}, {3.799113586080545*^9, 3.799113610336619*^9},
3.799113982242723*^9, 3.799119936454014*^9, 3.799121028438673*^9,
3.799121682590557*^9, 3.799121965014926*^9, 3.799122085602578*^9, {
3.7991222752168427`*^9, 3.799122302495413*^9}, 3.8000746749992437`*^9,
3.800088876294077*^9, 3.800148382767583*^9, {3.800539328605701*^9,
3.8005393463625402`*^9}, 3.800610204606415*^9, 3.800768510807893*^9,
3.800768557054263*^9, 3.8007983952296457`*^9, 3.800836721589457*^9,
3.801289306841797*^9, 3.801723120829804*^9, 3.8018052338776903`*^9,
3.80190122034149*^9},
CellLabel-> "Out[130]=",ExpressionUUID->"5a59dd1d-dcee-4d7a-954f-f98d277d06f5"]
}, Open ]],
Cell[CellGroupData[{
Cell[BoxData[{
RowBox[{"mysxscaserh", "=",
RowBox[{"Reverse", "[",
RowBox[{"Select", "[",
RowBox[{
RowBox[{"FileNames", "[",
RowBox[{"\"\<*Lev6/rh*\>\"", ",", "mysxscase", ",", "2"}], "]"}], ",",
RowBox[{
RowBox[{"Not", "@",
RowBox[{"StringMatchQ", "[",
RowBox[{"#", ",", "\"\<*raw*\>\""}], "]"}]}], "&"}]}], "]"}],
"]"}]}], "\[IndentingNewLine]",
RowBox[{"mysxscasemetafile", "=",
RowBox[{
RowBox[{"FileNames", "[",
RowBox[{"\"\<metadata.txt\>\"", ",", "mysxscase", ",", "2"}], "]"}], "[",
RowBox[{"[", "1", "]"}], "]"}]}], "\n",
RowBox[{
RowBox[{"metadata", "=",
RowBox[{"SXSMetaFilesToRules", "[", "mysxscasemetafile", "]"}]}],
";"}]}], "Code",
CellChangeTimes->{
3.71714565326084*^9, {3.717216492175115*^9, 3.71721650261952*^9}, {
3.717223876329924*^9, 3.717223876601008*^9}, {3.717233944513124*^9,
3.717233958472711*^9}, {3.77000148377827*^9, 3.770001483895461*^9}, {
3.770089594441771*^9, 3.7700896091594048`*^9}, {3.770091181515603*^9,
3.770091181662176*^9}, {3.770091214426573*^9, 3.770091214585648*^9}, {
3.770092858390904*^9, 3.770092858508835*^9}, {3.770096730341798*^9,
3.770096740949168*^9}, {3.772428168100793*^9, 3.772428170608094*^9}, {
3.772428408521916*^9, 3.7724284091136923`*^9}, {3.772430667886035*^9,
3.772430670003887*^9}, {3.772885188569561*^9, 3.772885188844552*^9}, {
3.77294957315034*^9, 3.772949625242341*^9}, {3.780998797734757*^9,
3.780998797920843*^9}, {3.7991136519052963`*^9, 3.799113707640517*^9}},
CellLabel->
"In[131]:=",ExpressionUUID->"5584846d-d985-4ced-98a1-ac35c86181ea"],
Cell[BoxData[
RowBox[{"{", \
"\<\"/Users/xisco/SXS/BBH_SKS_d14.3_q1.22_sA_0_0_0.330_sB_0_0_-0.440/Lev6/\
rhOverM_Asymptotic_GeometricUnits_CoM.h5\"\>", "}"}]], "Output",
CellChangeTimes->{{3.717216498380938*^9, 3.7172165029748363`*^9},
3.7172194183959084`*^9, 3.717219451757556*^9, 3.7172210443770657`*^9,
3.717221284233238*^9, 3.717223885556622*^9, 3.71723083183225*^9,
3.7172308731861277`*^9, {3.717233930309966*^9, 3.717233959115399*^9},
3.717234265666395*^9, 3.717389810564118*^9, 3.717397796938533*^9,
3.717409087835437*^9, 3.719042267675703*^9, 3.719725716500965*^9,
3.720497982239339*^9, 3.720601796148543*^9, 3.7620579313020773`*^9,
3.7621713078981047`*^9, 3.762185875336357*^9, 3.762230696430828*^9,
3.762579690340506*^9, 3.762579774070689*^9, 3.76260245990764*^9,
3.762603825455809*^9, 3.762681789861004*^9, 3.762681826567923*^9,
3.764310836001659*^9, 3.76439221450559*^9, 3.764406818443306*^9,
3.7644071506895514`*^9, 3.7644169407193203`*^9, 3.764471034222904*^9,
3.7645023104911823`*^9, 3.7645131854137583`*^9, 3.764567390841213*^9,
3.7685680107885113`*^9, 3.768572298355291*^9, 3.768645827647299*^9,
3.768649839712215*^9, 3.76870184589828*^9, 3.76870653364054*^9,
3.7688823539170313`*^9, 3.768909643980805*^9, 3.769135241869529*^9,
3.7691371202865047`*^9, 3.769138853846279*^9, 3.7691738611067*^9,
3.7692283084744673`*^9, 3.769488530350202*^9, 3.769488774949139*^9,
3.769754030864167*^9, 3.769849332079545*^9, 3.769927143513033*^9,
3.769929847048765*^9, 3.7699349547247353`*^9, {3.7700014592708473`*^9,
3.770001484137087*^9}, {3.770089589019122*^9, 3.770089614454945*^9},
3.770091150926076*^9, {3.7700911818740463`*^9, 3.770091214764297*^9},
3.770092858828313*^9, {3.770096724350412*^9, 3.770096741200189*^9},
3.770175903073946*^9, 3.770355811577338*^9, 3.770368573479495*^9, 3.77043101246592*^9, 3.770452812583289*^9, 3.771218939593539*^9,
3.771673231116873*^9, 3.771907031257121*^9, 3.771907438848817*^9,
3.771915621797638*^9, 3.77219083949329*^9, 3.77219651074804*^9,
3.7722478558236303`*^9, 3.772256628919487*^9, 3.7722590823521147`*^9,
3.772348455019628*^9, 3.772367801149004*^9, 3.7723708907016478`*^9, {
3.7724200004504337`*^9, 3.7724200042661133`*^9}, {3.7724281623433*^9,
3.7724281708790503`*^9}, 3.772428409531721*^9, {3.7724306652135067`*^9,
3.772430670282818*^9}, {3.772440455361487*^9, 3.772440460876882*^9},
3.772529271686309*^9, 3.772875914079034*^9, 3.772878050591577*^9,
3.772885194256702*^9, 3.7729428221831827`*^9, 3.772943402612479*^9,
3.772943541918996*^9, 3.77294405673494*^9, {3.772949569772958*^9,
3.7729496256799297`*^9}, 3.7729504628386717`*^9, 3.77295051555145*^9,
3.7730492856798477`*^9, {3.773049628186563*^9, 3.773049630622552*^9},
3.773144075790704*^9, 3.773460623262436*^9, 3.7734639987477903`*^9,
3.773471312015882*^9, 3.7734788840985126`*^9, 3.773479639599619*^9,
3.773481967713728*^9, 3.773645277444871*^9, 3.773646804195562*^9,
3.77398624805871*^9, 3.776510090600584*^9, 3.7765790822290382`*^9,
3.776662520362281*^9, 3.776668741495263*^9, 3.776669989208765*^9,
3.7766754714003077`*^9, 3.776676366016034*^9, 3.776677554881569*^9,
3.776678156846538*^9, 3.776679728131668*^9, 3.776679799421623*^9,
3.776679922424872*^9, 3.7766814664267282`*^9, 3.776762885737376*^9,
3.777024814804093*^9, 3.777117695893361*^9, 3.777184967307158*^9,
3.77762564240553*^9, 3.7776907196567163`*^9, 3.777725375935665*^9,
3.777786888657611*^9, 3.777953025450624*^9, 3.7779643276193657`*^9,
3.778220763522934*^9, 3.778241132862837*^9, 3.778242215351343*^9,
3.778389851692855*^9, 3.778489688754407*^9, 3.778834393505157*^9,
3.778913036673313*^9, 3.778918648505211*^9, 3.778999681983601*^9,
3.779014946982914*^9, 3.779016358626405*^9, 3.779091498506896*^9,
3.779109063467824*^9, 3.779438042912994*^9, 3.779439035500759*^9,
3.7794461298809433`*^9, 3.77951277696731*^9, 3.779544706262*^9,
3.7795459144103613`*^9, 3.7796202555125847`*^9, 3.779967327387167*^9,
3.780050387061942*^9, 3.780061151710363*^9, 3.780115034431876*^9,
3.7801154413407516`*^9, {3.78099764910747*^9, 3.7809976642917347`*^9},
3.780997931988578*^9, 3.7824639802216053`*^9, 3.783947167922887*^9,
3.785049267615994*^9, 3.785052848818983*^9, 3.785057076455508*^9,
3.785159356717268*^9, 3.7885250658656263`*^9, 3.788588588281761*^9,
3.789016374469722*^9, 3.789101163461742*^9, 3.797132624634801*^9,
3.79713502596594*^9, 3.799113592904676*^9, {3.7991136528442707`*^9,
3.799113708088469*^9}, 3.7991139845245934`*^9, 3.799119938366549*^9,
3.799121028590159*^9, 3.7991216828410063`*^9, 3.7991219650292177`*^9,
3.7991220856243143`*^9, {3.79912227538988*^9, 3.79912230259875*^9},
3.8000746751750097`*^9, 3.800088876448193*^9, 3.800148382940559*^9, {
3.8005393286214533`*^9, 3.800539348433199*^9}, 3.8006102046741123`*^9,
3.8007685109300833`*^9, 3.8007685576592083`*^9, 3.800798395327581*^9,
3.8008367217300673`*^9, 3.801289306918448*^9, 3.801723120900467*^9,
3.8018052341161327`*^9, 3.801901223143807*^9},
CellLabel->
"Out[131]=",ExpressionUUID->"41029220-aae2-4aaa-a49d-ef2c6d004855"],
Cell[BoxData["\<\"/Users/xisco/SXS/BBH_SKS_d14.3_q1.22_sA_0_0_0.330_sB_0_0_-0.\
440/Lev6/metadata.txt\"\>"], "Output",
CellChangeTimes->{{3.717216498380938*^9, 3.7172165029748363`*^9},
3.7172194183959084`*^9, 3.717219451757556*^9, 3.7172210443770657`*^9,
3.717221284233238*^9, 3.717223885556622*^9, 3.71723083183225*^9,
3.7172308731861277`*^9, {3.717233930309966*^9, 3.717233959115399*^9},
3.717234265666395*^9, 3.717389810564118*^9, 3.717397796938533*^9,
3.717409087835437*^9, 3.719042267675703*^9, 3.719725716500965*^9,
3.720497982239339*^9, 3.720601796148543*^9, 3.7620579313020773`*^9,
3.7621713078981047`*^9, 3.762185875336357*^9, 3.762230696430828*^9,
3.762579690340506*^9, 3.762579774070689*^9, 3.76260245990764*^9,
3.762603825455809*^9, 3.762681789861004*^9, 3.762681826567923*^9,
3.764310836001659*^9, 3.76439221450559*^9, 3.764406818443306*^9,
3.7644071506895514`*^9, 3.7644169407193203`*^9, 3.764471034222904*^9,
3.7645023104911823`*^9, 3.7645131854137583`*^9, 3.764567390841213*^9,
3.7685680107885113`*^9, 3.768572298355291*^9, 3.768645827647299*^9,
3.768649839712215*^9, 3.76870184589828*^9, 3.76870653364054*^9,
3.7688823539170313`*^9, 3.768909643980805*^9, 3.769135241869529*^9,
3.7691371202865047`*^9, 3.769138853846279*^9, 3.7691738611067*^9,
3.7692283084744673`*^9, 3.769488530350202*^9, 3.769488774949139*^9,
3.769754030864167*^9, 3.769849332079545*^9, 3.769927143513033*^9, 3.769929847048765*^9, 3.7699349547247353`*^9, {3.7700014592708473`*^9,
3.770001484137087*^9}, {3.770089589019122*^9, 3.770089614454945*^9},
3.770091150926076*^9, {3.7700911818740463`*^9, 3.770091214764297*^9},
3.770092858828313*^9, {3.770096724350412*^9, 3.770096741200189*^9},
3.770175903073946*^9, 3.770355811577338*^9, 3.770368573479495*^9,
3.77043101246592*^9, 3.770452812583289*^9, 3.771218939593539*^9,
3.771673231116873*^9, 3.771907031257121*^9, 3.771907438848817*^9,
3.771915621797638*^9, 3.77219083949329*^9, 3.77219651074804*^9,
3.7722478558236303`*^9, 3.772256628919487*^9, 3.7722590823521147`*^9,
3.772348455019628*^9, 3.772367801149004*^9, 3.7723708907016478`*^9, {
3.7724200004504337`*^9, 3.7724200042661133`*^9}, {3.7724281623433*^9,
3.7724281708790503`*^9}, 3.772428409531721*^9, {3.7724306652135067`*^9,
3.772430670282818*^9}, {3.772440455361487*^9, 3.772440460876882*^9},
3.772529271686309*^9, 3.772875914079034*^9, 3.772878050591577*^9,
3.772885194256702*^9, 3.7729428221831827`*^9, 3.772943402612479*^9,
3.772943541918996*^9, 3.77294405673494*^9, {3.772949569772958*^9,
3.7729496256799297`*^9}, 3.7729504628386717`*^9, 3.77295051555145*^9,
3.7730492856798477`*^9, {3.773049628186563*^9, 3.773049630622552*^9},
3.773144075790704*^9, 3.773460623262436*^9, 3.7734639987477903`*^9,
3.773471312015882*^9, 3.7734788840985126`*^9, 3.773479639599619*^9,
3.773481967713728*^9, 3.773645277444871*^9, 3.773646804195562*^9,
3.77398624805871*^9, 3.776510090600584*^9, 3.7765790822290382`*^9,
3.776662520362281*^9, 3.776668741495263*^9, 3.776669989208765*^9,
3.7766754714003077`*^9, 3.776676366016034*^9, 3.776677554881569*^9,
3.776678156846538*^9, 3.776679728131668*^9, 3.776679799421623*^9,
3.776679922424872*^9, 3.7766814664267282`*^9, 3.776762885737376*^9,
3.777024814804093*^9, 3.777117695893361*^9, 3.777184967307158*^9,
3.77762564240553*^9, 3.7776907196567163`*^9, 3.777725375935665*^9,
3.777786888657611*^9, 3.777953025450624*^9, 3.7779643276193657`*^9,
3.778220763522934*^9, 3.778241132862837*^9, 3.778242215351343*^9,
3.778389851692855*^9, 3.778489688754407*^9, 3.778834393505157*^9,
3.778913036673313*^9, 3.778918648505211*^9, 3.778999681983601*^9,
3.779014946982914*^9, 3.779016358626405*^9, 3.779091498506896*^9,
3.779109063467824*^9, 3.779438042912994*^9, 3.779439035500759*^9,
3.7794461298809433`*^9, 3.77951277696731*^9, 3.779544706262*^9,
3.7795459144103613`*^9, 3.7796202555125847`*^9, 3.779967327387167*^9,
3.780050387061942*^9, 3.780061151710363*^9, 3.780115034431876*^9,
3.7801154413407516`*^9, {3.78099764910747*^9, 3.7809976642917347`*^9},
3.780997931988578*^9, 3.7824639802216053`*^9, 3.783947167922887*^9,
3.785049267615994*^9, 3.785052848818983*^9, 3.785057076455508*^9,
3.785159356717268*^9, 3.7885250658656263`*^9, 3.788588588281761*^9,
3.789016374469722*^9, 3.789101163461742*^9, 3.797132624634801*^9,
3.79713502596594*^9, 3.799113592904676*^9, {3.7991136528442707`*^9,
3.799113708088469*^9}, 3.7991139845245934`*^9, 3.799119938366549*^9,
3.799121028590159*^9, 3.7991216828410063`*^9, 3.7991219650292177`*^9,
3.7991220856243143`*^9, {3.79912227538988*^9, 3.79912230259875*^9},
3.8000746751750097`*^9, 3.800088876448193*^9, 3.800148382940559*^9, {
3.8005393286214533`*^9, 3.800539348433199*^9}, 3.8006102046741123`*^9,
3.8007685109300833`*^9, 3.8007685576592083`*^9, 3.800798395327581*^9,
3.8008367217300673`*^9, 3.801289306918448*^9, 3.801723120900467*^9,
3.8018052341161327`*^9, 3.801901223145651*^9},
CellLabel->
"Out[132]=",ExpressionUUID->"93c7ae4a-68f8-4aff-99b3-3854a0d267da"]
}, Open ]],
Cell[CellGroupData[{
Cell[BoxData[{
RowBox[{"Print", "[",
RowBox[{"\"\<mass1 = \>\"", ",",
RowBox[{"mass1", "=",
RowBox[{
RowBox[{"(",
RowBox[{"\"\<initial-mass1\>\"", "/.", "metadata"}], ")"}], "[",
RowBox[{"[", "1", "]"}], "]"}]}]}], "]"}], "\[IndentingNewLine]",
RowBox[{"Print", "[",
RowBox[{"\"\<mass2 = \>\"", ",",
RowBox[{"mass2", "=",
RowBox[{
RowBox[{"(", RowBox[{"\"\<initial-mass2\>\"", "/.", "metadata"}], ")"}], "[",
RowBox[{"[", "1", "]"}], "]"}]}]}], "]"}], "\[IndentingNewLine]",
RowBox[{"Print", "[",
RowBox[{"\"\<\[Chi]1 = \>\"", ",",
RowBox[{"\[Chi]1", "=",
RowBox[{"Chop", "[",
RowBox[{"(",
RowBox[{
RowBox[{"(",
RowBox[{"\"\<initial-dimensionless-spin1\>\"", "/.", "metadata"}],
")"}], "[",
RowBox[{"[",
RowBox[{"-", "1"}], "]"}], "]"}], ")"}], "]"}]}]}],
"]"}], "\[IndentingNewLine]",
RowBox[{"Print", "[",
RowBox[{"\"\<\[Chi]2 = \>\"", ",",
RowBox[{"\[Chi]2", "=",
RowBox[{"Chop", "[",
RowBox[{"(",
RowBox[{
RowBox[{"(",
RowBox[{"\"\<initial-dimensionless-spin2\>\"", "/.", "metadata"}],
")"}], "[",
RowBox[{"[",
RowBox[{"-", "1"}], "]"}], "]"}], ")"}], "]"}]}]}],
"]"}], "\[IndentingNewLine]",
RowBox[{"Print", "[",
RowBox[{"\"\<m1/m2 = \>\"", ",",
RowBox[{"q", "=",
RowBox[{"Max", "[",
RowBox[{"{",
RowBox[{
RowBox[{"mass1", "/", "mass2"}], ",", "1"}], "}"}], "]"}]}]}],
"]"}], "\[IndentingNewLine]",
RowBox[{"Print", "[",
RowBox[{"\"\<m1\[CenterDot]m2 = \>\"", ",",
RowBox[{"\[Eta]", "=",
RowBox[{
RowBox[{"q", "/",
RowBox[{
RowBox[{"(",
RowBox[{"1", "+", "q"}], ")"}], "^", "2"}]}], "*", "1."}]}]}],
"]"}], "\n",
RowBox[{"Print", "[",
RowBox[{"\"\<af = \>\"", ",",
RowBox[{"af", "=",
RowBox[{
RowBox[{"(",
RowBox[{"\"\<remnant-dimensionless-spin\>\"", "/.", "metadata"}], ")"}],
"[",
RowBox[{"[",
RowBox[{"-", "1"}], "]"}], "]"}]}]}], "]"}], "\n",
RowBox[{"Print", "[",
RowBox[{"\"\<mf = \>\"", ",",
RowBox[{"mf", "=",
RowBox[{
RowBox[{"(",
RowBox[{"\"\<remnant-mass\>\"", "/.", "metadata"}], ")"}], "[",
RowBox[{"[", "1", "]"}], "]"}]}]}], "]"}], "\n",
RowBox[{"Print", "[",
RowBox[{"\"\<Tag = \>\"", ",",
RowBox[{"tag", "=",
RowBox[{
RowBox[{"(",
RowBox[{"\"\<alternative-names\>\"", "/.", "metadata"}], ")"}], "[",
RowBox[{"[", "2", "]"}], "]"}]}]}], "]"}]}], "Code",
CellChangeTimes->{{3.7171458048227043`*^9, 3.7171458352702723`*^9}, {
3.71721654881397*^9, 3.717216557507278*^9}, {3.717223896217493*^9,
3.717223976111596*^9}, {3.762057934955058*^9, 3.762057944429266*^9}, {
3.762681836996138*^9, 3.762681839233232*^9}, {3.7626820639192266`*^9, 3.7626820648200083`*^9}, {3.7700014602000523`*^9, 3.770001460770048*^9}, {
3.770355810960787*^9, 3.77035581334439*^9}, {3.779514788245801*^9,
3.779514823386899*^9}, 3.779514858783835*^9, {3.78002425250508*^9,
3.7800242983466043`*^9}, {3.7800243862787457`*^9, 3.780024392075529*^9}, {
3.7850570816980143`*^9, 3.7850570936459417`*^9}, {3.7850571258554573`*^9,
3.785057158276864*^9}},
CellLabel->
"In[134]:=",ExpressionUUID->"6e14f11c-a2d0-4ce1-82ab-1a2a53c717dd"],
Cell[CellGroupData[{
Cell[BoxData[
InterpretationBox[
RowBox[{"\<\"mass1 = \"\>", "\[InvisibleSpace]", "0.5497953688701543`"}],
SequenceForm["mass1 = ", 0.5497953688701543],
Editable->False]], "Print",
CellChangeTimes->{{3.717223923734548*^9, 3.717223928214951*^9}, {
3.717223971239379*^9, 3.7172239766293097`*^9}, {3.717233931322109*^9,
3.7172339602600603`*^9}, 3.717234266025383*^9, 3.717389811826817*^9,
3.717397798635911*^9, 3.717409088225647*^9, 3.719042269449914*^9,
3.71972571685752*^9, 3.720497982809791*^9, 3.7206017968494253`*^9,
3.762057944992012*^9, 3.762171308612159*^9, 3.762185878959035*^9,
3.762230704000032*^9, 3.76257969223776*^9, 3.762579775163062*^9,
3.76260246063175*^9, 3.762603825586153*^9, {3.7626818315079517`*^9,
3.76268183964472*^9}, {3.762682055451223*^9, 3.762682065259396*^9},
3.764310838838461*^9, 3.76439221561644*^9, 3.764406820576201*^9,
3.7644071524893913`*^9, 3.764416941641247*^9, 3.764471035389357*^9,
3.764502311266172*^9, 3.76451318673353*^9, 3.764567390943874*^9,
3.768568011438259*^9, 3.768572299139023*^9, 3.768645828428453*^9,
3.7686498406439943`*^9, 3.7687018467692337`*^9, 3.7687065345918627`*^9,
3.768882355269024*^9, 3.76890964428026*^9, 3.769135242785296*^9,
3.769137120470574*^9, 3.769138854628717*^9, 3.769173861479599*^9,
3.769228309277001*^9, 3.769488536636314*^9, 3.769488775493741*^9,
3.7697540311309013`*^9, 3.7698493329786863`*^9, 3.769927143788989*^9,
3.7699298474833517`*^9, 3.7699349643369827`*^9, {3.7700014615163918`*^9,
3.770001486342461*^9}, 3.770089615045055*^9, 3.770091151999776*^9,
3.7700911832234373`*^9, 3.77009121560664*^9, 3.770092860611184*^9, {
3.770096733677506*^9, 3.770096743552308*^9}, 3.770175904771923*^9,
3.770355813747511*^9, 3.770368575883328*^9, 3.7704310126844397`*^9,
3.7704528133156967`*^9, 3.771218941512669*^9, 3.771673232190048*^9,
3.771907031907504*^9, 3.771907439029517*^9, 3.771915621974977*^9,
3.7721908396758137`*^9, 3.77219651143331*^9, 3.772247856888829*^9,
3.7722566296994762`*^9, 3.772259084152637*^9, 3.772348455386256*^9,
3.772367801749072*^9, 3.772370891811442*^9, {3.77242000146596*^9,
3.7724200050229883`*^9}, 3.7724281731485243`*^9, 3.772428410081149*^9,
3.772430677028174*^9, 3.772440462784586*^9, 3.772529273853674*^9,
3.7728759142457733`*^9, 3.7728780519980917`*^9, 3.772885195029847*^9,
3.772942823379051*^9, 3.7729434028029613`*^9, 3.772943543318931*^9,
3.7729440577300873`*^9, 3.772949628763401*^9, 3.772950463819756*^9,
3.772950515768664*^9, 3.773049286443142*^9, 3.7730496313722057`*^9,
3.773144076011161*^9, 3.773460624127591*^9, 3.7734639993643*^9,
3.7734713127079906`*^9, 3.773478884641347*^9, 3.773479639752248*^9,
3.7734819682595043`*^9, 3.773645277661384*^9, 3.773646804937044*^9,
3.773986250904626*^9, 3.7765100953056383`*^9, 3.776579083431468*^9,
3.776662521396657*^9, 3.7766687415984173`*^9, 3.776669989743709*^9,
3.7766754714989023`*^9, 3.7766763661356173`*^9, 3.776677554983851*^9,
3.776678156945221*^9, 3.776679730384227*^9, 3.7766797995218353`*^9,
3.776679922527975*^9, 3.776681466538108*^9, 3.776762885901209*^9,
3.7770248152256727`*^9, 3.77711769599905*^9, 3.7771849674192142`*^9,
3.77762564250346*^9, 3.777690720738963*^9, 3.7777253762606087`*^9,
3.777786888754923*^9, 3.777953025556863*^9, 3.777964327722644*^9,
3.778220763638605*^9, 3.778241133006884*^9, 3.77824221545049*^9,
3.778389851806761*^9, 3.778489688892495*^9, 3.778834393767021*^9,
3.778913036796289*^9, 3.7789186486061497`*^9, 3.778999683599246*^9,
3.779014947266182*^9, 3.779016358891302*^9, 3.779091498784801*^9,
3.779109063617063*^9, 3.77943804301024*^9, 3.7794390356224527`*^9,
3.779446129998694*^9, 3.779512777100968*^9, 3.779514823682631*^9,
3.779514859562222*^9, 3.779544709953924*^9, 3.7795459145467453`*^9,
3.779620255613166*^9, 3.779967327511076*^9, {3.780024275082299*^9,
3.7800242987111483`*^9}, {3.7800243789563026`*^9, 3.780024392712452*^9}, 3.78005038726173*^9, 3.7800611518316393`*^9, 3.780115034568583*^9,
3.7801154414361067`*^9, {3.780997650057397*^9, 3.78099766439631*^9},
3.7824639813777742`*^9, 3.7839471685460453`*^9, 3.78504926862815*^9,
3.785052849502256*^9, 3.785057077348769*^9, {3.785057141637362*^9,
3.7850571599033833`*^9}, 3.785159356911799*^9, 3.7885250659535437`*^9,
3.7885885882949257`*^9, 3.789016374575654*^9, 3.7891011635298862`*^9,
3.797132624730732*^9, 3.7971350260463753`*^9, 3.799113710120328*^9,
3.7991139854914227`*^9, 3.799119939271639*^9, 3.799121028607518*^9,
3.79912168285885*^9, 3.799121965114169*^9, 3.7991220861947727`*^9, {
3.7991222754707747`*^9, 3.799122302610426*^9}, 3.800074675182968*^9,
3.8000888764634867`*^9, 3.8001483829496517`*^9, {3.800539328667665*^9,
3.80053935006597*^9}, 3.8006102047411127`*^9, 3.800768511041994*^9,
3.8007685584485817`*^9, 3.800798395343419*^9, 3.80083672174466*^9,
3.801289306970142*^9, 3.801723120960483*^9, 3.801805234128213*^9,
3.80190123947233*^9},
CellLabel->
"During evaluation of \
In[134]:=",ExpressionUUID->"b0646428-ecbf-4825-87c0-5f5a9a25a494"],
Cell[BoxData[
InterpretationBox[
RowBox[{"\<\"mass2 = \"\>", "\[InvisibleSpace]", "0.450204644189637`"}],
SequenceForm["mass2 = ", 0.450204644189637],
Editable->False]], "Print",
CellChangeTimes->{{3.717223923734548*^9, 3.717223928214951*^9}, {
3.717223971239379*^9, 3.7172239766293097`*^9}, {3.717233931322109*^9,
3.7172339602600603`*^9}, 3.717234266025383*^9, 3.717389811826817*^9,
3.717397798635911*^9, 3.717409088225647*^9, 3.719042269449914*^9,
3.71972571685752*^9, 3.720497982809791*^9, 3.7206017968494253`*^9,
3.762057944992012*^9, 3.762171308612159*^9, 3.762185878959035*^9,
3.762230704000032*^9, 3.76257969223776*^9, 3.762579775163062*^9,
3.76260246063175*^9, 3.762603825586153*^9, {3.7626818315079517`*^9,
3.76268183964472*^9}, {3.762682055451223*^9, 3.762682065259396*^9},
3.764310838838461*^9, 3.76439221561644*^9, 3.764406820576201*^9,
3.7644071524893913`*^9, 3.764416941641247*^9, 3.764471035389357*^9,
3.764502311266172*^9, 3.76451318673353*^9, 3.764567390943874*^9,
3.768568011438259*^9, 3.768572299139023*^9, 3.768645828428453*^9,
3.7686498406439943`*^9, 3.7687018467692337`*^9, 3.7687065345918627`*^9,
3.768882355269024*^9, 3.76890964428026*^9, 3.769135242785296*^9,
3.769137120470574*^9, 3.769138854628717*^9, 3.769173861479599*^9,
3.769228309277001*^9, 3.769488536636314*^9, 3.769488775493741*^9,
3.7697540311309013`*^9, 3.7698493329786863`*^9, 3.769927143788989*^9,
3.7699298474833517`*^9, 3.7699349643369827`*^9, {3.7700014615163918`*^9,
3.770001486342461*^9}, 3.770089615045055*^9, 3.770091151999776*^9,
3.7700911832234373`*^9, 3.77009121560664*^9, 3.770092860611184*^9, {
3.770096733677506*^9, 3.770096743552308*^9}, 3.770175904771923*^9,
3.770355813747511*^9, 3.770368575883328*^9, 3.7704310126844397`*^9,
3.7704528133156967`*^9, 3.771218941512669*^9, 3.771673232190048*^9,
3.771907031907504*^9, 3.771907439029517*^9, 3.771915621974977*^9,
3.7721908396758137`*^9, 3.77219651143331*^9, 3.772247856888829*^9,
3.7722566296994762`*^9, 3.772259084152637*^9, 3.772348455386256*^9,
3.772367801749072*^9, 3.772370891811442*^9, {3.77242000146596*^9,
3.7724200050229883`*^9}, 3.7724281731485243`*^9, 3.772428410081149*^9,
3.772430677028174*^9, 3.772440462784586*^9, 3.772529273853674*^9,
3.7728759142457733`*^9, 3.7728780519980917`*^9, 3.772885195029847*^9,
3.772942823379051*^9, 3.7729434028029613`*^9, 3.772943543318931*^9,
3.7729440577300873`*^9, 3.772949628763401*^9, 3.772950463819756*^9,
3.772950515768664*^9, 3.773049286443142*^9, 3.7730496313722057`*^9,
3.773144076011161*^9, 3.773460624127591*^9, 3.7734639993643*^9,
3.7734713127079906`*^9, 3.773478884641347*^9, 3.773479639752248*^9,
3.7734819682595043`*^9, 3.773645277661384*^9, 3.773646804937044*^9,
3.773986250904626*^9, 3.7765100953056383`*^9, 3.776579083431468*^9,
3.776662521396657*^9, 3.7766687415984173`*^9, 3.776669989743709*^9,
3.7766754714989023`*^9, 3.7766763661356173`*^9, 3.776677554983851*^9,
3.776678156945221*^9, 3.776679730384227*^9, 3.7766797995218353`*^9,
3.776679922527975*^9, 3.776681466538108*^9, 3.776762885901209*^9,
3.7770248152256727`*^9, 3.77711769599905*^9, 3.7771849674192142`*^9,
3.77762564250346*^9, 3.777690720738963*^9, 3.7777253762606087`*^9,
3.777786888754923*^9, 3.777953025556863*^9, 3.777964327722644*^9,
3.778220763638605*^9, 3.778241133006884*^9, 3.77824221545049*^9, 3.778389851806761*^9, 3.778489688892495*^9, 3.778834393767021*^9,
3.778913036796289*^9, 3.7789186486061497`*^9, 3.778999683599246*^9,
3.779014947266182*^9, 3.779016358891302*^9, 3.779091498784801*^9,
3.779109063617063*^9, 3.77943804301024*^9, 3.7794390356224527`*^9,
3.779446129998694*^9, 3.779512777100968*^9, 3.779514823682631*^9,
3.779514859562222*^9, 3.779544709953924*^9, 3.7795459145467453`*^9,
3.779620255613166*^9, 3.779967327511076*^9, {3.780024275082299*^9,
3.7800242987111483`*^9}, {3.7800243789563026`*^9, 3.780024392712452*^9},
3.78005038726173*^9, 3.7800611518316393`*^9, 3.780115034568583*^9,
3.7801154414361067`*^9, {3.780997650057397*^9, 3.78099766439631*^9},
3.7824639813777742`*^9, 3.7839471685460453`*^9, 3.78504926862815*^9,
3.785052849502256*^9, 3.785057077348769*^9, {3.785057141637362*^9,
3.7850571599033833`*^9}, 3.785159356911799*^9, 3.7885250659535437`*^9,
3.7885885882949257`*^9, 3.789016374575654*^9, 3.7891011635298862`*^9,
3.797132624730732*^9, 3.7971350260463753`*^9, 3.799113710120328*^9,
3.7991139854914227`*^9, 3.799119939271639*^9, 3.799121028607518*^9,
3.79912168285885*^9, 3.799121965114169*^9, 3.7991220861947727`*^9, {
3.7991222754707747`*^9, 3.799122302610426*^9}, 3.800074675182968*^9,
3.8000888764634867`*^9, 3.8001483829496517`*^9, {3.800539328667665*^9,
3.80053935006597*^9}, 3.8006102047411127`*^9, 3.800768511041994*^9,
3.8007685584485817`*^9, 3.800798395343419*^9, 3.80083672174466*^9,
3.801289306970142*^9, 3.801723120960483*^9, 3.801805234128213*^9,
3.80190123947443*^9},
CellLabel->
"During evaluation of \
In[134]:=",ExpressionUUID->"a5cf05bc-2cbc-4f46-a364-fa765f897969"],
Cell[BoxData[
InterpretationBox[
RowBox[{"\<\"\[Chi]1 = \"\>", "\[InvisibleSpace]", "0.330000004123`"}],
SequenceForm["\[Chi]1 = ", 0.330000004123],
Editable->False]], "Print",
CellChangeTimes->{{3.717223923734548*^9, 3.717223928214951*^9}, {
3.717223971239379*^9, 3.7172239766293097`*^9}, {3.717233931322109*^9,
3.7172339602600603`*^9}, 3.717234266025383*^9, 3.717389811826817*^9,
3.717397798635911*^9, 3.717409088225647*^9, 3.719042269449914*^9,
3.71972571685752*^9, 3.720497982809791*^9, 3.7206017968494253`*^9,
3.762057944992012*^9, 3.762171308612159*^9, 3.762185878959035*^9,
3.762230704000032*^9, 3.76257969223776*^9, 3.762579775163062*^9,
3.76260246063175*^9, 3.762603825586153*^9, {3.7626818315079517`*^9,
3.76268183964472*^9}, {3.762682055451223*^9, 3.762682065259396*^9},
3.764310838838461*^9, 3.76439221561644*^9, 3.764406820576201*^9,
3.7644071524893913`*^9, 3.764416941641247*^9, 3.764471035389357*^9,
3.764502311266172*^9, 3.76451318673353*^9, 3.764567390943874*^9,
3.768568011438259*^9, 3.768572299139023*^9, 3.768645828428453*^9,
3.7686498406439943`*^9, 3.7687018467692337`*^9, 3.7687065345918627`*^9,
3.768882355269024*^9, 3.76890964428026*^9, 3.769135242785296*^9,
3.769137120470574*^9, 3.769138854628717*^9, 3.769173861479599*^9,
3.769228309277001*^9, 3.769488536636314*^9, 3.769488775493741*^9,
3.7697540311309013`*^9, 3.7698493329786863`*^9, 3.769927143788989*^9,
3.7699298474833517`*^9, 3.7699349643369827`*^9, {3.7700014615163918`*^9,
3.770001486342461*^9}, 3.770089615045055*^9, 3.770091151999776*^9,
3.7700911832234373`*^9, 3.77009121560664*^9, 3.770092860611184*^9, {
3.770096733677506*^9, 3.770096743552308*^9}, 3.770175904771923*^9,
3.770355813747511*^9, 3.770368575883328*^9, 3.7704310126844397`*^9,
3.7704528133156967`*^9, 3.771218941512669*^9, 3.771673232190048*^9,
3.771907031907504*^9, 3.771907439029517*^9, 3.771915621974977*^9,
3.7721908396758137`*^9, 3.77219651143331*^9, 3.772247856888829*^9,
3.7722566296994762`*^9, 3.772259084152637*^9, 3.772348455386256*^9,
3.772367801749072*^9, 3.772370891811442*^9, {3.77242000146596*^9,
3.7724200050229883`*^9}, 3.7724281731485243`*^9, 3.772428410081149*^9,
3.772430677028174*^9, 3.772440462784586*^9, 3.772529273853674*^9,
3.7728759142457733`*^9, 3.7728780519980917`*^9, 3.772885195029847*^9,
3.772942823379051*^9, 3.7729434028029613`*^9, 3.772943543318931*^9,
3.7729440577300873`*^9, 3.772949628763401*^9, 3.772950463819756*^9,
3.772950515768664*^9, 3.773049286443142*^9, 3.7730496313722057`*^9,
3.773144076011161*^9, 3.773460624127591*^9, 3.7734639993643*^9,
3.7734713127079906`*^9, 3.773478884641347*^9, 3.773479639752248*^9,
3.7734819682595043`*^9, 3.773645277661384*^9, 3.773646804937044*^9,
3.773986250904626*^9, 3.7765100953056383`*^9, 3.776579083431468*^9, 3.776662521396657*^9, 3.7766687415984173`*^9, 3.776669989743709*^9,
3.7766754714989023`*^9, 3.7766763661356173`*^9, 3.776677554983851*^9,
3.776678156945221*^9, 3.776679730384227*^9, 3.7766797995218353`*^9,
3.776679922527975*^9, 3.776681466538108*^9, 3.776762885901209*^9,
3.7770248152256727`*^9, 3.77711769599905*^9, 3.7771849674192142`*^9,
3.77762564250346*^9, 3.777690720738963*^9, 3.7777253762606087`*^9,
3.777786888754923*^9, 3.777953025556863*^9, 3.777964327722644*^9,
3.778220763638605*^9, 3.778241133006884*^9, 3.77824221545049*^9,
3.778389851806761*^9, 3.778489688892495*^9, 3.778834393767021*^9,
3.778913036796289*^9, 3.7789186486061497`*^9, 3.778999683599246*^9,
3.779014947266182*^9, 3.779016358891302*^9, 3.779091498784801*^9,
3.779109063617063*^9, 3.77943804301024*^9, 3.7794390356224527`*^9,
3.779446129998694*^9, 3.779512777100968*^9, 3.779514823682631*^9,
3.779514859562222*^9, 3.779544709953924*^9, 3.7795459145467453`*^9,
3.779620255613166*^9, 3.779967327511076*^9, {3.780024275082299*^9,
3.7800242987111483`*^9}, {3.7800243789563026`*^9, 3.780024392712452*^9},
3.78005038726173*^9, 3.7800611518316393`*^9, 3.780115034568583*^9,
3.7801154414361067`*^9, {3.780997650057397*^9, 3.78099766439631*^9},
3.7824639813777742`*^9, 3.7839471685460453`*^9, 3.78504926862815*^9,
3.785052849502256*^9, 3.785057077348769*^9, {3.785057141637362*^9,
3.7850571599033833`*^9}, 3.785159356911799*^9, 3.7885250659535437`*^9,
3.7885885882949257`*^9, 3.789016374575654*^9, 3.7891011635298862`*^9,
3.797132624730732*^9, 3.7971350260463753`*^9, 3.799113710120328*^9,
3.7991139854914227`*^9, 3.799119939271639*^9, 3.799121028607518*^9,
3.79912168285885*^9, 3.799121965114169*^9, 3.7991220861947727`*^9, {
3.7991222754707747`*^9, 3.799122302610426*^9}, 3.800074675182968*^9,
3.8000888764634867`*^9, 3.8001483829496517`*^9, {3.800539328667665*^9,
3.80053935006597*^9}, 3.8006102047411127`*^9, 3.800768511041994*^9,
3.8007685584485817`*^9, 3.800798395343419*^9, 3.80083672174466*^9,
3.801289306970142*^9, 3.801723120960483*^9, 3.801805234128213*^9,
3.8019012394828453`*^9},
CellLabel->
"During evaluation of \
In[134]:=",ExpressionUUID->"21ae58be-bb6b-4910-95e4-c68634b6ddf2"],
Cell[BoxData[
InterpretationBox[
RowBox[{"\<\"\[Chi]2 = \"\>", "\[InvisibleSpace]",
RowBox[{"-", "0.440000004931`"}]}],
SequenceForm["\[Chi]2 = ", -0.440000004931],
Editable->False]], "Print",
CellChangeTimes->{{3.717223923734548*^9, 3.717223928214951*^9}, {
3.717223971239379*^9, 3.7172239766293097`*^9}, {3.717233931322109*^9,
3.7172339602600603`*^9}, 3.717234266025383*^9, 3.717389811826817*^9,
3.717397798635911*^9, 3.717409088225647*^9, 3.719042269449914*^9,
3.71972571685752*^9, 3.720497982809791*^9, 3.7206017968494253`*^9,
3.762057944992012*^9, 3.762171308612159*^9, 3.762185878959035*^9,
3.762230704000032*^9, 3.76257969223776*^9, 3.762579775163062*^9,
3.76260246063175*^9, 3.762603825586153*^9, {3.7626818315079517`*^9,
3.76268183964472*^9}, {3.762682055451223*^9, 3.762682065259396*^9},
3.764310838838461*^9, 3.76439221561644*^9, 3.764406820576201*^9,
3.7644071524893913`*^9, 3.764416941641247*^9, 3.764471035389357*^9,
3.764502311266172*^9, 3.76451318673353*^9, 3.764567390943874*^9,
3.768568011438259*^9, 3.768572299139023*^9, 3.768645828428453*^9,
3.7686498406439943`*^9, 3.7687018467692337`*^9, 3.7687065345918627`*^9,
3.768882355269024*^9, 3.76890964428026*^9, 3.769135242785296*^9,
3.769137120470574*^9, 3.769138854628717*^9, 3.769173861479599*^9,
3.769228309277001*^9, 3.769488536636314*^9, 3.769488775493741*^9,
3.7697540311309013`*^9, 3.7698493329786863`*^9, 3.769927143788989*^9,
3.7699298474833517`*^9, 3.7699349643369827`*^9, {3.7700014615163918`*^9,
3.770001486342461*^9}, 3.770089615045055*^9, 3.770091151999776*^9,
3.7700911832234373`*^9, 3.77009121560664*^9, 3.770092860611184*^9, {
3.770096733677506*^9, 3.770096743552308*^9}, 3.770175904771923*^9,
3.770355813747511*^9, 3.770368575883328*^9, 3.7704310126844397`*^9,
3.7704528133156967`*^9, 3.771218941512669*^9, 3.771673232190048*^9,
3.771907031907504*^9, 3.771907439029517*^9, 3.771915621974977*^9,
3.7721908396758137`*^9, 3.77219651143331*^9, 3.772247856888829*^9,
3.7722566296994762`*^9, 3.772259084152637*^9, 3.772348455386256*^9,
3.772367801749072*^9, 3.772370891811442*^9, {3.77242000146596*^9,
3.7724200050229883`*^9}, 3.7724281731485243`*^9, 3.772428410081149*^9, 3.772430677028174*^9, 3.772440462784586*^9, 3.772529273853674*^9,
3.7728759142457733`*^9, 3.7728780519980917`*^9, 3.772885195029847*^9,
3.772942823379051*^9, 3.7729434028029613`*^9, 3.772943543318931*^9,
3.7729440577300873`*^9, 3.772949628763401*^9, 3.772950463819756*^9,
3.772950515768664*^9, 3.773049286443142*^9, 3.7730496313722057`*^9,
3.773144076011161*^9, 3.773460624127591*^9, 3.7734639993643*^9,
3.7734713127079906`*^9, 3.773478884641347*^9, 3.773479639752248*^9,
3.7734819682595043`*^9, 3.773645277661384*^9, 3.773646804937044*^9,
3.773986250904626*^9, 3.7765100953056383`*^9, 3.776579083431468*^9,
3.776662521396657*^9, 3.7766687415984173`*^9, 3.776669989743709*^9,
3.7766754714989023`*^9, 3.7766763661356173`*^9, 3.776677554983851*^9,
3.776678156945221*^9, 3.776679730384227*^9, 3.7766797995218353`*^9,
3.776679922527975*^9, 3.776681466538108*^9, 3.776762885901209*^9,
3.7770248152256727`*^9, 3.77711769599905*^9, 3.7771849674192142`*^9,
3.77762564250346*^9, 3.777690720738963*^9, 3.7777253762606087`*^9,
3.777786888754923*^9, 3.777953025556863*^9, 3.777964327722644*^9,
3.778220763638605*^9, 3.778241133006884*^9, 3.77824221545049*^9,
3.778389851806761*^9, 3.778489688892495*^9, 3.778834393767021*^9,
3.778913036796289*^9, 3.7789186486061497`*^9, 3.778999683599246*^9,
3.779014947266182*^9, 3.779016358891302*^9, 3.779091498784801*^9,
3.779109063617063*^9, 3.77943804301024*^9, 3.7794390356224527`*^9,
3.779446129998694*^9, 3.779512777100968*^9, 3.779514823682631*^9,
3.779514859562222*^9, 3.779544709953924*^9, 3.7795459145467453`*^9,
3.779620255613166*^9, 3.779967327511076*^9, {3.780024275082299*^9,
3.7800242987111483`*^9}, {3.7800243789563026`*^9, 3.780024392712452*^9},
3.78005038726173*^9, 3.7800611518316393`*^9, 3.780115034568583*^9,
3.7801154414361067`*^9, {3.780997650057397*^9, 3.78099766439631*^9},
3.7824639813777742`*^9, 3.7839471685460453`*^9, 3.78504926862815*^9,
3.785052849502256*^9, 3.785057077348769*^9, {3.785057141637362*^9,
3.7850571599033833`*^9}, 3.785159356911799*^9, 3.7885250659535437`*^9,
3.7885885882949257`*^9, 3.789016374575654*^9, 3.7891011635298862`*^9,
3.797132624730732*^9, 3.7971350260463753`*^9, 3.799113710120328*^9,
3.7991139854914227`*^9, 3.799119939271639*^9, 3.799121028607518*^9,
3.79912168285885*^9, 3.799121965114169*^9, 3.7991220861947727`*^9, {
3.7991222754707747`*^9, 3.799122302610426*^9}, 3.800074675182968*^9,
3.8000888764634867`*^9, 3.8001483829496517`*^9, {3.800539328667665*^9,
3.80053935006597*^9}, 3.8006102047411127`*^9, 3.800768511041994*^9,
3.8007685584485817`*^9, 3.800798395343419*^9, 3.80083672174466*^9,
3.801289306970142*^9, 3.801723120960483*^9, 3.801805234128213*^9,
3.8019012394853697`*^9},
CellLabel->
"During evaluation of \
In[134]:=",ExpressionUUID->"5f689c14-b703-4074-a06a-efd763c742b7"],
Cell[BoxData[
InterpretationBox[
RowBox[{"\<\"m1/m2 = \"\>", "\[InvisibleSpace]", "1.2212121220112675`"}],
SequenceForm["m1/m2 = ", 1.2212121220112675`],
Editable->False]], "Print",
CellChangeTimes->{{3.717223923734548*^9, 3.717223928214951*^9}, {
3.717223971239379*^9, 3.7172239766293097`*^9}, {3.717233931322109*^9,
3.7172339602600603`*^9}, 3.717234266025383*^9, 3.717389811826817*^9,
3.717397798635911*^9, 3.717409088225647*^9, 3.719042269449914*^9,
3.71972571685752*^9, 3.720497982809791*^9, 3.7206017968494253`*^9,
3.762057944992012*^9, 3.762171308612159*^9, 3.762185878959035*^9,
3.762230704000032*^9, 3.76257969223776*^9, 3.762579775163062*^9,
3.76260246063175*^9, 3.762603825586153*^9, {3.7626818315079517`*^9,
3.76268183964472*^9}, {3.762682055451223*^9, 3.762682065259396*^9},
3.764310838838461*^9, 3.76439221561644*^9, 3.764406820576201*^9,
3.7644071524893913`*^9, 3.764416941641247*^9, 3.764471035389357*^9,
3.764502311266172*^9, 3.76451318673353*^9, 3.764567390943874*^9,
3.768568011438259*^9, 3.768572299139023*^9, 3.768645828428453*^9,
3.7686498406439943`*^9, 3.7687018467692337`*^9, 3.7687065345918627`*^9,
3.768882355269024*^9, 3.76890964428026*^9, 3.769135242785296*^9,
3.769137120470574*^9, 3.769138854628717*^9, 3.769173861479599*^9,
3.769228309277001*^9, 3.769488536636314*^9, 3.769488775493741*^9,
3.7697540311309013`*^9, 3.7698493329786863`*^9, 3.769927143788989*^9,
3.7699298474833517`*^9, 3.7699349643369827`*^9, {3.7700014615163918`*^9,
3.770001486342461*^9}, 3.770089615045055*^9, 3.770091151999776*^9,
3.7700911832234373`*^9, 3.77009121560664*^9, 3.770092860611184*^9, { 3.770096733677506*^9, 3.770096743552308*^9}, 3.770175904771923*^9,
3.770355813747511*^9, 3.770368575883328*^9, 3.7704310126844397`*^9,
3.7704528133156967`*^9, 3.771218941512669*^9, 3.771673232190048*^9,
3.771907031907504*^9, 3.771907439029517*^9, 3.771915621974977*^9,
3.7721908396758137`*^9, 3.77219651143331*^9, 3.772247856888829*^9,
3.7722566296994762`*^9, 3.772259084152637*^9, 3.772348455386256*^9,
3.772367801749072*^9, 3.772370891811442*^9, {3.77242000146596*^9,
3.7724200050229883`*^9}, 3.7724281731485243`*^9, 3.772428410081149*^9,
3.772430677028174*^9, 3.772440462784586*^9, 3.772529273853674*^9,
3.7728759142457733`*^9, 3.7728780519980917`*^9, 3.772885195029847*^9,
3.772942823379051*^9, 3.7729434028029613`*^9, 3.772943543318931*^9,
3.7729440577300873`*^9, 3.772949628763401*^9, 3.772950463819756*^9,
3.772950515768664*^9, 3.773049286443142*^9, 3.7730496313722057`*^9,
3.773144076011161*^9, 3.773460624127591*^9, 3.7734639993643*^9,
3.7734713127079906`*^9, 3.773478884641347*^9, 3.773479639752248*^9,
3.7734819682595043`*^9, 3.773645277661384*^9, 3.773646804937044*^9,
3.773986250904626*^9, 3.7765100953056383`*^9, 3.776579083431468*^9,
3.776662521396657*^9, 3.7766687415984173`*^9, 3.776669989743709*^9,
3.7766754714989023`*^9, 3.7766763661356173`*^9, 3.776677554983851*^9,
3.776678156945221*^9, 3.776679730384227*^9, 3.7766797995218353`*^9,
3.776679922527975*^9, 3.776681466538108*^9, 3.776762885901209*^9,
3.7770248152256727`*^9, 3.77711769599905*^9, 3.7771849674192142`*^9,
3.77762564250346*^9, 3.777690720738963*^9, 3.7777253762606087`*^9,
3.777786888754923*^9, 3.777953025556863*^9, 3.777964327722644*^9,
3.778220763638605*^9, 3.778241133006884*^9, 3.77824221545049*^9,
3.778389851806761*^9, 3.778489688892495*^9, 3.778834393767021*^9,
3.778913036796289*^9, 3.7789186486061497`*^9, 3.778999683599246*^9,
3.779014947266182*^9, 3.779016358891302*^9, 3.779091498784801*^9,
3.779109063617063*^9, 3.77943804301024*^9, 3.7794390356224527`*^9,
3.779446129998694*^9, 3.779512777100968*^9, 3.779514823682631*^9,
3.779514859562222*^9, 3.779544709953924*^9, 3.7795459145467453`*^9,
3.779620255613166*^9, 3.779967327511076*^9, {3.780024275082299*^9,
3.7800242987111483`*^9}, {3.7800243789563026`*^9, 3.780024392712452*^9},
3.78005038726173*^9, 3.7800611518316393`*^9, 3.780115034568583*^9,
3.7801154414361067`*^9, {3.780997650057397*^9, 3.78099766439631*^9},
3.7824639813777742`*^9, 3.7839471685460453`*^9, 3.78504926862815*^9,
3.785052849502256*^9, 3.785057077348769*^9, {3.785057141637362*^9,
3.7850571599033833`*^9}, 3.785159356911799*^9, 3.7885250659535437`*^9,
3.7885885882949257`*^9, 3.789016374575654*^9, 3.7891011635298862`*^9,
3.797132624730732*^9, 3.7971350260463753`*^9, 3.799113710120328*^9,
3.7991139854914227`*^9, 3.799119939271639*^9, 3.799121028607518*^9,
3.79912168285885*^9, 3.799121965114169*^9, 3.7991220861947727`*^9, {
3.7991222754707747`*^9, 3.799122302610426*^9}, 3.800074675182968*^9,
3.8000888764634867`*^9, 3.8001483829496517`*^9, {3.800539328667665*^9,
3.80053935006597*^9}, 3.8006102047411127`*^9, 3.800768511041994*^9,
3.8007685584485817`*^9, 3.800798395343419*^9, 3.80083672174466*^9,
3.801289306970142*^9, 3.801723120960483*^9, 3.801805234128213*^9,
3.801901239487548*^9},
CellLabel->
"During evaluation of \
In[134]:=",ExpressionUUID->"3f8c41fc-6da7-4446-ba56-ed613b63b498"],
Cell[BoxData[
InterpretationBox[
RowBox[{"\<\"m1\[CenterDot]m2 = \"\>", "\[InvisibleSpace]",
"0.24752042195416796`"}],
SequenceForm["m1\[CenterDot]m2 = ", 0.24752042195416796`],
Editable->False]], "Print",
CellChangeTimes->{{3.717223923734548*^9, 3.717223928214951*^9}, {
3.717223971239379*^9, 3.7172239766293097`*^9}, {3.717233931322109*^9,
3.7172339602600603`*^9}, 3.717234266025383*^9, 3.717389811826817*^9,
3.717397798635911*^9, 3.717409088225647*^9, 3.719042269449914*^9,
3.71972571685752*^9, 3.720497982809791*^9, 3.7206017968494253`*^9,
3.762057944992012*^9, 3.762171308612159*^9, 3.762185878959035*^9,
3.762230704000032*^9, 3.76257969223776*^9, 3.762579775163062*^9,
3.76260246063175*^9, 3.762603825586153*^9, {3.7626818315079517`*^9,
3.76268183964472*^9}, {3.762682055451223*^9, 3.762682065259396*^9},
3.764310838838461*^9, 3.76439221561644*^9, 3.764406820576201*^9,
3.7644071524893913`*^9, 3.764416941641247*^9, 3.764471035389357*^9,
3.764502311266172*^9, 3.76451318673353*^9, 3.764567390943874*^9, 3.768568011438259*^9, 3.768572299139023*^9, 3.768645828428453*^9,
3.7686498406439943`*^9, 3.7687018467692337`*^9, 3.7687065345918627`*^9,
3.768882355269024*^9, 3.76890964428026*^9, 3.769135242785296*^9,
3.769137120470574*^9, 3.769138854628717*^9, 3.769173861479599*^9,
3.769228309277001*^9, 3.769488536636314*^9, 3.769488775493741*^9,
3.7697540311309013`*^9, 3.7698493329786863`*^9, 3.769927143788989*^9,
3.7699298474833517`*^9, 3.7699349643369827`*^9, {3.7700014615163918`*^9,
3.770001486342461*^9}, 3.770089615045055*^9, 3.770091151999776*^9,
3.7700911832234373`*^9, 3.77009121560664*^9, 3.770092860611184*^9, {
3.770096733677506*^9, 3.770096743552308*^9}, 3.770175904771923*^9,
3.770355813747511*^9, 3.770368575883328*^9, 3.7704310126844397`*^9,
3.7704528133156967`*^9, 3.771218941512669*^9, 3.771673232190048*^9,
3.771907031907504*^9, 3.771907439029517*^9, 3.771915621974977*^9,
3.7721908396758137`*^9, 3.77219651143331*^9, 3.772247856888829*^9,
3.7722566296994762`*^9, 3.772259084152637*^9, 3.772348455386256*^9,
3.772367801749072*^9, 3.772370891811442*^9, {3.77242000146596*^9,
3.7724200050229883`*^9}, 3.7724281731485243`*^9, 3.772428410081149*^9,
3.772430677028174*^9, 3.772440462784586*^9, 3.772529273853674*^9,
3.7728759142457733`*^9, 3.7728780519980917`*^9, 3.772885195029847*^9,
3.772942823379051*^9, 3.7729434028029613`*^9, 3.772943543318931*^9,
3.7729440577300873`*^9, 3.772949628763401*^9, 3.772950463819756*^9,
3.772950515768664*^9, 3.773049286443142*^9, 3.7730496313722057`*^9,
3.773144076011161*^9, 3.773460624127591*^9, 3.7734639993643*^9,
3.7734713127079906`*^9, 3.773478884641347*^9, 3.773479639752248*^9,
3.7734819682595043`*^9, 3.773645277661384*^9, 3.773646804937044*^9,
3.773986250904626*^9, 3.7765100953056383`*^9, 3.776579083431468*^9,
3.776662521396657*^9, 3.7766687415984173`*^9, 3.776669989743709*^9,
3.7766754714989023`*^9, 3.7766763661356173`*^9, 3.776677554983851*^9,
3.776678156945221*^9, 3.776679730384227*^9, 3.7766797995218353`*^9,
3.776679922527975*^9, 3.776681466538108*^9, 3.776762885901209*^9,
3.7770248152256727`*^9, 3.77711769599905*^9, 3.7771849674192142`*^9,
3.77762564250346*^9, 3.777690720738963*^9, 3.7777253762606087`*^9,
3.777786888754923*^9, 3.777953025556863*^9, 3.777964327722644*^9,
3.778220763638605*^9, 3.778241133006884*^9, 3.77824221545049*^9,
3.778389851806761*^9, 3.778489688892495*^9, 3.778834393767021*^9,
3.778913036796289*^9, 3.7789186486061497`*^9, 3.778999683599246*^9,
3.779014947266182*^9, 3.779016358891302*^9, 3.779091498784801*^9,
3.779109063617063*^9, 3.77943804301024*^9, 3.7794390356224527`*^9,
3.779446129998694*^9, 3.779512777100968*^9, 3.779514823682631*^9,
3.779514859562222*^9, 3.779544709953924*^9, 3.7795459145467453`*^9,
3.779620255613166*^9, 3.779967327511076*^9, {3.780024275082299*^9,
3.7800242987111483`*^9}, {3.7800243789563026`*^9, 3.780024392712452*^9},
3.78005038726173*^9, 3.7800611518316393`*^9, 3.780115034568583*^9,
3.7801154414361067`*^9, {3.780997650057397*^9, 3.78099766439631*^9},
3.7824639813777742`*^9, 3.7839471685460453`*^9, 3.78504926862815*^9,
3.785052849502256*^9, 3.785057077348769*^9, {3.785057141637362*^9,
3.7850571599033833`*^9}, 3.785159356911799*^9, 3.7885250659535437`*^9,
3.7885885882949257`*^9, 3.789016374575654*^9, 3.7891011635298862`*^9,
3.797132624730732*^9, 3.7971350260463753`*^9, 3.799113710120328*^9,
3.7991139854914227`*^9, 3.799119939271639*^9, 3.799121028607518*^9,
3.79912168285885*^9, 3.799121965114169*^9, 3.7991220861947727`*^9, {
3.7991222754707747`*^9, 3.799122302610426*^9}, 3.800074675182968*^9,
3.8000888764634867`*^9, 3.8001483829496517`*^9, {3.800539328667665*^9,
3.80053935006597*^9}, 3.8006102047411127`*^9, 3.800768511041994*^9,
3.8007685584485817`*^9, 3.800798395343419*^9, 3.80083672174466*^9,
3.801289306970142*^9, 3.801723120960483*^9, 3.801805234128213*^9,
3.80190123948921*^9},
CellLabel->
"During evaluation of \
In[134]:=",ExpressionUUID->"699de0d3-51d3-456d-b993-57ce112918b1"],
Cell[BoxData[
InterpretationBox[
RowBox[{"\<\"af = \"\>", "\[InvisibleSpace]", "0.692085186818`"}],
SequenceForm["af = ", 0.692085186818],
Editable->False]], "Print",
CellChangeTimes->{{3.717223923734548*^9, 3.717223928214951*^9}, {
3.717223971239379*^9, 3.7172239766293097`*^9}, {3.717233931322109*^9,
3.7172339602600603`*^9}, 3.717234266025383*^9, 3.717389811826817*^9,
3.717397798635911*^9, 3.717409088225647*^9, 3.719042269449914*^9, 3.71972571685752*^9, 3.720497982809791*^9, 3.7206017968494253`*^9,
3.762057944992012*^9, 3.762171308612159*^9, 3.762185878959035*^9,
3.762230704000032*^9, 3.76257969223776*^9, 3.762579775163062*^9,
3.76260246063175*^9, 3.762603825586153*^9, {3.7626818315079517`*^9,
3.76268183964472*^9}, {3.762682055451223*^9, 3.762682065259396*^9},
3.764310838838461*^9, 3.76439221561644*^9, 3.764406820576201*^9,
3.7644071524893913`*^9, 3.764416941641247*^9, 3.764471035389357*^9,
3.764502311266172*^9, 3.76451318673353*^9, 3.764567390943874*^9,
3.768568011438259*^9, 3.768572299139023*^9, 3.768645828428453*^9,
3.7686498406439943`*^9, 3.7687018467692337`*^9, 3.7687065345918627`*^9,
3.768882355269024*^9, 3.76890964428026*^9, 3.769135242785296*^9,
3.769137120470574*^9, 3.769138854628717*^9, 3.769173861479599*^9,
3.769228309277001*^9, 3.769488536636314*^9, 3.769488775493741*^9,
3.7697540311309013`*^9, 3.7698493329786863`*^9, 3.769927143788989*^9,
3.7699298474833517`*^9, 3.7699349643369827`*^9, {3.7700014615163918`*^9,
3.770001486342461*^9}, 3.770089615045055*^9, 3.770091151999776*^9,
3.7700911832234373`*^9, 3.77009121560664*^9, 3.770092860611184*^9, {
3.770096733677506*^9, 3.770096743552308*^9}, 3.770175904771923*^9,
3.770355813747511*^9, 3.770368575883328*^9, 3.7704310126844397`*^9,
3.7704528133156967`*^9, 3.771218941512669*^9, 3.771673232190048*^9,
3.771907031907504*^9, 3.771907439029517*^9, 3.771915621974977*^9,
3.7721908396758137`*^9, 3.77219651143331*^9, 3.772247856888829*^9,
3.7722566296994762`*^9, 3.772259084152637*^9, 3.772348455386256*^9,
3.772367801749072*^9, 3.772370891811442*^9, {3.77242000146596*^9,
3.7724200050229883`*^9}, 3.7724281731485243`*^9, 3.772428410081149*^9,
3.772430677028174*^9, 3.772440462784586*^9, 3.772529273853674*^9,
3.7728759142457733`*^9, 3.7728780519980917`*^9, 3.772885195029847*^9,
3.772942823379051*^9, 3.7729434028029613`*^9, 3.772943543318931*^9,
3.7729440577300873`*^9, 3.772949628763401*^9, 3.772950463819756*^9,
3.772950515768664*^9, 3.773049286443142*^9, 3.7730496313722057`*^9,
3.773144076011161*^9, 3.773460624127591*^9, 3.7734639993643*^9,
3.7734713127079906`*^9, 3.773478884641347*^9, 3.773479639752248*^9,
3.7734819682595043`*^9, 3.773645277661384*^9, 3.773646804937044*^9,
3.773986250904626*^9, 3.7765100953056383`*^9, 3.776579083431468*^9,
3.776662521396657*^9, 3.7766687415984173`*^9, 3.776669989743709*^9,
3.7766754714989023`*^9, 3.7766763661356173`*^9, 3.776677554983851*^9,
3.776678156945221*^9, 3.776679730384227*^9, 3.7766797995218353`*^9,
3.776679922527975*^9, 3.776681466538108*^9, 3.776762885901209*^9,
3.7770248152256727`*^9, 3.77711769599905*^9, 3.7771849674192142`*^9,
3.77762564250346*^9, 3.777690720738963*^9, 3.7777253762606087`*^9,
3.777786888754923*^9, 3.777953025556863*^9, 3.777964327722644*^9,
3.778220763638605*^9, 3.778241133006884*^9, 3.77824221545049*^9,
3.778389851806761*^9, 3.778489688892495*^9, 3.778834393767021*^9,
3.778913036796289*^9, 3.7789186486061497`*^9, 3.778999683599246*^9,
3.779014947266182*^9, 3.779016358891302*^9, 3.779091498784801*^9,
3.779109063617063*^9, 3.77943804301024*^9, 3.7794390356224527`*^9,
3.779446129998694*^9, 3.779512777100968*^9, 3.779514823682631*^9,
3.779514859562222*^9, 3.779544709953924*^9, 3.7795459145467453`*^9,
3.779620255613166*^9, 3.779967327511076*^9, {3.780024275082299*^9,
3.7800242987111483`*^9}, {3.7800243789563026`*^9, 3.780024392712452*^9},
3.78005038726173*^9, 3.7800611518316393`*^9, 3.780115034568583*^9,
3.7801154414361067`*^9, {3.780997650057397*^9, 3.78099766439631*^9},
3.7824639813777742`*^9, 3.7839471685460453`*^9, 3.78504926862815*^9,
3.785052849502256*^9, 3.785057077348769*^9, {3.785057141637362*^9,
3.7850571599033833`*^9}, 3.785159356911799*^9, 3.7885250659535437`*^9,
3.7885885882949257`*^9, 3.789016374575654*^9, 3.7891011635298862`*^9,
3.797132624730732*^9, 3.7971350260463753`*^9, 3.799113710120328*^9,
3.7991139854914227`*^9, 3.799119939271639*^9, 3.799121028607518*^9,
3.79912168285885*^9, 3.799121965114169*^9, 3.7991220861947727`*^9, {
3.7991222754707747`*^9, 3.799122302610426*^9}, 3.800074675182968*^9,
3.8000888764634867`*^9, 3.8001483829496517`*^9, {3.800539328667665*^9,
3.80053935006597*^9}, 3.8006102047411127`*^9, 3.800768511041994*^9,
3.8007685584485817`*^9, 3.800798395343419*^9, 3.80083672174466*^9,
3.801289306970142*^9, 3.801723120960483*^9, 3.801805234128213*^9,
3.8019012394907827`*^9},
CellLabel->
"During evaluation of \
In[134]:=",ExpressionUUID->"b53e756f-7ced-4d44-9816-7a477e3725bc"],
Cell[BoxData[ InterpretationBox[
RowBox[{"\<\"mf = \"\>", "\[InvisibleSpace]", "0.952032939704`"}],
SequenceForm["mf = ", 0.952032939704],
Editable->False]], "Print",
CellChangeTimes->{{3.717223923734548*^9, 3.717223928214951*^9}, {
3.717223971239379*^9, 3.7172239766293097`*^9}, {3.717233931322109*^9,
3.7172339602600603`*^9}, 3.717234266025383*^9, 3.717389811826817*^9,
3.717397798635911*^9, 3.717409088225647*^9, 3.719042269449914*^9,
3.71972571685752*^9, 3.720497982809791*^9, 3.7206017968494253`*^9,
3.762057944992012*^9, 3.762171308612159*^9, 3.762185878959035*^9,
3.762230704000032*^9, 3.76257969223776*^9, 3.762579775163062*^9,
3.76260246063175*^9, 3.762603825586153*^9, {3.7626818315079517`*^9,
3.76268183964472*^9}, {3.762682055451223*^9, 3.762682065259396*^9},
3.764310838838461*^9, 3.76439221561644*^9, 3.764406820576201*^9,
3.7644071524893913`*^9, 3.764416941641247*^9, 3.764471035389357*^9,
3.764502311266172*^9, 3.76451318673353*^9, 3.764567390943874*^9,
3.768568011438259*^9, 3.768572299139023*^9, 3.768645828428453*^9,
3.7686498406439943`*^9, 3.7687018467692337`*^9, 3.7687065345918627`*^9,
3.768882355269024*^9, 3.76890964428026*^9, 3.769135242785296*^9,
3.769137120470574*^9, 3.769138854628717*^9, 3.769173861479599*^9,
3.769228309277001*^9, 3.769488536636314*^9, 3.769488775493741*^9,
3.7697540311309013`*^9, 3.7698493329786863`*^9, 3.769927143788989*^9,
3.7699298474833517`*^9, 3.7699349643369827`*^9, {3.7700014615163918`*^9,
3.770001486342461*^9}, 3.770089615045055*^9, 3.770091151999776*^9,
3.7700911832234373`*^9, 3.77009121560664*^9, 3.770092860611184*^9, {
3.770096733677506*^9, 3.770096743552308*^9}, 3.770175904771923*^9,
3.770355813747511*^9, 3.770368575883328*^9, 3.7704310126844397`*^9,
3.7704528133156967`*^9, 3.771218941512669*^9, 3.771673232190048*^9,
3.771907031907504*^9, 3.771907439029517*^9, 3.771915621974977*^9,
3.7721908396758137`*^9, 3.77219651143331*^9, 3.772247856888829*^9,
3.7722566296994762`*^9, 3.772259084152637*^9, 3.772348455386256*^9,
3.772367801749072*^9, 3.772370891811442*^9, {3.77242000146596*^9,
3.7724200050229883`*^9}, 3.7724281731485243`*^9, 3.772428410081149*^9,
3.772430677028174*^9, 3.772440462784586*^9, 3.772529273853674*^9,
3.7728759142457733`*^9, 3.7728780519980917`*^9, 3.772885195029847*^9,
3.772942823379051*^9, 3.7729434028029613`*^9, 3.772943543318931*^9,
3.7729440577300873`*^9, 3.772949628763401*^9, 3.772950463819756*^9,
3.772950515768664*^9, 3.773049286443142*^9, 3.7730496313722057`*^9,
3.773144076011161*^9, 3.773460624127591*^9, 3.7734639993643*^9,
3.7734713127079906`*^9, 3.773478884641347*^9, 3.773479639752248*^9,
3.7734819682595043`*^9, 3.773645277661384*^9, 3.773646804937044*^9,
3.773986250904626*^9, 3.7765100953056383`*^9, 3.776579083431468*^9,
3.776662521396657*^9, 3.7766687415984173`*^9, 3.776669989743709*^9,
3.7766754714989023`*^9, 3.7766763661356173`*^9, 3.776677554983851*^9,
3.776678156945221*^9, 3.776679730384227*^9, 3.7766797995218353`*^9,
3.776679922527975*^9, 3.776681466538108*^9, 3.776762885901209*^9,
3.7770248152256727`*^9, 3.77711769599905*^9, 3.7771849674192142`*^9,
3.77762564250346*^9, 3.777690720738963*^9, 3.7777253762606087`*^9,
3.777786888754923*^9, 3.777953025556863*^9, 3.777964327722644*^9,
3.778220763638605*^9, 3.778241133006884*^9, 3.77824221545049*^9,
3.778389851806761*^9, 3.778489688892495*^9, 3.778834393767021*^9,
3.778913036796289*^9, 3.7789186486061497`*^9, 3.778999683599246*^9,
3.779014947266182*^9, 3.779016358891302*^9, 3.779091498784801*^9,
3.779109063617063*^9, 3.77943804301024*^9, 3.7794390356224527`*^9,
3.779446129998694*^9, 3.779512777100968*^9, 3.779514823682631*^9,
3.779514859562222*^9, 3.779544709953924*^9, 3.7795459145467453`*^9,
3.779620255613166*^9, 3.779967327511076*^9, {3.780024275082299*^9,
3.7800242987111483`*^9}, {3.7800243789563026`*^9, 3.780024392712452*^9},
3.78005038726173*^9, 3.7800611518316393`*^9, 3.780115034568583*^9,
3.7801154414361067`*^9, {3.780997650057397*^9, 3.78099766439631*^9},
3.7824639813777742`*^9, 3.7839471685460453`*^9, 3.78504926862815*^9,
3.785052849502256*^9, 3.785057077348769*^9, {3.785057141637362*^9,
3.7850571599033833`*^9}, 3.785159356911799*^9, 3.7885250659535437`*^9,
3.7885885882949257`*^9, 3.789016374575654*^9, 3.7891011635298862`*^9,
3.797132624730732*^9, 3.7971350260463753`*^9, 3.799113710120328*^9,
3.7991139854914227`*^9, 3.799119939271639*^9, 3.799121028607518*^9,
3.79912168285885*^9, 3.799121965114169*^9, 3.7991220861947727`*^9, {
3.7991222754707747`*^9, 3.799122302610426*^9}, 3.800074675182968*^9,
3.8000888764634867`*^9, 3.8001483829496517`*^9, {3.800539328667665*^9,
3.80053935006597*^9}, 3.8006102047411127`*^9, 3.800768511041994*^9, 3.8007685584485817`*^9, 3.800798395343419*^9, 3.80083672174466*^9,
3.801289306970142*^9, 3.801723120960483*^9, 3.801805234128213*^9,
3.801901239492372*^9},
CellLabel->
"During evaluation of \
In[134]:=",ExpressionUUID->"92ca4415-ed1b-43c0-b8b4-efbe20d529f7"],
Cell[BoxData[
InterpretationBox[
RowBox[{"\<\"Tag = \"\>", "\[InvisibleSpace]", "\<\"SXS:BBH:0305\"\>"}],
SequenceForm["Tag = ", "SXS:BBH:0305"],
Editable->False]], "Print",
CellChangeTimes->{{3.717223923734548*^9, 3.717223928214951*^9}, {
3.717223971239379*^9, 3.7172239766293097`*^9}, {3.717233931322109*^9,
3.7172339602600603`*^9}, 3.717234266025383*^9, 3.717389811826817*^9,
3.717397798635911*^9, 3.717409088225647*^9, 3.719042269449914*^9,
3.71972571685752*^9, 3.720497982809791*^9, 3.7206017968494253`*^9,
3.762057944992012*^9, 3.762171308612159*^9, 3.762185878959035*^9,
3.762230704000032*^9, 3.76257969223776*^9, 3.762579775163062*^9,
3.76260246063175*^9, 3.762603825586153*^9, {3.7626818315079517`*^9,
3.76268183964472*^9}, {3.762682055451223*^9, 3.762682065259396*^9},
3.764310838838461*^9, 3.76439221561644*^9, 3.764406820576201*^9,
3.7644071524893913`*^9, 3.764416941641247*^9, 3.764471035389357*^9,
3.764502311266172*^9, 3.76451318673353*^9, 3.764567390943874*^9,
3.768568011438259*^9, 3.768572299139023*^9, 3.768645828428453*^9,
3.7686498406439943`*^9, 3.7687018467692337`*^9, 3.7687065345918627`*^9,
3.768882355269024*^9, 3.76890964428026*^9, 3.769135242785296*^9,
3.769137120470574*^9, 3.769138854628717*^9, 3.769173861479599*^9,
3.769228309277001*^9, 3.769488536636314*^9, 3.769488775493741*^9,
3.7697540311309013`*^9, 3.7698493329786863`*^9, 3.769927143788989*^9,
3.7699298474833517`*^9, 3.7699349643369827`*^9, {3.7700014615163918`*^9,
3.770001486342461*^9}, 3.770089615045055*^9, 3.770091151999776*^9,
3.7700911832234373`*^9, 3.77009121560664*^9, 3.770092860611184*^9, {
3.770096733677506*^9, 3.770096743552308*^9}, 3.770175904771923*^9,
3.770355813747511*^9, 3.770368575883328*^9, 3.7704310126844397`*^9,
3.7704528133156967`*^9, 3.771218941512669*^9, 3.771673232190048*^9,
3.771907031907504*^9, 3.771907439029517*^9, 3.771915621974977*^9,
3.7721908396758137`*^9, 3.77219651143331*^9, 3.772247856888829*^9,
3.7722566296994762`*^9, 3.772259084152637*^9, 3.772348455386256*^9,
3.772367801749072*^9, 3.772370891811442*^9, {3.77242000146596*^9,
3.7724200050229883`*^9}, 3.7724281731485243`*^9, 3.772428410081149*^9,
3.772430677028174*^9, 3.772440462784586*^9, 3.772529273853674*^9,
3.7728759142457733`*^9, 3.7728780519980917`*^9, 3.772885195029847*^9,
3.772942823379051*^9, 3.7729434028029613`*^9, 3.772943543318931*^9,
3.7729440577300873`*^9, 3.772949628763401*^9, 3.772950463819756*^9,
3.772950515768664*^9, 3.773049286443142*^9, 3.7730496313722057`*^9,
3.773144076011161*^9, 3.773460624127591*^9, 3.7734639993643*^9,
3.7734713127079906`*^9, 3.773478884641347*^9, 3.773479639752248*^9,
3.7734819682595043`*^9, 3.773645277661384*^9, 3.773646804937044*^9,
3.773986250904626*^9, 3.7765100953056383`*^9, 3.776579083431468*^9,
3.776662521396657*^9, 3.7766687415984173`*^9, 3.776669989743709*^9,
3.7766754714989023`*^9, 3.7766763661356173`*^9, 3.776677554983851*^9,
3.776678156945221*^9, 3.776679730384227*^9, 3.7766797995218353`*^9,
3.776679922527975*^9, 3.776681466538108*^9, 3.776762885901209*^9,
3.7770248152256727`*^9, 3.77711769599905*^9, 3.7771849674192142`*^9,
3.77762564250346*^9, 3.777690720738963*^9, 3.7777253762606087`*^9,
3.777786888754923*^9, 3.777953025556863*^9, 3.777964327722644*^9,
3.778220763638605*^9, 3.778241133006884*^9, 3.77824221545049*^9,
3.778389851806761*^9, 3.778489688892495*^9, 3.778834393767021*^9,
3.778913036796289*^9, 3.7789186486061497`*^9, 3.778999683599246*^9,
3.779014947266182*^9, 3.779016358891302*^9, 3.779091498784801*^9,
3.779109063617063*^9, 3.77943804301024*^9, 3.7794390356224527`*^9,
3.779446129998694*^9, 3.779512777100968*^9, 3.779514823682631*^9,
3.779514859562222*^9, 3.779544709953924*^9, 3.7795459145467453`*^9,
3.779620255613166*^9, 3.779967327511076*^9, {3.780024275082299*^9,
3.7800242987111483`*^9}, {3.7800243789563026`*^9, 3.780024392712452*^9},
3.78005038726173*^9, 3.7800611518316393`*^9, 3.780115034568583*^9,
3.7801154414361067`*^9, {3.780997650057397*^9, 3.78099766439631*^9},
3.7824639813777742`*^9, 3.7839471685460453`*^9, 3.78504926862815*^9,
3.785052849502256*^9, 3.785057077348769*^9, {3.785057141637362*^9, 3.7850571599033833`*^9}, 3.785159356911799*^9, 3.7885250659535437`*^9,
3.7885885882949257`*^9, 3.789016374575654*^9, 3.7891011635298862`*^9,
3.797132624730732*^9, 3.7971350260463753`*^9, 3.799113710120328*^9,
3.7991139854914227`*^9, 3.799119939271639*^9, 3.799121028607518*^9,
3.79912168285885*^9, 3.799121965114169*^9, 3.7991220861947727`*^9, {
3.7991222754707747`*^9, 3.799122302610426*^9}, 3.800074675182968*^9,
3.8000888764634867`*^9, 3.8001483829496517`*^9, {3.800539328667665*^9,
3.80053935006597*^9}, 3.8006102047411127`*^9, 3.800768511041994*^9,
3.8007685584485817`*^9, 3.800798395343419*^9, 3.80083672174466*^9,
3.801289306970142*^9, 3.801723120960483*^9, 3.801805234128213*^9,
3.8019012394938803`*^9},
CellLabel->
"During evaluation of \
In[134]:=",ExpressionUUID->"c39e3e62-f655-4ce9-9bde-77a8f501bd95"]
}, Open ]]
}, Open ]],
Cell[BoxData[
RowBox[{
RowBox[{"sxsrhs", "=",
RowBox[{"Flatten", "[",
RowBox[{
RowBox[{
RowBox[{
RowBox[{"Conjugate", "/@",
RowBox[{"GetAsymptoticMultiMode", "[",
RowBox[{"#", ",", "2", ",", "modes", ",",
RowBox[{"\"\<ReSample\>\"", "\[Rule]", "True"}]}], "]"}]}], "&"}], "/@",
"mysxscaserh"}], ",", "1"}], "]"}]}], ";"}]], "Code",
CellChangeTimes->{{3.7171635521596403`*^9, 3.717163558141595*^9}, {
3.7172166219975853`*^9, 3.717216627320222*^9}, {3.7172194728512707`*^9,
3.717219490830843*^9}, {3.717223993502838*^9, 3.717223996533388*^9}, {
3.72060202502007*^9, 3.720602032786186*^9}, {3.764404927143031*^9,
3.764404936617689*^9}, {3.764405546249731*^9, 3.764405548934586*^9}, {
3.764406129145973*^9, 3.764406130113237*^9}, {3.764406181485956*^9,
3.764406184078907*^9}, {3.770355814410892*^9, 3.7703558160102386`*^9}, {
3.772949641678318*^9, 3.772949677949027*^9}, {3.773653255970372*^9,
3.773653259449843*^9}, {3.7736533620280743`*^9, 3.773653364920361*^9}, {
3.773986356744377*^9, 3.7739863583243027`*^9}, {3.7991140470420923`*^9,
3.7991140508638773`*^9}},
CellLabel->
"In[143]:=",ExpressionUUID->"96195c64-1c6d-4895-96a0-d93b22b41c3c"]
}, Closed]]
}, Open ]]
}, Open ]],
Cell[CellGroupData[{
Cell["Results for the paper", "Section",
CellChangeTimes->{{3.717216037964575*^9, 3.717216039170602*^9}, {
3.762579673099522*^9, 3.762579673682219*^9}, {3.76369974318791*^9,
3.763699749547122*^9}, {3.7850573300501738`*^9,
3.785057336126431*^9}},ExpressionUUID->"55e918b2-ca1f-498b-a377-\
bd7458fe203c"],
Cell[CellGroupData[{
Cell["Vary the fundamental frequency", "Subsubsection",
CellChangeTimes->{{3.799115057201908*^9, 3.799115066138445*^9}, {
3.7991152472676363`*^9, 3.799115260025752*^9}, {3.799115347398109*^9,
3.799115355434946*^9}},ExpressionUUID->"d264caca-e5b5-4c88-941a-\
93010756d554"],
Cell[BoxData[{
RowBox[{"t0", "=",
RowBox[{
RowBox[{"TimeOfMaximum", "[",
RowBox[{"sxsrhs", "[",
RowBox[{"[", "1", "]"}], "]"}], "]"}], "+", "tshift"}]}], "\n",
RowBox[{ RowBox[{"data", "=",
RowBox[{"Select", "[",
RowBox[{
RowBox[{"sxsrhs", "[",
RowBox[{"[", "1", "]"}], "]"}], ",",
RowBox[{
RowBox[{
RowBox[{"#", "[",
RowBox[{"[", "1", "]"}], "]"}], "\[GreaterEqual]", " ", "t0"}],
"&"}]}], "]"}]}], ";"}]}], "Code",
CellChangeTimes->{{3.801901370029584*^9, 3.801901406233098*^9}, {
3.801901893503097*^9, 3.801901895231626*^9}, {3.801905966678508*^9,
3.801905966881785*^9}, {3.8019061500176573`*^9, 3.8019061832460957`*^9}, {
3.80190671422969*^9, 3.8019067208765078`*^9}, {3.80190675541617*^9,
3.801906761681551*^9}, {3.801996962552052*^9,
3.801996962743515*^9}},ExpressionUUID->"3c15dded-a57b-4656-b83e-\
c264f7e27e40"],
Cell[BoxData[
RowBox[{
RowBox[{"tab", "=",
RowBox[{"Table", "[", "\n",
RowBox[{
RowBox[{
RowBox[{"ansatz", "=",
RowBox[{"OvertoneModelV2", "[",
RowBox[{"ntones", ",",
RowBox[{"{",
RowBox[{"\[Eta]", ",", "\[Chi]1", ",", "\[Chi]2"}], "}"}], ",",
"t0", ",",
RowBox[{"\"\<Fit\[Alpha]\>\"", "\[Rule]",
RowBox[{"{", "}"}]}], ",",
RowBox[{"\"\<Vary\[Omega]\>\"", "\[Rule]", "True"}], ",",
RowBox[{"\"\<\[Omega]val\>\"", "\[Rule]",
RowBox[{"{",
RowBox[{
RowBox[{"-", "\[Omega]fact"}], ",", "\[Omega]fact"}], "}"}]}],
",",
RowBox[{"\"\<Export_\[Omega]val\>\"", "\[Rule]", "True"}], ",",
RowBox[{"\"\<ModesData\>\"", "->", "Modedata"}]}], "]"}]}], ";",
"\[IndentingNewLine]",
RowBox[{"cfit", "=",
RowBox[{"NonlinearModelFit", "[",
RowBox[{"data", ",",
RowBox[{"ansatz", "[",
RowBox[{"[", "1", "]"}], "]"}], ",",
RowBox[{"{",
RowBox[{
"x0", ",", "x1", ",", "x2", ",", "x3", ",", "x4", ",", "x5", ",",
"x6", ",", "x7"}], "}"}], ",", "t"}], "]"}]}], ";",
"\[IndentingNewLine]",
RowBox[{"cfitd", "=",
RowBox[{"Transpose", "[",
RowBox[{"{",
RowBox[{
RowBox[{"data", "[",
RowBox[{"[",
RowBox[{"All", ",", "1"}], "]"}], "]"}], ",",
RowBox[{
RowBox[{"Normal", "[", "cfit", "]"}], "/.",
RowBox[{"t", "->",
RowBox[{"data", "[",
RowBox[{"[",
RowBox[{"All", ",", "1"}], "]"}], "]"}]}]}]}], "}"}], "]"}]}],
";", "\n",
RowBox[{"{",
RowBox[{
RowBox[{"ansatz", "[",
RowBox[{"[", "2", "]"}], "]"}], ",",
RowBox[{ RowBox[{"{",
RowBox[{
"x0", ",", "x1", ",", "x2", ",", "x3", ",", "x4", ",", "x5", ",",
"x6", ",", "x7"}], "}"}], "/.",
RowBox[{"cfit", "[", "\"\<BestFitParameters\>\"", "]"}]}], ",",
RowBox[{"1", "-",
RowBox[{"EasyMatchT", "[",
RowBox[{"data", ",", "cfitd", ",", "t0", ",",
RowBox[{"t0", "+", "90"}]}], "]"}]}]}], "}"}]}], ",",
RowBox[{"{",
RowBox[{"j", ",", "npoints"}], "}"}]}], "]"}]}], ";"}]], "Code",
CellChangeTimes->{
3.8019012567182837`*^9, {3.801901306368492*^9, 3.801901347644397*^9},
3.801901410887135*^9, {3.80190187722987*^9, 3.801901879607212*^9}, {
3.801905987882124*^9, 3.801905989699266*^9}, {3.8019062871664743`*^9,
3.801906312482821*^9}, {3.8019063425225697`*^9, 3.801906429559026*^9}, {
3.801906653354439*^9, 3.801906654839933*^9}, {3.801906698689931*^9,
3.8019067056553392`*^9}, {3.801906766302006*^9, 3.801906766497856*^9}, {
3.801906809977264*^9, 3.8019068263340263`*^9}, {3.801906889494577*^9,
3.8019069354645443`*^9}, {3.801906976150055*^9, 3.801906980794434*^9}, {
3.8019970228363132`*^9,
3.8019970236951714`*^9}},ExpressionUUID->"b7995799-c38a-4b15-8646-\
daf87c3c734d"],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"Export", "[",
RowBox[{
RowBox[{"rootpath", "<>", "\"\<plots/freq_bfitpars_mmatch.dat\>\""}], ",",
"tab"}], "]"}]], "Code",
CellChangeTimes->{{3.8019065069435463`*^9, 3.80190652588253*^9}, {
3.801906563232386*^9, 3.80190658024755*^9}, {3.801906639024638*^9,
3.801906645087222*^9}, {3.801906945654804*^9, 3.801906958171398*^9}},
CellLabel->
"In[221]:=",ExpressionUUID->"599f1833-d91e-4999-b904-ddbb9f32dd32"],
Cell[BoxData["\<\"~/git/rdstackingproject/plots/freq_bfitpars_mmatch.dat\"\>"]\
, "Output",
CellChangeTimes->{3.801906632797943*^9, 3.801906958786167*^9},
CellLabel->
"Out[221]=",ExpressionUUID->"a4b2021d-c620-4905-999e-ab0338c84924"]
}, Open ]]
}, Open ]]
}, Open ]]
}, Open ]]
},
WindowSize->{Full, Full},
WindowMargins->{{436, Automatic}, {Automatic, 50}},
FrontEndVersion->"12.0 for Mac OS X x86 (64-bit) (April 8, 2019)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)
(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[CellGroupData[{
Cell[580, 22, 168, 3, 98, "Title",ExpressionUUID->"b53761d4-1b24-4b86-b93e-9c44d4b29d65"],
Cell[751, 27, 250, 5, 30, "Input",ExpressionUUID->"980e79f6-9c07-4072-bd1e-300ac946e8c3"],
Cell[CellGroupData[{
Cell[1026, 36, 201, 4, 67, "Section",ExpressionUUID->"fe9e23be-b684-46cc-9b54-957e00ada2e3"],
Cell[CellGroupData[{
Cell[1252, 44, 232, 4, 54, "Subsection",ExpressionUUID->"e952a5b9-bbfe-43a9-9668-8567ca76bf42"],Cell[1487, 50, 554, 10, 52, "Code",ExpressionUUID->"ded0ca57-00ba-4008-906d-ee5041d06f91"],
Cell[2044, 62, 712, 13, 91, "Code",ExpressionUUID->"df18942a-15e3-4e8b-880e-923fca9e18d7"],
Cell[2759, 77, 590, 17, 129, "Code",ExpressionUUID->"d0e5d8fb-ddc3-4ce5-93a5-2a83bdf621d4"]
}, Closed]],
Cell[CellGroupData[{
Cell[3386, 99, 187, 3, 38, "Subsection",ExpressionUUID->"46366112-0a9f-4a64-97d3-72f08bceb964"],
Cell[CellGroupData[{
Cell[3598, 106, 1861, 37, 186, "Code",ExpressionUUID->"2c37e1b9-5cb9-4980-a108-163d386515de"],
Cell[5462, 145, 208, 4, 34, "Output",ExpressionUUID->"04cad810-7ff2-4fa4-88b1-a3cb2e1eaf60"],
Cell[5673, 151, 240, 4, 34, "Output",ExpressionUUID->"3a282303-63af-4a85-87c0-523e68dad528"]
}, Open ]]
}, Closed]],
Cell[CellGroupData[{
Cell[5962, 161, 162, 3, 38, "Subsection",ExpressionUUID->"2b160476-27ad-417e-ba6b-aba2f29dfd9c"],
Cell[CellGroupData[{
Cell[6149, 168, 1113, 27, 148, "Code",ExpressionUUID->"91d0d186-3694-4095-b0ad-f5476cf3486b"],
Cell[7265, 197, 1782, 27, 34, "Output",ExpressionUUID->"4661e259-0c1d-4e90-b4af-b76e6ee6e71c"],
Cell[9050, 226, 27811, 724, 34, "Output",ExpressionUUID->"d642568b-341d-4984-ae2d-30a8eeb30577"]
}, Open ]]
}, Closed]],
Cell[CellGroupData[{
Cell[36910, 956, 165, 3, 38, "Subsection",ExpressionUUID->"44e4b7cd-d026-4079-8fd8-3f63cf14cfa7"],
Cell[37078, 961, 4114, 106, 338, "Code",ExpressionUUID->"48958061-652d-46b3-9aca-2ec5e658ec2a"],
Cell[41195, 1069, 3650, 91, 281, "Code",ExpressionUUID->"0204b879-5f51-44f5-8a4f-95cc490da72f"]
}, Closed]],
Cell[CellGroupData[{
Cell[44882, 1165, 157, 3, 38, "Subsection",ExpressionUUID->"8488942b-dc04-47e9-8ca6-119712f58bc1"],
Cell[CellGroupData[{
Cell[45064, 1172, 163, 3, 45, "Subsubsection",ExpressionUUID->"f0a39a1f-5153-4be0-b03d-fd24806eca3d"],
Cell[45230, 1177, 1075, 29, 52, "Input",ExpressionUUID->"7694c13e-a765-44ce-9adc-09e59459c3d2"],
Cell[CellGroupData[{
Cell[46330, 1210, 1237, 31, 72, "Code",ExpressionUUID->"763f62ae-036b-47b4-86e4-d1b0e24fc92b"],
Cell[47570, 1243, 1110, 27, 34, "Output",ExpressionUUID->"81efd3bd-daab-4062-8b73-d2a7435e97df"]
}, Open ]],
Cell[CellGroupData[{
Cell[48717, 1275, 2162, 50, 72, "Code",ExpressionUUID->"37d6369f-02c3-4af4-91b6-0490e703d64a"],
Cell[50882, 1327, 1981, 31, 56, "Output",ExpressionUUID->"1fd520ec-5fee-43d1-967d-a3a48f981509"]
}, Open ]],
Cell[52878, 1361, 416, 10, 72, "Code",ExpressionUUID->"827b1be5-400a-4c10-97cf-df0471536b52"]
}, Closed]],
Cell[CellGroupData[{
Cell[53331, 1376, 168, 3, 37, "Subsubsection",ExpressionUUID->"e69b91c9-6c5f-47bb-9540-07f42335cb1a"],
Cell[CellGroupData[{
Cell[53524, 1383, 777, 14, 52, "Code",ExpressionUUID->"b361f1b5-e01e-4264-bf04-6b665f472fd9"],
Cell[54304, 1399, 4850, 71, 34, "Output",ExpressionUUID->"5a59dd1d-dcee-4d7a-954f-f98d277d06f5"]
}, Open ]],
Cell[CellGroupData[{
Cell[59191, 1475, 1643, 35, 91, "Code",ExpressionUUID->"5584846d-d985-4ced-98a1-ac35c86181ea"],
Cell[60837, 1512, 5241, 75, 34, "Output",ExpressionUUID->"41029220-aae2-4aaa-a49d-ef2c6d004855"],
Cell[66081, 1589, 5189, 73, 34, "Output",ExpressionUUID->"93c7ae4a-68f8-4aff-99b3-3854a0d267da"]
}, Open ]],
Cell[CellGroupData[{
Cell[71307, 1667, 3350, 90, 205, "Code",ExpressionUUID->"6e14f11c-a2d0-4ce1-82ab-1a2a53c717dd"],
Cell[CellGroupData[{
Cell[74682, 1761, 5190, 76, 24, "Print",ExpressionUUID->"b0646428-ecbf-4825-87c0-5f5a9a25a494"],
Cell[79875, 1839, 5188, 76, 24, "Print",ExpressionUUID->"a5cf05bc-2cbc-4f46-a364-fa765f897969"],
Cell[85066, 1917, 5189, 76, 24, "Print",ExpressionUUID->"21ae58be-bb6b-4910-95e4-c68634b6ddf2"],
Cell[90258, 1995, 5209, 77, 24, "Print",ExpressionUUID->"5f689c14-b703-4074-a06a-efd763c742b7"],
Cell[95470, 2074, 5192, 76, 24, "Print",ExpressionUUID->"3f8c41fc-6da7-4446-ba56-ed613b63b498"],
Cell[100665, 2152, 5219, 77, 24, "Print",ExpressionUUID->"699de0d3-51d3-456d-b993-57ce112918b1"],
Cell[105887, 2231, 5179, 76, 24, "Print",ExpressionUUID->"b53e756f-7ced-4d44-9816-7a477e3725bc"],
Cell[111069, 2309, 5177, 76, 24, "Print",ExpressionUUID->"92ca4415-ed1b-43c0-b8b4-efbe20d529f7"],
Cell[116249, 2387, 5186, 76, 24, "Print",ExpressionUUID->"c39e3e62-f655-4ce9-9bde-77a8f501bd95"]
}, Open ]]
}, Open ]],
Cell[121462, 2467, 1236, 24, 52, "Code",ExpressionUUID->"96195c64-1c6d-4895-96a0-d93b22b41c3c"]
}, Closed]]
}, Open ]]
}, Open ]],
Cell[CellGroupData[{Cell[122759, 2498, 311, 5, 67, "Section",ExpressionUUID->"55e918b2-ca1f-498b-a377-bd7458fe203c"],
Cell[CellGroupData[{
Cell[123095, 2507, 278, 4, 45, "Subsubsection",ExpressionUUID->"d264caca-e5b5-4c88-941a-93010756d554"],
Cell[123376, 2513, 912, 23, 72, "Code",ExpressionUUID->"3c15dded-a57b-4656-b83e-c264f7e27e40"],
Cell[124291, 2538, 3099, 74, 129, "Code",ExpressionUUID->"b7995799-c38a-4b15-8646-daf87c3c734d"],
Cell[CellGroupData[{
Cell[127415, 2616, 449, 9, 52, "Code",ExpressionUUID->"599f1833-d91e-4999-b904-ddbb9f32dd32"],
Cell[127867, 2627, 238, 4, 34, "Output",ExpressionUUID->"a4b2021d-c620-4905-999e-ab0338c84924"]
}, Open ]]
}, Open ]]
}, Open ]]
}, Open ]]
}
]
*)