% Extract F and OD experimental data
evol=readExperiment3;

%% FOD and OD 
fname = 'Exp1_YFP_data.mat';
% MEFL/cell factor
MEFL_factor = 0; %273.7 GFP/ 121.42 sfGFP
% Efficiency
ef = 0.48390*0.7657; 
% Excitation & Emission efficiency
    % GFP: 0.6606 & 0.5088
    % sfGFP: 0.8967 & 0.6333;
    % YFP: 0.48390*0.7657

% Extract FOD and OD matrices
plotExperiment_yb;

% Data reshape (matrices name changes depending on the experiment). In this
% case, we're reshaping YFP matrices
FOD = [estimation.FOD_pITU71_0  estimation.FOD_pITU72_0  estimation.FOD_pITU73_0...
    estimation.FOD_pITU81_0  estimation.FOD_pITU82_0 estimation.FOD_pITU83_0...
    estimation.FOD_pITU201_0]/ef;   
MOL = FOD*MEFL_factor;
OD = [estimation.OD_pITU71_0  estimation.OD_pITU72_0  estimation.OD_pITU73_0...
    estimation.OD_pITU81_0  estimation.OD_pITU82_0 estimation.OD_pITU83_0...
    estimation.OD_pITU201_0]/ef; 
time = estimation.Time_min;

% Save data
if MEFL_factor == 0
    save(fname,'estimation','FOD','OD','time');
else 
    save(fname,'estimation','FOD','OD','MOL','time');
end