function reportECG(ECG_w, detailLevel, report_mode, win_lengths, report_format, filename)

% This function creates a report of the signals handled by the ECG wrapper
% object ECG_w. The report includes several views of the signals, from a
% wide to a narrow scale. Some aspects of the report can be configured as
% the detail degree, the length of each time scale and the report format.
% 
% Arguments:
%     
%     +ECG_w: [ECGwrapper object] REQUIRED. The signal handler.
% 
%     +detailLevel: [char] OPTIONAL. The report detail level:
%                   'HighDetail', 'MediumDetail', 'LowDetail'. A higher
%                   detail level means report the whole recording at every
%                   time resolution defined in "win_lengths". High
%                   resolution also means larger reports. Default 'LowDetail'. 
% 
%     +report_mode: [char] OPTIONAL. Information from other tasks like QRS
%                   detection/delineation/classification added to the
%                   signals in case available mode. Possible values are:
%                   'full' 'ECG only' 'QRS detection' 'Wave delineation'
%                   'Heartbeat classification' 
%                   Default "ECG only". 
% 
%     +win_lengths: [numeric] OPTIONAL. The amount and size (in seconds) of
%                   each scale length present in the report. Default 
%                   [60*60 30*60 60 7] it means 1 hour - 30 min - 1 min
%                   and 7 seconds.
% 
%     +report_format: [char] OPTIONAL. The report format of the document.
%                  Default PDF.
% 
%     +filename: [char] OPTIONAL. The report filename. Default same folder
%                of the original recording handled by ECG_w.
% 
% 
% Author: Mariano Llamedo Soria (llamedom at {electron.frba.utn.edu.ar; unizar.es}
% Version: 0.1 beta
% Last update: 17/6/2014
% Birthdate  : 15/8/2012



%% constants

cKnownModes = {'full' 'ECG only' 'QRS detection' 'Wave delineation' 'Heartbeat classification' };
cKnownFormats = {'pdf' 'eps' 'png' 'tiff' 'jpg' 'bmp' };
cDetails = {'HighDetail' 'MediumDetail' 'LowDetail'};

%% argument parsing

if( nargin < 5 || ~any(strcmpi(cKnownFormats, report_format)) )
    report_format = 'pdf';
end

heasig = ECG_w.ECG_header;

if( nargin < 3 || isempty(intersect(cKnownModes, report_mode)) )
    report_mode = 'ECG only';
end

file_suffix = rowvec(colvec(char(report_mode)'));
prev_length = length(file_suffix);
file_suffix = strrep(file_suffix, '  ', ' ');
while( length(file_suffix) ~= prev_length  )
    prev_length = length(file_suffix);
    file_suffix = strrep(file_suffix, '  ', ' ');
end
file_suffix = ['_' file_suffix];

if( any(strcmpi(report_mode, 'full')) )
    report_mode = cKnownModes(2:end);
end

if( nargin < 6 || isempty(filename) )

    [~, file_recname ]= fileparts(ECG_w.recording_name);
    report_path = ECG_w.output_path;
    
    if( isfield(heasig, 'recname') )
        filename = [report_path heasig.recname file_suffix '.' report_format];
    else
        filename = [report_path file_recname file_suffix '.' report_format];
    end

end

if( nargin < 2 || ~any(strcmpi(cDetails, detailLevel)) )
    detailLevel = 'LowDetail'; % seconds
end

if( nargin < 4 || isempty(win_lengths) || ~isnumeric(win_lengths) )
    win_lengths = [60*60 30*60 60 7]; % seconds
end

if( ~isobject(ECG_w) || ~isa(ECG_w, 'ECGwrapper') )
    error('reportECG:BadArg', 'Need an ECGwrapper object to work.')
end

init_ghostscript();

win_lengths = bsxfun( @min, win_lengths, floor(heasig.nsamp/heasig.freq) );
win_lengths = sort(unique(win_lengths), 'descend');

% Activate the progress_struct bar.
pb = progress_bar('Report ECG');

pb.Loops2Do = length(win_lengths);

ii = 1;
fig_hdl = [];

for win_length = win_lengths
    
    pb.start_loop();
    
    if( (win_length * heasig.freq) > heasig.nsamp )
        win_length = round(heasig.nsamp / heasig.freq);
    end
    
    switch detailLevel
        case 'LowDetail'
            max_cant_figs = 10;
        case 'MediumDetail'
            max_cant_figs = 30;
        case 'HighDetail'
            max_cant_figs = realmax;
    end
    
    cant_figs = min(max_cant_figs, ceil(heasig.nsamp / win_length / heasig.freq));
    starts = linspace(0, round(heasig.nsamp / heasig.freq) - win_length, cant_figs);

    str_aux = ['Reporting for window size: ' Seconds2HMS(win_length) ];
    pb1 = pb.AddInerLoop( str_aux );
    pb1.Loops2Do = cant_figs;
    
    jj = 1;
    
    for start_time = starts
    
        pb1.start_loop();

        fig_hdl = prepare_fig_hdl( fig_hdl );
        
        
        eval_str = 'plot_ecg_strip(ECG_w, ''PrettyPrint'', true, ''Figure_handle'', fig_hdl, ''Start_time'', start_time, ''End_time'', start_time + win_length ';
        
        if( any(strcmpi(report_mode, 'QRS detection')) )
            eval_str = [ eval_str ',''QRS_locations'', ECG_w' ];
        end
            
        if( any(strcmpi(report_mode, 'Wave delineation')) )
            eval_str = [ eval_str ',''ECG_delineation_single_lead'', ECG_w, ''ECG_delineation_multilead'', ECG_w' ];
        end
            
        if( any(strcmpi(report_mode, 'Heartbeat classification')) )
            eval_str = [ eval_str ',''Heartbeat_classification'', ECG_w' ];
        end

        eval_str = [ eval_str ');'];
        
        eval(eval_str);
        
        set(fig_hdl, 'Visible', 'off');
        
        str_aux2 = [ str_aux ' - Start: ' Seconds2HMS(start_time) ];
        pb1.checkpoint(str_aux2);
        
        switch( report_format )
            case {'pdf' 'tiff' }
                
                if( jj == 1 )
%                     export_fig(filename, '-nocrop', '-bookmark', Seconds2HMS(win_length), '-append', ['-' report_format], fig_hdl);
                    export_fig(filename, '-nocrop', '-append', ['-' report_format], fig_hdl);
                else
%                     export_fig(filename, '-nocrop', '-bookmark', ['    ' Seconds2HMS(start_time)], '-append', ['-' report_format], fig_hdl);
                    export_fig(filename, '-nocrop', '-append', ['-' report_format], fig_hdl);
                end
                
            otherwise
                [aux_path, aux_filename] = fileparts(filename);
                aux_filename = [aux_path filesep aux_filename];
                aux_filename = sprintf( '%s_%04d.%s', aux_filename, ii, report_format);
                export_fig(aux_filename, '-nocrop', ['-' report_format], fig_hdl);
        end
        
        pb1.checkpoint(str_aux2);

        ii = ii + 1;
        jj = jj + 1;
        
        pb1.end_loop();
        
    end
    
    pb.end_loop();
    
end

clear pb


function fig_hdl = prepare_fig_hdl( fig_hdl )

    paper_size = [0 0 29.7 21.0 ];

    if( ishandle(fig_hdl) )
        close(fig_hdl)
    end
    fig_hdl = figure('Visible', 'off');
    set(fig_hdl, 'NumberTitle', 'off');
    set(fig_hdl, 'PaperUnits', 'centimeters');
    set(fig_hdl, 'PaperType', 'A4');
    set(fig_hdl, 'PaperOrientation', 'landscape');
    set(fig_hdl, 'PaperPosition', paper_size);
    set(fig_hdl, 'Units', 'centimeters');
    set(fig_hdl, 'Position', paper_size);