2014-02-04
Adrien Peyrache
Recording from the anterior thalamus, and post-subiculum, left hemisphere.
Post-subiculum:
Electrode Groups: 1-4
4 shank probe, Neuronexus Buz32 design. Shank #1 is the most medial, #4 the most lateral
Probe is perpendicular to the midline (coronal plane), tilted by 10 degrees (tips pointed away from midline), mounted on a movable drive.
Insertion coordinates: AP -4.25mm, ML 1.2-1.9mm
Approx depth from surface: 2.26mm
Anterior Thalamus:
Electrode groups: 5-8.
4 shank probe, Neuronexus Buz32 design. Shank #1 (5th electrode group) is the most lateral, #4 (8th electrode group) the most medial
Probe is perpendicular to the midline (coronal plane), mounted on a movable drive.
Insertion coordinates: AP -0.7mm, ML 0.5-1.2mm
Approx depth from surface: 2.72mm
16
64
20000
20
1000
0
1250
16
19
18
17
20
28
24
26
21
22
25
30
27
23
29
31
8
11
0
7
9
5
1
3
13
14
15
12
2
10
4
6
48
51
50
49
52
60
56
58
53
54
57
62
59
55
61
63
40
43
32
39
41
37
33
35
45
46
47
44
34
42
36
38
16
19
18
17
20
28
24
26
32
16
3
21
22
25
30
27
23
29
31
32
16
3
8
11
0
7
9
5
1
3
32
16
3
13
14
15
12
2
10
4
6
32
16
3
48
51
50
49
52
60
56
58
32
16
3
53
54
57
62
59
55
61
63
32
16
3
40
43
32
39
41
37
33
35
32
16
3
45
46
47
44
34
42
36
38
32
16
3
2
2
2
3
2
4
2
5
--
--
--
--
--
4
2
4
3
4
4
4
5
4
6
4
7
4
8
4
9
4
10
4
11
4
12
4
13
5
2
5
3
5
4
5
5
5
6
------------
------------
------------
------------
------------
5
7
------------
------------
------------
------------
------------
5
8
----
----
----
----
----
5
9
--
--
--
--
--
5
10
--
--
--
--
--
5
11
--
--
--
--
--
0.2
32
16
0
#2dfffc
#2dfffc
#2dfffc
0
0
1
#2dfffc
#2dfffc
#2dfffc
1
0
2
#ebff34
#ebff34
#ebff34
2
0
3
#2dfffc
#2dfffc
#2dfffc
3
0
4
#ebff34
#ebff34
#ebff34
4
0
5
#2dfffc
#2dfffc
#2dfffc
5
0
6
#ebff34
#ebff34
#ebff34
6
0
7
#2dfffc
#2dfffc
#2dfffc
7
0
8
#2dfffc
#2dfffc
#2dfffc
8
0
9
#2dfffc
#2dfffc
#2dfffc
9
0
10
#ebff34
#ebff34
#ebff34
10
0
11
#2dfffc
#2dfffc
#2dfffc
11
0
12
#ebff34
#ebff34
#ebff34
12
0
13
#ebff34
#ebff34
#ebff34
13
0
14
#ebff34
#ebff34
#ebff34
14
0
15
#ebff34
#ebff34
#ebff34
15
0
16
#0080ff
#0080ff
#0080ff
16
0
17
#0080ff
#0080ff
#0080ff
17
0
18
#0080ff
#0080ff
#0080ff
18
0
19
#0080ff
#0080ff
#0080ff
19
0
20
#0080ff
#0080ff
#0080ff
20
0
21
#ff265c
#ff265c
#ff265c
21
0
22
#ff265c
#ff265c
#ff265c
22
0
23
#ff265c
#ff265c
#ff265c
23
0
24
#0080ff
#0080ff
#0080ff
24
0
25
#ff265c
#ff265c
#ff265c
25
0
26
#0080ff
#0080ff
#0080ff
26
0
27
#ff265c
#ff265c
#ff265c
27
0
28
#0080ff
#0080ff
#0080ff
28
0
29
#ff265c
#ff265c
#ff265c
29
0
30
#ff265c
#ff265c
#ff265c
30
0
31
#ff265c
#ff265c
#ff265c
31
0
32
#55d2ff
#55d2ff
#55d2ff
32
0
33
#55d2ff
#55d2ff
#55d2ff
33
0
34
#ffb06f
#ffb06f
#ffb06f
34
0
35
#55d2ff
#55d2ff
#55d2ff
35
0
36
#ffb06f
#ffb06f
#ffb06f
36
0
37
#55d2ff
#55d2ff
#55d2ff
37
0
38
#ffb06f
#ffb06f
#ffb06f
38
0
39
#55d2ff
#55d2ff
#55d2ff
39
0
40
#55d2ff
#55d2ff
#55d2ff
40
0
41
#55d2ff
#55d2ff
#55d2ff
41
0
42
#ffb06f
#ffb06f
#ffb06f
42
0
43
#55d2ff
#55d2ff
#55d2ff
43
0
44
#ffb06f
#ffb06f
#ffb06f
44
0
45
#ffb06f
#ffb06f
#ffb06f
45
0
46
#ffb06f
#ffb06f
#ffb06f
46
0
47
#ffb06f
#ffb06f
#ffb06f
47
0
48
#14ff37
#14ff37
#14ff37
48
0
49
#14ff37
#14ff37
#14ff37
49
0
50
#14ff37
#14ff37
#14ff37
50
0
51
#14ff37
#14ff37
#14ff37
51
0
52
#14ff37
#14ff37
#14ff37
52
0
53
#c56fff
#c56fff
#c56fff
53
0
54
#c56fff
#c56fff
#c56fff
54
0
55
#c56fff
#c56fff
#c56fff
55
0
56
#14ff37
#14ff37
#14ff37
56
0
57
#c56fff
#c56fff
#c56fff
57
0
58
#14ff37
#14ff37
#14ff37
58
0
59
#c56fff
#c56fff
#c56fff
59
0
60
#14ff37
#14ff37
#14ff37
60
0
61
#c56fff
#c56fff
#c56fff
61
0
62
#c56fff
#c56fff
#c56fff
62
0
63
#c56fff
#c56fff
#c56fff
63
0
process_calcThreshold
executeScriptBool
1
Mandatory
rmsWindowLength
2
Mandatory
useage:
process_calcThreshold fileBaseName electrodeNumber
This script calculates the spike detection threshold and puts it in the file fileBaseName.threshold.electrodeNumber
parameters:
- executeScriptBool determines whether the script should be run with the next batch process.
- rmsWindowLength represents the root mean square (rms) integration window for detection.
process_detectSpikes
executeScriptBool
1
Mandatory
thresholdFactor
7
Mandatory
refractoryPeriod
8
Mandatory
rmsWindowLength
2
Mandatory
peakWindowLength
10
Mandatory
windowPeakAlign
5
Mandatory
spikesFromDatBool
1
Mandatory
This script does the following:
* detects the spikes (creates the .res.# file),
* aligns the peaks,
* gets the waveforms (creates the .spk.# file),
The variables used by the different programs called in this script depend on the sampling rate. The default is for a sampling rate of 10khz.
All the script parameters are multiplifaction factor of 10khz except the thresholdFactor.
parameters:
- executeScriptBool determines whether the script should be run with the next batch process.
- refractoryPeriod represents the after detection, how many samples to skip before detecting a new spike if one is detected.
- rmsWindowLength represents the root mean square (rms) integration window for detection.
- peakWindowLength represents the window (number of samples to test) to search for spike peak, seams to be used for eigen calculation or finding of the true spike peak.
- windowPeakAlign represents the middle of the peakWindowLength, where to take the spike peak in peakWindowLength.
- thresholdFactor is a factor to compute the detection threshold:
threshold = thresholdFactor * baseline
process_findCommonThreshold
executeScriptBool
1
Mandatory
maxThreshold
Mandatory
useage:
process_findCommonThreshold fileBaseName [fileBaseName ...]
for each electrodeGroup in fileBaseName.xml
the value in fileBaseName.threshold.electrodeGroup is replaced by the minimum threshold value across passed fileBaseNames.
parameters:
- executeScriptBool determines whether the script should be run with the next batch process.
- maxThreshold default = 500
process_merge
executeScriptBool
1
Mandatory
mergeBaseName
Mouse25-140204
Mandatory
process_merge
usage:
process_merge -n $mergeBaseName $fileBaseNames
merges .spk .res etc for multi-file clustering.
parameters:
- executeScriptBool determines whether the script should be run with the next batch process.
- mergeBaseName designates the fileBaseName of the merged files.
*** Requirements ***
* calls process_pca_light : to do pca without presence of .fil file
* calls pproc
* each fileBaseName.* must be in a directory fileBaseName - which can be performed by process_dirFiles
process_mhipass
executeScriptBool
1
Mandatory
frequency
800
Mandatory
libraryPath
/u16/local/matlab6.5/extern/lib/glnx86:/u16/local/matlab6.5/bin/glnx86:/u16/local/matlab6.5/sys/os/glnx86
Mandatory
This script does the high-pass filtering for subsequent spike extraction. The .fil file is created by the program called in this script.
The frequency parameter is the lowest frequency cut off.
The highest frequency is computed based on the sampling rate as 0.95 % of the Nyquist frequency:
sampling rate * 0.5 * 0.95
The program called in this script is currently a matlab program, therefore this part of the process relies on Matlab until this program is rewritten.
parameters:
- executeScriptBool determines whether the script should be run with the next batch process.
- libraryPath is the path for Matlab.
process_mlopass
executeScriptBool
1
Mandatory
libraryPath
/u16/local/matlab6.5/extern/lib/glnx86:/u16/local/matlab6.5/bin/glnx86:/u16/local/matlab6.5/sys/os/glnx86
Mandatory
This script does the low-pass filtering (for local field potentials). The .eeg file is created by the program called in this script.
It calls a resampling program which first apply an anti aliasing low pass filter to the data and then samples down to the local field potential sampling rate provided in the parameter file.
The resampling program is currently a matlab program, therefore this part of the process relies on Matlab until this program is rewritten.
parameters:
- executeScriptBool determines whether the script should be run with the next batch process.
- libraryPath is the path for Matlab.
process_multi_start
cpFirstXmlBool
1
Mandatory
rmFilFileBool
1
Mandatory
This script starts eeg and spike data processing on a list of input fileBaseNames
This script is called by:
process_multi_start fileBaseName [fileBaseNames]
Parameters:
- rmFilFileBool
1 = remove the high pass filtered .fil file to save space
- cpFirstXmlBool
1 = copy the xml file from the first fileBaseName to all other fileBaseNames to save time
* Note that ideally the common threshold would be for single channels rather than single electrodes (may be a problem of process_baseline)
process_mvFiles2Dirs
executeScriptBool
1
Mandatory
This script organizes files in preparation for process_merge by moving fileBaseName.* into a directory fileBaseName
useage:
process_mvFiles2Dirs fileBaseName
parameters:
- executeScriptBool determines whether the script should be run with the next batch process.
process_pca
executeScriptBool
1
Mandatory
beforePeak
4
Mandatory
afterPeak
4
Mandatory
nSamplesInPCA
16
Mandatory
This script performs the principle component analysis (PCA ). It computes the variance (creates the .m1m2.# file) and creates the .fet.# and .mm.# files.
The variables used to compute the PCA depend on the sampling rate. The default is for a sampling rate of 10khz.
The parameters in the parameter file are multiplifaction factor of 10khz.
parameters:
- executeScriptBool determines whether the script should be run with the next batch process.
- beforePeak represents the number of samples before peak to be used for the PCA. There is one value for all the electrodes.
- afterPeak represents the number of samples after peak to be used for the PCA. There is one value for all the electrodes.
- nSamplesInPCA represents the number of samples used for PCA. There is one value for all the electrodes.
updateParFile
executeScriptBool
1
Mandatory
samplesInWaveform
16
Mandatory
peakSample
8
Mandatory
This script updates the parameter file.
It adds for each spike group, the number of samples in each waveform, the sample index of the peak in a waveform and the number of feaures to use to do the PCA.
The variables used by the different programs depend on the sampling rate. The default is for a sampling rate of 10khz.
All the script parameters are multiplifaction factor of 10khz.
parameters:
- executeScriptBool determines whether the script should be run with the next batch process.
- The parameter samplesInWaveform represents the number of samples in each waveform .
- The parameter peakSample represents the sample index of the peak in a waveform.