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

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

load('GFW/SEY.mat')
FF=Spectrum(3:11)./Integral(3:11)*1e-24;
eFF=Err(3:11)./Integral(3:11)*1e-24;

Freq=300./(6*(1-cos((62:7:118)/180*pi)));
% figure
% hold on
% errorbar(Freq, FF/9,eFF/9)


figure

cd('./KK')
for c=1:100
    
FF2=6*eFF.*rand(size(eFF))+FF-3*eFF;    
[RSP,ISP]=SRf(flip(Freq),flip(FF2)/9);
 PhaseH=-(imag(hilbert(log(abs(RSP)))));
RPH=real((ifft(RSP.*exp(1i.*(PhaseH)))));
[mx,ps]=max(RPH);

RPH= [zeros(1,150-ps), RPH(1:end-(150-ps))];


% plot(RSP,'-r')
hold on
plot((1:length(RPH))*1e-9/length(RPH)*1e+12,RPH,'Color',[0.8 0.8 0.8])
RP(c,:)=RPH(:);
 FWHM(c)=sum(RPH>max(RPH)/2)*1e-9/length(RPH)

end
plot((1:length(RPH))*1e-9/length(RPH)*1e+12,mean(RP),'k','linew',2)
 set(gca,'fontsize',16)
 xlabel('Time, [ps]')
 title('Reconstructed profile')
ylabel('Amplitude, [a.u.]')
xlim([0, 50])
grid on

cd('..')
mean(FWHM)
std(FWHM)