PEAK_F=zeros(1,21);
PEAK_T=zeros(1,21);
cells=zeros(1,21);

for i=1:1:21

PATH=(i-1)*5;

%%
load(sprintf('Net_burst_extKo_IN_path_%d.mat',PATH));    
    
T=5000;         % first duration
t_old=1;
dT=T;           % shift every second dT
a=0;

peak_T=10;          % threshould for the peak value on the spectrum
                    % former value, 35 for in vitro
                    % another threshould

peak_t=zeros;        % timings of the peaks
peak_v=zeros;        % peak value
peak_f=zeros;        % peak frequency location

for t=(1000/dt):length(time)  % start analysis after 1000ms
   
    if t*dt==T
    
    a=a+1;         % counter for the sweeps
    % figure(a);
    [A,B,C,D]=spect_peak(LFP(t_old:t),dt/1000,50);
            
    if D>peak_T
    peak_t(a)=t*dt;
    peak_v(a)=D;        % peak value
    peak_f(a)=C;        % peak frequency location
    else
        peak_t(a)=0;
        peak_v(a)=0;
        peak_f(a)=0;
    end
    
    if  t_old==1
        t_old=0;
    end
    
    %title(sprintf('%d - %d time s',t_old*dt/1000, t*dt/1000));
    %ylim([0 6]);
    
    T=T+dT;         % time for the next sweep 
    t_old=t;        % moment of the previous sweep
    
    end
    
end


figure('units','normalized','outerposition',[0 0 1 1]); % full screen

path=PATH;
%ceil(length(Npath)/841*100);


subplot(3,1,1);
plot(peak_t,peak_f,'.','Markersize',10);
title(sprintf('KCC2(-) %d ',path));
xlim([0 time(end)]);
set(gca,'Fontsize',10);
xlabel('time, ms');
ylabel('Peak frequency, Hz');

subplot(3,1,2);
plot(peak_t,peak_v,'.','Markersize',10);
xlim([0 time(end)]);
set(gca,'Fontsize',10);
xlabel('time, ms');
ylabel('Peak amlitude');

subplot(3,1,3);
plot(time,LFP);
set(gca,'Fontsize',10);
xlabel('time, ms');
ylabel('LFP, \muV');
%%
%
%savefig(sprintf('KCC2(-)%d_LFP.fig',path));           % safe *.fig

saveas(gcf,sprintf('KCC2(-)%d_LFP.jpg',path),'jpg');  % save *.jpg

save(sprintf('KCC2(-)%d_LFP.mat',path),'peak_f','peak_t','peak_T','peak_v','path','T','time','LFP');

% RECORD PEAK POSITION AND TIME TO SEIZURE

if max(peak_f)>0      % remember peak and timing to the peak if there is one
idx = find(peak_f~=0, 1, 'first'); % index of the firs element == 0


if peak_f(idx)>=max(peak_f)
PEAK_F(i)=peak_f(idx);        % >= than the max peak
PEAK_T(i)=peak_t(idx);
else
PEAK_F(i)=max(peak_f);        % if peak is smaller than max - take max peak
PEAK_T(i)=peak_t(find(peak_f==max(peak_f)));
end


%ind=find(peak_f==max(peak_f)); % find the index of the maximal frequency peak, it is always decreasing

cells(i)=PATH;

else
    
PEAK_F(i)=0;
PEAK_T(i)=0;
cells(i)=PATH;

end    

%%
close
end

%}