pb_vReadExp.m 8.27 KB
Newer Older
1
function [block,cfg] = pb_vReadExp(cfg)
Jesse Heckman's avatar
Jesse Heckman committed
2
% PB_VREADEXP(expfile)
Jesse Heckman's avatar
Jesse Heckman committed
3
%
Jesse Heckman's avatar
Jesse Heckman committed
4
5
% PB_VREADEXP() reads experimental data from expfile and loads BLOCK and
% CFG parameters.
Jesse Heckman's avatar
Jesse Heckman committed
6
%
Jesse Heckman's avatar
Jesse Heckman committed
7
% See also PB_VRUNEXP
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
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
   cfg.comment = checkcomment(fid);
   cfg         = hread(fid,cfg);
   n           = cfg.Blocks;
   
   %% EXP
   
   block      = struct([]);
   block(1).signal                     = [];
   block(1).trial                      = [];

   bn = 0;
   while ~feof(fid)
      tline       = fgetl(fid);
      firstCell   = sscanf(tline,'%s',1);
      nchar				= length(firstCell);
      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);

            
         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);
         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
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
      
   
   end
   cfg		= pb_vLookup(cfg);
   
   %% EXTRA:  Azimuth Elevation
   for iBlck = 1:cfg.Blocks
      ntrials = length(block(iBlck).trial);
      for iTrl = 1:ntrials % for every trial
         s			= block(iBlck).trial(iTrl).stim;
         block(iBlck).trial(iTrl).nstim = numel(s); % number of stimuli per trial
         for stmIdx	= 1:block(iBlck).trial(iTrl).nstim % for every stimulus in a trial
            X			= block(iBlck).trial(iTrl).stim(stmIdx).X;
            Y			= block(iBlck).trial(iTrl).stim(stmIdx).Y;
            mod			= block(iBlck).trial(iTrl).stim(stmIdx).modality;
            if ~isempty(X) % for every stimulus that has an X and Y parameter, determine azimuth and elevation
               if cfg.Lab==1 % Hoop lab
                  if strcmpi(mod,'sky')
                     [Az,El] = hoopsky2azel(X,Y);
                  else
                     [Az,El]	= hoopXY2azel(X,Y);
                  end
               elseif ismember(cfg.Lab,[2 3]) % Sphere lab
                  channel = cfg.interpolant(X,Y);
                  Az		= cfg.lookup(channel+1,5);
                  El		= cfg.lookup(channel+1,6);
               elseif ismember(cfg.Lab,4) % SphereMinor lab
                  channel = cfg.interpolant(X,Y);
                  Az		= cfg.lookup(channel+1,4);
                  El		= cfg.lookup(channel+1,5);
               elseif ismember(cfg.Lab,5) % vPrime lab
                  channel = cfg.interpolant(X,Y);
                  Az		= cfg.lookup(channel+1,4);
                  El		= cfg.lookup(channel+1,5);
               end
               block(iBlck).trial(iTrl).stim(stmIdx).azimuth		= Az;
               block(iBlck).trial(iTrl).stim(stmIdx).elevation    = El;
            end
         end
      end
Jesse Heckman's avatar
Jesse Heckman committed
158
   end
159
   
Jesse Heckman's avatar
Jesse Heckman committed
160
   %% CHECK OUT
161
	fclose(fid);
Jesse Heckman's avatar
Jesse Heckman committed
162
163
end

Jesse Heckman's avatar
Jesse Heckman committed
164
165
166
167
168
169
170
171
172
173
174
175
176
function comment = checkcomment(fid)
   isComment	= true;
   cnt			= 0; % counter
   comment		= cell(1);
   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
177
   end
Jesse Heckman's avatar
Jesse Heckman committed
178
179
180
181
182
183
184
   fseek(fid,position,'bof');
end

function cfg = hread(fid,cfg)
   cnt = 0;
   isBody = false;
   header = cell(1);
Jesse Heckman's avatar
Jesse Heckman committed
185
   
Jesse Heckman's avatar
Jesse Heckman committed
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
   while ~isBody
      position	= ftell(fid);
      str      = fscanf(fid,'%s',1);
      isBody   = strcmp(str,'~~~');  
      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
205
206
end

Jesse Heckman's avatar
Jesse Heckman committed
207
function signal = readVest(line)
Jesse Heckman's avatar
Jesse Heckman committed
208
   types = {'none','sine','noise','turn'};
Jesse Heckman's avatar
Jesse Heckman committed
209
210
211
212
213
214
215
   n = str2num(erase(sscanf(line,'%s',2),sscanf(line,'%s',1)))+1;
   type  = types{n};
   signal.type = type;
   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)));
   end
Jesse Heckman's avatar
Jesse Heckman committed
216
end
Jesse Heckman's avatar
Jesse Heckman committed
217

Jesse Heckman's avatar
Jesse Heckman committed
218
219
220
221
222
223
224
% % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % 
%                                                           %
%       Part of Programmeer Beer Toolbox (PBToolbox)        %
%       Written by: Jesse J. Heckman (2018)                 %
%                                                           %
% % % % % % % % % % % % % % % % % % % % % % % % % % % % % % %