convert2nwbEcephysMatNnx/nwbSessionParams.m

% Session-specific NWB parameters
% Parameters that are different for different animals are:
%   sessionID
%   sessionDescription
%   sessionNotes
%   areas
%   endCh
%   probeInserted
%   electrodeName
%   nChannelsPerShank
%   electrodeCoordinates
%   electrodeImplantationType

% General info
sessionID = {'20190122191103'};
sessionDescription = {'awake restrained'};
sessionNotes = {'Probes advanced to 700 um (z-axis) into the left S1bf and 3250 um (z-axis) into the left VB.   Recording was done in an awake head-restraint animal.'};
endCh{1}{1} = [16]; % Corresponding probe end channels starting from the tip of the probe. Corresponding and previous end channels are used to work out probe channels that reside in the corresponding brain area.
endCh{1}{2} = [6 10 18 24 32];

for iSess = 1:numel(sessionID)
  sessionStartTime{iSess} = datetime(str2double(sessionID{iSess}(1:4)), str2double(sessionID{iSess}(5:6)), str2double(sessionID{iSess}(7:8)),...
    str2double(sessionID{iSess}(9:10)), str2double(sessionID{iSess}(11:12)), str2double(sessionID{iSess}(13:14))); %#ok<*SAGROW>
  
  % Probe #1 info
  ref = 1; % probe reference
  probeInserted{iSess}{ref} = true; % If the probe used at all
  if probeInserted{iSess}{ref} && ~isempty(endCh{iSess}{ref})
    electrodeName{iSess}{ref} = 'CM16LP-A1x16-Poly2-5mm-50s-177';
    electrodeDescription{iSess}{ref} = ['Single shank low density probe in position ' num2str(ref)];
    electrodeManufacturer{iSess}{ref} = 'Neuronexus';
    electrodeFolder{iSess}{ref} = [animalRawDataFolder filesep sessionID{iSess} '1'];
    if ~exist(electrodeFolder{iSess}{ref}, 'file') % in case a probe is missing
      electrodeFolder{iSess}{ref} = electrodeFolder{iSess}{ref}(1:end-1);
    end
    electrodeMap = [electrodeFolder{iSess}{ref} filesep 'forPRB_' electrodeName{iSess}{ref} '.mat'];
    nShanks{iSess}{ref} = 1;
    nChannelsPerShank{iSess}{ref} = 16;
    nCh{iSess}{ref} = nChannelsPerShank{iSess}{ref}*nShanks{iSess}{ref}; % total number of probe channels
    areas = {'S1'}; % brain areas that this probe spans
    electrodeLocation{iSess}{ref} = electrodeLocations(areas, endCh{iSess}{ref}, nCh{iSess}{ref}); % Brain area assigned to each recording channel.
    electrodeCoordinates{iSess}{ref} = [-1.8, -2.5, 0]; % Electrode insertion location on the cortical surface in Paxinos coords: AP (posterior negative), ML (left negative), DV (recording site position starting with the tip of the probe.
    electrodeCoordinates{iSess}{ref} = repmat(electrodeCoordinates{iSess}{ref}, nCh{iSess}{ref}, 1); % Coordinates of each probe recording channel (the probe rotation angle is not taken into account). Y coordinates are relative to the tip of the probe.
    load(electrodeMap, 'ycoords');
    if size(ycoords,1) == 1
      electrodeCoordinates{iSess}{ref}(:,3) = ycoords(1:nCh{iSess}{ref})'./1000;
    else
      electrodeCoordinates{iSess}{ref}(:,3) = ycoords(1:nCh{iSess}{ref})./1000;
    end
    electrodeLabel{iSess}{ref} = ['probe' num2str(ref)];
    electrodeImplantationType{iSess}{ref} = 'chronic';
  else % The case when the probe #1 is missing
    electrodeName{iSess}{ref} = []; %#ok<*UNRCH>
    electrodeDescription{iSess}{ref} = [];
    electrodeManufacturer{iSess}{ref} = [];
    electrodeFolder{iSess}{ref} = [];
    nShanks{iSess}{ref} = [];
    nChannelsPerShank{iSess}{ref} = [];
    nCh{iSess}{ref} = [];
    electrodeLocation{iSess}{ref} = [];
    electrodeCoordinates{iSess}{ref} = [];
    electrodeLabel{iSess}{ref} = [];
    electrodeImplantationType{iSess}{ref} = [];
  end
  
  % Probe #2 info
  ref = 2;
  probeInserted{iSess}{ref} = true;
  if probeInserted{iSess}{ref} && ~isempty(endCh{iSess}{ref})
    electrodeName{iSess}{ref} = 'CM32-A1x32-Edge-5mm-100-177';
    electrodeDescription{iSess}{ref} = ['Single shank low density probe in position ' num2str(ref)];
    electrodeManufacturer{iSess}{ref} = 'Neuronexus';
    electrodeFolder{iSess}{ref} = [animalRawDataFolder filesep sessionID{iSess} '26'];
    if ~exist(electrodeFolder{iSess}{ref},'file') || numel(dir(electrodeFolder{iSess}{ref})) <= 2 % in case the raw data folder name is shorter
      electrodeFolder{iSess}{ref} = electrodeFolder{iSess}{ref}(1:end-1);
    end
    if ~exist(electrodeFolder{iSess}{ref}, 'file') % in case a probe is missing
      electrodeFolder{iSess}{ref} = electrodeFolder{iSess}{ref}(1:end-1);
    end
    electrodeMap = [electrodeFolder{iSess}{ref} filesep 'forPRB_' electrodeName{iSess}{ref} '.mat'];
    nShanks{iSess}{ref} = 1;
    nChannelsPerShank{iSess}{ref} = 32;
    nCh{iSess}{ref} = nChannelsPerShank{iSess}{ref}*nShanks{iSess}{ref};
    areas = {'Po','LP','DG','CA1','RSC'};
    electrodeLocation{iSess}{ref} = electrodeLocations(areas, endCh{iSess}{ref}, nCh{iSess}{ref});
    electrodeCoordinates{iSess}{ref} = [-1.8, -0.5, 0];
    electrodeCoordinates{iSess}{ref} = repmat(electrodeCoordinates{iSess}{ref}, nCh{iSess}{ref}, 1);
    load(electrodeMap, 'ycoords');
    if size(ycoords,1) == 1
      electrodeCoordinates{iSess}{ref}(:,3) = ycoords(1:nCh{iSess}{ref})'./1000;
    else
      electrodeCoordinates{iSess}{ref}(:,3) = ycoords(1:nCh{iSess}{ref})./1000;
    end
    electrodeLabel{iSess}{ref} = ['probe' num2str(ref)];
    electrodeImplantationType{iSess}{ref} = 'chronic';
  else % The case when the probe #2 is missing
    electrodeName{iSess}{ref} = []; %#ok<*UNRCH>
    electrodeDescription{iSess}{ref} = [];
    electrodeManufacturer{iSess}{ref} = [];
    electrodeFolder{iSess}{ref} = [];
    nShanks{iSess}{ref} = [];
    nChannelsPerShank{iSess}{ref} = [];
    nCh{iSess}{ref} = [];
    electrodeLocation{iSess}{ref} = [];
    electrodeCoordinates{iSess}{ref} = [];
    electrodeLabel{iSess}{ref} = [];
    electrodeImplantationType{iSess}{ref} = [];
  end
end



%% Local functions
function electrodeLocationVec = electrodeLocations(areas, endCh, nCh)
% electrodeLocationVec = electrodeLocations(areas, endCh, nCh)
% Function assigns brain area to a channel

for iCh = 1:nCh
  areaInd = find(endCh >= iCh,1);
  electrodeLocationVec{iCh} = areas{areaInd};
end
end