%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%PREPROCESSING Taco ERF
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
cfg = [];
cfg.channel = {'MEG' 'STIM'};
cfg.dataset = 'TacoBakker_MIVI_20050201_01.ds';
cfg.artfctdef.eog.sgn ={'hEOG' 'vEOG'};
cfg.fsample=1200;
cfg.datatype='continuous';

% Trial definition
cfg.trialfun='trials_extract';
cfg.trialdef.eventtype  = '?';
cfg.trialdef.stim = [21:25 31:35 41:45 51:55];
cfg.trialdef.resp = [11 12];
cfg.trialdef.prestim = 0.2; % trial starts 0.2 sec before stimulus onset
cfg.trialdef.poststim = 0.6; % trial ends 0.6 sec after stim onset
%cfg.trialdef.preresp = 0.2; % trial ends 0.2 sec before response onset
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Artifact options
cfg.rejecteog         = 'yes'
cfg.artfctdef.eog.feedback = 'yes'
cfg.rejectmuscle      = 'yes'
cfg.artfctdef.muscle.feedback = 'yes'
cfg.rejectjump        = 'yes'
cfg.artfctdef.muscle.feedback = 'yes'


cfg.blc = 'yes';                % baseline correction
cfg.blcwindow = [-0.1 0];       % baseline correction window

cfg.lpfilter = 'yes'            % low-pass filter
cfg.lpfreq = 30;                % low-pass filter of 30 Hertz

cfg.padding = 0;                % trial is long enough


Taco_EF_stimlocked = preprocessing(cfg);  % puts the preprocessed data in data-structure

save Taco_EF_stimlocked Taco_EF_stimlocked
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%