sc_auto.m

function autoCall = sc_auto( firingRate, cData ) 
    vpp = cData.vpp;
    SNR = cData.SNR;
    ISIRatio = cData.ISIRatio;
    ISIViolations = cData.ISIViolations;
    IsoDist = cData.IsoDist;
    firing_std = cData.firingStd;
    
    SNR_high_thresh = 6;
    SNR_low_thresh = 2.5; %for low > SNR > high, check the isolation distance
    %SNR_multi_thresh = 2; %Too little signal to be sure of separation from background, useful if uniform firing
    ISIRatio_single_thresh = 0.005; %(max 0.5%)
    ISIRatio_ok_thresh = 0.015; 
    ISIViolations_thresh = 5;
    iso_thresh = 100;
    firingRate_thresh = 0.05;      %in Hz
    
    firing_std_thresh = 0.5;
    
    autoCall = "None";
    if firingRate < firingRate_thresh
        return;
    end
     
    if (SNR >= SNR_high_thresh) 
        % call a single if ISI conditions are met
        % check the ISI
        if ( (ISIRatio <= ISIRatio_single_thresh) || (ISIViolations <= ISIViolations_thresh) )
           autoCall = "single";
           return;
        elseif (ISIRatio <= ISIRatio_ok_thresh)
           autoCall = "ok";
        else
               if firing_std < 0.5*firing_std_thresh
                    autoCall = "multi";
               elseif firing_std < firing_std_thresh
                    autoCall = "multi/None";
               end      
        end
    
    elseif (SNR >= SNR_low_thresh)
        % if isolated, can be a single or ok
        if (isnan(IsoDist) || IsoDist > iso_thresh )
            if ( (ISIRatio <= ISIRatio_single_thresh) || (ISIViolations <= ISIViolations_thresh) )
               autoCall = "single/ok";  % whether unit is good enough to be single a function of human judged isolation
               return;
            elseif (ISIRatio <= ISIRatio_ok_thresh)
               autoCall = "ok"; % Ok
            else
               if firing_std < 0.5*firing_std_thresh
                    autoCall = "multi";
               elseif firing_std < firing_std_thresh
                    autoCall = "multi/None";
               end
            end
        else
            % if not isolated, call as a multi/ok if ISI is low enough
            % otherwise call as multi if unniform enough
            if ( (ISIRatio <= ISIRatio_ok_thresh) || (ISIViolations <= ISIViolations_thresh))
                autoCall = "ok";
            else
                if firing_std < 0.5*firing_std_thresh
                    autoCall = "multi";
                elseif firing_std < firing_std_thresh
                    autoCall = "multi/None";
                end
            end
        end
%     elseif (SNR >= SNR_multi_thresh)
%         % just check for stable firing. If "quite stable" very likely
%         % qualifies as a analyzable MUA; otherwise up to curator
%             if firing_std < 0.5*firing_std_thresh
%                 autoCall = "multi";
%             elseif firing_std < firing_std_thresh
%                 autoCall = "multi/None";
%             end
    end

    
end