ChronicGui/S150406_EventBatch.m at master · jamesjun/ChronicGui · GitHub

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
51
52
53
54
55
56
57
58
59
60
61
62
csAnimals = {'ANM282996', 'ANM286577', 'ANM286578', 'ANM292757', ...
    'ANM279097', 'ANM279094', 'ANM287075', 'ANM287074'};
fSave = 1;
tDuration = 900;
fParfor = 1;
nAnimals = numel(csAnimals);
viDay = [];
viShank = [];
fPlot = 0;
freqLim = [500 6000];
fMeanSubt = 2;
thresh = [];
fDiffPair = 0;
fMeasureEvent = 1;
vcFet = 'vpp';
vcPeak = 'vpp';

fCluster = 1;
funcFet = [];
SIGMA_FACTOR = 6;

fDebug = 1;

%%
% %60% spent on load time
if fDebug
    nAnimals = 1; tDuration = 200; fSave = 0; fParfor = 0;
    viDay = 20; viShank = 1; fPlot = 1;%debug
    fDiffPair = 0; fMeasureEvent = 1; fCluster = 1; fMeanSubt = 6; thresh = [];
end

% for fDiffPair = 0:1
for iAnimal = 1:nAnimals
    tic
    try
        vcAnimal = csAnimals{iAnimal};
        obj = Animal(vcAnimal, 'readDuration', [0 tDuration], ...
            'fParfor', fParfor, 'fOverwrite', 0);
        obj = obj.getFet('fCluster', fCluster, 'fPlot', fPlot, ...
            'fDiffPair', fDiffPair, 'readDuration', [0 tDuration], ...
            'fParfor', fParfor, ...
            'vcFet', vcFet, 'funcFet', funcFet, ...
            'viDay', viDay, 'viShank', viShank, 'maxAmp', 500, ...
            'fKeepNoiseClu', 0, 'freqLim', freqLim, ...
            'SIGMA_FACTOR', SIGMA_FACTOR, ...
            'fMeasureEvent', fMeasureEvent, 'vcPeak', vcPeak, ...
            'fMeanSubt', fMeanSubt, 'thresh', thresh, 'vcDist', 'euclidean');
        obj.plotBarClu();
        obj.plotEvents();
        if fSave
            eval(sprintf('%s = obj;', vcAnimal));
            clear obj;
            eval(sprintf('save D150407clu_%s_%d %s -v7.3;', ...
                vcAnimal, tDuration, vcAnimal));
            clear(vcAnimal);
        end
    catch
        fprintf('%s, %s\n', lasterr, vcAnimal);
    end
    toc
end
% end