clear all
close all

load('pyramidal_response.mat')
threshold = 20; % mV

figure()
hold on
plot(time, neuronal_data, 'color', [0.2 0.5 0.7], 'linewidth', 1., ...
    'displayname', 'membrane voltage')
plot(spikes, ones(size(spikes)) .* threshold, 'ro', 'markersize', 5, ...
    'displayname', 'spike times')
line([time(1) time(end)], [threshold threshold], 'linestyle', '--', ...
    'linewidth', 0.75, 'color', [0.5 0.5 0.9], 'displayname', 'threshold')

xlabel('time [s]', 'fontname', 'times', 'fontsize', 11)
ylabel('potential [mV]', 'fontname', 'times', 'fontsize', 11)
title('ELL pyramidal neuron', 'fontname', 'times', 'fontsize', 12)
ylim([0 35])
xlim([0 2.25])
box('off')
l = legend(gca,'show');
set(l,'location','northwest', 'fontsize', 7, 'linewidth', 1.);
set(gca, 'xminortick','on','yminortick','on')
set(gca, 'tickdir','out', 'linewidth', 1.5, 'fontname', 'times', ...
    'fontsize', 11)
set(gcf, 'paperunits', 'centimeters', 'papersize', [15 7.5]);
set(gcf, 'paperposition',[0.0 0.0 15, 7.5], 'color', 'white')

saveas(gcf, [pwd '/../lecture/images/spike_detection.pdf'], 'pdf')