pb_vReadExp.m 6.57 KB
Newer Older
Jesse Heckman's avatar
Jesse Heckman committed
1
2
function [block,cfg] = pb_vReadExp(expfile)
% 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
13
14
15

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

   %% INITIALIZE
   
   if ~pb_fexist(expfile); return; end
   fid = fopen(expfile,'r');
   
Jesse Heckman's avatar
Jesse Heckman committed
16
   %% HEADER
Jesse Heckman's avatar
Jesse Heckman committed
17
   
Jesse Heckman's avatar
Jesse Heckman committed
18
   cfg         = struct; 
Jesse Heckman's avatar
Jesse Heckman committed
19
   
Jesse Heckman's avatar
Jesse Heckman committed
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
   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

   end
Jesse Heckman's avatar
Jesse Heckman committed
121
122
123
124
   %% CHECK OUT
   fclose(fid);
end

Jesse Heckman's avatar
Jesse Heckman committed
125
126
127
128
129
130
131
132
133
134
135
136
137
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
138
   end
Jesse Heckman's avatar
Jesse Heckman committed
139
140
141
142
143
144
145
   fseek(fid,position,'bof');
end

function cfg = hread(fid,cfg)
   cnt = 0;
   isBody = false;
   header = cell(1);
Jesse Heckman's avatar
Jesse Heckman committed
146
   
Jesse Heckman's avatar
Jesse Heckman committed
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
   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
166
167
end

Jesse Heckman's avatar
Jesse Heckman committed
168
function signal = readVest(line)
Jesse Heckman's avatar
Jesse Heckman committed
169
   types = {'none','sine','noise','turn'};
Jesse Heckman's avatar
Jesse Heckman committed
170
171
172
173
174
175
176
   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
177
end
Jesse Heckman's avatar
Jesse Heckman committed
178

Jesse Heckman's avatar
Jesse Heckman committed
179
180
181
182
183
184
185
% % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % 
%                                                           %
%       Part of Programmeer Beer Toolbox (PBToolbox)        %
%       Written by: Jesse J. Heckman (2018)                 %
%                                                           %
% % % % % % % % % % % % % % % % % % % % % % % % % % % % % % %