Added contour plots of absorption.

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,