ScannerSyncMatlab - MRC CBU Imaging Wiki

Please enter your password of your account at the remote wiki below.
/!\ You should trust both wikis because the password could be read by the particular administrators.

Clear message
location: ScannerSyncMatlab

*****

NB - as of early 2015, the method of using scanner sync in matlab has changed. The new method uses Tibor's ScannerSynchClass rather than Rhodri's MRISync / ScannerSync software

*****

Using ScannerSync from Matlab - NEW METHOD

Open Matlab and check that the ScannerSynchClass is on the path, e.g. using:

which ScannerSynchClass

This should return the path to the ScannerSynchClass.m file. On the stim machines, this should be C:\Program Files\MATLAB\R2014a\toolbox\CBSU.

If you get an error message like "ScannerSynchClass not found", try:

addpath('C:\Program Files\MATLAB\R2014a\toolbox\CBSU')

Once the ScannerSynchClass is available, create an instance of the class:

SSO=ScannerSynchClass

Once the instance has been created, you can use the class properties and methods to collect pulses and button box responses:

== Properties (internal variables) ==
IsValid         = device set and operational
TR              = set a non-zero value (in seconds) for emulation mode (will not detect "real"          pulses)
Keys    = set a cell of key names for emulation mode. For key names look for KbName.m.
        N.B.: Requires PTB.
        N.B.: Suppress passing keypresses to MATLAB during the whole experiment

Clock   = interal clock (seconds past since the first scanner pulse or clock reset)
Synch   = current state of the scanner synch pulse
TimeOfLastPulse         = time (according to the internal clock) of the last pulse 
MeasuredTR      = estimated TR
SynchCount      = number of scanner synch pulses
MissedSynch     = number of missed scanner synch pulses

EmulSynch       = is scanner synch pulse emulated 
EmulButtons     = is button box emulated
%
Buttons                 = current state of the any button
LastButtonPress         = index/indices of the last button(s) pressed
TimeOfLastButtonPress   = time (according to the internal clock) of the last button press (any) 
BBoxTimeout     = set a non-Inf value (in seconds) to wait for button press only                               = set a negative value (in seconds) to wait even in case of response

== Methods (internal functions) ==
ScannerSynchClass       = constructor
delete  = destructor
ResetClock      = reset internal clock

ResetSynchCount         = reset scanner synch pulse counter
SetSynchReadoutTime(t)  = blocks scanner synch pulse readout after a pulse for 't' seconds
WaitForSynch    = wait until a scanner synch pulse arrives
CheckSynch(t)   = wait for a scanner synch pulse for 't' seconds or unitl a scanner synch pulse arrives (whichever first) and returns whether a scanner synch pulse was detected

        SetButtonReadoutTime(t)         = blocks individual button readout after a button press for 't' seconds (detection of other buttons is still possible) 
SetButtonBoxReadoutTime(t)      = blocks the whole button box readout after a button press for 't' seconds (detection of other buttons is also not possible) 
WaitForButtonPress                      = wait until a button is pressed
WaitForButtonRelease            = wait until a button is released



%
% USAGE
%
% Initialise:
%       SSO = ScannerSynchClass;
%       % SSO = ScannerSynchClass(1);   % emulate scanner synch pulse
%       % SSO = ScannerSynchClass(0,1); % emulate button box
%       % SSO = ScannerSynchClass(1,1); % emulate scanner synch pulse and button box
%
% Close:
%       SSO.delete;
%
% Example for scanner synch pulse #1: - Simple case
%       SSO.SetSynchReadoutTime(0.5);
%       SSO.TR = 2;                % allows detecting missing pulses
%       while SSO.SynchCount < 10  % polls 10 pulses
%       SSO.WaitForSynch;
%       fprintf('Pulse %d: %2.3f. Measured TR = %2.3fs\n',...
%           SSO.SynchCount,...
%           SSO.TimeOfLastPulse,...
%           SSO.MeasuredTR);
%       end
%
% Example for scanner synch pulse #2 - Chance for missing pulse
%       SSO.SetSynchReadoutTime(0.5);
%       SSO.TR = 2;                    % allows detecting missing pulses
%       while SSO.SynchCount < 10      % until 10 pulses
%       WaitSecs(Randi(100)/1000); % in every 0-100 ms ...
%       if SSO.CheckSynch(0.01)    % ... waits for 10 ms for a pulse
%               fprintf('Pulse %d: %2.3f. Measured TR = %2.3fs. %d synch pulses has/have been missed\n',...
%               SSO.SynchCount,...
%               SSO.TimeOfLastPulse,...
%               SSO.MeasuredTR,...
%               SSO.MissedSynch);
%       end
%       end
%
% Example for buttons:
%       SSO.SetButtonReadoutTime(0.5);      % block individual buttons
%       % SSO.SetButtonBoxReadoutTime(0.5); % block the whole buttonbox
%   % SSO.Keys = {'f1','f2','f3','f4'}; % emulation Buttons #1-#4 with F1-F4
%       n = 0;
%   % SSO.BBoxTimeout = 1.5;            % Wait for button press for 1.5s
%   % SSO.BBoxTimeout = -1.5;           % Wait for button press for 1.5s even in case of response
%       SSO.ResetClock;
%       while n ~= 10                       % polls 10 button presses
%       SSO.WaitForButtonPress;         % Wait for any button to be pressed
%       SSO.WaitForButtonRelease;       % Wait for any button to be     released
%       % SSO.WaitForButtonPress([],2); % Wait for Button #2
%       % SSO.WaitForButtonPress(2);    % Wait for any button for 2s (overrides SSO.BBoxTimeout only for this event)
%       % SSO.WaitForButtonPress(-2);   % Wait for any button for 2s even in case of response (overrides SSO.BBoxTimeout only for this event)
%       % SSO.WaitForButtonPress(2,2);  % Wait for Button #2 for 2s (overrides SSO.BBoxTimeout only for this event)
%       % SSO.WaitForButtonPress(-2,2); % Wait for Button #2 for 2s even in case of response (overrides SSO.BBoxTimeout only for this event)
%       n = n + 1;
%       fprintf('At %2.3fs, ',SSO.Clock);
%       fprintf('Button %d ',SSO.LastButtonPress);
%       fprintf('pressed: %2.3fs\n',SSO.TimeOfLastButtonPress);
%       end

Using ScannerSync from Matlab - OLD METHOD

Make sure you have the ScannerSync control installed.

Somewhere near the top of your code, set up ScannerSync & initialise communication with the input-output board (will not work if you're using a machine without the board, e.g., not the mimic or stim delivery machines).

%% SCANNERSYNC
TR=1000; % TR in ms
numdummies=8;

% Create & initialise scanner sync object
objSS=actxserver('MRISync.ScannerSync');  %Create a scanner object
invoke(objSS,'Initialize','') %Initialise the Keithley board and object
invoke(objSS,'SetMSPerSample',2); %Set StartExperiment routines’ sampling

Now, wait for first pulse from scanner. Also, tell ScannerSync the approximate TR.

% ScannerSync - waits for the first pulse, resets timer to 0, sets TR
invoke(objSS,'StartExperiment', double(TR));

Wait for dummies (one less as StartExperiment command above will already have heard a pulse)

for countdown=(numdummies-1):-1:1
    count_text=sprintf('%d', countdown);
  % uncomment the next two if you're using the PsychToolbox
%    DrawFormattedText(window, count_text, 'center', 'center', white);
%    Screen(window, 'Flip');
    invoke(objSS,'SynchroniseExperiment',1,0); %  1= force wait for actual pulse; 0=return this many ms after pulse
end;

In your trial loop, you'll need to spend some of the time listening for pulses

    invoke(objSS,'CheckPulseSynchronyForTime', double(500)); % spend 500 ms listening for any pulses

You should dump out the timing of critical events and the measured TR

    picstarttime=invoke(objSS,'SSGetTimer');
    measuredTR=invoke(objSS,'GetMeasuredTR');

% ...add your own code to write these to your output file

Also, for some designs you may want to occasionally synchronise your trials to the scanner

    invoke(objSS,'SynchroniseExperiment',1,0); %  1= force wait for actual pulse; 0=return this many ms after pulse

Using pretend mode

If you'd like to test code on a machine without the card, you may use a feature of ScannerSync called "pretend mode". To do this:

(1) When in pretend mode, don't do "Initialise", and instead issue a "SetPretendMode" command. So change the first block to something like

pretendmode=1;
% Create & initialise scanner sync object
objSS=actxserver('MRISync.ScannerSync');  %Create a scanner object
if (~pretendmode)
  invoke(objSS,'Initialize','') %Initialise the Keithley board and object
else
  invoke(objSS,'SetPretendMode',1)
end;
invoke(objSS,'SetMSPerSample',2); %Set StartExperiment routines’ sampling

(2) You'll then see a warning message when you get to the "StartExperiment" command, to which you should click "OK". If you're using the PsychToolbox and have already set up the screen, make sure that you only use the HideCursor command after the StartExperiment command, otherwise you won't be able to click OK.

If you would like to use the button box in the practice room of the Scanner building, you will need this code instead to make the button box work:

pretendmode=1;
% Create & initialise scanner sync object
objSS=actxserver('MRISync.ScannerSync');  %Create a scanner object
if (~pretendmode)
  invoke(objSS,'Initialize','') %Initialise the Keithley board and object
else
  invoke(objSS,'SetPretendModeExtended',1,0)
  invoke(objSS,'Initialize','') %Initialise the Keithley board and object
end;
invoke(objSS,'SetMSPerSample',2); %Set StartExperiment routines’ sampling

Collecting Responses

To collect responses in matlab using the button box add a code similar to this:

 tic
 resp=0;
 gotresp=false;
   while toc<rt_window
      if (~gotresp)
         resp=bitand(30,invoke(objSS,'GetResponse'));
         if (resp==samekey || resp==diffkey)
            gotresp=true;
         end;
       end;
    end;

This will collect the responses made by the subject during a pre-defined time interval (rt_window). You can also need to define which buttons the subject should press (in this example they are called samekey and diffkey). The codes returned by pressing the buttons are 28, 26, 22, 14 (i.e. when the right hand is placed on the button box and the subject used his index finger the code returned will be 28).

The GetResponse() function will by default only return first four button presses. If there is a need to use more than 4 buttons one could use GetResponseExtended() or ReadPIOValue(). The example below uses the latter function and ignores responses 254 and 255, which relate to no-response events.

 tic
 resp=0;
 gotresp=false;
 while toc<rt_window
    if (~gotresp)
       resp=invoke(objSS,'ReadPIOValue');
       if isempty(intersect([255,254],resp))
          gotresp=true;
       end;
    end;
 end;

Collecting Response Times

To collect RTs you could try the following:

resp=0;
gotresp=false;
respstarttime=invoke(objSS,'SSGetTimer');
tic
   while toc<rt_window
      if (~gotresp)
         resp=bitand(30,invoke(objSS,'GetResponse'));
         if (resp==samekey || resp==diffkey)
            rt=invoke(objSS,'SSGetTimer')-respstarttime; %or alternatively you can set rt=toc
            gotresp=true;
         end;
       end;
    end;