MAP

% function test2toneMovie

%

% Demonstration of two-tone suppression in the AN using a

%   F1 tone suppressed by a 1.5*F1 kHz tone

%

global  ANprobRateOutput DRNLoutput

dbstop if error

restorePath=path;

addpath (['..' filesep 'MAP'],    ['..' filesep 'wavFileStore'], ...

    ['..' filesep 'utilities'])

figure(5), clf

disp('')

primaryToneFrequency=8000;

probedBs=10 :10: 50;

nProbeLevels=length(probedBs);

% disp(['F1 F1 level= ' num2str([primaryToneFrequency probedB])])

suppressorLevels=0:10:80;

lowestBF=250; 	highestBF= 8000; 	numChannels=21;

BFlist=round(logspace(log10(lowestBF), log10(highestBF), numChannels));

% fix so that probe is middle range and is at a BF

BFlist=BFlist*primaryToneFrequency/2000;

suppressorFrequencies=BFlist;

BFchannel=find(round(BFlist)==round(primaryToneFrequency));

sampleRate= 100000; % Hz (higher sample rate needed for BF>8000 Hz)

dt=1/sampleRate; % seconds

duration=.050;		      % seconds

    tone2Duration=duration/2; % s

startSilenceDuration=0.010;

startSilence=...

    zeros(1,startSilenceDuration*sampleRate);

suppressorStartsPTR=...

    round((startSilenceDuration+duration/2)*sampleRate);

probedBCount=0;

for probedB=probedBs

    probedBCount=probedBCount+1;

peakResponse= zeros(length(suppressorLevels),length(suppressorFrequencies));

suppFreqCount=0;

for suppressorFrequency=suppressorFrequencies

    suppFreqCount=suppFreqCount+1;

    suppressorLevelCount=0;

for suppressorDB=suppressorLevels

suppressorLevelCount=suppressorLevelCount+1;

    % primary + suppressor

    primaryLevelDB=probedB;

    suppressorLevelDB=suppressorDB;

% primary BF tone

    time1=dt: dt: duration;

    amp=10^(primaryLevelDB/20)*28e-6;

    inputSignal=amp*sin(2*pi*primaryToneFrequency*time1);

    rampDuration=.005; rampTime=dt:dt:rampDuration;

    ramp=[0.5*(1+cos(2*pi*rampTime/(2*rampDuration)+pi)) ones(1,length(time1)-length(rampTime))];

    inputSignal=inputSignal.*ramp;

    inputSignal=inputSignal.*fliplr(ramp);

    % suppressor

    time2= dt: dt: tone2Duration;

    % B: tone on

    amp=10^(suppressorLevelDB/20)*28e-6;

    inputSignal2=amp*sin(2*pi*suppressorFrequency*time2);

    rampDuration=.005; rampTime=dt:dt:rampDuration;

    ramp=[0.5*(1+cos(2*pi*rampTime/(2*rampDuration)+pi)) ones(1,length(time2)-length(rampTime))];

    inputSignal2=inputSignal2.*ramp;

    inputSignal2=inputSignal2.*fliplr(ramp);

    % A: initial silence (delay to suppressor)

    silence=zeros(1,length(time2));

    inputSignal2=[silence inputSignal2];

    % add tone and suppressor components

    inputSignal=inputSignal+inputSignal2;

    inputSignal=...

        [startSilence inputSignal ];

    % run MAP

    MAPparamsName='Normal';

    AN_spikesOrProbability='probability';

    paramChanges={'IHCpreSynapseParams.tauCa=80e-6;'};

    MAP1_14(inputSignal, sampleRate, BFlist, ...

        MAPparamsName, AN_spikesOrProbability, paramChanges);

    response=ANprobRateOutput;

peakResponse(suppressorLevelCount,suppFreqCount)=...

    mean(response(BFchannel,suppressorStartsPTR:end));

disp(['F2 level= ', num2str([suppressorFrequency suppressorDB ...

    peakResponse(suppressorLevelCount,suppFreqCount)])])

%% make movie

    figure(6), subplot(4,1,2)

    axis tight

    set(gca,'nextplot','replacechildren');

    plot(dt:dt:dt*length(inputSignal), inputSignal, 'k')

    ylim([-.5 .5])

    figure(6), subplot(2,1,2)

    PSTHbinWidth=0.010;

    PSTH= UTIL_PSTHmakerb(...

        response, dt, PSTHbinWidth);

    [nY nX]=size(PSTH);

    time=PSTHbinWidth*(0:nX-1);

    surf(time, BFlist, PSTH)

    zlim([0 500]),

    caxis([0 500])

    shading interp

    colormap(jet)

    set(gca, 'yScale','log')

    xlim([0 max(time)+dt])

    ylim([0 max(BFlist)])

    view([0 90]) % view([-8 68])

    title([' primary + suppressor '])

end

end

%% plot matrix

peakResponse=peakResponse/peakResponse(1,1);

figure(5), subplot(2,nProbeLevels, probedBCount)

surf(suppressorFrequencies,suppressorLevels,peakResponse)

shading interp

view([0 90])

set(gca,'Xscale','log')

xlim([min(suppressorFrequencies) max(suppressorFrequencies)])

set(gca,'Xtick', [1000  4000],'xticklabel',{'1000', '4000'})

title([num2str(primaryToneFrequency) ' Hz'])

subplot(2,nProbeLevels,nProbeLevels+probedBCount)

contourf(suppressorFrequencies,suppressorLevels,peakResponse, ...

    [.1:.1:.9 0.99] )

set(gca,'Xscale','log')

set(gca,'Xtick', [1000  4000],'xticklabel',{'1000', '4000'})

set(gcf, 'name',['primaryToneFrequency= ' num2str(primaryToneFrequency)])

title([num2str(probedB) ' dB'])

end

path=restorePath;