%% 2022-04-19
% Load in healthy hemo recipient data

function [supName,outDir,expName, expFile, expDir, irange, irangeby2, cutoff] = loadHealthyHemoTrxData
    supName = "Healthy recipients";
    expFolder = uigetdir('Select directory with healthy M recipient transfusion data');
    outDir = uigetdir('Select output directory');

    expName{1} = "Hemo_CM_A";
    expFile{1} = '02-28-2020-08-29-13__Circadian_CsWF-Hemo6_M_Injected_1-128_Day1_margoConvert_denoise.mat';
    expDir{1} = strcat(expFolder,'\02-28-2020-08-29-13__Circadian_CsWF-Hemo6_M_Injected_1-128_Day1\');
    irange{1} = 1:64;
    irangeby2{1} = 1:2:64; 
    cutoff{1} = exp(-1*(0.3+0.18)/2); %% based on speed histogram found above

    % Hemo 5 data
    expName{2} = "Hemo_CM_B";
    expFile{2} = '01-30-2020-08-16-51__Circadian_CsWF-Recipients_M_Injected_1-128_Day1_margoConvert_denoise.mat';
    expDir{2} = strcat(expFolder,'01-30-2020-08-16-51__Circadian_CsWF-Recipients_M_Injected_1-128_Day1\');
    irange{2} = 1:64;
    irangeby2{2} = 1:2:64; 
    cutoff{2} = exp(-1*(0.299+0.182)/2); %% based on speed histogram found above

    %Hemo 4 data
    expName{3} = "Hemo_CM_C";
    expFile{3} = '12-12-2019-08-42-57__Circadian_CsWF-Hemo4_1-96_Day1_margoConvert_denoise.mat';
    expDir{3} = strcat(expFolder,'12-12-2019-08-42-57__Circadian_CsWF-Hemo4_1-96_Day1\');
    irange{3} = [1:26,33:64];
    irangeby2{3} = [1:2:26,33:2:64]; 
    cutoff{3} = exp(-1*(0.376+0.264)/2);     
end