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