== ***** == == 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