rmeddis@28: function errormsg=nextStimulus(handles)
rmeddis@28: % Handles everything concerned with the stimulus presentation
rmeddis@28: %  called from startNewRun in subjGUI
rmeddis@28: 
rmeddis@28: global experiment stimulusParameters withinRuns betweenRuns
rmeddis@28: experiment.status='presentingStimulus';
rmeddis@28: errormsg='';
rmeddis@28: 
rmeddis@28: % interrupt by 'stop' button
rmeddis@28: if experiment.stop
rmeddis@28:     disp('******** experiment manually stopped  *****************')
rmeddis@28:     experiment.status= 'waitingForStart';
rmeddis@28:     addToMsg('manually stopped',1)
rmeddis@28:     return
rmeddis@28: end
rmeddis@28: 
rmeddis@28: % -----------------------------------------choose catch trials at random
rmeddis@28: % catch trials are for subject threshold measurements only
rmeddis@28: % this is the only place where withinRuns.catchTrial is set
rmeddis@28: if experiment.allowCatchTrials
rmeddis@28:     if withinRuns.trialNumber==1;
rmeddis@28:         % first trial is never a catch trial
rmeddis@28:         withinRuns.catchTrial=0;
rmeddis@28:         withinRuns.catchTrialCount=0;   % reset count on first trial
rmeddis@28:     elseif withinRuns.trialNumber==2 ...
rmeddis@28:             && withinRuns.catchTrialCount==0
rmeddis@28:         % second trial is always a catch trial
rmeddis@28:         withinRuns.catchTrial=1;
rmeddis@28:         withinRuns.catchTrialCount=1;   % this must be the first
rmeddis@28:     elseif withinRuns.thisIsRepeatTrial
rmeddis@28:         % for requested repeats do not change catch trial status
rmeddis@28:         withinRuns.thisIsRepeatTrial=0; % reset toggle
rmeddis@28:     else
rmeddis@28:         % choose whether or not to have a catch trial
rmeddis@28:         R=rand;
rmeddis@28:         if R<stimulusParameters.catchTrialRate
rmeddis@28:             % catch trial
rmeddis@28:             withinRuns.catchTrial=1;
rmeddis@28:             addToMsg('Catch Trial',1)
rmeddis@28:             withinRuns.catchTrialCount=withinRuns.catchTrialCount+1;
rmeddis@28:         else
rmeddis@28:             % not a catch trial
rmeddis@28:             withinRuns.catchTrial=0;
rmeddis@28:         end
rmeddis@28:     end
rmeddis@28: else
rmeddis@28:     % no catch trials for statistical evaluations or 2AIFC or (poss) MAP
rmeddis@28:     withinRuns.catchTrial=0;
rmeddis@28: end
rmeddis@28: 
rmeddis@28: %------------ during stimulus presentation show appropriate button images
rmeddis@28: switch experiment.ear
rmeddis@28:     case {'statsModelLogistic', 'statsModelRareEvent',...
rmeddis@28:             'MAPmodel', 'MAPmodelMultiCh', 'MAPmodelSingleCh'}
rmeddis@28:         % no buttons shown
rmeddis@28:     otherwise
rmeddis@28:         switch experiment.threshEstMethod
rmeddis@28:             case {'2I2AFC++', '2I2AFC+++'}
rmeddis@28:                 %Except for 2I2AFC
rmeddis@28:                 % For 2I2AFC the buttons on the screen ab initio
rmeddis@28:                 set(handles.frame1,'visible','off')
rmeddis@28:                 set(handles.pushbuttonGO,'visible','off')
rmeddis@28:                 set(handles.pushbuttoNotSure,'visible','off')
rmeddis@28:                 set(handles.pushbuttonWrongButton,'visible','off')
rmeddis@28:                 set(handles.pushbutton3,'visible','off')
rmeddis@28:                 set(handles.pushbutton0,'visible','off')
rmeddis@28:                 set(handles.pushbutton1,'visible','on')
rmeddis@28:                 set(handles.pushbutton2,'visible','on')
rmeddis@28:                 drawnow
rmeddis@28:             otherwise
rmeddis@28:                 % i.e. single interval/ maxLikelihood
rmeddis@28:                 set(handles.frame1,'backgroundColor','w')
rmeddis@28:                 set(handles.frame1,'visible','off')
rmeddis@28:                 set(handles.pushbuttoNotSure,'visible','off')
rmeddis@28:                 set(handles.pushbuttonWrongButton,'visible','off')
rmeddis@28:                 set(handles.pushbutton3,'visible','off')
rmeddis@28:                 set(handles.pushbutton2,'visible','off')
rmeddis@28:                 set(handles.pushbutton1,'visible','off')
rmeddis@28:                 set(handles.pushbutton0,'visible','off')
rmeddis@28:                 pause(.1)
rmeddis@28:         end
rmeddis@28: end
rmeddis@28: 
rmeddis@28: set(handles.textMSG,'BackgroundColor','w', 'ForegroundColor', 'b')
rmeddis@28: 
rmeddis@28: % Now the serious business of crafting and presenting the stimulus
rmeddis@28: errormsg= stimulusMakeAndPlay (handles);
rmeddis@28: 
rmeddis@28: if ~isempty(errormsg)
rmeddis@28:     % e.g. clipping. subjGUI will service the error
rmeddis@28:     return
rmeddis@28: end
rmeddis@28: 
rmeddis@28: % after playing the stimulus, reset the subjectGUI
rmeddis@28: switch experiment.ear
rmeddis@28:     case {'statsModelLogistic', 'statsModelRareEvent',...
rmeddis@28:             'MAPmodel', 'MAPmodelMultiCh', 'MAPmodelSingleCh'}
rmeddis@28:         % no changes required if model used
rmeddis@28:         % NB these changes do occur is 'MAPmodelListen' is selected
rmeddis@28:     otherwise
rmeddis@28:         switch experiment.threshEstMethod
rmeddis@28:             case {'2I2AFC++', '2I2AFC+++'}
rmeddis@28:                 % buttons already visible
rmeddis@28:             otherwise
rmeddis@28:                 % single interval now make buttons visible
rmeddis@28:                 set(handles.frame1,'visible','on')
rmeddis@28:                 set(handles.pushbuttoNotSure,'visible','on')
rmeddis@28:                 % set(handles.pushbuttonWrongButton,'visible','on')
rmeddis@28:                 set(handles.pushbutton0,'visible','on')
rmeddis@28:                 set(handles.pushbutton1,'visible','on')
rmeddis@28:                 set(handles.pushbutton2,'visible','on')
rmeddis@28:                 set(handles.pushbutton3,'visible','on')
rmeddis@28:         end
rmeddis@28: end
rmeddis@28: 
rmeddis@28: switch experiment.paradigm
rmeddis@28:     case 'SRT'
rmeddis@28:         set(handles.frame1,'backgroundColor','w')
rmeddis@28:         set(handles.frame1,'visible','off')
rmeddis@28:         set(handles.pushbuttoNotSure,'visible','off')
rmeddis@28:         set(handles.pushbuttonWrongButton,'visible','off')
rmeddis@28:         set(handles.pushbutton3,'visible','off')
rmeddis@28:         set(handles.pushbutton2,'visible','off')
rmeddis@28:         set(handles.pushbutton1,'visible','off')
rmeddis@28:         set(handles.pushbutton0,'visible','off')
rmeddis@28:         set(handles.editdigitInput,'visible','on')
rmeddis@28:         set(handles.editdigitInput,'string',[])
rmeddis@28:         pause(.2)
rmeddis@28:         uicontrol(handles.editdigitInput)
rmeddis@28:         
rmeddis@28:     otherwise
rmeddis@28:         set(handles.editdigitInput,'visible','off')
rmeddis@28: end
rmeddis@28: 
rmeddis@28: 
rmeddis@28: experiment.status='waitingForResponse';
rmeddis@28: % home again
rmeddis@28: 
rmeddis@28: % ------------------------------------------------------------------------------------------stimulusMakeAndPlay
rmeddis@28: function errormsg=stimulusMakeAndPlay (handles)
rmeddis@28: global experiment stimulusParameters betweenRuns withinRuns expGUIhandles  audio
rmeddis@28: % creates the stimulus and plays it; there are two stimuli; cue and test
rmeddis@28: %  called from nextStimulus
rmeddis@28: 
rmeddis@28: errormsg='';
rmeddis@28: 
rmeddis@28: % first post the subjects instructions on subjGUI
rmeddis@28: set(handles.textMSG,'string', stimulusParameters.subjectText)
rmeddis@28: 
rmeddis@28: % select the new levels of the between runs variables
rmeddis@28: num=betweenRuns.runNumber;
rmeddis@28: cmd=(['stimulusParameters.' betweenRuns.variableName1 '= ' ...
rmeddis@28:     num2str(betweenRuns.var1Sequence(num)) ';']);
rmeddis@28: % e.g.  stimulusParameters.targetFrequency= 1000;
rmeddis@28: eval(cmd);
rmeddis@28: 
rmeddis@28: cmd=(['stimulusParameters.' betweenRuns.variableName2 '= ' ...
rmeddis@28:     num2str(betweenRuns.var2Sequence(num)) ';']);
rmeddis@28: % e.g.  stimulusParameters.targetDuration= 0.1;
rmeddis@28: eval(cmd);
rmeddis@28: 
rmeddis@28: switch experiment.paradigm
rmeddis@28:     % target level may vary between runs
rmeddis@28:     case {'trainingIFMC', 'TMC','TMC_16ms', 'TMC - ELP', 'IFMC','IFMC_8ms','IFMC_16ms'}
rmeddis@28:         idx=floor(num/length(betweenRuns.variableList1)-0.01)+1;
rmeddis@28:         cmd=(['stimulusParameters.targetLevel = ' ...
rmeddis@28:             num2str(stimulusParameters.targetLevels(idx)) ';']);
rmeddis@28:         eval(cmd);
rmeddis@28:         if withinRuns.trialNumber==1
rmeddis@28:             disp(['targetLevel=' num2str(stimulusParameters.targetLevel)])
rmeddis@28:         end
rmeddis@28: end
rmeddis@28: 
rmeddis@28: 
rmeddis@28: % for more readable code use shorter variable names;
rmeddis@28: % NB these may change below; these are only the starting values
rmeddis@28: 
rmeddis@28: targetType=        stimulusParameters.targetType;
rmeddis@28: targetDuration=    stimulusParameters.targetDuration;
rmeddis@28: targetLevel=       stimulusParameters.targetLevel;
rmeddis@28: targetFrequency=   stimulusParameters.targetFrequency;
rmeddis@28: 
rmeddis@28: maskerType=        stimulusParameters.maskerType;
rmeddis@28: maskerDuration=    stimulusParameters.maskerDuration;
rmeddis@28: maskerLevel=       stimulusParameters.maskerLevel;
rmeddis@28: maskerRelativeFrequency=  stimulusParameters.maskerRelativeFrequency;
rmeddis@28: maskerFrequency=   maskerRelativeFrequency*targetFrequency;
rmeddis@28: 
rmeddis@28: gapDuration=       stimulusParameters.gapDuration;
rmeddis@28: 
rmeddis@28: rampDuration=      stimulusParameters.rampDuration;
rmeddis@28: AFCsilenceDuration=stimulusParameters.AFCsilenceDuration; % 2I2AFC gap
rmeddis@28: backgroundLevel=   stimulusParameters.backgroundLevel;
rmeddis@28: 
rmeddis@28: % Set level of within runs variable
rmeddis@28: % this is the first change to one of the values shown above
rmeddis@28: cmd=[stimulusParameters.WRVname '= withinRuns.variableValue;' ];
rmeddis@28: % e.g.: maskerLevel= withinRuns.variableValue;
rmeddis@28: eval(cmd);
rmeddis@28: 
rmeddis@28: % cue and test stimuli are identical except for a single difference
rmeddis@28: % depending on the paradigm
rmeddis@28: cueTestDifference= stimulusParameters.cueTestDifference;
rmeddis@28: %  cue characteristics before adding cue differences
rmeddis@28: cueTargetLevel=targetLevel;
rmeddis@28: cueMaskerFrequency=maskerFrequency;
rmeddis@28: cueMaskerDuration=maskerDuration;
rmeddis@28: cueMaskerLevel=maskerLevel;
rmeddis@28: cueTargetFrequency=targetFrequency;
rmeddis@28: cueGapDuration=gapDuration;
rmeddis@28: 
rmeddis@28: % ----------------------------paradigm sensitive cue and masker settings
rmeddis@28: % switch off unwanted components and base cue on target values
rmeddis@28: % for catch trials switch off the target
rmeddis@28: switch experiment.paradigm
rmeddis@28:     % OHIO is a temporary special arrangement
rmeddis@28:     case{'OHIOrand','OHIOspect','OHIOtemp','OHIOspectemp','OHIOabs'}
rmeddis@28:         % these values must be set in MAPparamsOHIO
rmeddis@28:         targetFrequency=experiment.OHIOfrequencies;
rmeddis@28:         numOHIOtones=stimulusParameters.numOHIOtones;
rmeddis@28:         maskerType='OHIO';
rmeddis@28:         targetType='OHIO';
rmeddis@28:         
rmeddis@28:         % globalStimParams.beginSilences says when each tone begins
rmeddis@28:         % targetFrequency specifies the frequency of each tone
rmeddis@28:         % targetLevel is the level of each tone
rmeddis@28:         switch experiment.paradigm
rmeddis@28:             case 'OHIOabs'
rmeddis@28:                 % only one tone
rmeddis@28:                 targetFrequency=targetFrequency(numOHIOtones);
rmeddis@28:                 globalStimParams.beginSilences=0.01;
rmeddis@28:                 targetDuration=0.01;
rmeddis@28:             case 'OHIOrand'
rmeddis@28:                 % select random frequencies from the list
rmeddis@28:                 x=rand(12,1); [x idx]=sort(x);
rmeddis@28:                 targetFrequency=targetFrequency(idx(1:numOHIOtones));
rmeddis@28:                 thresholds=experiment.OHIOthresholds;
rmeddis@28:                 targetLevel=thresholds(idx(1:numOHIOtones))+targetLevel;
rmeddis@28:                 globalStimParams.beginSilences=0.01:0.02:...
rmeddis@28:                     0.01+0.02*(numOHIOtones-1);
rmeddis@28:                 targetDuration=numOHIOtones*0.02;
rmeddis@28:             case 'OHIOspect'
rmeddis@28:                 % only one tone with multiple frequencies
rmeddis@28:                 targetFrequency=fliplr(targetFrequency);
rmeddis@28:                 targetFrequency=targetFrequency(1:numOHIOtones);
rmeddis@28:                 thresholds=experiment.OHIOthresholds;
rmeddis@28:                 thresholds=fliplr(thresholds);
rmeddis@28:                 targetLevel=thresholds(1:numOHIOtones)+targetLevel;
rmeddis@28:                 globalStimParams.beginSilences=...
rmeddis@28:                     repmat(0.01, 1, numOHIOtones);
rmeddis@28:                 targetDuration=0.02;
rmeddis@28:             case 'OHIOspectemp'
rmeddis@28:                 % only one tone with multiple frequencies
rmeddis@28:                 targetFrequency=fliplr(targetFrequency);
rmeddis@28:                 targetFrequency=targetFrequency(1:numOHIOtones);
rmeddis@28:                 thresholds=experiment.OHIOthresholds;
rmeddis@28:                 thresholds=fliplr(thresholds);
rmeddis@28:                 targetLevel=thresholds(1:numOHIOtones)+targetLevel;
rmeddis@28:                 globalStimParams.beginSilences=0.01:0.02:...
rmeddis@28:                     0.01+0.02*(numOHIOtones-1);
rmeddis@28:                 targetDuration=numOHIOtones*0.02;
rmeddis@28:             case 'OHIOtemp'
rmeddis@28:                 % use only one tone repeatedly
rmeddis@28:                 tonesToUse=10;
rmeddis@28:                 targetFrequency=targetFrequency(tonesToUse);
rmeddis@28:                 targetFrequency=repmat(targetFrequency,1,numOHIOtones);
rmeddis@28:                 thresholds=experiment.OHIOthresholds;
rmeddis@28:                 thresholds=thresholds(tonesToUse);
rmeddis@28:                 targetLevel=repmat(thresholds,1,numOHIOtones)+targetLevel;
rmeddis@28:                 globalStimParams.beginSilences=0.01:0.02:...
rmeddis@28:                     0.01+0.02*(numOHIOtones-1);
rmeddis@28:                 targetDuration=numOHIOtones*0.02;
rmeddis@28:         end
rmeddis@28:         % still in OHIO
rmeddis@28:         % Dummy values to make things work although no masker or cue used
rmeddis@28:         % target values have changed and this affects the cue values
rmeddis@28:         cueMaskerLevel=targetLevel;
rmeddis@28:         cueTargetLevel=targetLevel;
rmeddis@28:         cueTargetFrequency=targetFrequency;
rmeddis@28:         cueMaskerFrequency=targetFrequency;
rmeddis@28:         maskerFrequency=targetFrequency;
rmeddis@28:         maskerLevel=targetLevel;
rmeddis@28:         disp(['OHIO frequencies= ' num2str(targetFrequency)])
rmeddis@28: end
rmeddis@28: 
rmeddis@28: % --- set cueTarget level according to assessment method
rmeddis@28: % cue-test difference applies only with singleInterval
rmeddis@28: switch experiment.threshEstMethod
rmeddis@28:     case {'2I2AFC++', '2I2AFC+++'}
rmeddis@28:         % For 2IFC the cue stimulus (masker + probe) is the 'no' window
rmeddis@28:         % and the target stimulus (masker+probe) is the 'yes' window
rmeddis@28:         % the order of presentation is decided at the last minute.
rmeddis@28:         cueTargetLevel=-100;    % the target is never in the 'no' window
rmeddis@28:         cueMaskerLevel=maskerLevel; % masker level is the same in both
rmeddis@28:     otherwise
rmeddis@28:         % 'single interval' or max likelihood
rmeddis@28:         switch experiment.paradigm
rmeddis@28:             % cue target is more audible
rmeddis@28:             case {'training','absThreshold', 'absThreshold_8',  ...
rmeddis@28:                     'TENtest', 'threshold_duration','discomfort',...
rmeddis@28:                     'overShoot','overShootB','overShootMB1', ...
rmeddis@28:                     'overShootMB2', 'OHIO','OHIOabs','OHIOspect'}
rmeddis@28:                 cueTargetLevel=targetLevel+cueTestDifference;
rmeddis@28:                 
rmeddis@28:             case {'forwardMasking','forwardMaskingD','trainingIFMC', ...
rmeddis@28:                     'TMC','TMC_16ms', 'TMC - ELP', 'IFMC','IFMC_8ms', 'FMreProbe'}
rmeddis@28:                 % cue masker is weaker to make target more audible
rmeddis@28:                 cueMaskerLevel=maskerLevel-cueTestDifference;
rmeddis@28:         end
rmeddis@28: end
rmeddis@28: 
rmeddis@28: % ----------------------------- catch trial
rmeddis@28: if withinRuns.catchTrial
rmeddis@28:     targetLevel=-100;	% no target
rmeddis@28: end
rmeddis@28: 
rmeddis@28: % ----------------------------- calibration of sound output
rmeddis@28: calibrationCorrectiondB=stimulusParameters.calibrationdB;
rmeddis@28: if calibrationCorrectiondB<-50
rmeddis@28:     if maskerFrequency==targetFrequency
rmeddis@28:         load 'calibrationFile'  % calibrationFrequency calibrationAttenutation
rmeddis@28:         idx=find(calibrationFrequency==targetFrequency);
rmeddis@28:         if isempty(idx)
rmeddis@28:             error('Calibration bty file; frequency not found')
rmeddis@28:         else
rmeddis@28:             calibrationCorrectiondB=calibrationAttenutation(idx)
rmeddis@28:         end
rmeddis@28:     else
rmeddis@28:         error('calibration by file requested but masker frequency is not the same as target')
rmeddis@28:     end
rmeddis@28: end
rmeddis@28: 
rmeddis@28: 
rmeddis@28: % -------------------------------------- Checks on excessive signal level
rmeddis@28: 
rmeddis@28: % clipping is relevant only for soundcard use (not modelling)
rmeddis@28: switch experiment.ear
rmeddis@28:     case {'left', 'right', 'diotic',...
rmeddis@28:             'dichotic', 'dioticLeft', 'dichoticRight'}
rmeddis@28:         experiment.headphonesUsed=1;
rmeddis@28:     otherwise
rmeddis@28:         experiment.headphonesUsed=0;
rmeddis@28: end
rmeddis@28: 
rmeddis@28: % NB calibration *reduces* the level of the soundCard output
rmeddis@28: switch experiment.ear
rmeddis@28:     case {'left', 'right', 'diotic',...
rmeddis@28:             'dichotic', 'dioticLeft', 'dichoticRight'}
rmeddis@28:         clippingLevel=91+calibrationCorrectiondB;
rmeddis@28:         soundCardMinimum=clippingLevel-20*log10(2^24);
rmeddis@28:     otherwise
rmeddis@28:         clippingLevel=inf;
rmeddis@28:         soundCardMinimum=-inf;
rmeddis@28: end
rmeddis@28: 
rmeddis@28: % Check for extreme WRV values and abort if necessary
rmeddis@28: % WRVname specifies the value that changes from trial to trial
rmeddis@28: withinRuns.forceThreshold=[];
rmeddis@28: switch stimulusParameters.WRVname
rmeddis@28:     % check for extreme values. Note that one of the tones might be switched off
rmeddis@28:     case 'maskerLevel'
rmeddis@28:         upperLevel=stimulusParameters.WRVlimits(2);
rmeddis@28:         lowerLevel=stimulusParameters.WRVlimits(1);
rmeddis@28:         if max(maskerLevel, cueMaskerLevel)> upperLevel
rmeddis@28:             errormsg=['Level(' num2str(max(maskerLevel,cueMaskerLevel)) ...
rmeddis@28:                 ') is too high ***'];
rmeddis@28:             withinRuns.forceThreshold=upperLevel;
rmeddis@28:             withinRuns.forceThreshold=NaN;
rmeddis@28:             return
rmeddis@28:         end
rmeddis@28:         if max(maskerLevel, cueMaskerLevel)< lowerLevel
rmeddis@28:             errormsg=['Level(' num2str(max(maskerLevel,cueMaskerLevel)) ...
rmeddis@28:                 ') is too low ***'];
rmeddis@28:             withinRuns.forceThreshold=lowerLevel;
rmeddis@28:             withinRuns.forceThreshold=NaN;
rmeddis@28:             return
rmeddis@28:         end
rmeddis@28:         
rmeddis@28:         if max(maskerLevel, cueMaskerLevel)> clippingLevel
rmeddis@28:             errormsg=['Level(' num2str(max(maskerLevel,cueMaskerLevel)) ...
rmeddis@28:                 ') is clipping ***'];
rmeddis@28:             withinRuns.forceThreshold=clippingLevel;
rmeddis@28:             withinRuns.forceThreshold=NaN;
rmeddis@28:             return
rmeddis@28:         end
rmeddis@28:         
rmeddis@28:     case 'targetLevel'
rmeddis@28:         upperLevel=stimulusParameters.WRVlimits(2);
rmeddis@28:         lowerLevel=stimulusParameters.WRVlimits(1);
rmeddis@28:         if ~withinRuns.catchTrial
rmeddis@28:             if max(targetLevel, cueTargetLevel)> upperLevel
rmeddis@28:                 errormsg=['target level (' ...
rmeddis@28:                     num2str(max(targetLevel, cueTargetLevel)) ...
rmeddis@28:                     ')  is too high ***'];
rmeddis@28:                 withinRuns.forceThreshold=upperLevel;
rmeddis@28:                 withinRuns.forceThreshold=NaN;
rmeddis@28:                 return
rmeddis@28:             end
rmeddis@28:             if max(targetLevel, cueTargetLevel)< lowerLevel
rmeddis@28:                 errormsg=['target level (' ...
rmeddis@28:                     num2str(max(targetLevel, cueTargetLevel)) ...
rmeddis@28:                     ')  is too low ***'];
rmeddis@28:                 withinRuns.forceThreshold=lowerLevel;
rmeddis@28:                 withinRuns.forceThreshold=NaN;
rmeddis@28:                 return
rmeddis@28:             end
rmeddis@28:             if max(targetLevel, cueTargetLevel)> clippingLevel
rmeddis@28:                 errormsg=['target level (' ...
rmeddis@28:                     num2str(max(targetLevel, cueTargetLevel)) ...
rmeddis@28:                     ')  is clipping ***'];
rmeddis@28:                 withinRuns.forceThreshold=upperLevel;
rmeddis@28:                 withinRuns.forceThreshold=NaN;
rmeddis@28:                 return
rmeddis@28:             end
rmeddis@28:         end
rmeddis@28:     case 'maskerDuration'
rmeddis@28:         % this is odd! but harmless
rmeddis@28:         if max(maskerDuration, cueMaskerDuration)> ...
rmeddis@28:                 stimulusParameters.WRVlimits(2)
rmeddis@28:             errormsg=['maskerDuration (' ...
rmeddis@28:                 num2str(max(maskerDuration, cueMaskerDuration))...
rmeddis@28:                 ') is too long ***'];
rmeddis@28:             withinRuns.forceThreshold=stimulusParameters.WRVlimits(2);
rmeddis@28:             withinRuns.forceThreshold=NaN;
rmeddis@28:             return
rmeddis@28:         end
rmeddis@28:         
rmeddis@28:         if min(maskerDuration, cueMaskerDuration)...
rmeddis@28:                 < stimulusParameters.WRVlimits(1)
rmeddis@28:             errormsg=['maskerDuration (' num2str(maskerLevel) ...
rmeddis@28:                 ') too short ***'];
rmeddis@28:             withinRuns.forceThreshold=stimulusParameters.WRVlimits(1);
rmeddis@28:             withinRuns.forceThreshold=NaN;
rmeddis@28:             return
rmeddis@28:         end
rmeddis@28:         
rmeddis@28:         % legacy programming
rmeddis@28:     case 'gapDuration'
rmeddis@28:         if gapDuration<0
rmeddis@28:             errormsg=['gapDuration (' num2str(gapDuration) ...
rmeddis@28:                 ') is less than zero  ***'];
rmeddis@28:             return
rmeddis@28:         end
rmeddis@28:         
rmeddis@28:     case 'maskerFrequency'
rmeddis@28:         switch experiment.paradigm
rmeddis@28:             case 'bandwidth'
rmeddis@28:                 frequency=maskerFrequency';
rmeddis@28:                 if stimulusParameters.WRVstep<0
rmeddis@28:                     lowerLevel=stimulusParameters.targetFrequency;
rmeddis@28:                     upperLevel=stimulusParameters.targetFrequency*2;
rmeddis@28:                 else
rmeddis@28:                     lowerLevel=stimulusParameters.targetFrequency/3;
rmeddis@28:                     upperLevel=stimulusParameters.targetFrequency;
rmeddis@28:                 end
rmeddis@28:                 
rmeddis@28:                 if frequency(1)>upperLevel || frequency(1)<lowerLevel
rmeddis@28:                     errormsg=['frequency out of range: ' ...
rmeddis@28:                         num2str(frequency)];
rmeddis@28:                     withinRuns.forceThreshold=frequency;
rmeddis@28:                     return
rmeddis@28:                 end
rmeddis@28:             otherwise
rmeddis@28:         end
rmeddis@28:         
rmeddis@28:     case 'maskerRelativeFrequency'
rmeddis@28:         if  maskerRelativeFrequency<stimulusParameters.WRVlimits(1)
rmeddis@28:             errormsg=['masker frequency (' ...
rmeddis@28:                 num2str(frequencyDifference) ...
rmeddis@28:                 ')  is outside WRV limits ***'];
rmeddis@28:             withinRuns.forceThreshold=stimulusParameters.WRVlimits(1) ;
rmeddis@28:             return
rmeddis@28:         end
rmeddis@28:         if   maskerRelativeFrequency>stimulusParameters.WRVlimits(2)
rmeddis@28:             errormsg=['masker frequency (' ...
rmeddis@28:                 num2str(frequencyDifference) ...
rmeddis@28:                 ')  is outside WRV limits ***'];
rmeddis@28:             withinRuns.forceThreshold=stimulusParameters.WRVlimits(2) ;
rmeddis@28:             return
rmeddis@28:         end
rmeddis@28:         
rmeddis@28: end
rmeddis@28: 
rmeddis@28: % --------------------------------Ear ----------------------------------
rmeddis@28: globalStimParams.ears='specified';
rmeddis@28: % ear: 1=left, 2=right
rmeddis@28: switch experiment.ear
rmeddis@28:     case 'left'
rmeddis@28:         maskerEar=1;
rmeddis@28:         targetEar=1;
rmeddis@28:     case 'right'
rmeddis@28:         maskerEar=2;
rmeddis@28:         targetEar=2;
rmeddis@28:     case 'dichoticLeft'
rmeddis@28:         maskerEar=2;
rmeddis@28:         targetEar=1;
rmeddis@28:     case 'dichoticRight'
rmeddis@28:         maskerEar=1;
rmeddis@28:         targetEar=2;
rmeddis@28:     case 'diotic'
rmeddis@28:         maskerEar=1;
rmeddis@28:         targetEar=1;
rmeddis@28:         globalStimParams.ears='diotic';
rmeddis@28:     case {'MAPmodel', 'MAPmodelMultiCh', 'MAPmodelSingleCh', 'MAPmodelListen',...
rmeddis@28:             'statsModelLogistic', 'statsModelRareEvent'}
rmeddis@28:         maskerEar=1;
rmeddis@28:         targetEar=1;
rmeddis@28: end
rmeddis@28: 
rmeddis@28: backgroundType=stimulusParameters.backgroundType;
rmeddis@28: switch stimulusParameters.backgroundType
rmeddis@28:     case {'noiseDich', 'pinkNoiseDich'}
rmeddis@28:         %     case 'Dich'
rmeddis@28:         % dich means put the background in the ear opposite to the target
rmeddis@28:         backgroundType=backgroundType(1:end-4);
rmeddis@28:         switch targetEar
rmeddis@28:             case 1
rmeddis@28:                 backgroundEar=2;
rmeddis@28:             case 2
rmeddis@28:                 backgroundEar=1;
rmeddis@28:         end
rmeddis@28:     otherwise
rmeddis@28:         %             case {'none','noise', 'pinkNoise', 'TEN','babble'}
rmeddis@28:         backgroundEar=targetEar;
rmeddis@28: end
rmeddis@28: 
rmeddis@28: % ------------------------------- Make Stimulus -------------------
rmeddis@28: % single interval up/down plays cue then target stimulus
rmeddis@28: % 2IFC uses cue stimulus as interval with no target
rmeddis@28: globalStimParams.FS=stimulusParameters.sampleRate;
rmeddis@28: dt=1/stimulusParameters.sampleRate;
rmeddis@28: globalStimParams.dt=dt;
rmeddis@28: stimulusParameters.dt=dt; % for use later
rmeddis@28: 
rmeddis@28: 
rmeddis@28: 
rmeddis@28: globalStimParams.audioOutCorrection=10^(calibrationCorrectiondB/20);
rmeddis@28: % the output will be reduced by this amount in stimulusCreate
rmeddis@28: % i.e.  audio=audio/globalStimParams.audioOutCorrection
rmeddis@28: % A 91 dB level will yield a peak amp of 1 for calibration=0
rmeddis@28: % A 91 dB level will yield a peak amp of 0.4467 for calibration=7
rmeddis@28: % A 98 dB level will yield a peak amp of 1 for calibration=7
rmeddis@28: 
rmeddis@28: precedingSilence=stimulusParameters.stimulusDelay;
rmeddis@28: % all stimuli have 20 ms terminal silence.
rmeddis@28: % this is clearance for modelling late-ringing targets
rmeddis@28: terminalSilence=.03;
rmeddis@28: 
rmeddis@28: % Now compute overall duration of the stimulus
rmeddis@28: % note that all endsilence values are set to -1
rmeddis@28: %  so that they will fill with terminal silence as required to make
rmeddis@28: %  components equal in length
rmeddis@28: % We need to find the longest possible duration
rmeddis@28: duration(1)=precedingSilence+maskerDuration+cueGapDuration...
rmeddis@28:     +targetDuration+terminalSilence;
rmeddis@28: duration(2)=precedingSilence+maskerDuration+gapDuration...
rmeddis@28:     +targetDuration+ terminalSilence;
rmeddis@28: % If the gap is negative we need to ignore it when estimating total length
rmeddis@28: duration(3)=precedingSilence+maskerDuration+ terminalSilence;
rmeddis@28: globalStimParams.overallDuration=max(duration);
rmeddis@28: globalStimParams.nSignalPoints=...
rmeddis@28:     round(globalStimParams.overallDuration/globalStimParams.dt);
rmeddis@28: 
rmeddis@28: %           ----------------------------------------------cue stimulus
rmeddis@28: % cue masker
rmeddis@28: componentNo=1;
rmeddis@28: precedingSilence=stimulusParameters.stimulusDelay;
rmeddis@28: stimComponents(maskerEar,componentNo).type=maskerType;
rmeddis@28: stimComponents(maskerEar,componentNo).toneDuration=cueMaskerDuration;
rmeddis@28: stimComponents(maskerEar,componentNo).frequencies=cueMaskerFrequency;
rmeddis@28: stimComponents(maskerEar,componentNo).amplitudesdB=cueMaskerLevel;
rmeddis@28: stimComponents(maskerEar,componentNo).beginSilence=precedingSilence;
rmeddis@28: stimComponents(maskerEar,componentNo).endSilence=-1;
rmeddis@28: stimComponents(maskerEar,componentNo).AMfrequency=0;
rmeddis@28: stimComponents(maskerEar,componentNo).AMdepth=0;
rmeddis@28: if rampDuration<maskerDuration
rmeddis@28:     % ramps must be shorter than the signal
rmeddis@28:     stimComponents(maskerEar,componentNo).rampOnDur=rampDuration;
rmeddis@28:     stimComponents(maskerEar,componentNo).rampOffDur=rampDuration;
rmeddis@28: else
rmeddis@28:     % or squeeze the ramp in
rmeddis@28:     stimComponents(maskerEar,componentNo).rampOnDur=maskerDuration/2;
rmeddis@28:     stimComponents(maskerEar,componentNo).rampOffDur=maskerDuration/2;
rmeddis@28: end
rmeddis@28: stimComponents(maskerEar,componentNo).phases=...
rmeddis@28:     stimulusParameters.maskerPhase;
rmeddis@28: stimComponents(maskerEar,componentNo).niterations=0;    % for IRN only
rmeddis@28: % stimComponents(targetEar,componentNo)
rmeddis@28: 
rmeddis@28: % cue target
rmeddis@28: componentNo=2;
rmeddis@28: precedingSilence=precedingSilence + maskerDuration+cueGapDuration;
rmeddis@28: stimComponents(targetEar,componentNo).type=targetType;
rmeddis@28: stimComponents(targetEar,componentNo).toneDuration=targetDuration;
rmeddis@28: stimComponents(targetEar,componentNo).frequencies=cueTargetFrequency;
rmeddis@28: stimComponents(targetEar,componentNo).amplitudesdB=cueTargetLevel;
rmeddis@28: stimComponents(targetEar,componentNo).beginSilence=precedingSilence;
rmeddis@28: stimComponents(targetEar,componentNo).endSilence=-1;
rmeddis@28: stimComponents(targetEar,componentNo).AMfrequency=0;
rmeddis@28: stimComponents(targetEar,componentNo).AMdepth=0;
rmeddis@28: if rampDuration<targetDuration
rmeddis@28:     % ramps must be shorter than the signal
rmeddis@28:     stimComponents(targetEar,componentNo).rampOnDur=rampDuration;
rmeddis@28:     stimComponents(targetEar,componentNo).rampOffDur=rampDuration;
rmeddis@28: else
rmeddis@28:     stimComponents(targetEar,componentNo).rampOnDur=0;
rmeddis@28:     stimComponents(targetEar,componentNo).rampOffDur=0;
rmeddis@28: end
rmeddis@28: stimComponents(targetEar,componentNo).phases=...
rmeddis@28:     stimulusParameters.targetPhase;
rmeddis@28: % stimComponents(targetEar,componentNo)
rmeddis@28: 
rmeddis@28: % background same ear as target
rmeddis@28: componentNo=3;
rmeddis@28: stimComponents(backgroundEar,componentNo).type=backgroundType;
rmeddis@28: switch backgroundType
rmeddis@28:     case 'TEN'
rmeddis@28:         fileName=['..' filesep '..' filesep ...
rmeddis@28:             'multithresholdResources' filesep ...
rmeddis@28:             'backgrounds and maskers'...
rmeddis@28:             filesep 'ten.wav'];
rmeddis@28:         [tenNoise, FS]=wavread(fileName);
rmeddis@28:         tenNoise=resample(tenNoise, globalStimParams.FS, FS);
rmeddis@28:         stimComponents(backgroundEar,componentNo).type='file';
rmeddis@28:         stimComponents(backgroundEar,componentNo).stimulus=tenNoise';
rmeddis@28: end
rmeddis@28: stimComponents(backgroundEar,componentNo).toneDuration=...
rmeddis@28:     globalStimParams.overallDuration;
rmeddis@28: stimComponents(backgroundEar,componentNo).amplitudesdB=backgroundLevel;
rmeddis@28: stimComponents(backgroundEar,componentNo).beginSilence=0;
rmeddis@28: stimComponents(backgroundEar,componentNo).endSilence=-1;
rmeddis@28: stimComponents(backgroundEar,componentNo).AMfrequency=0;
rmeddis@28: stimComponents(backgroundEar,componentNo).AMdepth=0;
rmeddis@28: stimComponents(backgroundEar,componentNo).rampOnDur=rampDuration;
rmeddis@28: stimComponents(backgroundEar,componentNo).rampOffDur=rampDuration;
rmeddis@28: 
rmeddis@28: [cueStimulus, errormsg]=...
rmeddis@28:     stimulusCreate(globalStimParams, stimComponents, 0);
rmeddis@28: if ~isempty(errormsg)     % e.g. limits exceeded
rmeddis@28:     errormsg
rmeddis@28:     return
rmeddis@28: end
rmeddis@28: 
rmeddis@28: %               ------------------------------------------ test stimulus
rmeddis@28: % masker
rmeddis@28: componentNo=1;
rmeddis@28: precedingSilence=stimulusParameters.stimulusDelay;
rmeddis@28: stimComponents(maskerEar,componentNo).type=maskerType;
rmeddis@28: stimComponents(maskerEar,componentNo).toneDuration=maskerDuration;
rmeddis@28: stimComponents(maskerEar,componentNo).frequencies=maskerFrequency;
rmeddis@28: stimComponents(maskerEar,componentNo).amplitudesdB=maskerLevel;
rmeddis@28: stimComponents(maskerEar,componentNo).beginSilence=precedingSilence;
rmeddis@28: stimComponents(maskerEar,componentNo).endSilence=-1;
rmeddis@28: stimComponents(maskerEar,componentNo).AMfrequency=0;
rmeddis@28: stimComponents(maskerEar,componentNo).AMdepth=0;
rmeddis@28: if rampDuration<maskerDuration
rmeddis@28:     % ramps must be shorter than the signal
rmeddis@28:     stimComponents(maskerEar,componentNo).rampOnDur=rampDuration;
rmeddis@28:     stimComponents(maskerEar,componentNo).rampOffDur=rampDuration;
rmeddis@28: else
rmeddis@28:     stimComponents(maskerEar,componentNo).rampOnDur=maskerDuration/2;
rmeddis@28:     stimComponents(maskerEar,componentNo).rampOffDur=maskerDuration/2;
rmeddis@28: end
rmeddis@28: stimComponents(maskerEar,componentNo).phases=...
rmeddis@28:     stimulusParameters.maskerPhase;
rmeddis@28: stimComponents(maskerEar,componentNo).niterations=0;    % for IRN only
rmeddis@28: 
rmeddis@28: % target
rmeddis@28: componentNo=2;
rmeddis@28: targetDelay=precedingSilence+ maskerDuration+ gapDuration;
rmeddis@28: stimComponents(targetEar,componentNo).type=targetType;
rmeddis@28: stimComponents(targetEar,componentNo).toneDuration=targetDuration;
rmeddis@28: stimComponents(targetEar,componentNo).frequencies=targetFrequency;
rmeddis@28: stimComponents(targetEar,componentNo).amplitudesdB=targetLevel;
rmeddis@28: stimComponents(targetEar,componentNo).beginSilence=targetDelay;
rmeddis@28: stimComponents(targetEar,componentNo).endSilence=-1;
rmeddis@28: stimComponents(targetEar,componentNo).AMfrequency=0;
rmeddis@28: stimComponents(targetEar,componentNo).AMdepth=0;
rmeddis@28: if rampDuration<targetDuration
rmeddis@28:     % ramps must be shorter than the signal
rmeddis@28:     stimComponents(targetEar,componentNo).rampOnDur=rampDuration;
rmeddis@28:     stimComponents(targetEar,componentNo).rampOffDur=rampDuration;
rmeddis@28: else
rmeddis@28:     stimComponents(targetEar,componentNo).rampOnDur=0;
rmeddis@28:     stimComponents(targetEar,componentNo).rampOffDur=0;
rmeddis@28: end
rmeddis@28: stimComponents(targetEar,componentNo).phases=stimulusParameters.targetPhase;
rmeddis@28: % stimComponents(targetEar,componentNo)
rmeddis@28: 
rmeddis@28: % background same ear as target
rmeddis@28: componentNo=3;
rmeddis@28: stimComponents(backgroundEar,componentNo).type=backgroundType;
rmeddis@28: switch backgroundType
rmeddis@28:     case 'TEN'
rmeddis@28:         fileName=['..' filesep '..' filesep ...
rmeddis@28:             'multithresholdResources' filesep ...
rmeddis@28:             'backgrounds and maskers'...
rmeddis@28:             filesep 'ten.wav'];
rmeddis@28:         [tenNoise, FS]=wavread(fileName);
rmeddis@28:         
rmeddis@28:         tenNoise=resample(tenNoise, globalStimParams.FS, FS);
rmeddis@28:         stimComponents(backgroundEar,componentNo).type='file';
rmeddis@28:         stimComponents(backgroundEar,componentNo).stimulus=tenNoise';
rmeddis@28: end
rmeddis@28: stimComponents(backgroundEar,componentNo).toneDuration=...
rmeddis@28:     globalStimParams.overallDuration;
rmeddis@28: stimComponents(backgroundEar,componentNo).amplitudesdB=backgroundLevel;
rmeddis@28: stimComponents(backgroundEar,componentNo).beginSilence=0;
rmeddis@28: stimComponents(backgroundEar,componentNo).endSilence=-1;
rmeddis@28: stimComponents(backgroundEar,componentNo).rampOnDur=rampDuration;
rmeddis@28: stimComponents(backgroundEar,componentNo).rampOffDur=rampDuration;
rmeddis@28: stimComponents(backgroundEar,componentNo).AMfrequency=0;
rmeddis@28: stimComponents(backgroundEar,componentNo).AMdepth=0;
rmeddis@28: 
rmeddis@28: % timings used when evaluating MAP peripheral model
rmeddis@28: % this is the Slope during which spikes are counted
rmeddis@28: switch experiment.paradigm
rmeddis@28:     case 'gapDetection'
rmeddis@28:         % gap is the 'target' in this case
rmeddis@28:         stimulusParameters.testTargetBegins=...
rmeddis@28:             stimulusParameters.stimulusDelay...
rmeddis@28:             +stimulusParameters.maskerDuration;
rmeddis@28:         stimulusParameters.testTargetEnds=...
rmeddis@28:             stimulusParameters.testTargetBegins+withinRuns.variableValue;
rmeddis@28:         %     case 'SRT'
rmeddis@28:         %         set(handles.editdigitInput,'visible','off')
rmeddis@28:     otherwise
rmeddis@28:         stimulusParameters.testTargetBegins=targetDelay;
rmeddis@28:         stimulusParameters.testTargetEnds=targetDelay+targetDuration;
rmeddis@28: end
rmeddis@28: 
rmeddis@28: % ------------------------------------------------------------- play!
rmeddis@28: % Create and play stimulus (as required by different paradigms)
rmeddis@28: switch experiment.ear
rmeddis@28:     case {'statsModelLogistic', 'statsModelRareEvent'}
rmeddis@28:         audio=[0;0];            % no need to compute stimulus
rmeddis@28:         
rmeddis@28:     otherwise                   % create the stimulus
rmeddis@28:         [targetStimulus, errormsg]= ...
rmeddis@28:             stimulusCreate(globalStimParams, stimComponents, 0);
rmeddis@28:         
rmeddis@28:         if ~isempty(errormsg)   % e.g. limits exceeded
rmeddis@28:             errormsg
rmeddis@28:             return
rmeddis@28:         end
rmeddis@28:         
rmeddis@28:         switch experiment.ear
rmeddis@28:             case {'MAPmodel' , 'MAPmodelMultiCh', 'MAPmodelSingleCh', 'MAPmodelListen'}
rmeddis@28:                 % model requires no calibration correction;
rmeddis@28:                 %  signal is already in Pascals
rmeddis@28:                 globalStimParams.audioOutCorrection=1;
rmeddis@28:                 % use only the targetStimulus for the MAP model
rmeddis@28:                 audio=targetStimulus;
rmeddis@28:                 
rmeddis@28:             otherwise           % left, right diotic dichotic
rmeddis@28:                 if stimulusParameters.includeCue
rmeddis@28:                     audio= [cueStimulus;  targetStimulus];
rmeddis@28:                 else            % no cue
rmeddis@28:                     audio=targetStimulus;
rmeddis@28:                 end
rmeddis@28:         end
rmeddis@28:         
rmeddis@28:         % playtime
rmeddis@28:         % order of the cue and test stimuli varies for 2AFC
rmeddis@28:         switch experiment.threshEstMethod
rmeddis@28:             case {'2I2AFC++', '2I2AFC+++'}
rmeddis@28:                 % intervening silence (currently none; masking delay serves this purpose)
rmeddis@28:                 IAFCinterveningSilence=zeros(round(AFCsilenceDuration/dt),2);
rmeddis@28:                 if rand>0.5				% put test stimulus first
rmeddis@28:                     stimulusParameters.testTargetBegins=targetDelay ;
rmeddis@28:                     stimulusParameters.testTargetEnds= ...
rmeddis@28:                         targetDelay+targetDuration;
rmeddis@28:                     stimulusParameters.testNonTargetBegins=...
rmeddis@28:                         length(cueStimulus)*dt ...
rmeddis@28:                         + AFCsilenceDuration +targetDelay ;
rmeddis@28:                     stimulusParameters.testNonTargetEnds=...
rmeddis@28:                         length(cueStimulus)*dt ...
rmeddis@28:                         + AFCsilenceDuration+targetDelay+targetDuration;
rmeddis@28:                     
rmeddis@28:                     set(handles.pushbutton1,'backgroundcolor','r'), drawnow
rmeddis@28:                     y=audioplayer(targetStimulus, globalStimParams.FS, 24);
rmeddis@28:                     playblocking(y)
rmeddis@28:                     set(handles.pushbutton1,'backgroundcolor',...
rmeddis@28:                         get(0,'defaultUicontrolBackgroundColor')), drawnow
rmeddis@28:                     y=audioplayer(IAFCinterveningSilence, ...
rmeddis@28:                         globalStimParams.FS, 24);
rmeddis@28:                     playblocking(y)
rmeddis@28:                     set(handles.pushbutton2,'backgroundcolor','r'), drawnow
rmeddis@28:                     y=audioplayer(cueStimulus, globalStimParams.FS, 24);
rmeddis@28:                     playblocking(y)
rmeddis@28:                     set(handles.pushbutton2,'backgroundcolor',...
rmeddis@28:                         get(0,'defaultUicontrolBackgroundColor')), drawnow
rmeddis@28:                     withinRuns.stimulusOrder='targetFirst';
rmeddis@28:                     audio= [targetStimulus; IAFCinterveningSilence; ...
rmeddis@28:                         cueStimulus];   % for plotting purposes later
rmeddis@28:                     
rmeddis@28:                 else					% put test stimulus second
rmeddis@28:                     stimulusParameters.testTargetBegins=...
rmeddis@28:                         length(cueStimulus)*dt ...
rmeddis@28:                         + AFCsilenceDuration +targetDelay ;
rmeddis@28:                     stimulusParameters.testTargetEnds=...
rmeddis@28:                         length(cueStimulus)*dt ...
rmeddis@28:                         + AFCsilenceDuration+targetDelay+targetDuration;
rmeddis@28:                     stimulusParameters.testNonTargetBegins=targetDelay ;
rmeddis@28:                     stimulusParameters.testNonTargetEnds=...
rmeddis@28:                         targetDelay+targetDuration;
rmeddis@28:                     
rmeddis@28:                     set(handles.pushbutton1,'backgroundcolor','r'),drawnow
rmeddis@28:                     y=audioplayer(cueStimulus, globalStimParams.FS, 24);
rmeddis@28:                     playblocking(y)
rmeddis@28:                     set(handles.pushbutton1,'backgroundcolor',...
rmeddis@28:                         get(0,'defaultUicontrolBackgroundColor')), drawnow
rmeddis@28:                     y=audioplayer(IAFCinterveningSilence, ...
rmeddis@28:                         globalStimParams.FS, 24);
rmeddis@28:                     playblocking(y)
rmeddis@28:                     set(handles.pushbutton2,'backgroundcolor','r'), drawnow
rmeddis@28:                     y=audioplayer(targetStimulus, globalStimParams.FS, 24);
rmeddis@28:                     playblocking(y)
rmeddis@28:                     set(handles.pushbutton2,'backgroundcolor',...
rmeddis@28:                         get(0,'defaultUicontrolBackgroundColor')), drawnow
rmeddis@28:                     withinRuns.stimulusOrder='targetSecond';
rmeddis@28:                     audio= [cueStimulus; IAFCinterveningSilence; ...
rmeddis@28:                         targetStimulus];    % for plotting purposes later
rmeddis@28:                 end
rmeddis@28:             otherwise                       % singleInterval
rmeddis@28:                 if strcmp(experiment.ear,'MAPmodel') ...
rmeddis@28:                         || strcmp(experiment.ear,'MAPmodelMultiCh') ...
rmeddis@28:                         || strcmp(experiment.ear,'MAPmodelSingleCh') ...
rmeddis@28:                         ||strcmp(experiment.ear,'MAPmodelListen')
rmeddis@28:                     % don't play for MAPmodel
rmeddis@28:                     switch experiment.ear
rmeddis@28:                         % except on special request
rmeddis@28:                         case {'MAPmodelListen'}
rmeddis@28:                             y=audioplayer(audio, globalStimParams.FS, 24);
rmeddis@28:                             playblocking(y) % suspends operations until completed
rmeddis@28:                     end
rmeddis@28:                 else
rmeddis@28:                     y=audioplayer(audio, globalStimParams.FS, 24);
rmeddis@28:                     playblocking(y)
rmeddis@28:                 end	  %   if experiment.ear
rmeddis@28:         end	% switch experiment.threshEstMethod
rmeddis@28: end % switch experiment.ear
rmeddis@28: 
rmeddis@28: 
rmeddis@28: % switch experiment.ear
rmeddis@28: %     case	{'MAPmodel', 'MAPmodelListen',  'MAPmodelMultiCh','MAPmodelSingleCh'}
rmeddis@28: %         % save audio for later reference or for input to MAP model
rmeddis@28: %         wavwrite(audio/max(audio), globalStimParams.FS,32,'stimulus')
rmeddis@28: % end
rmeddis@28: 
rmeddis@28: % Panel 1
rmeddis@28: % graphical presentation of the stimulus
rmeddis@28: % NB shown *after* the stimulus has been presented
rmeddis@28: axes(expGUIhandles.axes1), cla
rmeddis@28: % plot is HW rectified and plotted as dB re 28e-6
rmeddis@28: % calibration is ignored
rmeddis@28: t=dt:dt:dt*length(audio);
rmeddis@28: plot(t,stimulusParameters.calibrationdB+20*log10((abs(audio)+1e-10)/28e-6))
rmeddis@28: % set(gca,'xtick',[])
rmeddis@28: ylim([-20 100])
rmeddis@28: ylabel('stimulus (dB SPL)')
rmeddis@28: xlim([0 t(end)])
rmeddis@28: grid on
rmeddis@28: header=[betweenRuns.variableName1  ': ' ...
rmeddis@28:     num2str(betweenRuns.var1Sequence(betweenRuns.runNumber))];
rmeddis@28: header=[header '      ' num2str(...
rmeddis@28:     betweenRuns.var2Sequence(betweenRuns.runNumber)) ':' ...
rmeddis@28:     betweenRuns.variableName2 ];
rmeddis@28: title(header)
rmeddis@28: