forked from MouseLand/Kilosort
-
Notifications
You must be signed in to change notification settings - Fork 1
/
Copy pathStandardConfig_MOVEME.m
50 lines (38 loc) · 2.19 KB
/
StandardConfig_MOVEME.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
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
ops.chanMap = 'D:\GitHub\KiloSort2\configFiles\neuropixPhase3A_kilosortChanMap.mat';
% ops.chanMap = 1:ops.Nchan; % treated as linear probe if no chanMap file
% sample rate
ops.fs = 30000;
% frequency for high pass filtering (150)
ops.fshigh = 150;
% minimum firing rate on a "good" channel (0 to skip)
ops.minfr_goodchannels = 0.1;
% threshold on projections (like in Kilosort1, can be different for last pass like [10 4])
ops.Th = [10 4];
% how important is the amplitude penalty (like in Kilosort1, 0 means not used, 10 is average, 50 is a lot)
ops.lam = 10;
% splitting a cluster at the end requires at least this much isolation for each sub-cluster (max = 1)
ops.AUCsplit = 0.9;
% minimum spike rate (Hz), if a cluster falls below this for too long it gets removed
ops.minFR = 1/50;
% number of samples to average over (annealed from first to second value)
ops.momentum = [20 400];
% spatial constant in um for computing residual variance of spike
ops.sigmaMask = 30;
% threshold crossings for pre-clustering (in PCA projection space)
ops.ThPre = 8;
%% danger, changing these settings can lead to fatal errors
% options for determining PCs
ops.spkTh = -6; % spike threshold in standard deviations (-6)
ops.reorder = 1; % whether to reorder batches for drift correction.
ops.nskip = 25; % how many batches to skip for determining spike PCs
ops.GPU = 1; % has to be 1, no CPU version yet, sorry
% ops.Nfilt = 1024; % max number of clusters
ops.nfilt_factor = 4; % max number of clusters per good channel (even temporary ones)
ops.ntbuff = 64; % samples of symmetrical buffer for whitening and spike detection
ops.NT = 64*1024+ ops.ntbuff; % must be multiple of 32 + ntbuff. This is the batch size (try decreasing if out of memory).
ops.whiteningRange = 32; % number of channels to use for whitening each channel
ops.nSkipCov = 25; % compute whitening matrix from every N-th batch
ops.scaleproc = 200; % int16 scaling of whitened data
ops.nPCs = 3; % how many PCs to project the spikes into
ops.useRAM = 0; % not yet available
%%