fidlist = {'2011-04-26_pulse-sweep-raw_TTX_wSALPA','2011-04-26_pulse-sweep-raw_TTX_woSALPA',...
    '2011-04-26_pulse-sweep-raw_preTTX_wSALPA','2011-04-26_pulse-sweep-raw_preTTX_woSALPA',...
    '2011-04-26_pulse-sweep-raw_postTTX_wSALPA','2011-04-26_pulse-sweep-raw_postTTX_woSALPA'};

for i = 1:6
    x = loadspike([fidlist{i} '.spk']);
    c = x.channel(x.time<60);
    i_badspk = find(c==1 | c==8 | c==31 | c==33 | c==41 | c==57 | c==63 | c==64);
    n_badspk = length(i_badspk);
    n{i} = length(c) - n_badspk;
    subplot(6,1,i);
    plot(x.time,x.channel,'.','Color',[0 0 0]); axis([0 60 0 65]); hold on;
    plot(x.time(i_badspk),x.channel(i_badspk),'.','Color',[0.7 0.7 0.7]);
end

% n{1}; % falsely detected spikes w/ SALPA
% n{2}; % falsely detected spikes w/o SALPA
% n{3}; % total detected spikes w/ SALPA before adding TTX
% n{4}; % total detected spikes w/o SALPA before adding TTX
% n{5}; % total detected spikes w/ SALPA after wash
% n{6}; % total detected spikes w/o SALPA after wash

dDR_TTX = (n{2}-n{1})/n{2}
dDR_pre = (n{4}-n{3})/n{4}
dDR_post = (n{6}-n{5})/n{6}

