rmeddis@0: function testRP(BFs,MAPparamsName,paramChanges)
rmeddis@0: % testIHC used for IHC I/O function
rmeddis@0: 
rmeddis@0: global experiment method inputStimulusParams
rmeddis@0: global stimulusParameters IHC_VResp_VivoParams IHC_cilia_RPParams
rmeddis@0: savePath=path;
rmeddis@0: addpath (['..' filesep 'utilities'],['..' filesep 'MAP'])
rmeddis@0: dbstop if error
rmeddis@0: 
rmeddis@0: figure(4), clf,
rmeddis@0: set (gcf, 'name', ['IHC'])
rmeddis@0: % set(gcf,'position',[613   354   360   322])
rmeddis@0: drawColors='rgbkmcy';
rmeddis@0: drawnow
rmeddis@0: 
rmeddis@0: if nargin<3, paramChanges=[]; end
rmeddis@0: 
rmeddis@0: levels=-20:10:100;
rmeddis@0: nLevels=length(levels);
rmeddis@0: toneDuration=.05;
rmeddis@0: silenceDuration=.01;
rmeddis@0: sampleRate=50000;
rmeddis@0: dt=1/sampleRate;
rmeddis@0: 
rmeddis@0: allIHC_RP_peak=[];
rmeddis@0: allIHC_RP_dc=[];
rmeddis@0: 
rmeddis@0: changes{1}=[];
rmeddis@0: % changes{1}='IHC_cilia_RPParams.Et=	0.100;';
rmeddis@0: % changes{2}='IHC_cilia_RPParams.Et=	0.070;';
rmeddis@0: 
rmeddis@0: for changeNo=1:length(changes)
rmeddis@0:     if isempty(changes{1})
rmeddis@0:         paramChanges=[];
rmeddis@0:     else
rmeddis@0:         paramChanges{1}=changes{changeNo};
rmeddis@0:     end
rmeddis@0: for BFno=1:length(BFs)
rmeddis@0:     BF=BFs(BFno);
rmeddis@0:     targetFrequency=BF;
rmeddis@0:     % OR
rmeddis@0:     %Patuzzi and Sellick test (see ELP & AEM, 2006)
rmeddis@0:     % targetFrequency=100;
rmeddis@0: 
rmeddis@0:     IHC_RP_peak=zeros(nLevels,1);
rmeddis@0:     IHC_RP_min=zeros(nLevels,1);
rmeddis@0:     IHC_RP_dc=zeros(nLevels,1);
rmeddis@0: 
rmeddis@0:     time=dt:dt:toneDuration;
rmeddis@0: 
rmeddis@0:     rampDuration=0.004;
rmeddis@0:     rampTime=dt:dt:rampDuration;
rmeddis@0:     ramp=[0.5*(1+cos(2*pi*rampTime/(2*rampDuration)+pi)) ...
rmeddis@0:         ones(1,length(time)-length(rampTime))];
rmeddis@0:     ramp=ramp.*fliplr(ramp);
rmeddis@0: 
rmeddis@0:     silence=zeros(1,round(silenceDuration/dt));
rmeddis@0: 
rmeddis@0:     toneStartptr=length(silence)+1;
rmeddis@0:     toneMidptr=toneStartptr+round(toneDuration/(2*dt)) -1;
rmeddis@0:     toneEndptr=toneStartptr+round(toneDuration/dt) -1;
rmeddis@0: 
rmeddis@0:     levelNo=0;
rmeddis@0:     for leveldB=levels
rmeddis@0:         levelNo=levelNo+1;
rmeddis@0:         % replicate at all levels
rmeddis@0:         amp=28e-6*10^(leveldB/20);
rmeddis@0: 
rmeddis@0:         %% create signal (leveldB/ frequency)
rmeddis@0:         inputSignal=amp*sin(2*pi*targetFrequency'*time);
rmeddis@0:         inputSignal= ramp.*inputSignal;
rmeddis@0:         inputSignal=[silence inputSignal silence];
rmeddis@0:         inputStimulusParams.sampleRate=1/dt;
rmeddis@0: %         global IHC_ciliaParams
rmeddis@0: 
rmeddis@0:         %% disable efferent for fast processing
rmeddis@0:         method.DRNLSave=1;
rmeddis@0:         method.IHC_cilia_RPSave=1;
rmeddis@0:         method.IHCpreSynapseSave=1;
rmeddis@0:         method.IHC_cilia_RPSave=1;
rmeddis@0:         method.segmentDuration=-1;
rmeddis@0:         moduleSequence=1:4;
rmeddis@0: 
rmeddis@0:         %% run the model
rmeddis@0:         global  DRNLoutput IHC_cilia_output IHCrestingCiliaCond IHCrestingV...
rmeddis@0:             IHCoutput
rmeddis@0:         AN_spikesOrProbability='probability';
rmeddis@0: 
rmeddis@0:         MAP1_14(inputSignal, sampleRate, BF, ...
rmeddis@0:             MAPparamsName, AN_spikesOrProbability, paramChanges);
rmeddis@0: 
rmeddis@0:         % DRNL
rmeddis@0:         DRNLoutput=DRNLoutput;
rmeddis@0:         DRNL_peak(levelNo,1)=max(DRNLoutput(toneMidptr:toneEndptr));
rmeddis@0:         DRNL_min(levelNo,1)=min(DRNLoutput(toneMidptr:toneEndptr));
rmeddis@0:         DRNL_dc(levelNo,1)=mean(DRNLoutput(toneMidptr:toneEndptr));
rmeddis@0: 
rmeddis@0:         % cilia
rmeddis@0:         IHC_ciliaData=IHC_cilia_output;
rmeddis@0:         IHC_ciliaData=IHC_ciliaData;
rmeddis@0:         IHC_cilia_peak(levelNo,1)=...
rmeddis@0:             max(IHC_ciliaData(toneMidptr:toneEndptr));
rmeddis@0:         IHC_cilia_min(levelNo,1)=...
rmeddis@0:             min(IHC_ciliaData(toneMidptr:toneEndptr));
rmeddis@0:         IHC_cilia_dc(levelNo,1)=...
rmeddis@0:             mean(IHC_ciliaData(toneMidptr:toneEndptr));
rmeddis@0: 
rmeddis@0:         % RP
rmeddis@0:         IHC_RPData=IHCoutput;
rmeddis@0:         IHC_RPData=IHC_RPData;
rmeddis@0:         IHC_RP_peak(levelNo,1)=...
rmeddis@0:             max(IHC_RPData(toneMidptr:toneEndptr));
rmeddis@0:         IHC_RP_min(levelNo,1)=...
rmeddis@0:             min(IHC_RPData(toneMidptr:toneEndptr));
rmeddis@0:         IHC_RP_dc(levelNo,1)=...
rmeddis@0:             mean(IHC_RPData(toneMidptr:toneEndptr));
rmeddis@0: 
rmeddis@0:     end % level
rmeddis@0: 
rmeddis@0: 
rmeddis@0:     disp(['parameter file was: ' MAPparamsName])
rmeddis@0:     fprintf('\n')
rmeddis@0: 
rmeddis@0:     %% plot DRNL
rmeddis@0:     subplot(2,2,1)
rmeddis@0: %     referenceDisp= 9e-9*1000/targetFrequency;
rmeddis@0: %     plot(levels,20*log10(DRNL_peak/referenceDisp), drawColors(BFno), ...
rmeddis@0: %         'linewidth',2), hold on
rmeddis@0:     referenceDisp=10e-9;
rmeddis@0:     plot(levels,20*log10(DRNL_peak/referenceDisp), drawColors(BFno), ...
rmeddis@0:         'linewidth',2), hold on
rmeddis@0:     title([' DRNL peak:  ' num2str(BFs) ' Hz'])
rmeddis@0:     ylabel ('log10DRNL(m)'), xlabel('dB SPL')
rmeddis@0:     xlim([min(levels) max(levels)]), ylim([-10 50])
rmeddis@0:     grid on
rmeddis@0: 
rmeddis@0:     %% plot cilia displacement
rmeddis@0:     figure(4)
rmeddis@0:     subplot(2,2,2)
rmeddis@0:     restingIHC_cilia=IHCrestingCiliaCond;
rmeddis@0:     plot(levels, IHC_cilia_peak,'k', 'linewidth',2), hold on
rmeddis@0:     plot(levels, IHC_cilia_min,'r', 'linewidth',2)
rmeddis@0:     hold on,
rmeddis@0:     plot([min(levels) max(levels)], ...
rmeddis@0:         [restingIHC_cilia restingIHC_cilia], 'g')
rmeddis@0:     title(' IHC apical cond.')
rmeddis@0:     ylabel ('IHCcilia(conductance)'), xlabel('dB SPL')
rmeddis@0:     xlim([min(levels) max(levels)])
rmeddis@0:     grid on
rmeddis@0: 
rmeddis@0:     %% plot receptor potentials
rmeddis@0:     figure(4)
rmeddis@0:     subplot(2,2,3)
rmeddis@0:     % RP I/O function min and max
rmeddis@0:     restingRP=IHC_RP_peak(1);
rmeddis@0:     restingRP=IHCrestingV;
rmeddis@0:     toPlot= [fliplr(IHC_RP_min(:,1)') IHC_RP_peak(:,1)'];
rmeddis@0:     microPa=   28e-6*10.^(levels/20);
rmeddis@0:     microPa=[-fliplr(microPa) microPa];
rmeddis@0:     plot(microPa,toPlot, drawColors(BFno), 'linewidth',2)
rmeddis@0:     % ylim([0 300])
rmeddis@0: 
rmeddis@0:     %% Dallos and Harris data
rmeddis@0:     dallosx=[-0.9	-0.1	-0.001	0.001	0.01	0.9];
rmeddis@0:     dallosy=[-8	-7.8	-6.5	11	16.5	22]/1000 + restingRP;
rmeddis@0:     hold on, plot(dallosx,dallosy, 'o')
rmeddis@0:     plot([-1 1], [restingRP restingRP], 'r')
rmeddis@0:     title(' Dallos(86) data at 800 Hz')
rmeddis@0:     ylabel ('receptor potential(V)'), xlabel('Pa')
rmeddis@0:     ylim([-0.08 -0.02]), xlim([-1 1])
rmeddis@0:     grid on
rmeddis@0: 
rmeddis@0:     %% RP I/O function min and max
rmeddis@0:     figure(4)
rmeddis@0:     subplot(2,2,4)
rmeddis@0:     restingRP=IHC_RP_peak(1);
rmeddis@0:     peakRP=max(IHC_RP_peak);
rmeddis@0:     plot(levels, IHC_RP_peak,drawColors(BFno), 'linewidth',2)
rmeddis@0:     hold on
rmeddis@0:     plot(levels, IHC_RP_dc, [drawColors(BFno) ':'], 'linewidth',2)
rmeddis@0:     hold on,
rmeddis@0:     plot([min(levels) max(levels)], [restingRP restingRP], 'r')
rmeddis@0:     xlim([min(levels) max(levels)])
rmeddis@0:     grid on
rmeddis@0:     title(['Et= ' num2str(IHC_cilia_RPParams.Et) ':  RP (AC- / DC:'])
rmeddis@0:     ylabel ('RP(V)'), xlabel('dB SPL')
rmeddis@0:     ylim([-0.08 -0.02])
rmeddis@0:     allIHC_RP_peak=[allIHC_RP_peak IHC_RP_peak];
rmeddis@0:     allIHC_RP_dc=[allIHC_RP_dc IHC_RP_dc];
rmeddis@0: 
rmeddis@0:     fprintf('level\t peak\t DC\n')
rmeddis@0:     UTIL_printTabTable([levels' IHC_RP_peak IHC_RP_dc])
rmeddis@0:     % disp(['restingIHC_cilia= ' num2str(restingIHC_cilia)])
rmeddis@0:     fprintf('peakRP= \t%6.3f', peakRP)
rmeddis@0:     fprintf('\nrestingRP= \t%6.3f', restingRP)
rmeddis@0:     fprintf('\ndifference= \t%6.3f\n', (peakRP-restingRP))
rmeddis@0:     drawnow
rmeddis@0: end
rmeddis@0: end
rmeddis@0: % UTIL_showStruct(IHC_VResp_VivoParams, 'IHC_VResp_VivoParams')
rmeddis@0: UTIL_showStruct(IHC_cilia_RPParams, 'IHC_cilia_RPParams')
rmeddis@0:     fprintf('level\t peak\n')
rmeddis@0:     UTIL_printTabTable([levels' allIHC_RP_peak])
rmeddis@0:     fprintf('level\t DC\n')
rmeddis@0:     UTIL_printTabTable([levels' allIHC_RP_dc])
rmeddis@0: 
rmeddis@0: path(savePath);
rmeddis@0: disp(paramChanges)