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

clear all; close all; clc;
date=20170524;
load data2405.mat ;
cut=(key_value(:,8)>1e+3)&(abs(data(:,11,2))<100)&(abs(data(:,12,2))<50);

position=key_value(:,4)*7.94e-6;
hour=rem(key_value(:,6),date);
diap(1,:)=(hour>143200)&(hour<144000);
diap(2,:)=(hour>144000)&(hour<145000);
diap(3,:)=(hour>145000)&(hour<145600);
diap(4,:)=(hour>145600)&(hour<150300);
diap(5,:)=(hour>150300)&(hour<151000);
diap(6,:)=(hour>151000)&(hour<151700);
diap(7,:)=(hour>151700)&(hour<152400);


fromd=-14;
bkgsh=1;
for detector=1:11
for dd=1:7
cut2=cut&diap(dd,:)';
x=position(cut2)-48.58-3.46;
tmpE=key_value(cut2,8);
Es(dd)=mean(tmpE(abs(tmpE-mean(tmpE))<std(tmpE)));
Ess(dd)=std(tmpE(abs(tmpE-mean(tmpE))<std(tmpE)));

y=data(cut2,detector,2)./Es(dd);

fb(dd,detector) = polyfit(x(x<(fromd+bkgsh)),y(x<(fromd+bkgsh)),0);

%[fs, gofs] = fit(x(x>fromd),y(x>fromd),'exp1');
[fs, gofs] = fit(x((x>fromd)&y>fb(dd,detector)),y((x>fromd)&y>fb(dd,detector)),'exp1','Lower',[0 0.015]);
chi(dd,detector) = gofs.adjrsquare;
lamb(dd,detector)=fs.b;
amp(dd,detector)=fs.a;
% figure
% plot(fs,x,y)
tmp1=confint(fs);
tmp2=(tmp1(1:2:3)-tmp1(2:2:4))/4;
Ierrp(dd,detector,:)=tmp2(:);

tmp1=confint(fs);
tmp2=(tmp1(1:2:3)-tmp1(2:2:4))/4;
Ierrp(dd,detector,:)=tmp2(:);


end 
end

Bphase=[2.4 2.2 2.0 1.8 1.6 1.4 1.2];
dd=3;
cut2=cut&diap(dd,:)';
x=position(cut2)-48.58-3.46;
xfits=fromd:1:(-7);
xfitb=(-22):1:fromd;

figure 
hold on
detector=7;
y=data(cut2,detector,2)./mean(key_value(cut2,8))*1e+3;
semilogy(x,y,'^r','linewidth',2)
semilogy(xfits,amp(dd,detector)*exp(xfits*lamb(dd,detector))*1e+3,'r','linewidth',2)
semilogy(xfitb,repmat(fb(dd,detector),size(xfitb))*1e+3,'Color',[0.8 0.8 0.8],'linewidth',2)
set(gca, 'YScale', 'log')
detector=5;
y=data(cut2,detector,2)./mean(key_value(cut2,8))*1e+3;
semilogy(x,y,'>k','linewidth',2)
semilogy(xfits,amp(dd,detector)*exp(xfits*lamb(dd,detector))*1e+3,'k','linewidth',2)
semilogy(xfitb,repmat(fb(dd,detector),size(xfitb))*1e+3,'Color',[0.8 0.8 0.8],'linewidth',2)

detector=9;
y=data(cut2,detector,2)./mean(key_value(cut2,8))*1e+3;
semilogy(x,y,'<b','linewidth',2)
semilogy(xfits,amp(dd,detector)*exp(xfits*lamb(dd,detector))*1e+3,'b','linewidth',2)
semilogy(xfitb,repmat(fb(dd,detector),size(xfitb))*1e+3,'Color',[0.8 0.8 0.8],'linewidth',2)

set(gca,'fontsize',16)
legend('76^o','Signal','Bkgnd','90^o','Signal','Bkgnd', '104^o','Signal','Bkgnd')
xlabel('Position, [mm]')
ylabel('Amplitude, [a.u.]')
grid on 
title('Buncher power = 2 [a.u.]')
xlim([-22 0])
% ylim([6 21])



figure
errorbar(Bphase,Es/1000,Ess/1000,'-*r','linewidth',2)
set(gca,'fontsize',16)
xlabel('Buncher power, [a.u.]')
ylabel('Electron intensity, [a.u.]')
grid on 




figure
chi(chi<0)=0;
[hc hc]=contourf(Bphase(:), 48:7:118, chi',40)
set(hc,'LineStyle','none');
colorbar
set(gca,'fontsize',16)
 ylabel('Angle, [deg.]')
xlabel('Buncher power, [a.u.]')
title('R^2')

figure
[hc hc]=contourf(Bphase(:), 48:7:118, 2./lamb',40)
set(hc,'LineStyle','none');
colorbar
set(gca,'fontsize',16)
 ylabel('Angle, [deg.]')
xlabel('Buncher power, [a.u.]')
title('Evanescent wavelength, [mm]')

%3d diap is bad

figure
hold on 
errorbar(Bphase',2./lamb(:,4),2./lamb(:,4).^2.*Ierrp(:,4,2),'-.^b','linewidth',2)
errorbar(Bphase',2./lamb(:,7),2./lamb(:,7).^2.*Ierrp(:,7,2),'-.*r','linewidth',2)
errorbar(Bphase',2./lamb(:,10),2./lamb(:,10).^2.*Ierrp(:,10,2),'-.ok','linewidth',2)
set(gca, 'YScale', 'log')
set(gca,'fontsize',16)
xlabel('Buncher power, [a.u.]')
ylabel('Evanescent wavelength, [mm]')
grid on 
gam=((38:1:44)+0.2)/0.511+1;
bet=sqrt(1-1./gam.^2);
phi=2/180*pi;
plot(Bphase',6*(1./bet-cos((69)/180*pi)).*bet.*gam/(2*pi)./(sqrt(1+(bet.*gam*sin((69)/180*pi)*sin(phi)).^2)),'-b','linewidth',2)
plot(Bphase',6*(1./bet-cos((90)/180*pi)).*bet.*gam/(2*pi)./(sqrt(1+(bet.*gam*sin((90)/180*pi)*sin(phi)).^2)),'-r','linewidth',2)
plot(Bphase',6*(1./bet-cos((111)/180*pi)).*bet.*gam/(2*pi)./(sqrt(1+(bet.*gam*sin((111)/180*pi)*sin(phi)).^2)),'-k','linewidth',2)
legend('D. 69^o','D. 90^o', 'D. 111^o','C. 69^o','C. 90^o', 'C. 111^o')

figure
hold on
plot(Bphase',amp(:,4),'-^b','linewidth',2)
plot(Bphase',amp(:,7),'-*r','linewidth',2)
plot(Bphase',amp(:,10),'-ok','linewidth',2)

set(gca,'fontsize',16)
legend('69^o','90^o', '111^o')
xlabel('Buncher power, [a.u.]')
ylabel('Amplitude, [a.u.]')
grid on 



figure
[hc hc]=contourf(Bphase(:), 48:7:118, fb',40)
  set(hc,'LineStyle','none');
colorbar
set(gca,'fontsize',16)
 ylabel('Angle, [deg.]')
xlabel('Buncher power, [a.u.]')
title('Background, [a.u.]')
%  (hf,'-depsc','file.eps')

figure
[hc hc]=contourf(Bphase(:), 48:7:118, amp',40)
  set(hc,'LineStyle','none');
colorbar
set(gca,'fontsize',16)
 ylabel('Angle, [deg.]')
xlabel('Buncher power, [a.u.]')
title('Amplitude, [a.u.]')
%  (hf,'-depsc','file.eps')




figure
hold on
 
 
gam=44.2/0.511+1;
bet=sqrt(1-1/gam^2);
phi=2/180*pi;
plot(48:7:118,6*(1/bet-cos((48:7:118)/180*pi))*bet*gam/(2*pi)./(sqrt(1+(bet*gam*sin((48:7:118)/180*pi)*sin(phi)).^2)),'m','linewidth',2)
phi=1/180*pi;
plot(48:7:118,6*(1/bet-cos((48:7:118)/180*pi))*bet*gam/(2*pi)./(sqrt(1+(bet*gam*sin((48:7:118)/180*pi)*sin(phi)).^2)),'g','linewidth',2)
phi=0/180*pi;
plot(48:7:118,6*(1/bet-cos((48:7:118)/180*pi))*bet*gam/(2*pi)./(sqrt(1+(bet*gam*sin((48:7:118)/180*pi)*sin(phi)).^2)),'b','linewidth',2)
  
W=1./(Ierrp(:,:,2)).^2;
chilevel=0;
for c=1:11
    LL(c)=sum(lamb((chi(:,c)>chilevel),c).*W((chi(:,c)>chilevel),c))/sum(W((chi(:,c)>chilevel),c));
    LE(c)=sqrt(1/sum(W((chi(:,c)>chilevel),c)));
end
errorbar(48:7:118,2*(LL.^-1),2./(LL.^2).*LE,'-..k','linewidth',2)

 grid on
 set(gca,'fontsize',16)
 xlabel('Angle, [deg.]')
 ylabel('Evanescent wavelength, [mm]')
 legend('\phi=2^o','\phi=1^o','\phi=0^o','data','Location','northwest')
 separation=-10;

figure
hold on

bp=1;
Spectrum=amp(bp,:).*exp(lamb(bp,:)*separation);
errorbar(48:7:118,Spectrum,Spectrum.*sqrt((Ierrp(bp,:,1)./amp(bp,:)).^2+(separation*Ierrp(bp,:,2)).^2),'k','linewidth',2)
bp=3;
Spectrum=amp(bp,:).*exp(lamb(bp,:)*separation);
errorbar(48:7:118,Spectrum,Spectrum.*sqrt((Ierrp(bp,:,1)./amp(bp,:)).^2+(separation*Ierrp(bp,:,2)).^2),'r','linewidth',2)
 bp=5;
Spectrum=amp(bp,:).*exp(lamb(bp,:)*separation);
errorbar(48:7:118,Spectrum,Spectrum.*sqrt((Ierrp(bp,:,1)./amp(bp,:)).^2+(separation*Ierrp(bp,:,2)).^2),'b','linewidth',2)
  
%  bp=7;
% Spectrum=amp(bp,:).*exp(lamb(bp,:)*separation);
% errorbar(48:7:118,Spectrum,Spectrum.*sqrt((Ierrp(bp,:,1)./amp(bp,:)).^2+(separation*Ierrp(bp,:,2)).^2),'g','linewidth',2)
%   
set(gca,'fontsize',16)
legend('BP=2.4','BP=2.0','BP=1.6')
xlabel('Angle, [deg.]')
ylabel('Amplitude, [a.u.]')
title('Spectrum')
grid on 
ylim([0 0.02])