comparison testPrograms/testFM.m @ 38:c2204b18f4a2 tip

End nov big change

 %
 
 
 global inputStimulusParams outerMiddleEarParams DRNLParams
 global IHC_VResp_VivoParams IHCpreSynapseParams  AN_IHCsynapseParams
-global  ANprobRateOutput  ANoutput ANtauCas  ANdt
+global  ANprobRateOutput  ANoutput ANtauCas  dtSpikes
 dbstop if error
 restorePath=path;
 addpath (['..' filesep 'MAP'], ['..' filesep 'utilities'], ...
     ['..' filesep 'parameterStore'],  ['..' filesep 'wavFileStore'],...
     ['..' filesep 'testPrograms'])

     paramsName=('Normal');
     AN_spikesOrProbability='spikes';
     paramChanges=[];
 else
     if nargin<3
-    paramChanges=[];
+        paramChanges=[];
     end
 end
 
 % masker and probe levels are relative to this threshold
 thresholdAtCF=10; % dB SPL

         time=dt: dt: gapDuration;
         gap=zeros(1,length(time));
 
         inputSignal=...
             [initialSilence maskerPa gap probe finalSilence];
-%         time=dt:dt:length(inputSignal)*dt;
-%         figure(99), plot(time,inputSignal)
+        %         time=dt:dt:length(inputSignal)*dt;
+        %         figure(99), plot(time,inputSignal)
 
         % **********************************  run MAP model
-%         showPlotsAndDetails=0;
-nChanges=length(paramChanges);
-paramChanges{nChanges+1}='AN_IHCsynapseParams.numFibers=	500;'; 
+        %         showPlotsAndDetails=0;
+        nChanges=length(paramChanges);
+        paramChanges{nChanges+1}='AN_IHCsynapseParams.numFibers=	500;';
         MAP1_14(inputSignal, 1/dt, targetFrequency, ...
             paramsName, AN_spikesOrProbability, paramChanges);
 
         if strcmp(AN_spikesOrProbability,'probability')
             [nFibers c]=size(ANprobRateOutput);
             nLSRfibers=nFibers/length(ANtauCas);
             ANresponse=ANprobRateOutput(end-nLSRfibers:end,:);
-            ANdt=dt; % no adjustment for spikes speedup
+            dtSpikes=dt; % no adjustment for spikes speedup
         else
             [nFibers c]=size(ANoutput);
             nLSRfibers=nFibers/length(ANtauCas);
             ANresponse=ANoutput(end-nLSRfibers:end,:);
         end
 
         ANresponse=sum(ANresponse);
-%         ANresponseTimes=ANdt:ANdt:length(ANresponse)*ANdt;
-%         figure(99), plot(ANresponseTimes,ANresponse)
+        %         ANresponseTimes=dtSpikes:dtSpikes:length(ANresponse)*dtSpikes;
+        %         figure(99), plot(ANresponseTimes,ANresponse)
 
         % analyse results
         probeStart=initialSilenceDuration+maskerDuration+gapDuration;
-        PSTHbinWidth=ANdt;
+        PSTHbinWidth=dtSpikes;
         responseDelay=0.005;
         probeTimes=probeStart+responseDelay:...
             PSTHbinWidth:probeStart+probeDuration+responseDelay;
         probeIDX=round(probeTimes/PSTHbinWidth);
         probePSTH=ANresponse(probeIDX);

         if showPSTHs
             nLevels=length(maskerLevels);
             nDurations=length(gapDurations);
             figure(8)
             PSTHbinWidth=1e-3;
-            PSTH=UTIL_PSTHmaker(ANresponse, ANdt, PSTHbinWidth);
-            PSTH=PSTH*ANdt/PSTHbinWidth;
+            PSTH=UTIL_PSTHmaker(ANresponse, dtSpikes, PSTHbinWidth);
+            PSTH=PSTH*dtSpikes/PSTHbinWidth;
             PSTHplotCount=PSTHplotCount+1;
             subplot(nLevels,nDurations,PSTHplotCount)
             PSTHtime=PSTHbinWidth:PSTHbinWidth:...
                 PSTHbinWidth*length(PSTH);
             if strcmp(AN_spikesOrProbability, 'spikes')
-            bar(PSTHtime, PSTH/PSTHbinWidth/nFibers)
-%                 ylim([0 500])
+                bar(PSTHtime, PSTH/PSTHbinWidth/nFibers)
+                %                 ylim([0 500])
             else
-            bar(PSTHtime, PSTH)
+                bar(PSTHtime, PSTH)
                 ylim([0 500])
             end
-%             xlim([0 longestSignalDuration])
+            xlim([0 longestSignalDuration])
             grid on
 
             if PSTHplotCount< (nLevels-1) * nDurations+1
                 set(gca,'xticklabel',[])
             end

 
             if PSTHplotCount<=nDurations
                 title([num2str(1000*gapDurations(PSTHplotCount)) 'ms'])
             end
 
-%         figure(99),            bar(PSTHtime, PSTH)
+            %         figure(99),            bar(PSTHtime, PSTH)
 
         end % showPSTHs
 
     end     % gapDurations duration
     summaryFiringRates=[summaryFiringRates allFiringRates'];