PlayMovieDriftBar5.m

% this script do not need the prepared movie text.
% Niki 2016/4/1
% add square grating on the bar to enhance response
% Niki 2016/6/20

clear
clc

Duration = 20;% seconds
NumSweeps = 10;
BarWidth = pi/45; % in radians
FlashFrequency = 5; % Hz
InterBlockRest = 20; % seconds

TriggerWidth = 60; % pixel size of Trigger Singal (a white square on the bottom right corner on the screen)
TriggerDuration = 1; % seconds

MoveDirection = [1,2,3,4];% [1:right, 2:bottom, 3:left, 4:top]
GratingCycles = [6,8,6,8];
BlocksForEachCondition = 3;

CameraFPS = 30;
DataFPS = 7.5;
PostStimDur = 8;
TotalDuration = ((Duration*NumSweeps+InterBlockRest)*length(MoveDirection)*BlocksForEachCondition)/60;
SingleStimDuration = (Duration*NumSweeps+PostStimDur)*1000;
NumCameraFrames = (Duration*NumSweeps+PostStimDur)*CameraFPS;
NumDataFrames = (Duration*NumSweeps+PostStimDur)*DataFPS;
NumPostStimFrames = PostStimDur*DataFPS;

stimSeq = Shuffle(repmat((MoveDirection)',1,BlocksForEachCondition));
s = screen;
s.view_distance = 250; % cm
BarWidthPix = round(s.pix_per_rad*BarWidth);
clear s

%%
try
    Screen('Preference', 'SkipSyncTests', 1);
    % This script calls Psychtoolbox commands available only in OpenGL-based
    % versions of the Psychtoolbox. (So far, the OS X Psychtoolbox is the
    % only OpenGL-base Psychtoolbox.)  The Psychtoolbox command AssertPsychOpenGL will issue
    % an error message if someone tries to execute this script on a computer without
    % an OpenGL Psychtoolbox
    AssertOpenGL;
    
    % Get the list of screens and choose the one with the highest screen number.
    % Screen 0 is, by definition, the display with the menu bar. Often when
    % two monitors are connected the one without the menu bar is used as
    % the stimulus display.  Chosing the display with the highest dislay number is
    % a best guess about where you want the stimulus displayed.
    screens=Screen('Screens');
    screenNumber=max(screens);
    
    % Find the color values which correspond to white and black: Usually
    % black is always 0 and white 255, but this rule is not true if one of
    % the high precision framebuffer modes is enabled via the
    % PsychImaging() commmand, so we query the true values via the
    % functions WhiteIndex and BlackIndex:
    white=WhiteIndex(screenNumber);
    black=BlackIndex(screenNumber);
    
    % Round gray to integral number, to avoid roundoff artifacts with some
    % graphics cards:
    gray=round((white+black)/2);
    
    % This makes sure that on floating point framebuffers we still get a
    % well defined gray. It isn't strictly neccessary in this demo:
    if gray == white
        gray=white / 2;
    end
    
    % Open a double buffered fullscreen window and select a black background
    % color:
    [w, screenRect]=Screen('OpenWindow',screenNumber, black);
    
    frameRate=Screen('FrameRate',screenNumber);
    % If MacOSX does not know the frame rate the 'FrameRate' will return 0.
    % That usually means we run on a flat panel with 60 Hz fixed refresh
    % rate:
    if frameRate == 0
        frameRate=60;
    end
    flashdurframes = frameRate/FlashFrequency/2;
    
    HideCursor;
    ListenChar(2); % suppress input from keybroad to appear in the command window or script
    
    while 1
        WaitSecs(0.2);
        [~, ~, keyCode] = KbCheck();
        if any(strcmp(KbName(keyCode), {'ESCAPE','Escape','ESC','esc'}))
            ListenChar(1)
            sca;
            return;
        elseif any(strcmp(KbName(keyCode), {'space','s'}))
            break
        end
    end
    
    % waitframes = 1 means: Redraw every monitor refresh. If your GPU is
    % not fast enough to do this, you can increment this to only redraw
    % every n'th refresh. All animation paramters will adapt to still
    % provide the proper grating. However, if you have a fine grating
    % drifting at a high speed, the refresh rate must exceed that
    % "effective" grating speed to avoid aliasing artifacts in time, i.e.,
    % to make sure to satisfy the constraints of the sampling theorem
    % (See Wikipedia: "Nyquist?Shannon sampling theorem" for a starter, if
    % you don't know what this means):
    waitframes = 1;
    % Query duration of one monitor refresh interval:
    ifi=Screen('GetFlipInterval', w);
    % Translate frames into seconds for screen update interval:
    waitduration = (waitframes - 0.5) * ifi;
    
    for i_block = 1:BlocksForEachCondition
        for j = 1:length(MoveDirection);
            BarRect = MakeDriftBarRect(Duration, BarWidthPix, stimSeq(j,i_block), screenRect, frameRate);
            BarRectSplit = SplitBarRect(BarRect, GratingCycles(stimSeq(j,i_block)), screenRect);
            
            TriggerFrames = TriggerDuration * frameRate;
            % right bottom:
            %TriggerRect = [screenRect([3,4])-TriggerWidth, screenRect([3,4])]';
            % right top:
            TriggerRect = [screenRect(3)-TriggerWidth, screenRect([2,3]), TriggerWidth]';
                        
            singleSweepFrames = size(BarRect,2);
            totalframes = singleSweepFrames*NumSweeps;
            GratingType = mod(ceil((1:totalframes)/flashdurframes),2)==0;
            % onset = nan(totalframes,1);
            frame_idx = mod((1:totalframes)-1,singleSweepFrames)+1;  
            % Draw the 1st frame:
            %Screen('Fillrect', w, gray);
            Screen('Fillrect', w, white, [squeeze(BarRectSplit(:,1,1:2:end)), TriggerRect]);
            Screen('Fillrect', w, black, squeeze(BarRectSplit(:,1,2:2:end)));
            % Perform initial Flip to sync us to the VBL and for getting an initial
            % VBL-Timestamp as timing baseline for our redraw loop:
            vbl=Screen('Flip', w);
            
            % Animation loop:
            for i=2:totalframes
                % Draw image:
                if i<=TriggerFrames
                    if GratingType(i)
                        Screen('Fillrect', w, white, [squeeze(BarRectSplit(:,frame_idx(i),1:2:end)), TriggerRect]);
%                         Screen('Fillrect', w, black, squeeze(BarRectSplit(:,frame_idx(i),2:2:end)));
                    else
%                         Screen('Fillrect', w, black, squeeze(BarRectSplit(:,frame_idx(i),1:2:end)));
                        Screen('Fillrect', w, white, [squeeze(BarRectSplit(:,frame_idx(i),2:2:end)), TriggerRect]);
                    end
                else
                    if GratingType(i)
                        Screen('Fillrect', w, white, squeeze(BarRectSplit(:,frame_idx(i),1:2:end)));
%                         Screen('Fillrect', w, black, squeeze(BarRectSplit(:,frame_idx(i),2:2:end)));
                    else
%                         Screen('Fillrect', w, black, squeeze(BarRectSplit(:,frame_idx(i),1:2:end)));
                        Screen('Fillrect', w, white, squeeze(BarRectSplit(:,frame_idx(i),2:2:end)));
                    end
                end
                % Screen('DrawText', w, [num2str(i_rep), num2str(frameRate)],10 );
                % Flip 'waitframes' monitor refresh intervals after last redraw.
                % Providing this 'when' timestamp allows for optimal timing
                % precision in stimulus onset, a stable animation framerate and at
                % the same time allows the built-in "skipped frames" detector to
                % work optimally and report skipped frames due to hardware
                % overload:
                %[vbl, onset(i)]= Screen('Flip', w, vbl + waitduration);
                vbl= Screen('Flip', w, vbl + waitduration);
                
                %[~, onset(i)]= Screen('Flip', w);
            end
            DrawFormattedText(w, 'Waiting for next block\n','center','center',white);
            Screen('Flip',w);
            WaitSecs(InterBlockRest);
        end
        
    end
    
    Priority(0);
    
    % Close window:
    Screen('CloseAll');
    ListenChar(1);
    ShowCursor;
    % difonset = diff(onset(:));
    % plot(difonset);
    % title(['total duration: ' num2str(range(onset(:)))])
    
    save(['StimSequcence',datestr(now,'yyyymmddTHHMMSS')], 'stimSeq')
catch
    
    Screen('CloseAll');
    Priority(0);
    ListenChar(1);
    ShowCursor;
    psychrethrow(psychlasterror);
    
end %try..catch..