basis_transformation;

fnames=glob('results/*.mat');

Nsteps=length(fnames)
B_field=zeros(1,Nsteps);
xi_linear=zeros(1,Nsteps);
xi_left=zeros(1,Nsteps);
xi_right=zeros(1,Nsteps);

% read the information from resulting files
for i=1:Nsteps
	d=load(fnames{i});
	B_field(i)=d.B_field;
	xi_linear(i)=d.xi_linear;
	xi_left(i)=d.xi_left;
	xi_right(i)=d.xi_right;
	Ep.linear(i)=d.E_field_pos_freq.linear;
	Ep.left(i)=d.E_field_pos_freq.left;
	Ep.right(i)=d.E_field_pos_freq.right;
endfor

Ep_out.linear=(1-xi_linear).*Ep.linear;
Ep_out.left=(1-xi_left).*Ep.left;
Ep_out.right=(1-xi_right).*Ep.right;

Ep_out.x= (Ep_out.left + Ep_out.right)/sqrt(2);
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(3);
plot(B_field, real(xi_left), B_field, real(xi_right));
legend('left','right');
title('xi-real');

figure(4);
plot(B_field, imag(xi_left), B_field, imag(xi_right));
legend('left','right');
title('xi-imag');

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.y).^2 );
legend('y');
title('Signal out');