1;

N_detun_steps=1000;
detuning_p_min=-200.0;
%detuning_p_max=-detuning_p_min;
detuning_p_max=1000;
detuning_freq=linspace(detuning_p_min,detuning_p_max,N_detun_steps);

gmg=.7; % gyro magnetic ration for ground level
zeeman_splitting=+0.000;
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/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;


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);