fid = 'sushi19308_12DIV_pre-drug_01';
t_end = loadspk([fid '.spk'],'last');
spk = loadspk([fid '.spk'],[0 t_end],0);
ss = SqueakySpk(fid,spk.meta.fs_Hz,1,spk);

%length(ss.time(ss.unit>0))/max(ss.time)/max(ss.unit)

avgFR_1hr = length(ss.time(ss.unit>0 & ss.time>max(ss.time)-1*3600))/(1*3600)/max(ss.unit)
avgFR_3hr = length(ss.time(ss.unit>0 & ss.time>max(ss.time)-3*3600))/(3*3600)/max(ss.unit)
avgFR_6hr = length(ss.time(ss.unit>0 & ss.time>max(ss.time)-6*3600))/(6*3600)/max(ss.unit)
avgFR_12hr = length(ss.time(ss.unit>0 & ss.time>max(ss.time)-12*3600))/(12*3600)/max(ss.unit)

%%
bin = max(ss.time(ss.unit>0))-12*3600:60:max(ss.time(ss.unit>0));
asdr = histc(ss.time(ss.unit>0),bin);
plot(bin,asdr/max(ss.unit)/60,'k-'); hold on;
plot([max(ss.time(ss.unit>0))-1*3600,max(ss.time(ss.unit>0))],[avgFR_1hr,avgFR_1hr],'-m','LineWidth',2);
plot([max(ss.time(ss.unit>0))-3*3600,max(ss.time(ss.unit>0))],[avgFR_3hr,avgFR_3hr],'-g','LineWidth',2);
plot([max(ss.time(ss.unit>0))-6*3600,max(ss.time(ss.unit>0))],[avgFR_6hr,avgFR_6hr],'-b','LineWidth',2);
plot([max(ss.time(ss.unit>0))-12*3600,max(ss.time(ss.unit>0))],[avgFR_12hr,avgFR_12hr],'-r','LineWidth',2);