2013-02-04
Adrien Peyrache / Peter Petersen
Recording from the anterior thalamus and the hippocampus, left hemisphere.
Anterior Thalamus:
Electrode groups: 1-8
8 shank probe, Neuronexus Buz64 design. Shank #1 is the most lateral, #8 the most medial
Probe is perpendicular to the midline (coronal plane), tilted by 15 degrees (tips pointed toward midline), mounted on a movable drive.
Insertion coordinates: AP -0.7mm, ML 0.5mm
Approx depth from surface: 3.2mm
Hippocampus:
Electrode group: 9
6 tungsten wires (fifty micrometer in diamater), one (at least) is in the pyramidal layer of CA1 as shown by the presence of Sharp-Wave/Ripples.
Insertion coordinates: AP -2.2mm, ML 1.3mm, depth: 1.7mm
16
70
20000
20
1000
0
1250
63
56
62
57
61
58
60
59
55
48
54
49
53
50
52
51
16
25
17
24
18
22
20
21
26
19
28
23
29
27
30
31
13
4
9
2
5
3
1
0
7
14
6
15
8
12
11
10
46
41
47
40
44
43
45
42
38
33
39
32
36
35
37
34
64
65
66
67
68
69
63
56
62
57
61
58
60
59
32
16
3
55
48
54
49
53
50
52
51
32
16
3
16
25
17
24
18
22
20
21
32
16
3
26
19
28
23
29
27
30
31
32
16
3
13
4
9
2
5
3
1
0
32
16
3
7
14
6
15
8
12
11
10
32
16
3
46
41
47
40
44
43
45
42
32
16
3
38
33
39
32
36
35
37
34
32
16
3
64
65
66
67
68
69
32
16
3
1
2
2
2
2
3
3
2
3
3
3
4
3
5
3
6
------
------
------
------
------
3
7
--
--
--
--
--
4
2
--
--
--
--
--
5
2
5
3
5
4
--
--
--
--
5
5
------
------
------
------
6
2
6
3
6
4
6
5
6
6
6
7
6
8
----
----
----
----
----
6
9
7
2
----
----
----
----
----
7
3
------
------
------
------
------
8
2
8
3
--------------
--------------
--------------
--------------
--------------
8
4
--
--
--
--
--
9
2
9
3
9
4
9
5
9
6
9
7
--------
--------
--------
--------
--------
9
8
----------
----------
----------
----------
----------
0.2
32
16
0
#b341ff
#b341ff
#b341ff
0
0
1
#b341ff
#b341ff
#b341ff
1
0
2
#b341ff
#b341ff
#b341ff
2
0
3
#b341ff
#b341ff
#b341ff
3
0
4
#b341ff
#b341ff
#b341ff
4
0
5
#b341ff
#b341ff
#b341ff
5
0
6
#69ff64
#69ff64
#69ff64
6
0
7
#69ff64
#69ff64
#69ff64
7
0
8
#69ff64
#69ff64
#69ff64
8
0
9
#b341ff
#b341ff
#b341ff
9
0
10
#69ff64
#69ff64
#69ff64
10
0
11
#69ff64
#69ff64
#69ff64
11
0
12
#69ff64
#69ff64
#69ff64
12
0
13
#b341ff
#b341ff
#b341ff
13
0
14
#69ff64
#69ff64
#69ff64
14
0
15
#69ff64
#69ff64
#69ff64
15
0
16
#50b0ff
#50b0ff
#50b0ff
16
0
17
#50b0ff
#50b0ff
#50b0ff
17
0
18
#50b0ff
#50b0ff
#50b0ff
18
0
19
#ffe46d
#ffe46d
#ffe46d
19
0
20
#50b0ff
#50b0ff
#50b0ff
20
0
21
#50b0ff
#50b0ff
#50b0ff
21
0
22
#50b0ff
#50b0ff
#50b0ff
22
0
23
#ffe46d
#ffe46d
#ffe46d
23
0
24
#50b0ff
#50b0ff
#50b0ff
24
0
25
#50b0ff
#50b0ff
#50b0ff
25
0
26
#ffe46d
#ffe46d
#ffe46d
26
0
27
#ffe46d
#ffe46d
#ffe46d
27
0
28
#ffe46d
#ffe46d
#ffe46d
28
0
29
#ffe46d
#ffe46d
#ffe46d
29
0
30
#ffe46d
#ffe46d
#ffe46d
30
0
31
#ffe46d
#ffe46d
#ffe46d
31
0
32
#e4ffb4
#e4ffb4
#e4ffb4
32
0
33
#e4ffb4
#e4ffb4
#e4ffb4
33
0
34
#e4ffb4
#e4ffb4
#e4ffb4
34
0
35
#e4ffb4
#e4ffb4
#e4ffb4
35
0
36
#e4ffb4
#e4ffb4
#e4ffb4
36
0
37
#e4ffb4
#e4ffb4
#e4ffb4
37
0
38
#e4ffb4
#e4ffb4
#e4ffb4
38
0
39
#e4ffb4
#e4ffb4
#e4ffb4
39
0
40
#0d00ff
#0d00ff
#0d00ff
40
0
41
#0d00ff
#0d00ff
#0d00ff
41
0
42
#0d00ff
#0d00ff
#0d00ff
42
0
43
#0d00ff
#0d00ff
#0d00ff
43
0
44
#0d00ff
#0d00ff
#0d00ff
44
0
45
#0d00ff
#0d00ff
#0d00ff
45
0
46
#0d00ff
#0d00ff
#0d00ff
46
0
47
#0d00ff
#0d00ff
#0d00ff
47
0
48
#ff1df4
#ff1df4
#ff1df4
48
0
49
#ff1df4
#ff1df4
#ff1df4
49
0
50
#ff1df4
#ff1df4
#ff1df4
50
0
51
#ff1df4
#ff1df4
#ff1df4
51
0
52
#ff1df4
#ff1df4
#ff1df4
52
0
53
#ff1df4
#ff1df4
#ff1df4
53
0
54
#ff1df4
#ff1df4
#ff1df4
54
0
55
#ff1df4
#ff1df4
#ff1df4
55
0
56
#41ff67
#41ff67
#41ff67
56
0
57
#41ff67
#41ff67
#41ff67
57
0
58
#41ff67
#41ff67
#41ff67
58
0
59
#41ff67
#41ff67
#41ff67
59
0
60
#41ff67
#41ff67
#41ff67
60
0
61
#41ff67
#41ff67
#41ff67
61
0
62
#41ff67
#41ff67
#41ff67
62
0
63
#41ff67
#41ff67
#41ff67
63
0
64
#ff0000
#ff0000
#ff0000
64
0
65
#ff0000
#ff0000
#ff0000
65
0
66
#ff0000
#ff0000
#ff0000
66
0
67
#ff0000
#ff0000
#ff0000
67
0
68
#ff0000
#ff0000
#ff0000
68
0
69
#ff0000
#ff0000
#ff0000
69
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
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
Mouse17-130204
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.