source: ETALON/CLIO/Analysis/ProfileX.m @ 715

close all; clear all; clc
load FFX.mat
cd('./KK')

for c=[2:10]%, 14:18
    c
% for c=1:20
SP=flip(Freq(3:end));
SS(c,:)=flip(FF(c,3:end))/9;
% SSM(c)=max(SS);
[RSP,ISP(c,:)]=SRf(SP,SS(c,:));
% 

% figure (1)
% plot(SP,SS(c,:),'or')
% hold on
%  plot(RSP)

 PhaseH=-(imag(hilbert(log(abs(RSP)))));
RPH=real((ifft(RSP.*exp(1i.*(PhaseH)))));


% figure(2)
[mx,ps]=max(RPH)
RPH= [zeros(1,80-ps), RPH(1:end-(80-ps))]
% plot((1:length(RPH))*1e-9/length(RPH),RPH)

Gauss=@(mu,sig,A,x)A.*exp(-(x-mu).^2/(2*sig^2))
% hold on
%2.1e-12
Profile5=Gauss(ps*1e-9/length(RPH),2.1e-12,mx,(1:length(RPH))*1e-9/length(RPH));
% plot((1:length(RPH))*1e-9/length(RPH),Profile5,'r')
 FWHM(c-1)=sum(RPH>max(RPH)/2)*1e-9/length(RPH)
  FW01M(c-1)=sum(RPH>max(RPH)*0.1)*1e-9/length(RPH)
 FW09M(c-1)=sum(RPH>max(RPH)*0.9)*1e-9/length(RPH)

 
PR(c-1,:)=RPH;
sum(Profile5>max(Profile5)/2)*1e-9/length(Profile5)
% figure(1)
% plot(abs(fft(RPH)),'r')
% plot(abs(fft(Profile5./sum(Profile5))),'k')

 end
Bphase=sort([8.14, 8.09, 8.04, 7.94, 7.84, 7.74, 8.19, 8.24, 8.34, 8.44, 8.54, 8.64, 8.74, 8.84, 8.94, 9.04, 9.14, 9.24]);

% ds=1:20;
 figure
 hold on
plot(Bphase(2:10),FW01M(:)*1e+12,'->k','linewidth',2)
plot(Bphase(2:10),FWHM(:)*1e+12,'->b','linewidth',2)
plot(Bphase(2:10),FW09M(:)*1e+12,'->r','linewidth',2)

 set(gca,'fontsize',16)
  xlabel('\phi_B, [a.u.]')
 ylabel('Pulsewidth, [ps]')
 legend('FW0.1M','FWHM','FW0.9M','Orientation','horizontal')
 grid on
 title('Pulsewidth vs. Buncher phase')
 
figure
[hc hc]=contourf((1:length(RPH))*1e-9/length(RPH)*1e+12,Bphase(2:10),PR*500)
colorbar
  set(hc,'LineStyle','none');
xlim([0 50])
 set(gca,'fontsize',16)

 xlabel('Time, [ps]')
 ylabel('\phi_B, [a.u.]')
 title('Bunch profile')
 
 
 figure
 
 hold on
plot(SP,  SS(2,:),'ob','linewidth',2)
plot(0:(length(ISP)-1),  ISP(2,:),'-.b','linewidth',2)
plot(SP,  SS(4,:),'om','linewidth',2)
plot(0:(length(ISP)-1),  ISP(4,:),'-.m','linewidth',2)
plot(SP,  SS(6,:),'ok','linewidth',2)
plot(0:(length(ISP)-1),  ISP(6,:),'-.k','linewidth',2)
plot(SP,  SS(8,:),'og','linewidth',2)
plot(0:(length(ISP)-1),  ISP(8,:),'-.g','linewidth',2)
plot(SP,  SS(10,:),'or','linewidth',2)
plot(0:(length(ISP)-1),  ISP(10,:),'-.r','linewidth',2)

 set(gca,'fontsize',16)
 xlabel('Frequency, [GHz]')
 legend('\phi_B=7.84','','\phi_B=8.04','','\phi_B=8.14','','\phi_B=8.24','','\phi_B=8.44','')
 % 2 4 6 8 10
 title('Form Factor (BGS=10 mm)')
ylabel('Amplitude, [a.u.]')
xlim([0, 200])
grid on
 
  figure
 
 hold on
plot((0:(length(RPH)-1))*1e-9/length(RPH)*1e+12,  PR(1,:)*500,'-.b','linewidth',2)
plot((0:(length(RPH)-1))*1e-9/length(RPH)*1e+12,  PR(3,:)*500,'-.m','linewidth',2)
plot((0:(length(RPH)-1))*1e-9/length(RPH)*1e+12,  PR(5,:)*500,'-.k','linewidth',2)
plot((0:(length(RPH)-1))*1e-9/length(RPH)*1e+12,  PR(6,:)*500,'-.g','linewidth',2)
plot((0:(length(RPH)-1))*1e-9/length(RPH)*1e+12,  PR(9,:)*500,'-.r','linewidth',2)

 set(gca,'fontsize',16)
 xlabel('Time, [ps]')
 legend('\phi_B=7.84','\phi_B=8.04','\phi_B=8.14','\phi_B=8.24','\phi_B=8.44')
 % 2 4 6 8 10
 title('Reconstructed profile')
ylabel('Amplitude, [a.u.]')
xlim([0, 50])
grid on
 
  
  figure
 
 hold on
plot((0:(length(RPH)-1))*1e-9/length(RPH)*1e+12,  PR(5,:)*500,'-.b','linewidth',2)
plot((0:(length(RPH)-1))*1e-9/length(RPH)*1e+12,  PR(6,:)*500,'-.m','linewidth',2)
plot((0:(length(RPH)-1))*1e-9/length(RPH)*1e+12,  PR(7,:)*500,'-.k','linewidth',2)
plot((0:(length(RPH)-1))*1e-9/length(RPH)*1e+12,  PR(8,:)*500,'-.g','linewidth',2)
plot((0:(length(RPH)-1))*1e-9/length(RPH)*1e+12,  PR(9,:)*500,'-.r','linewidth',2)

 set(gca,'fontsize',16)
 xlabel('Time, [ps]')
 legend('\phi_B=8.14','\phi_B=8.19','\phi_B=8.24','\phi_B=8.34','\phi_B=8.44')
 % 2 4 6 8 10
 title('Reconstructed profile')
ylabel('Amplitude, [a.u.]')
xlim([0, 50])
grid on
 
 
 time =(0:(length(RPH)-1))*1e-9/length(RPH)*1e+9;%ns
 Profile=PR(6,:)*500;
 
 save P6.mat time Profile
 
 
 
 
 
 
 
 
 
 
 
 
 
 cd('..')