fid = 'octave358_13DIV_40um-cnqx-CL-stim';
s = loadspk([fid '.spk']);
b = loadbang([fid '.bang']);

%%
close all
figure
ss = SqueakySpk(fid,s.meta.fs_Hz,1,s);
subplot(311)
ss.ASDR(1)
subplot(312)
ss.ASDR(300)
subplot(313)
ss.ASDR(3600)
hold on
plot([0 max(ss.asdr.bin)],[1.4360 1.4360],'r')
plot(b.time,b.out1)

