source: ETALON/reconstruction/c123/together.m @ 751

close all
clear all

Np=3

[Profile,~]=PR_func(Np, 1,1);
aFT=abs(fft(Profile));
samp=[25.6459,   16.7967,   12.0000,    9.1188,    7.2608,    6.0000,    5.1123,    4.4709,    4.0000,    3.6523,  0.6795,...
              3.3974,    5.1292,    3.3593,    2.4000,    1.8238,    1.4522,    1.2000,    1.0225,    0.8942,    0.8000,    0.7305,...
              0.8549,    0.5599,    0.4000,    0.3040,    0.2420,    0.2000,    0.1704,    0.1490,    0.1333,    0.1217,    0.1132];
SFr=sort(samp);
SP=unique(round(2100.*SFr./max(25.6459)));%ACHTUNG!!! remove same points
SS=aFT(SP);

[RSP,ISP]=SRf(SP,SS);
rho=ISP(1:2000)';
angl2=(imag(hilbert(log(abs([rho; flip(rho)])))));
% angl2=-KK_phase(rho')';
angl=angl2(1:end/2);
t0=-0.2;
step=0.005;
profH=repmat(rho,[1,1000]).*cos(2*pi*repmat((1:2000)'*(25.6459/2100),[1,1000]).*repmat((t0+step*(0:999)),[length(rho),1])-repmat(angl,[1,1000]));
prof2=2*sum(profH);
lm=prof2>=0;
% plot((t0+step*(0:999)),prof2./sum(prof2(lm)),'linew',2)
plot((t0+step*(0:999)),prof2./max(prof2(lm)),'linew',2)

fileID = fopen(['../data/profile' num2str(Np) '.txt'],'r');
A=fscanf(fileID,'%f %f');
fclose(fileID);
prof=A(2:2:end);
hold on
lm=prof>=0;
plot((t0+step*(0:999)),prof./max(prof(lm)),'r','linew',2)
set(gca,'fontsize',14)
title('Profile')
xlabel('Time, [ps]')
legend('matlab','cauchy2')
grid on
ylabel('Intensity')
RPH=fliplr(real((ifft([rho; flip(rho)].*exp(1i.*(-angl2))))));
RPHIn=RPH(1:1000);
lm=RPHIn>=0;
% plot((t0+step*(0:999)),circshift(RPHIn./sum(RPHIn(lm)),[1 40]),'g','linew',2)
%   plot((t0+step*(0:999)),circshift(RPHIn./max(RPHIn(lm)),[1 1]),'g','linew',2)
%  plot(t0+step*(0:999),RPHIn,'g','linew',2)
% figure(123)
% plot(RPH)