comparison multithreshold 1.46/testAN.m @ 9:ecad0ea62b43

May27 mainly better parameters

-function testAN
+function vectorStrength=testAN(targetFrequency,BFlist, levels)
 % testIHC used either for IHC I/O function ...
 %  or receptive field (doReceptiveFields=1)
 
 global experiment method stimulusParameters
-global IHC_VResp_VivoParams IHCpreSynapseParams
+global IHC_VResp_VivoParams  IHC_cilia_RPParams IHCpreSynapseParams
 global AN_IHCsynapseParams
 % global saveMembranePotential MacGregorMultiParams
 dbstop if error
 
 addpath (['..' filesep 'MAP'], ['..' filesep 'utilities'], ...
     ['..' filesep 'parameterStore'],  ['..' filesep 'wavFileStore'],...
     ['..' filesep 'testPrograms'])
 
-levels=-10:10:80;
-% levels=80:10:90;
+if nargin<3
+    levels=-10:10:80;
+    % levels=80:10:90;
+end
+
 nLevels=length(levels);
 
 toneDuration=.2;
 rampDuration=0.002;
 silenceDuration=.02;
 localPSTHbinwidth=0.001;
 
 % Use only the first frequency in the GUI targetFrequency box to defineBF
-targetFrequency=stimulusParameters.targetFrequency(1);
-BFlist=targetFrequency;
+% targetFrequency=stimulusParameters.targetFrequency(1);
+% BFlist=targetFrequency;
 
 AN_HSRonset=zeros(nLevels,1);
 AN_HSRsaturated=zeros(nLevels,1);
 AN_LSRonset=zeros(nLevels,1);
 AN_LSRsaturated=zeros(nLevels,1);

 % ANoutput=zeros(200,200);
 
 figure(15), clf
 set(gcf,'position',[980   356   401   321])
 figure(5), clf
- set(gcf,'position', [980 34 400 295])
+set(gcf,'position', [980 34 400 295])
+set(gcf,'name',[num2str(BFlist), ' Hz']);
 drawnow
 
 levelNo=0;
 for leveldB=levels
     levelNo=levelNo+1;
 
-    % sample rate should be amultiple of the targetFrequency for PSTH below
+    %% sample rate should be amultiple of the targetFrequency for PSTH below
     sampleRate=50000;
+    sampleRate=20*targetFrequency;
+    if sampleRate<20000
+    sampleRate=round(40000/targetFrequency)*targetFrequency;
+    end
+    
+     
+    %% ananana
     dt=1/sampleRate;
     period=1/targetFrequency;
     dt=dt*(dt/period)*round(period/dt);
 
     fprintf('%4.0f\t', leveldB)

     MAPparamsName=experiment.name;
     showPlotsAndDetails=0;
 
     global ANoutput CNoutput ICoutput ICmembraneOutput tauCas
     global ARattenuation MOCattenuation
-    
-    MAP1_14(inputSignal, 1/dt, targetFrequency, ...
+
+    MAP1_14(inputSignal, 1/dt, BFlist, ...
         MAPparamsName, AN_spikesOrProbability);
 
     nTaus=length(tauCas);
 
     %LSR (same as HSR if no LSR fibers present)

 
     figure(5), subplot(2,2,2)
     hold off, bar(PSTHtime,PSTH, 'b')
     hold on,  bar(PSTHtime,PSTHLSR,'r')
     ylim([0 1000])
-xlim([0 length(PSTH)*localPSTHbinwidth])
-    
+    xlim([0 length(PSTH)*localPSTHbinwidth])
+
     % AN - CV
     %  CV is computed 5 times. Use the middle one (3) as most typical
     cvANHSR=  UTIL_CV(HSRspikes, dt);
 
     % AN - vector strength

     plot(time,ICmembraneOutput(2, 1:end),'k')
     ylim([-0.07 0])
     xlim([0 max(time)])
     title(['IC  ' num2str(leveldB,'%4.0f') 'dB'])
     drawnow
-    
+
     figure(5), subplot(2,2,1)
     plot(20*log10(MOC), 'k'),
     title(' MOC'), ylabel('dB attenuation')
     ylim([-30 0])
 
 
 end % level
 figure(5), subplot(2,2,1)
-    plot(levels,20*log10(MOC), 'k'),
-    title(' MOC'), ylabel('dB attenuation')
-    ylim([-30 0])
+plot(levels,20*log10(MOC), 'k'),
+title(' MOC'), ylabel('dB attenuation')
+ylim([-30 0])
 xlim([0 max(levels)])
 
 fprintf('\n')
 toneDuration=2;
 rampDuration=0.004;

 
 
 % ---------------------------------------------------- display parameters
 
 figure(15), clf
-set(gcf,'position',[1000   356   381   321])
 nRows=2; nCols=2;
 
 % AN rate - level ONSET functions
 subplot(nRows,nCols,1)
 plot(levels,AN_LSRonset,'ro'), hold on

 ylim([0 400])
 set(gca,'ytick',0:50:300)
 xlim([min(levels) max(levels)])
 set(gca,'xtick',[levels(1):20:levels(end)]), grid on
 
+
 ttl=['spont=' num2str(mean(ICHSRsaturated(1,:)),'%4.0f') ...
     '  sat=' num2str(mean(ICHSRsaturated(end,1)),'%4.0f')];
 title( ttl)
 xlabel('level dB SPL'), ylabel ('adapted rate (sp/s)')
 text(0, 350, 'IC', 'fontsize', 16)
 set(gcf,'name',' AN CN IC rate/level')
 
-UTIL_showStruct(IHCpreSynapseParams, 'IHCpreSynapseParams')
-UTIL_showStruct(AN_IHCsynapseParams, 'AN_IHCsynapseParams')
+peakVectorStrength=max(vectorStrength);
 
 fprintf('\n')
 disp('levels vectorStrength')
 fprintf('%3.0f \t %6.4f \n', [levels; vectorStrength'])
 fprintf('\n')
-fprintf('Phase locking, max vector strength= %6.4f\n\n',...
+fprintf('Phase locking, max vector strength=\t %6.4f\n\n',...
     max(vectorStrength))
 
 allData=[ levels'  AN_HSRonset AN_HSRsaturated...
     AN_LSRonset AN_LSRsaturated ...
     CNHSRsaturated CNLSRrate...

 fprintf('\n levels \tANHSR Onset \tANHSR adapted\tANLSR Onset \tANLSR adapted\tCNHSR\tCNLSR\tICHSR  \tICLSR \n');
 UTIL_printTabTable(round(allData))
 fprintf('VS (phase locking)= \t%6.4f\n\n',...
     max(vectorStrength))
 
+UTIL_showStruct(IHC_cilia_RPParams, 'IHC_cilia_RPParams')
+UTIL_showStruct(IHCpreSynapseParams, 'IHCpreSynapseParams')
+UTIL_showStruct(AN_IHCsynapseParams, 'AN_IHCsynapseParams')
+
+fprintf('\n')
+disp('levels vectorStrength')
+fprintf('%3.0f \t %6.4f \n', [levels; vectorStrength'])
+fprintf('\n')
+fprintf('Phase locking, max vector strength= \t%6.4f\n\n',...
+    max(vectorStrength))
+
+allData=[ levels'  AN_HSRonset AN_HSRsaturated...
+    AN_LSRonset AN_LSRsaturated ...
+    CNHSRsaturated CNLSRrate...
+    ICHSRsaturated ICLSRsaturated];
+fprintf('\n levels \tANHSR Onset \tANHSR adapted\tANLSR Onset \tANLSR adapted\tCNHSR\tCNLSR\tICHSR  \tICLSR \n');
+UTIL_printTabTable(round(allData))
+fprintf('VS (phase locking)= \t%6.4f\n\n',...
+    max(vectorStrength))
+