summaryrefslogtreecommitdiff
path: root/mathemathica_fwm/RbXMDSSetup1.nb
diff options
context:
space:
mode:
Diffstat (limited to 'mathemathica_fwm/RbXMDSSetup1.nb')
-rw-r--r--mathemathica_fwm/RbXMDSSetup1.nb215
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,
generated by cgit v1.2.3 (git 2.39.1) at 2024-11-16 07:22:13 +0000