source: ETALON/SPESO/analysis/analyseDataRange.m @ 751

%dirname='20140202'
%time_start_of_range='00-00-00';
%time_end_of_range='01-00-00';

function analyseDataRange(dirname,time_start_of_range,time_end_of_range)
    addSPESOpath
    basedir='/sps/sera/delerue/smith-purcell/SPESO/data/';
    cd(basedir)
    
base_save_file_name= [  dirname '/' 'data_' dirname '_from_' time_start_of_range '_to_' time_end_of_range ];
mat_save_file=[base_save_file_name '_save.mat' ];
%%% End of parameters
start_of_day=datenum(dirname,'yyyymmdd');
num_start_of_range=datenum([dirname ' ' time_start_of_range ],'yyyymmdd HH-MM-SS');
  num_end_of_range=datenum([dirname ' ' time_end_of_range ],'yyyymmdd HH-MM-SS');

%num_end_of_range-num_start_of_range

htmlOut=[ '<h1> Range summary : ' dirname ' from ' time_start_of_range ...
          ' until ' time_end_of_range ' </h1>'];
htmlOut=[ htmlOut '<A HREF=summary.html> Daily summary for ' dirname ' </A> <BR/>'];



txt=ls('--color=never','-1',[ dirname '/*_dec.txt.gz' ]);
txt=strrep(txt,char(32),char(10));
txt=[char(10) txt char(10) ]; %we are sure that the text starts and finishes with char(10)
txt=strrep(txt,[char(10) char(10) ], char(10));
txt=strrep(txt,[char(10) char(10) ], char(10));
allPos=strfind(txt,char(10));
nFiles=size(allPos,2)-1
iFile=0;
allData=0;
allSignals=0;
channelsToPlot=[ 8 10 12 13 14 ];


htmlOut=[ htmlOut ' <TABLE>'  ];
            
        htmlOut=[ htmlOut  '<TR><TD colspan=' num2str(length(channelsToPlot)) '> Average </TD></TR><TR>' ];
        for iloop=channelsToPlot
        htmlOut=[ htmlOut '<TD> <IMG width=400 SRC=../' base_save_file_name '_data_' num2str(iloop) '.png></TD>' ];
        end %for
            htmlOut=[ htmlOut ' </TR>'  ];

%%% Check the number of files to be analysed
iFileToAnalyse=0;
for iloop=2:size(allPos,2)    
  date_txt=txt(allPos(1,iloop-1)+21:allPos(1,iloop-1)+21+18);
    date_array(iloop-1)=datenum(date_txt,'yyyy-mm-dd_HH-MM-SS');
    if ((date_array(iloop-1)>=num_start_of_range)&&(date_array(iloop-1)<=num_end_of_range))
        iFileToAnalyse=iFileToAnalyse+1;
    end %if file OK
end %for each file

            matFileAlreadyDone=0;
    if (exist(mat_save_file)==2)
        disp(['Mat file exists ' mat_save_file ]);
        mat=load(mat_save_file);
        if (mat.iFile==iFileToAnalyse)
            matFileAlreadyDone=1;
        end
    else
        disp(['Mat file does not exist ' mat_save_file ]);
    end

%%% Analyse each file to be analysed
if (matFileAlreadyDone==0)

for iloop=2:size(allPos,2)    
    %    disp(['dir:' txt(allPos(1,iloop-1):allPos(1,iloop)) '/' ])
  date_txt=txt(allPos(1,iloop-1)+21:allPos(1,iloop-1)+21+18);
%    datenum(date_txt,'yyyy-mm-dd_HH-MM-SS')
%    datenum(date_txt,'yyyy-mm-dd_HH-MM-SS')-start_of_day
    date_array(iloop-1)=datenum(date_txt,'yyyy-mm-dd_HH-MM-SS');
    if ((date_array(iloop-1)>=num_start_of_range)&&(date_array(iloop-1)<=num_end_of_range))
        iFile=iFile+1;
        disp(['dir: ' num2str(iFile) ' => ' txt(allPos(1,iloop-1)+1:allPos(1,iloop)) ])        
        [data,signals]=analyseFile(txt(allPos(1,iloop-1)+1:allPos(1,iloop)));
        if(iFile==1) %first file to be analysed
        allData=data;
        allSignals=zeros(nFiles,size(signals,1),size(signals,2));
        else
            allData=allData+data;
        end
        allSignals(iFile,:,:)=signals;
        base_save_file_name_file=strrep(strrep(txt(allPos(1,iloop-1)+1:allPos(1,iloop)),'_dec.txt.gz',''),char(10),'');
            htmlOut=[ htmlOut ' <TR><TD><A HREF=../' base_save_file_name_file ...
                      '_dec.html> ' base_save_file_name_file ' </A></TD></TR><TR>'  ];
        for iloop=channelsToPlot
        htmlOut=[ htmlOut '<TD> <IMG width=400 SRC=../' base_save_file_name_file '_channel_' num2str(iloop) '.png></TD>' ];
                  
        end %for
            htmlOut=[ htmlOut ' </TR>'  ];
                  end %if file OK
end %for each file
            htmlOut=[ htmlOut ' </TABLE>'  ];
    save(mat_save_file,'iFile','dirname','time_start_of_range','time_end_of_range','allData','allSignals');

    rangeInfo=[ 'Data from ' dirname ' at ' time_start_of_range ' until ' time_end_of_range ];

%size(allPos)
figure(10041)
clf
hold on
plot(date_array-start_of_day)
line([ 0 size(allPos,2)-1 ],[num_start_of_range-start_of_day num_start_of_range-start_of_day],'Color','g')
line([ 0 size(allPos,2)-1 ],[num_end_of_range-start_of_day num_end_of_range-start_of_day],'Color','r')
hold off
ylabel('time for start of day (fraction of day)')
xlabel('file number')
    title({ rangeInfo; [ 'Time versus files - ' 'nFiles= ' num2str(iFile) ]});
    

print('-dpng',[base_save_file_name '_day.png' ]);
        htmlOut=[ htmlOut ' <IMG width=400 SRC=' base_save_file_name ...
                  '_day.png> <BR/>' ];

if (iFile>0) 
    plotRange(rangeInfo,allData,iFile,allSignals,base_save_file_name);
end    

        fid=fopen([base_save_file_name '_core.html' ],'w');
        fprintf(fid,'%s',htmlOut);
        fclose(fid);
end % if matFileAlreadyDone==0        

end %function