added files to calculate and plot psr combo

1 files changed, 37 insertions, 3 deletions

diff --git a/psr/psr_vs_detuning_combo.m b/psr/psr_vs_detuning_combo.m
index 2e8921d..1016a82 100644
--- a/psr/psr_vs_detuning_combo.m
+++ b/psr/psr_vs_detuning_combo.m
@@ -13,12 +13,46 @@ B_field=zeeman_splitting/gmg;
 %[psr_rad]=psr_vs_detuning(Ep, psi_el, B_field, theta, phi) 
 
 % phi is angle between linear polarization and axis x
-phi=pi*2/8;
+phi=pi/4;
 % theta is angle between lab z axis (light propagation direction) and magnetic field axis (z')
 theta=0;
 % psi_el is the ellipticity parameter (phase difference between left and right polarization)
 psi_el=-30/180*pi;
 
-Ep=0.2;
-[psr_rad]=psr_vs_detuning(detuning_freq, Ep, psi_el, B_field, theta, phi) ;
+
+figure(6);
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%5 
+% zero magnetic field,,  30 degree ellipticity
+zeeman_splitting=+0.000;
+B_field=zeeman_splitting/gmg;
+psi_el=30/180*pi;
+
+[psr_rad_tnEp_pos_el, psr_rad_tnEp_neg_el, psr_rad_smEp_pos_el, psr_rad_smEp_neg_el, psr_rad_lgEp_pos_el, psr_rad_lgEp_neg_el, psr_rad_grEp_pos_el, psr_rad_grEp_neg_el] =make_representative_psr_vs_detuning_for_given_B_and_psi_el(detuning_freq, B_field, psi_el, theta, phi);
+
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%5 
+% 0.1 G magnetic field,,  30 degree ellipticity
+
+zeeman_splitting=+0.070;
+B_field=zeeman_splitting/gmg;
+psi_el=30/180*pi;
+
+[psr_rad_tnEp_pos_el, psr_rad_tnEp_neg_el, psr_rad_smEp_pos_el, psr_rad_smEp_neg_el, psr_rad_lgEp_pos_el, psr_rad_lgEp_neg_el, psr_rad_grEp_pos_el, psr_rad_grEp_neg_el] =make_representative_psr_vs_detuning_for_given_B_and_psi_el(detuning_freq, B_field, psi_el, theta, phi);
+
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%5 
+% 0.0001 G magnetic field,,  30 degree ellipticity
+
+zeeman_splitting=+0.000070;
+B_field=zeeman_splitting/gmg;
+psi_el=30/180*pi;
+
+[psr_rad_tnEp_pos_el, psr_rad_tnEp_neg_el, psr_rad_smEp_pos_el, psr_rad_smEp_neg_el, psr_rad_lgEp_pos_el, psr_rad_lgEp_neg_el, psr_rad_grEp_pos_el, psr_rad_grEp_neg_el] =make_representative_psr_vs_detuning_for_given_B_and_psi_el(detuning_freq, B_field, psi_el, theta, phi);
+	
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%5 
+% 1.0 G magnetic field,,  30 degree ellipticity
+
+zeeman_splitting=+0.70;
+B_field=zeeman_splitting/gmg;
+psi_el=30/180*pi;
+
+[psr_rad_tnEp_pos_el, psr_rad_tnEp_neg_el, psr_rad_smEp_pos_el, psr_rad_smEp_neg_el, psr_rad_lgEp_pos_el, psr_rad_lgEp_neg_el, psr_rad_grEp_pos_el, psr_rad_grEp_neg_el] =make_representative_psr_vs_detuning_for_given_B_and_psi_el(detuning_freq, B_field, psi_el, theta, phi);