diff options
Diffstat (limited to 'mathemathica_fwm')
| -rw-r--r-- | mathemathica_fwm/RbXMDSSetup1.nb | 215 |
1 files changed, 165 insertions, 50 deletions
diff --git a/mathemathica_fwm/RbXMDSSetup1.nb b/mathemathica_fwm/RbXMDSSetup1.nb index 58a451b..aae1a2d 100644 --- a/mathemathica_fwm/RbXMDSSetup1.nb +++ b/mathemathica_fwm/RbXMDSSetup1.nb @@ -1692,7 +1692,7 @@ Cell[BoxData[{ RowBox[{"PopulationStyle", "\[Rule]", RowBox[{"PointSize", "[", ".018", "]"}]}]}], "]"}]}], "Input"], -Cell[BoxData[{ +Cell[BoxData[ RowBox[{ RowBox[{"table", "=", RowBox[{"Flatten", "[", @@ -1758,38 +1758,101 @@ Cell[BoxData[{ SuperscriptBox["10", RowBox[{"-", "4"}]]}]}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", - RowBox[{ - RowBox[{"{", - RowBox[{"e1m", ",", "e3m", ",", + RowBox[{"absorption", "=", + RowBox[{"dE", "[", + RowBox[{ + RowBox[{ + SubscriptBox["\[CapitalOmega]", + RowBox[{"2", ",", "p"}]], "/.", "params"}], ",", + RowBox[{"points", "[", + RowBox[{"[", + RowBox[{"2", ",", "2", ",", "1"}], "]"}], "]"}]}], "]"}]}], ";", + "\[IndentingNewLine]", + RowBox[{"slope", "=", + RowBox[{ + RowBox[{"(", RowBox[{ - RowBox[{"Subtract", "@@", - RowBox[{"(", + RowBox[{"d\[Phi]", "[", + RowBox[{ RowBox[{ - RowBox[{ - RowBox[{"d\[Phi]", "[", - RowBox[{ - RowBox[{ - SubscriptBox["\[CapitalOmega]", - RowBox[{"2", ",", "p"}]], "/.", "params"}], ",", "#"}], - "]"}], "&"}], "/@", - RowBox[{"points", "[", - RowBox[{"[", - RowBox[{"All", ",", "2", ",", "1"}], "]"}], "]"}]}], ")"}]}], - "/", - RowBox[{"Subtract", "@@", - RowBox[{"points", "[", - RowBox[{"[", - RowBox[{"All", ",", "1"}], "]"}], "]"}]}]}]}], "}"}], "/.", - "params"}]}], ",", + SubscriptBox["\[CapitalOmega]", + RowBox[{"2", ",", "p"}]], "/.", "params"}], ",", + RowBox[{"points", "[", + RowBox[{"[", + RowBox[{"2", ",", "2", ",", "1"}], "]"}], "]"}]}], "]"}], "-", + RowBox[{"d\[Phi]", "[", + RowBox[{ + RowBox[{ + SubscriptBox["\[CapitalOmega]", + RowBox[{"2", ",", "p"}]], "/.", "params"}], ",", + RowBox[{"points", "[", + RowBox[{"[", + RowBox[{"1", ",", "2", ",", "1"}], "]"}], "]"}]}], "]"}]}], + ")"}], "/", + RowBox[{"(", + RowBox[{ + RowBox[{"points", "[", + RowBox[{"[", + RowBox[{"2", ",", "1"}], "]"}], "]"}], "-", + RowBox[{"points", "[", + RowBox[{"[", + RowBox[{"1", ",", "1"}], "]"}], "]"}]}], ")"}]}]}], ";", + "\[IndentingNewLine]", + RowBox[{ + RowBox[{"{", + RowBox[{"e1m", ",", "e3m", ",", "absorption", ",", "slope"}], "}"}], + "/.", "params"}]}], ",", RowBox[{"{", RowBox[{"e1m", ",", RowBox[{"-", "1."}], ",", "3.", ",", ".2"}], "}"}], ",", RowBox[{"{", RowBox[{"e3m", ",", RowBox[{"-", "1."}], ",", "3.", ",", ".2"}], "}"}]}], "]"}], ",", - "1"}], "]"}]}], ";"}], "\n", + "1"}], "]"}]}], ";"}]], "Input"], + +Cell["\<\ +Plot fractional absorption of the electric field per absorption length.\ +\>", "Text"], + +Cell[BoxData[ RowBox[{"ListContourPlot", "[", - RowBox[{"table", ",", + RowBox[{ + RowBox[{"table", "[", + RowBox[{"[", + RowBox[{"All", ",", + RowBox[{"{", + RowBox[{"1", ",", "2", ",", "3"}], "}"}]}], "]"}], "]"}], ",", + RowBox[{"PlotRange", "\[Rule]", "All"}], ",", + RowBox[{"ContourLabels", "\[Rule]", "All"}], ",", + RowBox[{"FrameLabel", "\[Rule]", + RowBox[{"{", + RowBox[{ + "\"\<log(\!\(\*SubscriptBox[\(\[CapitalOmega]\), \(1, m\)]\) [Mrad/sec])\ +\>\"", ",", + "\"\<log(\!\(\*SubscriptBox[\(\[CapitalOmega]\), \(3, m\)]\) \ +[Mrad/sec])\>\""}], "}"}]}], ",", + RowBox[{"BaseStyle", "\[Rule]", "14"}], ",", + RowBox[{"Contours", "\[Rule]", + RowBox[{"Join", "[", + RowBox[{ + RowBox[{"{", + RowBox[{".001", ",", ".01"}], "}"}], ",", + RowBox[{"Range", "[", + RowBox[{".1", ",", "1", ",", ".1"}], "]"}]}], "]"}]}]}], + "]"}]], "Input"], + +Cell["\<\ +Plot slope of the phase shift in rad/(Mrad/s) per absorption length\ +\>", "Text"], + +Cell[BoxData[ + RowBox[{"ListContourPlot", "[", + RowBox[{ + RowBox[{"table", "[", + RowBox[{"[", + RowBox[{"All", ",", + RowBox[{"{", + RowBox[{"1", ",", "2", ",", "4"}], "}"}]}], "]"}], "]"}], ",", RowBox[{"PlotRange", "\[Rule]", "All"}], ",", RowBox[{"ContourLabels", "\[Rule]", "All"}], ",", RowBox[{"FrameLabel", "\[Rule]", @@ -1811,7 +1874,7 @@ Cell[BoxData[{ RowBox[{ RowBox[{"-", "4."}], ",", RowBox[{"-", "2."}], ",", "1"}], "]"}]}]}], "]"}]}]}], - "]"}]}], "Input"], + "]"}]], "Input"], Cell[BoxData[ RowBox[{ @@ -1879,29 +1942,50 @@ Cell[BoxData[ SuperscriptBox["10", RowBox[{"-", "4"}]]}]}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", - RowBox[{ - RowBox[{"{", - RowBox[{"e1m", ",", "e3m", ",", + RowBox[{"absorption", "=", + RowBox[{"dE", "[", + RowBox[{ RowBox[{ - RowBox[{"Subtract", "@@", - RowBox[{"(", + SubscriptBox["\[CapitalOmega]", + RowBox[{"2", ",", "p"}]], "/.", "params"}], ",", + RowBox[{"points", "[", + RowBox[{"[", + RowBox[{"2", ",", "2", ",", "1"}], "]"}], "]"}]}], "]"}]}], ";", + "\[IndentingNewLine]", + RowBox[{"slope", "=", + RowBox[{ + RowBox[{"(", + RowBox[{ + RowBox[{"d\[Phi]", "[", + RowBox[{ + RowBox[{ + SubscriptBox["\[CapitalOmega]", + RowBox[{"2", ",", "p"}]], "/.", "params"}], ",", + RowBox[{"points", "[", + RowBox[{"[", + RowBox[{"2", ",", "2", ",", "1"}], "]"}], "]"}]}], "]"}], "-", + RowBox[{"d\[Phi]", "[", + RowBox[{ RowBox[{ - RowBox[{ - RowBox[{"d\[Phi]", "[", - RowBox[{ - RowBox[{ - SubscriptBox["\[CapitalOmega]", - RowBox[{"2", ",", "p"}]], "/.", "params"}], ",", "#"}], - "]"}], "&"}], "/@", - RowBox[{"points", "[", - RowBox[{"[", - RowBox[{"All", ",", "2", ",", "1"}], "]"}], "]"}]}], ")"}]}], - "/", - RowBox[{"Subtract", "@@", - RowBox[{"points", "[", - RowBox[{"[", - RowBox[{"All", ",", "1"}], "]"}], "]"}]}]}]}], "}"}], "/.", - "params"}]}], ",", + SubscriptBox["\[CapitalOmega]", + RowBox[{"2", ",", "p"}]], "/.", "params"}], ",", + RowBox[{"points", "[", + RowBox[{"[", + RowBox[{"1", ",", "2", ",", "1"}], "]"}], "]"}]}], "]"}]}], + ")"}], "/", + RowBox[{"(", + RowBox[{ + RowBox[{"points", "[", + RowBox[{"[", + RowBox[{"2", ",", "1"}], "]"}], "]"}], "-", + RowBox[{"points", "[", + RowBox[{"[", + RowBox[{"1", ",", "1"}], "]"}], "]"}]}], ")"}]}]}], ";", + "\[IndentingNewLine]", + RowBox[{ + RowBox[{"{", + RowBox[{"e1m", ",", "e3m", ",", "absorption", ",", "slope"}], "}"}], + "/.", "params"}]}], ",", RowBox[{"{", RowBox[{"e1m", ",", RowBox[{"-", "1."}], ",", "3.", ",", ".2"}], "}"}], ",", @@ -1912,7 +1996,39 @@ Cell[BoxData[ Cell[BoxData[ RowBox[{"ListContourPlot", "[", - RowBox[{"table", ",", + RowBox[{ + RowBox[{"table", "[", + RowBox[{"[", + RowBox[{"All", ",", + RowBox[{"{", + RowBox[{"1", ",", "2", ",", "3"}], "}"}]}], "]"}], "]"}], ",", + RowBox[{"PlotRange", "\[Rule]", "All"}], ",", + RowBox[{"ContourLabels", "\[Rule]", "All"}], ",", + RowBox[{"FrameLabel", "\[Rule]", + RowBox[{"{", + RowBox[{ + "\"\<log(\!\(\*SubscriptBox[\(\[CapitalOmega]\), \(1, m\)]\) [Mrad/sec])\ +\>\"", ",", + "\"\<log(\!\(\*SubscriptBox[\(\[CapitalOmega]\), \(3, m\)]\) \ +[Mrad/sec])\>\""}], "}"}]}], ",", + RowBox[{"BaseStyle", "\[Rule]", "14"}], ",", + RowBox[{"Contours", "\[Rule]", + RowBox[{"Join", "[", + RowBox[{ + RowBox[{"{", + RowBox[{".001", ",", ".01"}], "}"}], ",", + RowBox[{"Range", "[", + RowBox[{".1", ",", "1", ",", ".1"}], "]"}]}], "]"}]}]}], + "]"}]], "Input"], + +Cell[BoxData[ + RowBox[{"ListContourPlot", "[", + RowBox[{ + RowBox[{"table", "[", + RowBox[{"[", + RowBox[{"All", ",", + RowBox[{"{", + RowBox[{"1", ",", "2", ",", "4"}], "}"}]}], "]"}], "]"}], ",", RowBox[{"PlotRange", "\[Rule]", "All"}], ",", RowBox[{"ContourLabels", "\[Rule]", "All"}], ",", RowBox[{"FrameLabel", "\[Rule]", @@ -2963,8 +3079,7 @@ Cell[BoxData[{ RowBox[{ SubscriptBox["\[CapitalOmega]", RowBox[{"j_", ",", "m"}]], "[", "z", "]"}], "\[Rule]", - RowBox[{"Em", "[", "j", "]"}]}]}], "}"}]}], "/.", "labelrep"}], "/.", - + RowBox[{"Em", "[", "j", "]"}]}]}], "}"}]}], "/.", "labelrep"}], "/.", RowBox[{ RowBox[{"Complex", "[", RowBox[{"0", ",", "a_"}], "]"}], "\[Rule]", @@ -3079,7 +3194,7 @@ Cell[BoxData[{ 3.56319258550938*^9}}] }, Open ]] }, -WindowSize->{838, 902}, +WindowSize->{998, 902}, WindowMargins->{{-2, Automatic}, {Automatic, 0}}, PrivateNotebookOptions->{"FileOutlineCache"->False}, ShowSelection->True, |