comparison multithreshold 1.46/subjGUI_MT.m @ 38:c2204b18f4a2 tip

End nov big change

             title([stimulusParameters.WRVname ' = ' ...
                 num2str(withinRuns.variableValue, '%5.1f')])
         otherwise
             title([stimulusParameters.WRVname ' = ' ...
                 num2str(withinRuns.variableValue, '%5.1f') ...
-                ';    TH= ' num2str(estThresh, '%5.1f')])
+                ';    TH= ' num2str(estThresh, '%5.1f') ' dB'])
     end
 end
 xlim([0 experiment.maxTrials+withinRuns.beginningOfPhase2]);
 ylim(stimulusParameters.WRVlimits)
+ylabel ('dB SPL')
 grid on
 
 % Panel 4: Summary of threshold estimates (not used here)
 % Estimates from previous runs are set in 'runCompleted'
 % It is only necessary to change title showing runs/trials remaining

     end
     hold on
 end
 xlimRM([ min(betweenRuns.variableList1) max(betweenRuns.variableList1) ])
 ylim(stimulusParameters.WRVlimits)
-ylabel('thresholds')
+ylabel('thresholds (dB)')
 xlabel(betweenRuns.variableName1)
 set(gca,'ytick', [0 20 40 60 80 100])
 try
     % problems if only one x value
     set(gca,'XTick', sort(betweenRuns.variableList1))

     if ispc, wavplay(y/100,fs, 'async'), else sound(y/100,fs, 'async'), end
 
     % update experimenter GUI
     addToMsg('Experiment completed.',1)
 
+    set(expGUIhandles.pushbuttonSave,'visible','on')
     printReport
     experiment.status='endOfExperiment';
     return
 else
     % No, hang on.

 % -------------------------------------------------------MAPmodel
 function [modelResponse, MacGregorResponse]=MAPmodel
 
 global experiment stimulusParameters audio withinRuns
 % global outerMiddleEarParams DRNLParams AN_IHCsynapseParams
-global ICoutput ANdt dt savedBFlist ANprobRateOutput expGUIhandles
+global ICoutput dtSpikes dt savedBFlist ANprobRateOutput expGUIhandles
 global paramChanges
 
 savePath=path;
 addpath(['..' filesep 'MAP'], ['..' filesep 'utilities'])
 modelResponse=[];

 if showPlotsAndDetails
     options.printModelParameters=0;
     options.showModelOutput=1;
     options.printFiringRates=1;
     options.showACF=0;
-    options.showEfferent=1;
+    options.showEfferent=0;
     options.surfProbability=0;
     showMapOptions.surfSpikes=0;
     UTIL_showMAP(options)
 end
 

 end
 % ---------------------------------------------------------- end model run
 
 if strcmp(AN_spikesOrProbability,'spikes')
     MacGregorResponse= sum(ICoutput,1);                 % use IC
-    dt=ANdt;
+    dt=dtSpikes;
     time=dt:dt:dt*length(MacGregorResponse);
 else
     % for one channel, ANprobResponse=ANprobRateOutput
     % for multi channel take strongest in any epoch
     nChannels=length(savedBFlist);

 
 % group delay on unit response - these values are iffy
 windowOnsetDelay= 0.004;
 windowOffsetDelay= 0.022; % long ringing time
 
-% now find the response of the MacGregor model during the target presentation + group delay
+% now find the response of the MacGregor model during 
+%  the target presentation + group delay
 switch experiment.threshEstMethod
     case {'2I2AFC++', '2I2AFC+++'}
         idx= time>stimulusParameters.testTargetBegins+windowOnsetDelay ...
             & time<stimulusParameters.testTargetEnds+windowOffsetDelay;
         nSpikesTrueWindow=sum(MacGregorResponse(:,idx));

         end
         disp(['level target dummy decision: ' ...
             num2str([withinRuns.variableValue nSpikesTrueWindow ...
             nSpikesFalseWindow  nSpikesDuringTarget], '%4.0f') ] )
     otherwise
-        % single interval
+        
+        % single interval up/down
         idx=find(time>stimulusParameters.testTargetBegins +windowOnsetDelay...
             & time<stimulusParameters.testTargetEnds+windowOffsetDelay);
         if strcmp(AN_spikesOrProbability,'spikes')
             nSpikesDuringTarget=sum(MacGregorResponse(:,idx));
             timeX=time(idx);

     fclose(serobj);
 catch
 end
 
 try
-    serobj = serial('COM4') ;           	% Creating serial port object now its connected to COM4 !!! button boxes in booths are connected to COM2
+    % !!! button boxes in booths are connected to COM2. User beware of
+    % connection port.
+    serobj = serial('COM2') ;           	% Creating serial port object
     serobj.Baudrate = 9600;           		% Set the baud rate at the specific value
     set(serobj, 'Parity', 'none') ;     	% Set parity as none
     set(serobj, 'Databits', 8) ;          	% set the number of data bits
     set(serobj, 'StopBits', 1) ;         	% set number of stop bits as 1
     set(serobj, 'Terminator', 'CR') ; 		% set the terminator value to carriage return