1 files changed, 70 insertions, 0 deletions
diff --git a/faraday/output_psr_results_vs_power.m b/faraday/output_psr_results_vs_power.m
new file mode 100644
index 0000000..2d28565
--- /dev/null
+++ b/faraday/output_psr_results_vs_power.m
@@ -0,0 +1,70 @@
+1;
+
+
+load '/tmp/xi_vs_power.mat' ;
+
+Er=(1+I*xi_right)*E_field_pos_freq.right;
+El=(1+I*xi_left) *E_field_pos_freq.left;
+
+Ex=(Er+El)/sqrt(2);
+Ey=I*(Er-El)/sqrt(2);
+
+%extra rotation to compensate rotation due to ellipticity
+% actually no need for it since x-polarization shifts by positive phase 
+% and y-pol by negative phase 
+%el_rot=0*psi_el;
+%Ex=cos(el_rot)*Ex-sin(el_rot)*Ey;
+%Ey=sin(el_rot)*Ex+cos(el_rot)*Ey;
+
+Ipos=(abs(Ey).^2)/2;
+Ineg=(abs(Ex).^2)/2;
+
+figure(1);
+hold off;
+plot(Ep.^2, real(xi_left-xi_right), '-');  
+title("differential real xi");
+xlabel("two photon detuning");
+
+figure(2);
+hold off;
+plot(Ep.^2, imag(xi_left-xi_right), '-');  
+title("differential imag xi");
+xlabel("two photon detuning");
+
+figure(3);
+hold off;
+plot(Ep.^2, imag(xi_left), '-', Ep, imag(xi_right), '-');  
+title("imag xi");
+xlabel("two photon detuning");
+
+figure(4);
+hold off;
+plot(Ep.^2, real(xi_left), '-', Ep.^2, real(xi_right), '-');  
+title("real xi");
+xlabel("two photon detuning");
+
+figure(5);
+hold off;
+%plot(Ep.^2, (Ipos-Ineg), '-');  
+semilogx(Ep.^2, (Ipos-Ineg), '-');  
+%semilogx(Ep.^2, (Ipos-Ineg)./(Ep.^2), '-');  
+title("BPD signal xi");
+xlabel("two photon detuning");
+
+%figure(1); 
+	%hold off;
+	%plot(detuning_freq, imag(xi_linear), '-1;linear;');  
+	%hold on;
+	%plot(detuning_freq, imag(xi_left),   '-2;left;');  
+	%plot(detuning_freq, imag(xi_right),  '-3;right;');  
+	%title("probe absorption");
+	%hold off;
+%figure(2); 
+	%hold off;
+	%plot(detuning_freq, real(xi_linear), '-1;linear;');  
+	%hold on;
+	%plot(detuning_freq, real(xi_left),   '-2;left;');  
+	%plot(detuning_freq, real(xi_right),  '-3;right;');  
+	%title("probe dispersion");
+	%hold off;
+