SS.ASDR(3600,[ceil(min(SS.time)/3600)*3600 ceil(max(SS.time)/3600)*3600])
SS.PlotUSDR([],1,[ceil(min(SS.time)/3600)*3600 ceil(max(SS.time)/3600)*3600],1,1); colormap([jet])
set(gca,'XTick',[-36:3:36]*3600/SS.asdr.dt,'XTickLabel',[-36:3:36],'TickDir','out')
xlim([-6 24])
set(gca,'FontSize',16)
xlabel('Time (hours)','FontSize',16);
ylabel('Unit #')

export_fig([SS.name 'USDR_bin1hr_logscale'],'-png','-transparent',gcf)


%%

numUnits = max(SS.unit);

FR = zeros(numUnits,1);

for i = 1:numUnits
    FR(i,1) = length(SS.time(SS.unit==i & SS.time<0 & SS.time>-3*3600))/(3*3600);
%     FR(i,2) = length(SS.time(SS.unit==i & SS.time>0 & SS.time<24*3600))/(24*3600);
end

FR = [nan; FR]

SS.usdr.usdr = SS.usdr.usdr./FR';

SS.PlotUSDR([],1,[ceil(min(SS.time)/3600)*3600 ceil(max(SS.time)/3600)*3600],1,0); colormap([jet])
set(gca,'XTick',[-36:3:36]*3600/SS.asdr.dt,'XTickLabel',[-36:3:36],'TickDir','out')
xlim([-6 24])
set(gca,'FontSize',16)
xlabel('Time (hours)','FontSize',16);
ylabel('Unit #')

% export_fig([SS.name 'USDR_bin1hr_logscale'],'-png','-transparent',gcf)