close all;
clear all;
trgnum=input('trigger number (as from table) ');
filename=sprintf('00%u000\\bat\\rate\\sw00%u000brtms.lc',trgnum, trgnum)
gfilename=strcat(filename,'.gz');
gunzip(gfilename);
info=fitsinfo(filename)
%info.PrimaryData.Keywords
for i=1:length(info.PrimaryData.Keywords)
    if(strcmp(info.PrimaryData.Keywords{i,1},'OBJECT'))
        obj=info.PrimaryData.Keywords{i,2};
    end
    if (strcmp(info.PrimaryData.Keywords{i,1},'TRIGTIME'))
        trgt=info.PrimaryData.Keywords{i,2};
    end
end
data=fitsread(filename,'bintable');
gzip(filename);
delete(filename);
plot(data{:,1}-trgt, data{:,2});
range(1)=input('input seconds before trigger: ')
range(2)=input('input seconds after trigger: ')
%range=600;
t=data{:,1}-trgt;
counts=data{:,2}; %4 counters, on SWIFT
clear data
f=1/(t(2)-t(1)); %sampling frequency;
range_idx=floor(range*f)
idx0=find(abs(t)<1/f, 1); %find index of trigger time within the precision
s_i=counts(idx0-range_idx(1):idx0+range_idx(2), 1);
tt=t(idx0-range_idx(1):idx0+range_idx(2));
[b,a]=butter(5,.1);% Set b and a for a 5th order butterworth filter at 1/10 fs.
y=filter(b,a,s_i);
figure;
plot(tt,s_i,'b',tt,y,'r');
%%%%%%%%Analisi con CAO sul segnale filtrato%%%%%%%%%%%%%%%%%
ss=signal (y);
a=amutual(ss,128);
laga=firstmin(a);
c = cao(ss,30,laga,10,500);
cmax=firstmax(c);
e1 = embed(ss, 3, laga);
figure
view(e1)
%%%%%%%%%%%%% Analisi con SSA %%%%%%%%%%%%%%%
L=floor(length(s_i)/2); %%%%%%%%%%% Lunghezza finestra per ssa : pari alla metà della lunghezza della serie temporale
%
I=[1:12,L]; %%%%%%%%%%%%%%%Punto delicato per ottenere curva di luce filtrata da noise
figure;
[y,sr,vr]=ssa(s_i,ceil(L/3),I);