diff options
| -rw-r--r-- | psr/output_psr_results_vs_detuning.m | 48 |
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'); |