%% Make heatmaps for tone responses
type = 'events'; % Manipulations, events
field = 'tone';
regions = {'ENTm','DG','CA3','CA1'};
colors = [.5 .3 1; 0 1 0; 1 1 0; 0 0 1];
minRate = 0;
minPrctile = 
printfig = true;

% Create in- and exclude structure
include = struct('brainRegion',regions);
%[include.labels] = deal(opsins);
exclude = struct('cell_type',repmat({'unknown'},1,4));
criteria = struct('firingRate',repmat({[0 Inf]},1,4));

% DG plot
f1 = figure('position',[200 20 200*4 350]);
plot_psth(batch_metrics,field,'type',type,'labels',regions,...
    'include',include,'exclude',exclude,...
    'criteria',criteria,'orientation','row',...
    'colors',colors,'minRate',minRate,'minPrctile',minPrctile);

if printfig
    pos = get(f1,'Position');
    set(f1,'PaperPositionMode','Auto','PaperUnits','points',...
        'PaperSize',[pos(3), pos(4)],'Renderer','Painters')
    print(f1,sprintf('DS2_response_DG.pdf'),'-dpdf','-r0')
end