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

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

for c=1:6%, 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)=sum(RPH>max(RPH)/2)*1e-9/length(RPH)
  FW01M(c)=sum(RPH>max(RPH)*0.1)*1e-9/length(RPH)
 FW09M(c)=sum(RPH>max(RPH)*0.9)*1e-9/length(RPH)

 
PR(c,:)=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=[2.4 2.2 2.0 1.8 1.6 1.4 1.2];

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

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

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

 set(gca,'fontsize',16)
 xlabel('Frequency, [GHz]')
 legend('BP=2.4','','BP=2.2','','BP=2.0','','BP=1.8','','BP=1.6','')
 % 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(2,:)*500,'-.m','linewidth',2)
plot((0:(length(RPH)-1))*1e-9/length(RPH)*1e+12,  PR(3,:)*500,'-.k','linewidth',2)
plot((0:(length(RPH)-1))*1e-9/length(RPH)*1e+12,  PR(4,:)*500,'-.g','linewidth',2)
plot((0:(length(RPH)-1))*1e-9/length(RPH)*1e+12,  PR(5,:)*500,'-.r','linewidth',2)
plot((0:(length(RPH)-1))*1e-9/length(RPH)*1e+12,  PR(6,:)*500,'-.y','linewidth',2)

 set(gca,'fontsize',16)
 xlabel('Time, [ps]')
 legend('BP=2.4','BP=2.2','BP=2.0','BP=1.8','BP=1.6','BP=1.4')
 % 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(3,:)*500;
 
 save P6BP.mat time Profile
 
 
 
 
 
 
 
 
 
 
 
 
 
 cd('..')