clear all; clc;

CV='cv1';
Do.MRI = false;
Do.EPHYS = true;
Do.COMBINE = true;

%% MRI ====================================================================
if Do.MRI
    monkeys = {'M02','M01'};
    
    models = {...
        ['linhrf_' CV '_dhrf'],...
            ['linhrf_' CV '_mhrf'],...
        ['linhrf_' CV '_dhrf_neggain'],...
            ['linhrf_' CV '_mhrf_neggain'],...
        ['doghrf_' CV '_dhrf'],...
            ['doghrf_' CV '_mhrf'],...
        ['csshrf_' CV '_dhrf'],...
            ['csshrf_' CV '_mhrf']};
    
    output = ['AllFits_MRI_' CV];
    
    ck_GetMRI_pRF(monkeys,models,output);
end

%% EPHYS ==================================================================
if Do.EPHYS
    monkeys = {'M03','M04'};
    
     models = {...
         ['linear_ephys_' CV],...
         ['linear_ephys_' CV '_neggain'],...
         ['dog_ephys_' CV],...
         ['css_ephys_' CV],...
         'classicRF'};
    output = ['AllFits_ephys_' CV];
    
    ck_GetEphys_pRF(monkeys,models,output);
end

%% COMBINE ================================================================
if Do.COMBINE
    homefld = pwd;
    cd ../../../
    SHARED_ROOT_FLD = pwd;
    cd(homefld)

    fitres_path = fullfile(SHARED_ROOT_FLD,'FitResults');  
    
    % load mri
    fprintf('Loading MRI\n');
    R_MRI = load(fullfile(fitres_path,'MRI','Combined',['AllFits_MRI_' CV]));
    D99 = R_MRI.D99;
    R_MRI = R_MRI.R;
   
    % load ephys
    fprintf('Loading EPHYS\n');
    R_EPHYS = load(fullfile(fitres_path,'EPHYS','Combined',['AllFits_ephys_' CV]));
    R_EPHYS = R_EPHYS.R;
    
    % save
    fprintf('Saving results from both modalities\n');
    save(fullfile(fitres_path,'MULTIMODAL',['AllFits_' CV]),...
        'R_MRI','R_EPHYS','D99');
end

%% Done 
fprintf('ALL DONE!\n');