(* 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 ]] } ] *)