MAP

function showMAP (options)

% defaults

dbstop if warning

    DRNLoutput IHC_cilia_output IHCrestingCiliaCond IHCrestingV...

    IHCoutput ANprobRateOutput ANoutput savePavailable tauCas  ...

    CNoutput  ICoutput ICmembraneOutput ICfiberTypeRates MOCattenuation

restorePath=path;

addpath ( ['..' filesep 'utilities'], ['..' filesep 'parameterStore'])

    options.showModelParameters=1;

    options.showModelOutput=1;

    options.printFiringRates=1;

    options.showEfferent=1;

end

if options.showModelParameters

        duration=size(TMoutput,2)*dt;

        disp(['AN: ' num2str(sum(sum(ANoutput(end-nHSRfibers+1:end,:)))/...

            (nHSRfibers*duration))])

        nCNneurons=size(CNoutput,1);

        nHSRCNneuronss=nCNneurons/nANfiberTypes;

        disp(['CN: ' num2str(sum(sum(CNoutput(end-nHSRCNneuronss+1:end,:)))...

            /(nHSRCNneuronss*duration))])

    else

        disp(['AN: ' num2str(mean(mean(ANprobRateOutput)))])

    end

    plotInstructions.figureNo=99;

    signalRMS=mean(savedInputSignal.^2)^0.5;

    signalRMSdb=20*log10(signalRMS/20e-6);

    % plot signal (1)

    plotInstructions.displaydt=dt;

    plotInstructions.numPlots=6;

    r=size(savedInputSignal,1);

    if r==1, savedInputSignal=savedInputSignal'; end

    UTIL_plotMatrix(savedInputSignal', plotInstructions);

    % stapes (2)

    plotInstructions.subPlotNo=2;

    plotInstructions.title= ['stapes displacement'];

    UTIL_plotMatrix(OMEoutput, plotInstructions);

    % DRNL (3)

    plotInstructions.subPlotNo=3;

    plotInstructions.title= ['BM displacement'];

    plotInstructions.yValues= savedBFlist;

    UTIL_plotMatrix(DRNLoutput, plotInstructions);

    switch saveAN_spikesOrProbability

        case 'spikes'

            % AN (4)

                plotInstructions.rasterDotSize=3;

            end

            UTIL_plotMatrix(ANoutput, plotInstructions);

            % CN (5)

            plotInstructions.displaydt=ANdt;

            plotInstructions.subPlotNo=5;

                plotInstructions.rasterDotSize=3;

            end

            UTIL_plotMatrix(CNoutput, plotInstructions);

            % IC (6)

            plotInstructions.displaydt=ANdt;

            plotInstructions.subPlotNo=6;

                plotInstructions.zValuesRange= [-.1 0];

                UTIL_plotMatrix(ICmembraneOutput, plotInstructions);

            end

        otherwise % probability (4-6)

            plotInstructions.displaydt=dt;

            plotInstructions.numPlots=2;

    plotInstructions.title= ['AR strength.  Signal level= ' ...

        num2str(signalRMSdb,'%4.0f') ' dB SPL'];

    UTIL_plotMatrix(20*log10(ARattenuation), plotInstructions);

    plotInstructions.subPlotNo=2;

    plotInstructions.yValues= savedBFlist;

    plotInstructions.yLabel= 'BF';

    colorbar

end

%% ACF plot if required

if options.showACF

    method.dt=dt;

    method.segmentNo=1;

    method.nonlinCF=savedBFlist;

    minPitch=	80; maxPitch=	4000; numPitches=100;    % specify lags

    pitches=10.^ linspace(log10(minPitch), log10(maxPitch),numPitches);

    pitches=fliplr(pitches);

    filteredSACFParams.acfTau=	.003;       % time constant of running ACF

    filteredSACFParams.lambda=	0.12;       % slower filter to smooth ACF

    filteredSACFParams.lambda=	0.01;       % slower filter to smooth ACF

    filteredSACFParams.plotACFs=0;          % special plot (see code)

    filteredSACFParams.plotFilteredSACF=0;  % 0 plots unfiltered ACFs

    filteredSACFParams.plotMoviePauses=.3;          % special plot (see code)

    filteredSACFParams.usePressnitzer=0; % attenuates ACF at  long lags

    filteredSACFParams.lagsProcedure=  'useAllLags';

    % filteredSACFParams.lagsProcedure=  'useBernsteinLagWeights';

    % filteredSACFParams.lagsProcedure=  'omitShortLags';

    filteredSACFParams.criterionForOmittingLags=3;

    filteredSACFParams.plotACFsInterval=200;

    if filteredSACFParams.plotACFs

        % plot original waveform on ACF plot

        figure(13), clf

        xlim([0 t(end)])

        title(['stimulus: ' num2str(signalRMSdb, '%4.0f') ' dB SPL']);

    end

    % plot original waveform on summary/smoothed ACF plot

    subplot(2,1,1)

    t=dt*(1:length(savedInputSignal));

    plot(t,savedInputSignal)

    xlim([0 t(end)])

    title(['stimulus: ' num2str(signalRMSdb, '%4.0f') ' dB SPL']);

    % compute ACF

    switch saveAN_spikesOrProbability

        case 'probability'

        otherwise

            inputToACF=ANoutput;

    end

    disp ('computing ACF...')

    [P, BFlist, sacf, boundaryValue] = ...

        filteredSACF(inputToACF, method, filteredSACFParams);

    disp(' ACF done.')

    % SACF

    subplot(2,1,2)

    imagesc(P)

    set(gca,'xtickLabel', round(100*t(tt))/100)

end