2 files changed, 36 insertions, 26 deletions

diff --git a/faraday_and_psr/builder_tasks.m b/faraday_and_psr/builder_tasks.m
index 78eea9e..3d81057 100755
--- a/faraday_and_psr/builder_tasks.m
+++ b/faraday_and_psr/builder_tasks.m
@@ -14,14 +14,13 @@ gmg=.7; % gyro magnetic ration for ground level
 %[psr_rad]=psr_vs_detuning(Ep, psi_el, B_field, theta, phi) 
 
 % phi is angle between linear polarization and axis x
-%phi=pi/4;
 phi=0;
 % 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=0*3/180*pi;
 
-Ep=sqrt(10.1);
+Ep=sqrt(100);
 
 % we are going to sweep B so parts related to sweep paramer(s)
 zeeman_splitting=+0.5;
diff --git a/faraday_and_psr/pp_tasks.m b/faraday_and_psr/pp_tasks.m
index 82bcdac..a3919ba 100644
--- a/faraday_and_psr/pp_tasks.m
+++ b/faraday_and_psr/pp_tasks.m
@@ -29,37 +29,48 @@ Ep_out.y= 1i*(Ep_out.left - Ep_out.right)/sqrt(2);
 
 xi_x=(xi_right+xi_left)/sqrt(2);
 xi_y=(1i*xi_right-1i*xi_left)/sqrt(2);
-figure(1);
-plot(B_field, real(xi_x), B_field, real(xi_y));
-legend('x','y');
-title('xi-real');
 
-figure(2);
-plot(B_field, imag(xi_x), B_field, imag(xi_y));
-legend('x','y');
-title('xi-imag');
+%figure(1);
+%plot(B_field, real(xi_x), B_field, real(xi_y));
+%legend('x','y');
+%title('xi-real');
 
-figure(3);
-plot(B_field, real(xi_left), B_field, real(xi_right));
-legend('left','right');
-title('xi-real');
+%figure(2);
+%plot(B_field, imag(xi_x), B_field, imag(xi_y));
+%legend('x','y');
+%title('xi-imag');
 
-figure(4);
-plot(B_field, imag(xi_left), B_field, imag(xi_right));
-legend('left','right');
-title('xi-imag');
+%figure(3);
+%plot(B_field, real(xi_left), B_field, real(xi_right));
+%legend('left','right');
+%title('xi-real');
 
-figure(5);
-plot(B_field, abs(Ep_out.x).^2, B_field, abs(Ep_out.y).^2 );
-legend('x','y');
-title('Signal out');
+%figure(4);
+%plot(B_field, imag(xi_left), B_field, imag(xi_right));
+%legend('left','right');
+%title('xi-imag');
 
-figure(6);
-plot(B_field, abs(Ep_out.left).^2, B_field, abs(Ep_out.right).^2 );
-legend('left','right');
-title('Signal out');
+%figure(5);
+%plot(B_field, abs(Ep_out.x).^2, B_field, abs(Ep_out.y).^2 );
+%legend('x','y');
+%title('Signal out');
+
+%figure(6);
+%plot(B_field, abs(Ep_out.left).^2, B_field, abs(Ep_out.right).^2 );
+%legend('left','right');
+%title('Signal out');
 
 figure(7);
+plot(B_field, abs(Ep_out.x).^2 );
+legend('x');
+title('Signal out');
+
+figure(8);
 plot(B_field, abs(Ep_out.y).^2 );
 legend('y');
 title('Signal out');
+
+figure(9);
+plot(B_field, abs(Ep_out.x).^2-abs(Ep_out.y).^2 );
+legend('x-y');
+title('Differential x-y out');