0001 function make_trial_meg(p)
0002
0003
0004
0005
0006
0007 disp(mfilename);
0008
0009
0010 parm.Pretrigger_ms = p.Pretrigger_ms;
0011 parm.Posttrigger_ms = p.Posttrigger_ms;
0012 parm.condition={'stim'};
0013 parm.status_ch={'432'};
0014 parm.trig_type='analog';
0015 parm.slope='low_to_high';
0016 parm.status_level=0.5;
0017
0018 plot_parm.mode = 1;
0019 plot_parm.NXmax = 25;
0020 plot_parm.NYmax = 5;
0021
0022 Nfile = length(p.file_list);
0023
0024
0025 for file = 1:Nfile
0026
0027 parm.data_file = fullfile(p.meg_dirname, p.filter_dirname, [p.file_list{file} '.meg.mat']);
0028
0029
0030 parm.trig_file = fullfile(p.meg_dirname, p.filter_dirname, [p.file_list{file} '.trig.mat']);
0031
0032
0033 vb_job_trial_onset(p.proj_root,parm);
0034 vb_plot_status(fullfile(p.proj_root, parm.trig_file), plot_parm);
0035 end
0036
0037
0038
0039 for run = 1:2
0040
0041 trig_file = fullfile(p.proj_root, p.meg_dirname, p.filter_dirname, ['v' num2str(run) '.trig.mat']);
0042 load(trig_file, 'cond_id', 'ix_trial', 'trig')
0043
0044
0045 cond_id = cond_id(2:end);
0046 ix_trial = ix_trial(:,2:end);
0047 trig = trig(2:end);
0048
0049
0050 vb_save(trig_file, 'cond_id', 'ix_trial', 'trig');
0051 end
0052
0053
0054 data_file = fullfile(p.proj_root, p.meg_dirname, p.filter_dirname, 'a1.meg.mat');
0055 meg_info = vb_load_measurement_info(data_file);
0056 samplingrate = meg_info.SampleFreq;
0057 for file = 1:Nfile
0058
0059 trig_file = fullfile(p.proj_root, p.meg_dirname, p.filter_dirname, [p.file_list{file} '.trig.mat']);
0060 load(trig_file, 'trig')
0061 iti_meg = diff(trig);
0062
0063
0064 trig_file = fullfile(p.proj_root, p.eeg_dirname, p.filter_dirname, [p.file_list{file} '.trig.mat']);
0065 load(trig_file,'trig')
0066 iti_eeg = diff(trig);
0067
0068
0069 dif = abs(iti_meg-iti_eeg);
0070 disp(['Max difference of ITI is ' num2str(max(dif)/samplingrate*1000) ' ms for ' p.file_list{file}])
0071 end
0072
0073
0074 for file = 1:Nfile
0075
0076 filtered_file = fullfile(p.proj_root, p.meg_dirname, p.filter_dirname, [p.file_list{file} '.meg.mat']);
0077 proc_spec.trig_file = fullfile(p.proj_root, p.meg_dirname, p.filter_dirname, [p.file_list{file} '.trig.mat']);
0078
0079
0080 trial_file = fullfile(p.proj_root, p.meg_dirname, p.trial_dirname, [p.file_list{file} '.meg.mat']);
0081
0082
0083 vb_msrmnt_make_trial_data(filtered_file, proc_spec, trial_file);
0084 end
0085
0086
0087
0088