ExperimentScript_ME_v3.m

% MATLAB SCRIPT FOR THE AUTOMATICITY/NON-AUTOMATICITY EXPERIMENT
% An experiment in which automatic and non-automatic movements are tested.
% For the automatic movements, participants need to perform a pre-learned
% 12 digit sequence during a finger tapping and a foot stomping task. Also
% a non-automatic newly learned 12 digit sequence will be tested including both
% tasks. The script is randomized in such a way that the participant either
% starts with the automatic or non-automatic tasks. Within these options, it
% is also randomized whether they start with the finger tapping or foot
% stomping task. Tasks will automatically proceed after each other. 
% Note that for the non-automatic tasks the participant is also tested on
% automaticity for the movement during a short dual task, this holds for
% each limb.

%Before starting the experiment, clear the workspace.
clear all

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%LSL SET UP
%LSL outlet sending events
%
% Create and load the lab streaming layer library
lib = lsl_loadlib();
%
% Make a new stream outlet.
% info = lsl_streaminfo([lib handle],[name],[type],[channelcount],[fs],[channelformat],[sourceid])
% > name = name of stream; describes device/product
% > type = content type of stream (EEG, Markers)
% > channelcount = nr of channels per sample
% > fs = samplking rate (Hz) as advertized by data source
% > channelformat = cf_float32, cf__double64, cf_string, cf_int32, cf_int16
% > sourceid = unique identifier for source or device, if available
info = lsl_streaminfo(lib,'AutovsNAuto','Markers',1,0.0,'cf_string','sdfwerr32432');
%
% Open an outlet for the data to run through.
outlet = lsl_outlet(info);
%
% Create marker id's
instructions = 'instructions';
instructions_nonauto = 'instructions_nonauto';
instructions_auto = 'instructions_auto';
rest='rest';
nonauto='nonauto';
auto='auto';
finger_auto='finger_auto';
finger_nonauto='finger_nonauto';
finger_nonautotest='finger_nonautotest';
foot_auto='foot_auto';
foot_nonauto='foot_nonauto';
foot_nonautotest='foot_nonautotest';

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%INITIALISATION

%Open Phsychtoolbox.
PsychDefaultSetup(2);
KbName('UnifyKeyNames'); %Links the key presses to the key board names
% ListenChar; % option A
KbQueueCreate; % option B
KbQueueStart; % option B

%Skip screen synchronization to prevent Pyshtoolbox for freezing
Screen('Preference', 'SkipSyncTests', 1);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%SETUP PARAMETERS

t1 = 20; %Resting period in seconds
t2 = 5;  %Random interval around the resting period time
t3 = 8.5; %Duration of a trial (tapping/stomping the sequence 1 time)
N_letters = 8; % 8 letters presented during a trial for the automaticity test
N_trials = 8; % number of trials for the automaticity test

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%LOAD METRONOME SOUNDS (PsychToolbox)
audio_dir='C:\Users\mtabo\Documents\TryOutScript\metronomesounds';
cd(audio_dir)
[WAVMetronome8.wave,WAVMetronome8.fs]       = audioread('Metronome8.wav');
[WAVMetronome600.wave,WAVMetronome600.fs]       = audioread('Metronome600.wav');
[WAVMetronome5.wave,WAVMetronome5.fs]       = audioread('Metronome5.wav');

% change rows<>columns
WAVMetronome8.wave = WAVMetronome8.wave';         WAVMetronome8.nrChan=2;
WAVMetronome600.wave = WAVMetronome600.wave';         WAVMetronome600.nrChan=2;
WAVMetronome5.wave = WAVMetronome5.wave';         WAVMetronome5.nrChan=2;

% CREATE AND FILL AUDIO BUFFER
% Initialize Sounddriver
% This routine loads the PsychPortAudio sound driver for high precision, low latency,
% multichannel sound playback and recording
% Call it at the beginning of your experiment script, optionally providing the
% 'reallyneedlowlatency'-flag set to 1 to push really hard for low latency
InitializePsychSound(1);

priority = 0;                       % 0 = better quality, increased latency; 1 = minimum latency
duration = 1;                       % number of repetitions of the wav-file
PsychPortAudio('Verbosity',1);      % verbosity = "wordiness" -> 1= print errors

% Get audio device
h_device = PsychPortAudio ('GetDevices');

% Open handle
h_Metronome8   = PsychPortAudio('Open', [], [], priority, WAVMetronome8.fs, WAVMetronome8.nrChan);
h_Metronome600   = PsychPortAudio('Open', [], [], priority, WAVMetronome600.fs, WAVMetronome600.nrChan);
h_Metronome5   = PsychPortAudio('Open', [], [], priority, WAVMetronome5.fs, WAVMetronome5.nrChan);

% Fill buffer
PsychPortAudio('FillBuffer', h_Metronome8, WAVMetronome8.wave);
PsychPortAudio('FillBuffer', h_Metronome600, WAVMetronome600.wave);
PsychPortAudio('FillBuffer', h_Metronome5, WAVMetronome5.wave);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%SAVE FILES IN FOLDER

fprintf('Select the project directory \n')
root_dir=uigetdir('C:\Users\mtabo\Documents\TryOutScript\', 'Select the project directory');

complete=0;
while complete==0
    sub_ID=input('What is the subject ID (2 digit number) \n', 's');
    sub=sprintf('sub-%s', sub_ID);
        rec_n=input('What is the number of the recording? \n');
        rec=sprintf('rec-%.2d', rec_n);
     
    inf=fprintf('\n root_dir = %s \n sub = %s \n rec = %s \n', root_dir, sub, rec);
    correct=input('Is the above information correct? (y/n) \n', 's');
    if strcmp(correct, 'y')
        complete=1;
    else
        continue
    end
end

% go to subject folder
sub_dir=fullfile(root_dir, sub);
if ~exist(sub_dir)
    mkdir(sub_dir)
end
cd(sub_dir)

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%RANDOMIZATION

%Create a vector to represent the two different options (1=non-automatic
%test, 2=automatic test).
order_experiment=[1,2];
%Randomize and determine the order of the two different tests. Either
%[1,2] or [2,1].
%order_experiment=stim(randperm(length(stim)));
%Save the order of the experiment
save('order_experiment.mat', 'order_experiment');

%Create a vector to represent the two different options for non-automatic
%tasks (3 = finger tapping, 4 = foot stomping).
order_nonauto=[3,4];
%Randomize and determine the order of the limb tasks. Either [3,4] or
%[4,3].
%order_nonauto=nonauto_exp(randperm(length(nonauto_exp)));
%Save the order of the non-automatic experiment
save('order_nonauto.mat', 'order_nonauto');

%Create a vector to represent the two different options for automatic tasks
%(5 = finger tapping, 6 = foot stomping).
order_auto=[5,6];
%Randomize and determine the order of the limb tasks. Either [5,6] or
%[6,5].
%order_auto=auto_exp(randperm(length(auto_exp)));
%Save the order of the automatic experiment
save('order_auto.mat', 'order_auto');

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% SCREEN PREPARATION

% Get the screen numbers.
screens = Screen('Screens');

% Select the external screen if it is present, else revert to the native
% screen
screenNumber = max(screens);

% Define black, white and grey
black = BlackIndex(screenNumber);
white = WhiteIndex(screenNumber);

% Open an on screen window and color it grey
[window, windowRect] = PsychImaging('OpenWindow', screenNumber, black);

% Get the size of the on screen window in pixels
% For help see: Screen WindowSize?
[screenXpixels, screenYpixels] = Screen('WindowSize', window);

% Get the centre coordinate of the window in pixels
% For help see: help RectCenter
[xCenter, yCenter] = RectCenter(windowRect);

% Query the frame duration
ifi = Screen('GetFlipInterval', window);

% Set up alpha-blending for smooth (anti-aliased) lines
Screen('BlendFunction', window, 'GL_SRC_ALPHA', 'GL_ONE_MINUS_SRC_ALPHA');

% Preparations for the fixation cross so we only need to do this once
fixCrossDimPix = 40; % Here we set the size of the arms of our fixation cross
xCoords = [-fixCrossDimPix fixCrossDimPix 0 0]; % Set the coordinates (these are all relative to zero we will let the drawing routine center the cross in the center of our monitor for us)
yCoords = [0 0 -fixCrossDimPix fixCrossDimPix];
allCoords = [xCoords; yCoords];
lineWidthPix = 4;% Set the line width for the fixation cross

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% START OF THE EXPERIMENT

%Empty structure for key presses -> use later again so it saves the key
%presses within this structure -> save at the end
presses_handauto=struct([]);
presses_handnonauto=struct([]);
presses_handnonautotest=struct([]);
letters_handnonautotest=struct([]);
letters_footnonautotest=struct([]);

%Instruction experiment
Screen('TextSize',window,25);
DrawFormattedText(window,'You will now start with the experiment. \n You will either start with the automatic tasks, performing the at home studied sequence, \n or with the non-automatic tasks, for which you will be presented with a new sequence, \n which you can study for 10 minutes for each limb. \n Note that for the non-automatic tasks you will also perform an automaticity test (dual task). \n This is the same test as you just did for the automatic (at home studied) sequence. \n\n Instructions will appear at the start of each new task. \n\n Press any key to continue and see with which test you start.','center', 'center', white);
vbl = Screen('Flip', window);
KbStrokeWait; %wait for response to terminate instructions

%Start the randomization loop between non-automatic (=1) and automatic (=2)
for i=order_experiment; %Either [1,2] or [2,1] -> determines the order of the tasks
    outlet.push_sample({instructions})
    
    % NON-AUTOMATICITY TASKS
    if i==1; %Start with the non-automatic tasks
        outlet.push_sample({nonauto})
        
        % Show instructions
        Screen('TextSize', window, 25);
        DrawFormattedText(window, 'NON-AUTOMATICITY TASK \n\n You will perform the finger tapping and foot stomping task for a new sequence. \n You will either start with the finger tapping or the foot stomping task. \n You will get 10 minutes practice time before each task starts. \n Press any key to continue.', 'center', 'center', white);
        vbl= Screen('Flip', window);
        KbStrokeWait; %wait for response to terminate instructions
        
        %Determine to start with either hand or foot task.
        for j=order_nonauto %Either [3,4] or [4,3] -> determines the order of the limbs
            outlet.push_sample({instructions_nonauto})
            
            % NON-AUTOMATIC FINGER TAPPING TASK
            if j==3;
                outlet.push_sample({finger_nonauto})
                
                %Non automatic finger tapping task instructions
                %Practice new sequence
                Screen('TextSize', window, 25);
                DrawFormattedText(window, 'You will now perform the non-automaticity FINGER tapping task. \n For the next 10 minutes you can practice a new sequence for the finger tapping task, \n the same way you practiced at home. \n After that we will start with the finger tapping task. \n Press any key to see the new sequence and start practicing.', 'center', 'center', white);
                vbl= Screen('Flip', window);
                KbStrokeWait; %wait for response to terminate instructions
                
                %Presenting the new (non-automatic) sequence on the screen
                Screen('TextSize', window, 50);
                DrawFormattedText(window, '2 1 2 3 2 1 3 2 4 2 4 1', 'center', 'center', white);
                vbl= Screen('Flip', window);
                PsychPortAudio('Start', h_Metronome600, 1, [], []); %Play metronome sound file (10 minutes)
                
                Screen('TextSize', window, 25);
                DrawFormattedText(window, 'The time to practice the new sequence is over. \n Press any key to continue to the finger tapping experiment.', 'center', 'center', white);
                vbl= Screen('Flip', window);
                KbStrokeWait; %wait for response to terminate instructions
                
                %Instructions non-automatic finger tapping task
                Screen('TextSize',window,25);
                DrawFormattedText(window, 'You will now perform the finger tapping task for the sequence you just practiced (non-automatic). \n There is a total of 10 trials, so performing the sequence 10 times. \n  In between each trial there is a rest period of 20 seconds. \n During this rest you will hear a metronome sound, tap the sequence according to this interval sound. \n The trials and rest periods are indicated with Trial and Rest showing on the screen. \n\n When the experiment starts you cannot talk anymore. \n Furthermore, it is important to stay still except for your right hand. \n If you press any key, the experiment starts. \n It will start with a rest period. \n Whenever Trial appears on the screen, you should start tapping the sequence. \n Press any key to start the finger tapping experiment.','center','center', white);
                vbl = Screen('Flip', window);
                KbStrokeWait; %wait for response to terminate instructions
                
                %Stimulus for finger tapping non-automatic sequence
                for l = 1:10 %Will perform 10 trials
                    presses_handnonauto(l).key=[];
                    %Rest period between each sequence 20-25 seconds
                    outlet.push_sample({rest})
                    % Fixation cross during rest
                    Screen('TextSize', window, 36);
                    Screen('TextSize',window,40);
                    DrawFormattedText(window, 'Rest \n\n\n\n', 'center', 'center', white)
                    Screen('DrawLines', window, allCoords,...
                        lineWidthPix, white, [xCenter yCenter], 2);
                    Screen('Flip', window);
                    PsychPortAudio('Start', h_Metronome8, 1, [], []); % Play metronome sound file (8 seconds)
                    WaitSecs(t1+randi(t2));
                    
                    %Fixation cross during finger tapping trial
                    Screen('TextSize', window, 36);
                    Screen('TextSize',window,40);
                    DrawFormattedText(window, 'Trial \n\n\n\n', 'center', 'center', white)
                    Screen('DrawLines', window, allCoords,...
                        lineWidthPix, white, [xCenter yCenter], 2);
                    Screen('Flip', window)
                    m=1;
                    
                    % Record key presses
                    start_timer=GetSecs;
                    while GetSecs-start_timer < t3
                        if length(presses_handnonauto(l).key)<12
                            [secs, keyCode, deltaSecs] = KbWait([],2);
                            if any(keyCode)
                                fprintf('the key value was: %s \n', KbName(find(keyCode)));
                                presses_handnonauto(l).key{m}=KbName(find(keyCode));
                                presses_handnonauto(l).secs(m)=GetSecs;
                                m=m+1;
                            elseif KbName('ESCAPE')
                            elseif GetSecs-start_timer >= t3
                                break
                            end
                        end
                    end
                    
                end
                
                
                % End of the non-automatic finger tapping task is reached.
                Screen('TextSize', window, 25);
                DrawFormattedText(window, 'This is the end of the non-automatic finger tapping task. \n Press any key to continue with the automaticity test for this task.', 'center', 'center', white);
                vbl= Screen('Flip', window);
                KbStrokeWait; %wait for response to terminate instructions
                
                
                %Instruction automaticity task finger tapping
                Screen('TextSize',window,25);
                DrawFormattedText(window, 'You will now perform the sequence you just learned during an automaticity test for the FINGER tapping task. \n  Letters will be shown on the screen (A,G,O,L) while you perform the task. \n The goal is to perform the sequence tapping while counting how many times G is presented. \n After each time you tapped the full sequence, you should tell us how many times G was presented. \n We will perform 8 trails. \n\n Note that during the tapping task you cannot talk. \n Try to keep your body movements as still as possible exept for the right hand. \n\n In between the trials you will see a fixation cross for 10 seconds. \n During the first 5 seconds you will hear a metronome sound. \n Tap the sequence on this rhythm, which is the same as you practiced before. \n\n We will start with a fixation cross on the screen for 10 seconds. \n After that the first trial will start automatically. \n So start tapping the sequence as soon as a letter on the screen appears. \n Press any key to continue and start the test.','center','center', white);
                vbl = Screen('Flip', window);
                KbStrokeWait; %wait for response to terminate instructions
                
                for r=1:N_trials
                    %Presentation of the letters on the screen (dual task). -> is random.
                    %Participant has to count the amount that G was presented.
                    Letterlist= {'A', 'G', 'O', 'L'};
                    letter_order=randi(length(Letterlist), 1, N_letters);
                    letters_handnonautotest(r).presented =Letterlist(letter_order);
                    
                    % always start with a 10 seconds fixation cross with 5 seconds of metronome
                    % sound
                    Screen('TextSize', window, 36);
                    Screen('DrawLines', window, allCoords,...
                        lineWidthPix, white, [xCenter yCenter], 2);
                    Screen('Flip', window);
                    PsychPortAudio('Start', h_Metronome5, 1, [], []); % Play metronome sound file
                    WaitSecs(10);
                    
                    %Presentation of random letters on the screen during the finger
                    %tapping test + recording of the key presses
                    outlet.push_sample({finger_nonautotest})
                    v=1; % first key press
                    %FlushEvents('keyDown'); % option A: clear all previous key presses from the list
                    KbQueueFlush; % option B: clear all previous key presses from the list
                    for w=1:N_letters
                        % present random letter
                        Screen('TextSize', window, 100);
                        DrawFormattedText(window, [cell2mat(letters_handnonautotest(r).presented(w))],'center','center', white);
                        vbl = Screen('Flip', window);
                        time_letter=rand(1)+0.5; %Speed with which the letters are presented = A randomized value between 0 and 1, + 0.5 sec
                        
                        % record key presses
                        start_timer=GetSecs;
                        while GetSecs-start_timer<time_letter
                            %                   if CharAvail % option A
                            %                     [ch, when]=GetChar;
                            %                     fprintf('the key value was: %s \n', ch);
                            %                     presses_handtest(j).key{m}=ch;
                            %                     presses_handtest(j).secs(m)=when.secs;% better than using string arrays
                            %                     m=m+1;
                            %                   end
                            [ pressed, firstPress, ~, lastPress, ~]=KbQueueCheck; % option B
                            if pressed
                                if isempty(find(firstPress~=lastPress)) % no key was pressed twice
                                    keys=KbName(find(firstPress)); % find the pressed keys
                                    [timing, idx]=sort(firstPress(find(firstPress))); % get timing of key presses in ascending order
                                    keys=keys(idx); % sort the pressed keys in ascending order
                                    key_n=length(keys); % number of pressed keys
                                    presses_handnonautotest(r).key(v:v+key_n-1)=keys;
                                    presses_handnonautotest(r).secs(v:v+key_n-1)=timing;
                                    v=v+key_n;
                                    KbQueueFlush;
                                else
                                    error('key was pressed twice') % if this error occurs we need to find a way to handle this
                                end
                            end
                        end
                        
                        % between each letter red fixation cross
                        Screen('DrawLines', window, allCoords,...
                            lineWidthPix, [1 0 0], [xCenter yCenter], 2);
                        Screen('Flip', window);
                        WaitSecs (0.2);
                    end
                    
                    % present white fixation cross for some seconds to show that
                    % trial is over
                    Screen('TextSize', window, 36);
                    Screen('DrawLines', window, allCoords,...
                        lineWidthPix, white, [xCenter yCenter], 2);
                    Screen('Flip', window);
                    WaitSecs(2);
                    
                    % show feedback
                    % ask how many G's were presented
                    Screen('TextSize',window,30);
                    DrawFormattedText(window, 'How many times was G presented? ','center','center', white);
                    vbl = Screen('Flip', window);
                    [secs, keyCode, deltaSecs]=KbWait;
                    letters_handnonautotest(r).reported_G=KbName(find(keyCode));
                    DrawFormattedText(window, ['Your answer: ' letters_handnonautotest(r).reported_G '\n Press any key to continue.'],'center','center', white);
                    vbl = Screen('Flip', window);
                    KbStrokeWait; %wait for response to terminate instruction
                    % Show if the reported number of Gs is correct
                    if str2num(letters_handnonautotest(r).reported_G)==sum(strcmp(letters_handnonautotest(r).presented, 'G'))
                        DrawFormattedText(window, ['This is correct. \n Press any key to continue.' ],'center','center', white);
                        vbl = Screen('Flip', window);
                    else
                        DrawFormattedText(window, ['This is incorrect. \n Press any key to continue.' ],'center','center', white);
                        vbl = Screen('Flip', window);
                    end
                    KbStrokeWait;
                    % Show if the tempo was correct. ! reflect if you want to show
                    % this, because we cannot show this result for the foot stomping
                    margin=0.1; % margin of error: think about what is most convenient
                    if all(abs(diff(presses_handnonautotest(r).secs)-1/1.5)<margin)
                        DrawFormattedText(window, ['Your tempo was ok. \n Press any key to continue.' ],'center','center', white);
                        vbl = Screen('Flip', window);
                    else
                        DrawFormattedText(window, ['Your tempo was not ok. \n Press any key to continue.' ],'center','center', white);
                        vbl = Screen('Flip', window);
                    end
                    KbStrokeWait;
                    DrawFormattedText(window, 'Press any key to continue with the next trail. \n Note that you will first start with a fixation cross again. \n Start tapping the sequence as soon as a letter on the screen appears.' ,'center','center', white);
                    vbl = Screen('Flip', window);
                    KbStrokeWait;
                end
                
                % After all trials completed, the end of the finger tapping task is
                % reached.
                Screen('TextSize',window,30);
                DrawFormattedText(window, 'This is the end of the automaticity test for the non-automatic finger tapping task. \n  Press any key to continue with the experiment.' ,'center','center', white);
                vbl = Screen('Flip', window);
                save('letters_handnonautotest.mat', 'letters_handnonautotest'); % save the letters that were presented and the reported number of g's
                save('presses_handnonautotest.mat', 'presses_handnonautotest'); %Save which keys where pressed during the experiment
                KbStrokeWait; %wait for response to terminate instructions
                
            end
            
            
            %NON-AUTOMATIC FOOT STOMPING TASK
            if j==4;
                outlet.push_sample({foot_nonauto})
                
                %Non automatic foot stomping task instructions
                %Practice new sequence
                Screen('TextSize', window, 25);
                DrawFormattedText(window, 'You will now perform the non-automaticity FOOT stomping task. \n For the next 10 minutes you can practice the new sequence for the foot stomping task. \n After that the experiment for the foot stomping task will start. \n\n\n Press any key to see the new sequence and start practicing.', 'center', 'center', white);
                vbl= Screen('Flip', window);
                KbStrokeWait; %wait for response to terminate instructions
                
                % Presenting the new (non-automatic) sequence on the screen
                Screen('TextSize', window, 50);
                DrawFormattedText(window, '2 1 2 3 2 1 3 2 4 2 4 1', 'center', 'center', white);
                vbl= Screen('Flip', window);
                PsychPortAudio('Start', h_Metronome600, 1, [], []); % Play metronome sound file
                
                Screen('TextSize', window, 25);
                DrawFormattedText(window, 'The time to practice the new sequence is over. \n Press any key to continue to the foot stomping experiment.', 'center', 'center', white);
                vbl= Screen('Flip', window);
                KbStrokeWait; %wait for response to terminate instructions
                
                %Instructions non-automatic foot stomping task
                Screen('TextSize',window,25);
                DrawFormattedText(window,'You will now perform the foot stomping experiment for the sequence you just practiced (non-automatic). \n There is a total of 10 trials, so performing the sequence 10 times. \n In between each trial there is a rest period of 20 seconds. \n During this rest you will hear a metronome sound, stomp the sequence according to this interval sound. \n The trials and rest periods are indicated with Trial and Rest showing on the screen. \n\n When the experiment starts you cannot talk anymore. \n Furthermore, it is important to stay still except for your right leg. \n If you press any key, the experiment starts right away. \n It will start with a rest period. \n Whenever Trial appears on the screen, you should start stomping the sequence. \n Press any key to start the foot stomping experiment.','center', 'center', white);
                vbl = Screen('Flip', window);
                KbStrokeWait; %wait for response to terminate instructions
                
                %Stimulus for foot stomping non-automatic sequence
                for n = 1:10 %Will perform 10 trials
                    
                    % Rest period 20-25 seconds
                    outlet.push_sample({rest})
                    Screen('TextSize', window, 36);
                    Screen('TextSize',window,40);
                    DrawFormattedText(window, 'Rest \n\n\n\n', 'center', 'center', white)
                    Screen('DrawLines', window, allCoords,...
                        lineWidthPix, white, [xCenter yCenter], 2);
                    Screen('Flip', window);
                    PsychPortAudio('Start', h_Metronome8, 1, [], []); % Play metronome sound file (8 seconds)
                    WaitSecs(t1+randi(t2))
                    
                    %Fixation cross
                    Screen('TextSize', window, 36);
                    Screen('TextSize',window,40);
                    DrawFormattedText(window, 'Trial \n\n\n\n', 'center', 'center', white)
                    Screen('DrawLines', window, allCoords,...
                        lineWidthPix, white, [xCenter yCenter], 2);
                    %draw on screen
                    Screen('Flip', window)
                    WaitSecs(t3)
                    
                end
                
                % End of the non-automatic foot stomping task
                Screen('TextSize', window, 25);
                DrawFormattedText(window, 'This is the end of the non-automatic foot stomping task. \n Press any key to continue with the automaticity test.', 'center', 'center', white);
                vbl= Screen('Flip', window);
                KbStrokeWait; %wait for response to terminate instructions
                
                
                % Instruction automaticity task foot stomping
                Screen('TextSize',window,25);
                DrawFormattedText(window, 'You will now perform the sequence you just learned for an automaticity test for the FOOT stomping task. \n  Letters will be shown on the screen (A,G,O,L) while you perform the task. \n The goal is to perform the sequence stomping while counting how many times G is presented. \n After each time you stomped the full sequence, you should tell us how many times G was presented. \n We will perform 8 trials. \n\n Note that during the stomping task you cannot talk. \n Try to keep your body movements as still as possible exept for your right leg. \n\n In between the trials you will see a fixation cross for 10 seconds. \n During the first 5 seconds you will hear a metronome sound. \n Stomp the sequence on this rhythm, which is the same as you studied at home. \n\n We will start with a fixation cross on the screen for 10 seconds. \n After that the first trial will start automatically. \n So start stomping the sequence as soon as a letter on the screen appears. \n Press any key to continue and start the test.','center','center', white);
                vbl = Screen('Flip', window);
                KbStrokeWait; %wait for response to terminate instructions
                
                for r=1:N_trials
                    %Presentation of the letters on the screen (dual task). -> is random.
                    %Participant has to count the amount that G was presented.
                    Letterlist= {'A', 'G', 'O', 'L'};
                    letter_order=randi(length(Letterlist), 1, N_letters);
                    letters_footnonautotest(r).presented =Letterlist(letter_order);
                    
                    % always start with a fixation cross and 5 seconds of metronome
                    % sound
                    Screen('TextSize', window, 36);
                    Screen('DrawLines', window, allCoords,...
                        lineWidthPix, white, [xCenter yCenter], 2);
                    Screen('Flip', window);
                    PsychPortAudio('Start', h_Metronome5, 1, [], []); % Play metronome sound file
                    WaitSecs(10);
                    
                    %Presentation of random letters on the screen during the foot
                    %stomping test
                    outlet.push_sample({foot_nonautotest})
                    for w=1:N_letters
                        % present random letter
                        Screen('TextSize', window, 100);
                        DrawFormattedText(window, [cell2mat(letters_footnonautotest(r).presented(w))],'center','center', white);
                        vbl = Screen('Flip', window);
                        time_letter=rand(1)+0.5; %Speed with which the letters are presented = A randomized value between 0 and 1, + 0.5 sec
                        WaitSecs(time_letter);
                        
                        % between each letter red fixation cross
                        Screen('DrawLines', window, allCoords,...
                            lineWidthPix, [1 0 0], [xCenter yCenter], 2);
                        Screen('Flip', window);
                        WaitSecs (0.2);
                    end
                    
                    % present white fixation cross for some seconds to show that
                    % trial is over
                    Screen('TextSize', window, 36);
                    Screen('DrawLines', window, allCoords,...
                        lineWidthPix, white, [xCenter yCenter], 2);
                    Screen('Flip', window);
                    WaitSecs(2);
                    
                    % show feedback
                    % ask how many G's were presented
                    Screen('TextSize',window,30);
                    DrawFormattedText(window, 'How many times was G presented? ','center','center', white);
                    vbl = Screen('Flip', window);
                    [secs, keyCode, deltaSecs]=KbWait;
                    letters_footnonautotest(r).reported_G=KbName(find(keyCode));
                    DrawFormattedText(window, ['Your answer: ' letters_footnonautotest(r).reported_G '\n Press any key to continue.'],'center','center', white);
                    vbl = Screen('Flip', window);
                    KbStrokeWait; %wait for response to terminate instruction
                    % Show if the reported number of Gs is correct
                    if str2num(letters_footnonautotest(r).reported_G)==sum(strcmp(letters_footnonautotest(r).presented, 'G'))
                        DrawFormattedText(window, ['This is correct. \n Press any key to continue.' ],'center','center', white);
                        vbl = Screen('Flip', window);
                    else
                        DrawFormattedText(window, ['This is incorrect. \n Press any key to continue.' ],'center','center', white);
                        vbl = Screen('Flip', window);
                    end
                    KbStrokeWait;
                    DrawFormattedText(window, 'Press any key to continue with the next trail. \n Note that you will first start with a fixation cross again. \n Start tapping the sequence as soon as a letter on the screen appears.' ,'center','center', white);
                    vbl = Screen('Flip', window);
                    KbStrokeWait;
                end
                
                % After all trials completed, the end of the foot stomping task is reached.
                Screen('TextSize',window,25);
                DrawFormattedText(window, 'End of the automaticity test for the non-automatic foot stomping task. \n Press any key to continue with the experiment.','center','center', white);
                vbl = Screen('Flip', window);
                save('letters_footnonautotest.mat', 'letters_footnonautotest'); % save the letters that were presented and the reported number of g's
                KbStrokeWait; %wait for response to terminate instructions
            end
        end
    end
    
    
    
    % AUTOMATICITY TASKS
    if i==2; %Start with the automatic tasks
        outlet.push_sample({auto})
        
        % Show instructions
        Screen('TextSize', window, 25);
        DrawFormattedText(window, 'AUTOMATICITY TASK \n\n You will perform the finger tapping and foot stomping task for the sequence you studied at home. \n You will either start with the finger tapping or foot stomping task. \n Press any key to continue.', 'center', 'center', white);
        vbl= Screen('Flip', window);
        KbStrokeWait; %wait for response to terminate instructions
        
        %Determine to start with either hand or foot task.
        for k=order_auto; %Either [5,6] or [6,5] -> determines the order of the limbs
            outlet.push_sample({instructions_auto})
            
            % Automatic finger tapping experiment
            if k==5;
                outlet.push_sample({finger_auto})
                %Automatic finger tapping task instructions
                Screen('TextSize',window,25);
                DrawFormattedText(window, 'You will now perform the FINGER tapping task for the sequence you studied at home (automatic). \n There is a total of 10 trials, so performing the sequence 10 times. \n  In between each trial there is a rest period of 20 seconds. \n During this rest you will hear a metronome sound, tap the sequence according to this interval sound. \n The trials and rest periods are indicated with Trial and Rest showing on the screen. \n\n When the experiment starts you cannot talk anymore. \n Furthermore, it is important to stay still except for your right hand. \n If you press any key, the experiment starts right away. \n It will start with a rest period. \n Whenever Trial appears on the screen, you should start tapping the sequence. \n Press any key to start the finger tapping experiment.','center','center', white);
                vbl = Screen('Flip', window);
                KbStrokeWait;
                
                %Stimulus for finger tapping automatic sequence
                for o = 1:10 % Will perform 10 trials
                    presses_handauto(o).key=[];
                    %Rest period 20-25 seconds
                    outlet.push_sample({rest})
                    Screen('TextSize', window, 36);
                    Screen('TextSize',window,40);
                    DrawFormattedText(window, 'Rest \n\n\n\n\n', 'center', 'center', white)
                    Screen('DrawLines', window, allCoords,...
                        lineWidthPix, white, [xCenter yCenter], 2);
                    Screen('Flip', window);
                    PsychPortAudio('Start', h_Metronome8, 1, [], []); % Play metronome sound file
                    WaitSecs(t1+randi(t2))
                    
                    % Fixation cross
                    Screen('TextSize', window, 36);
                    Screen('TextSize',window,40);
                    DrawFormattedText(window, 'Trial \n\n\n\n', 'center', 'center', white)
                    Screen('DrawLines', window, allCoords,...
                        lineWidthPix, white, [xCenter yCenter], 2);
                    %draw on screen
                    Screen('Flip', window)
                    p=1;
                    
                    % Record key presses
                    start_timer=GetSecs;
                    while GetSecs-start_timer < t3
                        if length(presses_handauto(o).key)<12
                            [secs, keyCode, deltaSecs] = KbWait([],2);
                            if any(keyCode)
                                fprintf('the key value was: %s \n', KbName(find(keyCode)));
                                presses_handauto(o).key{p}=KbName(find(keyCode));
                                presses_handauto(o).secs(p)=GetSecs;
                                p=p+1;
                            elseif KbName('ESCAPE')
                            elseif GetSecs-start_timer >= t3
                                break
                            end
                        end
                    end
                    
                end
                
                %End of the automatic finger tapping task
                Screen('TextSize', window, 25);
                DrawFormattedText(window, 'This is the end of the automatic finger tapping task. \n Press any key to continue with the rest of the experiment.', 'center', 'center', white);
                vbl= Screen('Flip', window);
                KbStrokeWait; %wait for response to terminate instructions
                
            end
            
            
            % Automatic foot stomping task
            if k==6;
                outlet.push_sample({foot_auto})
                
                % Automatic foot stomping instructions
                Screen('TextSize',window,25);
                DrawFormattedText(window,'You will now perform the FOOT stomping experiment for the sequence you studied at home (automatic). \n There is a total of 10 trials, so performing the sequence 10 times. \n In between each trial there is a rest period of 20 seconds. \n During this rest you will hear a metronome sound, stomp the sequence according to this interval sound. \n The trials and rest periods are indicated with Trial and Rest showing on the screen. \n\n When the experiment starts you cannot talk anymore. \n Furthermore, it is important to stay still except for your right leg. \n If you press any key, the experiment starts right away. \n It will start with a rest period. \n Whenever Trial appears on the screen, you should start stomping the sequence. \n Press any key to start the foot stomping experiment.', 'center', 'center', white);
                vbl = Screen('Flip', window);
                KbStrokeWait;
                
                %Stimulus for foot stomping automatic sequence
                for q = 1:10 % Will perform 10 trials
                    
                    % Rest period 20-25 seconds
                    outlet.push_sample({rest})
                    %fixation cross
                    Screen('TextSize', window, 36);
                    Screen('TextSize',window,40);
                    DrawFormattedText(window, 'Rest \n\n\n\n', 'center', 'center', white)
                    Screen('DrawLines', window, allCoords,...
                        lineWidthPix, white, [xCenter yCenter], 2);
                    Screen('Flip', window);
                    PsychPortAudio('Start', h_Metronome8, 1, [], []); % Play metronome sound file
                    WaitSecs(t1+randi(t2))
                    
                    %Fixation cross
                    Screen('TextSize', window, 36);
                    Screen('TextSize',window,40);
                    DrawFormattedText(window, 'Trial \n\n\n\n', 'center', 'center', white)
                    Screen('DrawLines', window, allCoords,...
                        lineWidthPix, white, [xCenter yCenter], 2);
                    Screen('Flip', window)
                    WaitSecs (t3)
                    
                end
                
                % End of the automatic foot stomping task.
                Screen('TextSize', window, 25);
                DrawFormattedText(window, 'This is the end of the automatic foot stomping task. \n Press any key to continue with the rest of the experiment.', 'center', 'center', white);
                vbl= Screen('Flip', window);
                KbStrokeWait; %wait for response to terminate instructions
                
            end
        end
    end
    
end

% Save data for each task and experiment
save('presses_handauto.mat', 'presses_handauto');
save('presses_handnonauto.mat', 'presses_handnonauto');

% End of the experiment, thank the participant
Screen('TextSize',window,30);
DrawFormattedText(window,'This is the end of the experiment, thank you for participating.', 'center', 'center', white);
vbl = Screen('Flip', window);

KbStrokeWait;
sca