load vivy21369C_10-12DIV.mat;
% load vivy21369D_10-12DIV.mat;
% load vivy21369F_10-12DIV.mat;

dt = 1;
bin = [ceil(min(time)/dt)*dt:dt:floor(max(time)/dt)*dt];
asdr = histc(time,bin);

figure;
set(gcf,'Position',[455 361 2105 384]);
hold off;
stairs(bin/3600,asdr/dt,'-k','LineWidth',0.5);
hold on;
y_lim = ylim;

plot([timestamps(1) timestamps(1)]/3600,[0 y_lim(2)],'--r');
xlim([-6 24]);
set(gca,'XTick',[-20:2:28],'TickDir','out','TickLength',[0.0025 0.025]);
box off;

ylabel('Array-wide firing rate (Hz)','FontSize',14);
xlabel('Time (hours)','FontSize',14);