source: ETALON/CLIO/Analysis/Data24May/Analyse2bis.m @ 719

clear all; close all; clc;
load analyze.mat

idet=7;
figure
hold on
plot(amp(:,idet),amp(:,1),'*b','linewidth',2)
plot(amp(:,idet),amp(:,2),'^g','linewidth',2)
plot(amp(:,idet),amp(:,3),'*g','linewidth',2)
plot(amp(:,idet),amp(:,4),'^b','linewidth',2)
plot(amp(:,idet),amp(:,5),'ob','linewidth',2)
 plot(amp(:,idet),amp(:,6),'^m','linewidth',2)
 plot(amp(:,idet),amp(:,7),'-or','linewidth',2)
 plot(amp(:,idet),amp(:,8),'*r','linewidth',2)
 plot(amp(:,idet),amp(:,9),'*k','linewidth',2)
 plot(amp(:,idet),amp(:,10),'ok','linewidth',2)
xlim([0.0 0.055])
ylim([0 0.15])
set(gca,'fontsize',16)
legend('48^o','55^o', '62^o','69^o','76^o', '83^o','90^o','97^o','104^o', '111^o')
xlabel('Amplite 90^o, [a.u.]')
ylabel('Amplitude, [a.u.]')
title('Amplitude correlation')
grid on 

cut=amp(:,idet)>0.025;
for c=1:11
fx=fit(amp(cut,idet),amp(cut,c),'poly1');
f(c,:)=[fx.p1 fx.p2];
tmp1=confint(fx);
tmp2=(tmp1(1:2:3)-tmp1(2:2:4))/4;
Er(c,:)=tmp2;
end

figure
errorbar(48:7:118,f(:,1),Er(:,1),'-.ok','linewidth',2)
set(gca,'fontsize',16)
xlabel('Angle, [deg.]')
ylabel('Amplitude, [a.u.]')
title('Slope')
grid on 

figure
errorbar(48:7:118,f(:,2),Er(:,2),'-.or','linewidth',2)
set(gca,'fontsize',16)
xlabel('Angle, [deg.]')
ylabel('Amplitude, [a.u.]')
title('Base')
grid on 

cut=amp(:,7)>0.025;

for idet=1:11
for c=1:11
fx=fit(amp(cut,idet),amp(cut,c),'poly1');
f2d(c,idet,:)=[fx.p1 fx.p2];
Cx(c,idet)=corr(amp(cut,idet),amp(cut,c));
tmp1=confint(fx);
tmp2=(tmp1(1:2:3)-tmp1(2:2:4))/4;
Er2d(c,idet,:)=tmp2;
end
end 
figure
contourf(48:7:118,48:7:118,f2d(:,:,1)')
set(gca,'fontsize',16)
xlabel('Angle, [deg.]')
ylabel('Angle ref., [deg.]')
title('Slope')
grid on 
colorbar
figure
contourf(48:7:118,48:7:118,f2d(:,:,2)')
colorbar
set(gca,'fontsize',16)
xlabel('Angle, [deg.]')
ylabel('Angle ref., [deg.]')
title('Base')
grid on 
colorbar

% figure
% contourf(48:7:118,48:7:118,Er2d(:,:,1)')
% colorbar
% figure
% contourf(48:7:118,48:7:118,Er2d(:,:,2)')
% colorbar
figure
contourf(48:7:118,48:7:118,Cx(:,:)')
colorbar
set(gca,'fontsize',16)
xlabel('Angle, [deg.]')
ylabel('Angle ref., [deg.]')
title('Correlation')
grid on 
colorbar