small improvements

1 files changed, 20 insertions, 28 deletions

diff --git a/psr/output_psr_results_vs_detuning.m b/psr/output_psr_results_vs_detuning.m
index c6120fc..cf71711 100644
--- a/psr/output_psr_results_vs_detuning.m
+++ b/psr/output_psr_results_vs_detuning.m
@@ -23,46 +23,38 @@ figure(1);
 hold off;
 plot(detuning_freq, real(xi_left-xi_right), '-');  
 title("differential real xi");
-xlabel("two photon detuning");
+xlabel("two photon detuning (MHz)");
 
 figure(2);
 hold off;
 plot(detuning_freq, imag(xi_left-xi_right), '-');  
 title("differential imag xi");
-xlabel("two photon detuning");
+xlabel("two photon detuning (MHz)");
 
 figure(3);
 hold off;
-plot(detuning_freq, imag(xi_left), '-', detuning_freq, imag(xi_right), '-');  
-title("imag xi");
-xlabel("two photon detuning");
+plot(detuning_freq, real(xi_left), '-', detuning_freq, real(xi_right), '-');  
+title("real xi");
+xlabel("two photon detuning (MHz)");
 
 figure(4);
 hold off;
-plot(detuning_freq, real(xi_left), '-', detuning_freq, real(xi_right), '-');  
-title("real xi");
-xlabel("two photon detuning");
+plot(detuning_freq, imag(xi_left), '-', detuning_freq, imag(xi_right), '-');  
+title("imag xi");
+xlabel("two photon detuning (MHz)");
 
 figure(5);
 hold off;
-plot(detuning_freq, (Ipos-Ineg), '-');  
-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;
+I_probe=E_field_probe^2;
+psr_rad=(Ipos-Ineg)/(2*I_probe);
+plot(detuning_freq, psr_rad, '-');  
+subt_str=sprintf("Laser Rabi freq normalized to upper state decay %.3f, ellipticity %.1f degree, \n B field ground level splitting %.3f MHz", I_probe, psi_el*180/pi, B_field); 
+title(cstrcat("BPD normilized PSR signal at F_g=2 to F_e=1,2.\n ",subt_str) );
+xlabel("two photon detuning (MHz)");
+ylabel("PSR (radians)");
+
+print("psr_vs_detuning.ps");
+
+fname= sprintf("psr_vs_detuning_Fg=2toFe=1,2_Ip=%.3f_el_%.1f_B=%.3fMHz.mat", I_probe, psi_el*180/pi,B_field); 
+save(fname,'detuning_freq', 'psr_rad');