pb_vReadExp.m 8.75 KB
Newer Older
1
function [block,cfg] = pb_vReadExp(cfg)
2
% PB_VREADEXP(CFG)
Jesse Heckman's avatar
Jesse Heckman committed
3
%
4
5
% PB_VREADEXP(CFG) reads experimental data from expfile and loads it in 
% 'block' and 'cfg'.
Jesse Heckman's avatar
Jesse Heckman committed
6
%
7
% See also PB_VPRIME, PB_VPRIMEGUI, PB_VRUNEXP, PB_GETBLOCK
Jesse Heckman's avatar
Jesse Heckman committed
8
9
10
11
12

% PBToolbox (2018): JJH: j.heckman@donders.ru.nl

   %% INITIALIZE
   
13
   expfile = cfg.expfname;
Jesse Heckman's avatar
Jesse Heckman committed
14
15
16
   if ~pb_fexist(expfile); return; end
   fid = fopen(expfile,'r');
   
Jesse Heckman's avatar
Jesse Heckman committed
17
   %% HEADER
Jesse Heckman's avatar
Jesse Heckman committed
18
   
Jesse Heckman's avatar
Jesse Heckman committed
19
20
21
22
23
24
   cfg.comment = checkcomment(fid);
   cfg         = hread(fid,cfg);
   n           = cfg.Blocks;
   
   %% EXP
   
25
26
27
   block                	= struct([]);
   block(1).signal        	= [];
   block(1).trial        	= [];
Jesse Heckman's avatar
Jesse Heckman committed
28
29
30

   bn = 0;
   while ~feof(fid)
31
32
33
      
      tline          = fgetl(fid);
      firstCell      = sscanf(tline,'%s',1);
Jesse Heckman's avatar
Jesse Heckman committed
34
      nchar				= length(firstCell);
35
      
Jesse Heckman's avatar
Jesse Heckman committed
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
      if strcmp(firstCell,'%'); firstCell = sscanf(tline,'%s',2); firstCell = firstCell(2:end-1); end

      switch upper(firstCell)
         case 'BLOCK'
            % Updates counts for block
            bn = bn+1; tn = 0; sn = 0;
            fgetl(fid); 
         
         case 'TRIAL'
            % Updates counts for trial
            tn = tn+1; sn = 0;
            
         case 'HOR'
            % Writes horizontal vestibular signal for block
            block(bn).signal.hor = readVest(tline); 
            
         case 'VER'
            % Writes vertical vestibular signal for block
            block(bn).signal.ver = readVest(tline); 
            
         case {'SND','SND1','SND2'}
            % Updates count and writes stimulus
            sn    = sn+1;
            par   = sscanf(tline(nchar+1:end),'%d%d%d%f%d%d',[6,1]);
            
            block(bn).trial(tn).stim(sn).modality		= 'sound';
            block(bn).trial(tn).stim(sn).X            = par(1);
            block(bn).trial(tn).stim(sn).Y            = par(2);
            block(bn).trial(tn).stim(sn).matfile		= ['snd' num2str(par(3),'%03i') '.mat']; % for sphere
            block(bn).trial(tn).stim(sn).wavfile		= ['snd' num2str(par(3),'%03i') '.wav']; % double, for hoop
            block(bn).trial(tn).stim(sn).parameters	= par(3); % for sphereMinor

            block(bn).trial(tn).stim(sn).matfile		= ['snd' num2str(par(3),'%03i') '.mat'];
            block(bn).trial(tn).stim(sn).wavfile		= ['snd' num2str(par(3),'%03i') '.wav'];	 % double

            block(bn).trial(tn).stim(sn).intensity    = par(4);
            block(bn).trial(tn).stim(sn).onevent		= par(5);
            block(bn).trial(tn).stim(sn).ondelay		= par(6);
            block(bn).trial(tn).stim(sn).offevent		= par(5); % default duration
            block(bn).trial(tn).stim(sn).offdelay		= par(6)+150; % default duration
            
         case {'LED','LED1','LED2'}
            % Updates count and writes stimulus
            sn    = sn+1;
            par   = sscanf(tline(nchar+1:end),'%d%d%d%d%d%d%d',[7,1]);
				
            block(bn).trial(tn).stim(sn).modality		= 'LED';
				block(bn).trial(tn).stim(sn).X            = par(1);
				block(bn).trial(tn).stim(sn).Y            = par(2);

				block(bn).trial(tn).stim(sn).intensity    = par(3);

				block(bn).trial(tn).stim(sn).intensity    = par(3);  % hoop: range 0-255, sphere range 1-50

				block(bn).trial(tn).stim(sn).onevent		= par(4);
				block(bn).trial(tn).stim(sn).ondelay		= par(5);
				block(bn).trial(tn).stim(sn).offevent		= par(6);
				block(bn).trial(tn).stim(sn).offdelay		= par(7);
94
 
Jesse Heckman's avatar
Jesse Heckman committed
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
         case 'TRG0'
            % Updates count and writes stimulus
            sn    = sn+1;
            par	= sscanf(tline(nchar+1:end),'%d%d%d%d%d',[5,1]);
            
            block(bn).trial(tn).stim(sn).modality		= 'trigger';
            if par(1) == 1
               block(bn).trial(tn).stim(sn).detect		= 'rise';
            elseif par==2
               block(bn).trial(tn).stim(sn).detect		= 'fall';
            end
            block(bn).trial(tn).stim(sn).channel		= par(2);
            block(bn).trial(tn).stim(sn).onevent		= par(3);
            block(bn).trial(tn).stim(sn).ondelay		= par(4);
            block(bn).trial(tn).stim(sn).event        = par(5);
110
            
Jesse Heckman's avatar
Jesse Heckman committed
111
112
113
114
115
116
117
118
119
         case 'ACQ'
            % Updates count and writes stimulus
            sn    = sn+1;
            par	= sscanf(tline(nchar+1:end),'%d%d',[2,1]); % could also be 3
            
            block(bn).trial(tn).stim(sn).modality		= 'data acquisition';
            block(bn).trial(tn).stim(sn).onevent		= par(1);
            block(bn).trial(tn).stim(sn).ondelay		= par(2);
      end
120
121
   end
   cfg		= pb_vLookup(cfg);
122

123
   for iBlck = 1:cfg.Blocks
124
      
125
126
      ntrials = length(block(iBlck).trial);
      for iTrl = 1:ntrials % for every trial
127
128
         
         s = block(iBlck).trial(iTrl).stim;
129
         block(iBlck).trial(iTrl).nstim = numel(s); % number of stimuli per trial
130
131
132
133
134
135
         for iStm	= 1:block(iBlck).trial(iTrl).nstim % for every stimulus in a trial
            
            X			= block(iBlck).trial(iTrl).stim(iStm).X;
            Y			= block(iBlck).trial(iTrl).stim(iStm).Y;
            mod			= block(iBlck).trial(iTrl).stim(iStm).modality;
            
136
137
            if ~isempty(X) % for every stimulus that has an X and Y parameter, determine azimuth and elevation
               if cfg.Lab==1 % Hoop lab
138
                  
139
140
141
142
143
                  if strcmpi(mod,'sky')
                     [Az,El] = hoopsky2azel(X,Y);
                  else
                     [Az,El]	= hoopXY2azel(X,Y);
                  end
144
                  
145
               elseif ismember(cfg.Lab,[2 3]) % Sphere lab
146
                  
147
148
149
                  channel = cfg.interpolant(X,Y);
                  Az		= cfg.lookup(channel+1,5);
                  El		= cfg.lookup(channel+1,6);
150
                  
151
               elseif ismember(cfg.Lab,4) % SphereMinor lab
152
                  
153
154
155
                  channel = cfg.interpolant(X,Y);
                  Az		= cfg.lookup(channel+1,4);
                  El		= cfg.lookup(channel+1,5);
156
                  
157
               elseif ismember(cfg.Lab,5) % vPrime lab
158
                  
159
160
161
                  channel = cfg.interpolant(X,Y);
                  Az		= cfg.lookup(channel+1,4);
                  El		= cfg.lookup(channel+1,5);
162
                  
163
               end
164
165
               block(iBlck).trial(iTrl).stim(iStm).azimuth        = Az;
               block(iBlck).trial(iTrl).stim(iStm).elevation      = El;
166
167
168
            end
         end
      end
Jesse Heckman's avatar
Jesse Heckman committed
169
   end
170
	fclose(fid);
Jesse Heckman's avatar
Jesse Heckman committed
171
172
end

Jesse Heckman's avatar
Jesse Heckman committed
173
function comment = checkcomment(fid)
174
175
   % reads comments from the expfile

Jesse Heckman's avatar
Jesse Heckman committed
176
177
178
   isComment	= true;
   cnt			= 0; % counter
   comment		= cell(1);
179
   
Jesse Heckman's avatar
Jesse Heckman committed
180
181
182
183
184
185
186
187
188
   while isComment % do this for every line that starts with '%'
      position	= ftell(fid); % find the position in the file
      str			= fscanf(fid,'%s',1); % read the string (moving the position in the file)
      commentline = fgetl(fid); % get the entire line (again repositioning)
      isComment	= strncmp(str,'%',1); % and check whether the first string of the line actually indicated a comment
      if ~isempty(commentline) && isComment
         cnt				= cnt+1;
         comment(cnt)	= {commentline};
      end
Jesse Heckman's avatar
Jesse Heckman committed
189
   end
Jesse Heckman's avatar
Jesse Heckman committed
190
191
192
193
   fseek(fid,position,'bof');
end

function cfg = hread(fid,cfg)
194
195
   % reads the header from the expfile

Jesse Heckman's avatar
Jesse Heckman committed
196
197
198
   cnt = 0;
   isBody = false;
   header = cell(1);
Jesse Heckman's avatar
Jesse Heckman committed
199
   
Jesse Heckman's avatar
Jesse Heckman committed
200
   while ~isBody
201
      
Jesse Heckman's avatar
Jesse Heckman committed
202
203
      position	= ftell(fid);
      str      = fscanf(fid,'%s',1);
204
205
      isBody   = strcmp(str,'~~~');
      
Jesse Heckman's avatar
Jesse Heckman committed
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
      if ~isBody
         fseek(fid,position,'bof');
         cnt			= cnt+1;
         header(cnt)	= {fscanf(fid,'%s',1)};
         switch lower(header{cnt})
            case 'iti' % Inter-trial interval
               cfg.(header{cnt})	= fscanf(fid,'%d %d',[2 1]); % 2 integers: minimum and maximum possible inter trial interval
            case 'motor' % Motor
               cfg.(header{cnt})	= fscanf(fid,'%s',1); % String: yes or no
            otherwise
               cfg.(header{cnt})	= fscanf(fid,'%d',1); % Integer
         end
         checkcomment(fid);
      end  
   end
Jesse Heckman's avatar
Jesse Heckman committed
221
222
end

Jesse Heckman's avatar
Jesse Heckman committed
223
function signal = readVest(line)
224
225
   % reads VS signal from the expfile
   
Jesse Heckman's avatar
Jesse Heckman committed
226
   types = {'none','sine','noise','turn'};
Jesse Heckman's avatar
Jesse Heckman committed
227
228
229
   n = str2num(erase(sscanf(line,'%s',2),sscanf(line,'%s',1)))+1;
   type  = types{n};
   signal.type = type;
230
   
Jesse Heckman's avatar
Jesse Heckman committed
231
232
233
   if ~strcmp(type,'none')
      signal.amplitude  = str2double(erase(sscanf(line,'%s',3),sscanf(line,'%s',2)));
      signal.duration   = str2double(erase(sscanf(line,'%s',4),sscanf(line,'%s',3)));
234
235
236
237
238
      signal.frequency   = str2double(erase(sscanf(line,'%s',5),sscanf(line,'%s',4)));
   else 
      signal.amplitude  = 0;
      signal.duration   = 0;
      signal.frequency  = 0.1;
Jesse Heckman's avatar
Jesse Heckman committed
239
   end
Jesse Heckman's avatar
Jesse Heckman committed
240
end
Jesse Heckman's avatar
Jesse Heckman committed
241

Jesse Heckman's avatar
Jesse Heckman committed
242
243
244
245
246
247
248
% % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % 
%                                                           %
%       Part of Programmeer Beer Toolbox (PBToolbox)        %
%       Written by: Jesse J. Heckman (2018)                 %
%                                                           %
% % % % % % % % % % % % % % % % % % % % % % % % % % % % % % %