forked from arokem/surround_suppression
-
Notifications
You must be signed in to change notification settings - Fork 0
/
ss_run.py
199 lines (168 loc) · 8.05 KB
/
ss_run.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
"""
Surround suppression experiment, based on Zenger-Landolt and Heeger (2003)
And on the Psychtoolbox version used to get the data in Yoon et al. (2009 and
2010).
"""
import gc
import wx
import numpy as np
from ss_classes import Params
from psychopy import core,event
import psychopy.monitors.calibTools as calib
"""
Insert Marker
"""
#Allow to use commands from C library(EDC02032013)
from ctypes import *#(EDC02032013)
import os#(EDC02032013)
#Load the ViewPoint library#(EDC02032013)
vpxDylib = '/usr/local/lib/libvpx_interapp.dylib'#(EDC02032013)
vpx=cdll.LoadLibrary(vpxDylib)#(EDC02032013)
#vpx = CDLL( vpxDylib );#(EDC02032013)#not necessary on Mac
#Change this to match ethernet server and port#(EDC02032013)
vpx.VPX_ConnectToViewPoint('152.79.52.61',5000)#function from the vpx.h file #(EDC02032013)
#Allows synchrously insert marker
vpx.VPX_SendCommand('dataFile_AsynchMarkerData No');#(EDC02032013)
#Send the marker from the list
#not sure if the list is useful any more if changed to detection task
marker_list=["1","2","3","4"];#(EDC02032013)
"""
"""
#This brings in all of the classes defined in ss_classes:
from ss_classes import *
from ss_tools import start_data_file
if __name__ == "__main__":
""" The main function. This actually runs the experiment """
#Initialize params from file:
params = Params()
#For some reason, if this call is inside ss_classes, just importing
#ss_classes starts an instance of the GUI, so we put it out here:
app = wx.App()
app.MainLoop()
params.set_by_gui()
f = start_data_file(params.subject)
#Start by saving in the parameter setting:
params.save(f)
#For now, assume that the and the annulus are going to have the same
#orientation:
params.target_ori = params.annulus_ori
#This initializes the window:
#This is the folder where we keep the calibration data:
calib.monitorFolder = './calibration/'# over-ride the usual setting of
# where monitors are stored
mon = calib.Monitor(params.monitor) #Get the monitor object and pass that
#as an argument to win:
win = visual.Window(monitor=mon,
screen=params.screen,
fullscr=params.fullscreen,
units=params.display_units)
win.setMouseVisible(False)
#Initialize all the instances of stimulus bits and pieces for reuse:
bank = StimulusBank(win,params)
#Make a trial list:
trial_list = make_trial_list(win,params)
#Initialize the staircase, depending on which task is performed
if params.task == 'Annulus':
message = """ Is there a target ?\n Press 1 for yes, 2 for no,\n Press any key to start"""
if params.surround_ori == params.annulus_ori:
staircaseA = Staircase(params.start_target_contrastA,
params.annulus_contrast/params.contrast_increments,
harder = -1, #For this task, higher values are
#actually harder => closer to the annulus
#value
ub=params.targetA_contrast_max,
lb=params.targetA_contrast_min)
staircaseB = Staircase(params.start_target_contrastB,
params.start_target_contrastB/params.contrast_increments,
harder = -1, #For this task, higher values are
#actually harder => closer to the annulus
#value
ub=params.targetB_contrast_max,
lb=params.targetB_contrast_min)
else:
staircaseA = Staircase(params.start_target_orthog_contrastA,
params.annulus_contrast/params.contrast_increments,
harder = -1, #For this task, higher values are
# #actually harder => closer to the annulus
#value
ub=params.targetA_contrast_max,
lb=params.targetA_contrast_min)
staircaseB = Staircase(params.start_target_contrastB,
params.start_target_contrastB/params.contrast_increments,
harder = -1, #For this task, higher values are
#actually harder => closer to the annulus
#value
ub=params.targetB_contrast_max,
lb=params.targetB_contrast_min)
if params._replay is None:
#The fixation target appears but has a constant contrast set to the
#starting point
other_contrast = np.ones(len(trial_list)) * params.fix_target_start
else:
#Replay a previous run
other_contrast = params._replay
elif params.task == 'Fixation':
message = """ Is there a target ?\n Press 1 for yes, 2 for no,\n Press any key to start"""
#Just one staircase:
staircaseA = staircaseB = Staircase(params.fix_target_start,
params.fix_target_start/params.contrast_increments,
harder = -1,
ub=params.fix_target_max,
lb=params.fix_target_min
)
if params._replay is None:
#The annulus target appears and has a constant contrast set to the
#starting point:
other_contrast = np.ones(len(trial_list)) * params.fix_target_start
else:
#Replay a previous run:
other_contrast = params._replay
#Send a message to the screen and wait for a subject keypress:
Text(win,text=message,height=0.7)()
win.flip()
#If this is in the scanner, we would want to wait for the ttl pulse right
#here:
if params.scanner:
Text(win,text='',keys=['t'])() #Assuming a TTL is a 't' key
Beginning(win,params,bank)
#Loop over the event list, while consuming each event, by calling it:
for trial_idx,this_trial in enumerate(trial_list):
#Send a marker to the EyeTracker
#theMarker='T'#(EDC02032013)
vpx.VPX_SendCommand('dataFile_InsertMarker S')#(EDC02032013)
vpx.VPX_SendCommand('say "Inserting Marker S"')#(EDC02032013)
trial_clock = core.Clock()
if this_trial.block_type == 'A':
this_trial.finalize_stim(params,bank,staircaseA,other_contrast[trial_idx])
elif this_trial.block_type == 'B':
this_trial.finalize_stim(params,bank,staircaseB,other_contrast[trial_idx])
this_trial.stimulus(params,this_trial.block_type)
if this_trial.block_type == 'A':
this_trial.finalize_fix(params,bank,staircaseA,other_contrast[trial_idx])
if this_trial.block_type == 'B':
this_trial.finalize_fix(params,bank,staircaseB,other_contrast[trial_idx])
this_trial.fixation(params,this_trial.block_type)
this_trial.response.finalize(correct_key = this_trial.correct_key,
file_name=f)
this_trial.response()
if this_trial.correct_key is not None:
this_trial.feedback.finalize(this_trial.response.correct)
this_trial.feedback()
#On the first trial, insert the header:
if trial_idx == 0:
this_trial.save(f,insert_header=True)
else:
#On other trials, just insert the data:
this_trial.save(f)
#update after saving:
if this_trial.block_type == 'A':
staircaseA.update(this_trial.response.correct)
elif this_trial.block_type == 'B':
staircaseB.update(this_trial.response.correct)
this_trial.wait_iti(trial_clock)
#Flush unwanted events that may still be hanging out:
event.clearEvents()
#dbg
# trial_clock.getTime()
f.close()
core.quit()