source: ETALON/SPESO/ana2015/getDataForAnalysis.m @ 538

function dataset=getDataForAnalysis(listDays)

global basedir;
global plots_basedir;
global imageSizes;
global lsColorOption;

define_dirs_analysis;


dateAnalysed=9.9999e9;
eventsTiming=[];
eventsChannelsOn=[];
eventsSignalValue=[];
eventsPeakValue=[];
eventsRMSValue=[];
beamInfoDataTimed=[];

%%% Merges the beam info data
for iloop=1:length(listDays)
    dateStr=listDays{iloop};
    mergedInfoFile=[ basedir '/'  dateStr '/' 'merged_info_'  dateStr '.mat'];
    mif=load(mergedInfoFile);
    
    dateStr
    length(mif.eventsTiming)
    
    if (mif.dateAnalysed<dateAnalysed)
        dateAnalysed=mif.dateAnalysed;
    end
    if (length(eventsChannelsOn)>0)
        eventsTiming=[ eventsTiming mif.eventsTiming ];
        eventsChannelsOn=[ eventsChannelsOn'  mif.eventsChannelsOn' ]';
        eventsSignalValue=[ eventsSignalValue'  mif.eventsSignalValue' ]';
        eventsPeakValue=[ eventsPeakValue' mif.eventsPeakValue' ]';
        eventsRMSValue=[ eventsRMSValue' mif.eventsRMSValue' ]';
        beamInfoDataTimed=[ beamInfoDataTimed' mif.beamInfoDataTimed' ]';        
    else
        eventsTiming=mif.eventsTiming;
        eventsChannelsOn=mif.eventsChannelsOn;
        eventsSignalValue=mif.eventsSignalValue;
        eventsPeakValue=mif.eventsPeakValue;
        eventsRMSValue=mif.eventsRMSValue;
        beamInfoDataTimed=mif.beamInfoDataTimed;
    end %this is the first time so we do not merge files
end % for each file


%%% To get some grating pos info
%    dataset.hourFraction=floor(rem(eventsTiming*24,1)*6)
    
dataset.dateAnalysed=dateAnalysed;
dataset.eventsTiming=eventsTiming;
dataset.eventsChannelsOn=eventsChannelsOn;
dataset.eventsSignalValue=eventsSignalValue;
dataset.eventsPeakValue=eventsPeakValue;
dataset.eventsRMSValue=eventsRMSValue;
dataset.beamInfoDataTimed=beamInfoDataTimed;

end %function