Validate Stim Names for STR

@@ -174,6 +174,9 @@ for i0 = 1:size(t,1)
             chMatchLog{end+1,1} = Sensors.label{i1};
             chMatchLog{end+1,1} = Sensors.label{i1};
             chMatchLog{end,2} = Name{iChanPos};
             chMatchLog{end,2} = Name{iChanPos};
             Sensors.label{i1} = Name{iChanPos};
             Sensors.label{i1} = Name{iChanPos};
+            if max(Sensors.chanpos(:))>1 && strcmp(Sensors.unit,'m')
+                Sensors.unit='mm';
+            end
         else
         else
             disp(['ImaGIN> WARNING: ' Sensors.label{i1} ' not assigned']);
             disp(['ImaGIN> WARNING: ' Sensors.label{i1} ' not assigned']);
             chNotFound{end+1} = Sensors.label{i1};
             chNotFound{end+1} = Sensors.label{i1};
@@ -215,16 +218,18 @@ for i0 = 1:size(t,1)
     end
     end
     
     
     % Match file: Correspondance between SEEG-CSV-LENA conventions
     % Match file: Correspondance between SEEG-CSV-LENA conventions
-    if isfield(S, 'FileTxtOut') && ~isempty(S.FileTxtOut)
-        try
-            fid = fopen(S.FileTxtOut,'w');
-            fprintf(fid,'SEEG,CSV\n');
-            for i = 1:size(chMatchLog,1)
-                fprintf(fid,'%s,%s\n', chMatchLog{i,1}, chMatchLog{i,2});
+    if exist('S','var')
+        if isfield(S, 'FileTxtOut') && ~isempty(S.FileTxtOut)
+            try
+                fid = fopen(S.FileTxtOut,'w');
+                fprintf(fid,'SEEG,CSV\n');
+                for i = 1:size(chMatchLog,1)
+                    fprintf(fid,'%s,%s\n', chMatchLog{i,1}, chMatchLog{i,2});
+                end
+                fclose(fid);
+            catch
+                disp('Log with matched channels names SEEG-CSV not saved.')
             end
             end
-            fclose(fid);
-        catch
-            disp('Log with matched channels names SEEG-CSV not saved.')
         end
         end
     end
     end
     
     

@@ -100,7 +100,7 @@ for i1=1:NName
         Z=ones(1,101)*Pos(1,3);
         Z=ones(1,101)*Pos(1,3);
     end
     end
     D=sqrt((X-Pos(1,1)).^2+(Y-Pos(1,2)).^2+(Z-Pos(1,3)).^2);
     D=sqrt((X-Pos(1,1)).^2+(Y-Pos(1,2)).^2+(Z-Pos(1,3)).^2);
-    NElec=ceil(max(D)/d);
+    NElec=max([ceil(max(D)/d) 1]);
     for i2=1:NElec
     for i2=1:NElec
         NElectrode=NElectrode+1;
         NElectrode=NElectrode+1;
         NameElectrode{NElectrode}=[NameElec{i1} num2str(i2)];
         NameElectrode{NElectrode}=[NameElec{i1} num2str(i2)];

@@ -261,6 +261,53 @@ if strcmpi(patientCode(5:end),'FRE')
         fprintf('\n \n MESSAGE: .. %s parameters updated ..::\n',patientCode); 
         fprintf('\n \n MESSAGE: .. %s parameters updated ..::\n',patientCode); 
     end
     end
 end
 end
+str_flag = 0;
+% STR patients have amplitude that varies
+if strcmpi(patientCode(5:end),'STR')
+    load('/gin/data/database/02-raw/stim_parameters-ftract-str.mat','stim_params')
+    Loc = find(ismember(stim_params.PCode, patientCode), 1);
+    if ~isempty(Loc)
+        str_flag = 1;
+        Frq = stim_params.Freq{Loc};
+%         Amp = stim_params.Ampl{Loc};
+        Pul = stim_params.Pulse{Loc};
+        for n = 1:length(KeepEvent)
+            xsub1 = regexp(Notes{KeepEvent(n)},rxp1,'match');
+            xsub2 = regexp(Notes{KeepEvent(n)},rxp2,'match');
+            xsub3 = regexp(Notes{KeepEvent(n)},rxp3,'match');
+            
+%             if ~isempty(xsub1) % insert stim Amplitude
+%                 Notes{KeepEvent(n)} = char(strrep(Notes{KeepEvent(n)},xsub1,Amp));
+%             else
+%                 Notes{KeepEvent(n)} = [Notes{KeepEvent(n)} '_' Amp];
+%             end
+%             
+%             if isempty(xsub1) % insert stim Amplitude
+%                 Notes{KeepEvent(n)} = [Notes{KeepEvent(n)} '_' Amp];
+%             end
+            
+            if ~isempty(xsub2) % insert stim Frequency
+                Notes{KeepEvent(n)} = char(strrep(Notes{KeepEvent(n)},xsub2,Frq));
+            else
+                Notes{KeepEvent(n)} = [Notes{KeepEvent(n)} '_' Frq];
+            end
+            
+            if ~isempty(xsub3) % insert stim Pulse
+                Notes{KeepEvent(n)} = char(strrep(Notes{KeepEvent(n)},xsub3,Pul));
+            else
+                Notes{KeepEvent(n)} = [Notes{KeepEvent(n)} '_' Pul];
+            end
+            Notes{KeepEvent(n)} = char(strrep(Notes{KeepEvent(n)},'_3_','_'));
+            evt(KeepEvent(n)).type = Notes{KeepEvent(n)};
+        end
+        D = events(D,1,evt);
+        D2 = clone(D, D.fnamedat, [D.nchannels D.nsamples D.ntrials]);
+        D2(:,:,:) = D(:,:,:);
+        save(D2);
+        fprintf('\n \n MESSAGE: .. %s parameters updated ..::\n',patientCode); 
+        %set_final_status('OK') 
+    end
+end
 val_flag = 0;
 val_flag = 0;
 if strcmpi(patientCode(5:end),'VAL')
 if strcmpi(patientCode(5:end),'VAL')
     load('/gin/data/database/02-raw/stim_parameters-ftract-val.mat','stim_params')
     load('/gin/data/database/02-raw/stim_parameters-ftract-val.mat','stim_params')
@@ -287,7 +334,7 @@ if strcmpi(patientCode(5:end),'VAL')
 end
 end
 
 
 
 
-if mil_flag == 0 || fre_flag == 0 || buc_flag == 0 || val_flag == 0
+if mil_flag == 0 || fre_flag == 0 || buc_flag == 0 || str_flag == 0 || val_flag == 0
     pVals = [];
     pVals = [];
     for k = 1:length(KeepEvent)
     for k = 1:length(KeepEvent)
         xsub3 = regexp(Notes{KeepEvent(k)},rxp3,'match');
         xsub3 = regexp(Notes{KeepEvent(k)},rxp3,'match');

@@ -279,15 +279,15 @@ if size(DD,1)>1
 end
 end
 
 
 for i1 = 1:size(DD,1)
 for i1 = 1:size(DD,1)
-%     try
-%         [D,TimeWindow,TimeWindowWidth] = ImaGIN_spm_eeg_tf_main(deblank(DD(i1,:)),S{i1});
-%     catch
-%         try
+    try
+        [D,TimeWindow,TimeWindowWidth] = ImaGIN_spm_eeg_tf_main(deblank(DD(i1,:)),S{i1});
+    catch
+        try
             [D,TimeWindow,TimeWindowWidth] = ImaGIN_spm_eeg_tf_main(deblank(DD(i1,:)),S);
             [D,TimeWindow,TimeWindowWidth] = ImaGIN_spm_eeg_tf_main(deblank(DD(i1,:)),S);
-%         catch
-%             [D,TimeWindow,TimeWindowWidth] = ImaGIN_spm_eeg_tf_main(deblank(DD(i1,:)));
-%         end
-%     end
+        catch
+            [D,TimeWindow,TimeWindowWidth] = ImaGIN_spm_eeg_tf_main(deblank(DD(i1,:)));
+        end
+    end
     if (i1==1) && ~exist('S', 'var')
     if (i1==1) && ~exist('S', 'var')
         S.Method=D.tf.Method;
         S.Method=D.tf.Method;
         switch D.tf.Method
         switch D.tf.Method