%% Plots multidimentional signal in mosaic style
% This function plots several subplots in the same figure in order to do a
% mosaic with the different leads available in ECG. Annotations can be
% provided to each or for all the mosaic.
%   
% Prototype
% 
%   [ ECG_hdl axes_hdl fig_hdl all_yranges ] = plot_ecg_mosaic( ECG, varargin )
% 
% Arguments:
% 
%     +ECG: [numeric or cell] REQUIRED
%           
%           [numeric]: signal matrix of dimension [sig_length sig_size
%           repetitions_size] where:
%             - sig_length: time length in samples
%             - sig_size: number of ECG leads or number of signals.
%             - repetitions_size: number of repetitions of the same
%             signals. Typically used when time-synchronized events, like
%             heartbeats.  
%           
%           [cell]: cell array of length repetitions_size, where each cell
%           is (probably a time alligned event) a signal of dimension
%           [sig_length sig_size]
% 
%     +QRS_locations: [numeric] OPTIONAL. Default values enclosed in ()
%                 Synchronization sample. In ECG context, this values are
%                 the QRS fiducial point. (empty)
% 
%     +WinSize: [numeric] OPTIONAL. Default values enclosed in ()
%           Width of the window around each fiducial point provided in
%           QRS_locations. (empty)
% 
%     +ECG_header: [struct] OPTIONAL. 
% 
%             Description of the ECG typically available in the
%             header. Structure with fields:
% 
%               -freq: Sampling rate in Hz. (1)
% 
%               -nsig: Number of ECG leads. (size(ECG,2))
% 
%               -nsamp: Number of ECG samples. (size(ECG,1))
% 
%               -adczero: ADC offset (e.g. 2^(adc_res-1) when
%               using unsigned integers). ( repmat(0, ECG_header.nsig , 1) ) 
% 
%               -adcgain: ADC gain in units/adc_sample
%               (typically uV/adc_unit). ( repmat(1, ECG_header.nsig , 1) )
% 
%               -units: Measurement units. ( repmat('uV', ECG_header.nsig , 1) )
% 
%               -desc: Signal description. ( num2str(colvec(1:ECG_header.nsig)) )
%                 
%     +MaxECGrange: [numeric or string] OPTIONAL. Force a vertial range in order to
%                     ease visual comparison of signals in the mosaic. 
%                   [string] 
%                     'max': force the maximum range to be the range
%                            for all mosaics.
%                     'min', 'mean', 'median': are also available options.
%                     'none': Each mosaic with a different range. (Default).
%     
%     +RowsCols: [numeric] OPTIONAL. Number of rows and columns of the
%                 mosaic. If ommited or if rows * cols ~= ECG_header.nsig, these values are
%                 automatically adapted to the best fit mosaic in relation to the
%                 aspect ratio of the screen.   
%     
%     +FigureHdl: [figure handle] OPTIONAL. Choose the figure to be produced
%                the mosaic. (handle produced by gcf)
% 
%     +ECG_delineation: [struct] OPTIONAL. Default values enclosed in ()
%               Annotation struct generated with wavedet. 
% 
%     +ECG_annotations: [cell] OPTIONAL. Default values enclosed in ()
%               Annotations to be included in the mosaic. The funcion
%               accepts 2 type of annotations: points and lines. 
% 
% Output:
% 
%     + ECG_hdl: handle to the plotted signals.
% 
%     + axes_hdl: handle to the axes.
% 
%     + fig_hdl: handle to fig.
% 
%     + all_yranges: vertical ranges of the plotted signals.
% 
% 
% Example:
% 
% This example makes a synchronized plot of several random events (+1 0 -1)
% ocurryng randomly with additive white gaussian noise.  
% win_size = 100; 
% sig_samp = 10000;
% sig_size = 12;
% event_size = 50;
% x = 0.1*randn(sig_samp,sig_size); 
% event_locations = randsample(win_size:sig_samp-win_size, event_size);
% x(event_locations-1,:) = x(event_locations-1,:) + 1;
% x(event_locations+1,:) = x(event_locations+1,:) - 1;
% x_packed = pack_signal(x, event_locations, win_size);    
% 
% figure(1)
% % estimation of the signal averaged event
% plot_ecg_mosaic( mean(x_packed,3) );
% figure(2)
% % visualization of all events. In this case previous pack_signal call is
% % not needed. 
% plot_ecg_mosaic(x, 'QRS_locations', event_locations, 'WinSize', win_size);
% 
% figure(3)
% % introducing several kind of marks to the plot
% 
% h_line = cell(sig_size,7);
% h2_line = cell(sig_size,7);
% v_line = cell(sig_size,7);
% v2_line = cell(sig_size,7);
% point = cell(sig_size,7);
% a_line = cell(sig_size,7);
% 
% h_line(:,1) = {'line'};
% h_line(:,2) = { [ { 'String'                     'LineStyle' 'LineWidth' 'Color' 'TextColor' }; ...
%                   { 'horizontal line text'       '--'         1.5         'r'     'r'        } ]'}; 
% h_line(1:sig_size, [6 7]) = num2cell( repmat(-0.5,sig_size,2) );
% 
% h2_line(:,1) = {'line'};
% h2_line(:,2) = { [ { 'String'           'LineStyle' 'LineWidth' 'Color' 'TextColor' }; ...
%                   { 'other h-line'      '--'         1.5         'm'     'm'        } ]'}; 
% h2_line(1:sig_size, 4:7) = num2cell( [ repmat(60,sig_size,1) repmat(70,sig_size,1) repmat(0.5,sig_size,2) ] );
% 
% v_line(:,1) = {'line'};
% v_line(:,2) = { [ { 'String'             'LineStyle' 'LineWidth' 'Color' 'TextColor' }; ...
%                   { 'vertical line text' '--'         1.5         'g'     'g'        } ]'}; 
% v_line(1:sig_size, [4 5]) = num2cell( repmat(20,sig_size,2) );
% 
% v2_line(:,1) = {'line'};
% v2_line(:,2) = { [ { 'String'      'LineStyle' 'LineWidth' 'Color' 'TextColor' }; ...
%                   { 'other v-line' '--'         1.5         'b'     'b'        } ]'}; 
% v2_line(1:sig_size, 4:7) = num2cell( [ repmat(80,sig_size,2) repmat(-0.8,sig_size,1) repmat(-0.3,sig_size,1) ] );
% 
% point(:,1) = {'point'};
% point(:,2) = { [ { 'String'    'Color'       'TextColor'   }; ...
%                  { 'one-point' [0.2 0.3 0.4] [0.2 0.3 0.4] } ]'}; 
% point(1:sig_size,4) = num2cell( repmat(50,sig_size,1) );
% 
% a_line(:,1) = {'line'};                    
% a_line(:,2) = { [ { 'String'    'LineStyle' 'LineWidth' 'Color' 'TextColor' }; ...
%                   { 'line text' '--'         1.5         'k'     'k'        } ]'}; 
% a_line(1:sig_size,4:7) = num2cell( [ repmat(30,sig_size,1) repmat(40,sig_size,1) repmat(0.5,sig_size,1) repmat(-0.5,sig_size,1) ] );
%        
% aux_anns = cat(3,h_line,v_line,h2_line,v2_line,point,a_line);
% 
% plot_ecg_mosaic(mean(x_packed,3), 'ECG_annotations', aux_anns );
% 
% See also plot_ecg_strip
% 
% Author: Mariano Llamedo Soria llamedom@electron.frba.utn.edu.ar
% Last update: 19/10/2014
% Birthdate  : 16/2/2012
% Copyright 2008-2015
function [ ECG_hdl axes_hdl fig_hdl all_yranges ] = plot_ecg_mosaic( ECG, varargin )

    %% constants

    cWaveAnnotations = { 'Pon' 'P' 'Poff' 'QRSon' 'Q' 'R' 'S' 'QRSoff' 'Ton' 'T' 'Toff' };
    waveMode_idx = 1;
    lWaveAnnotations = length(cWaveAnnotations);
    
    cModesMaxYrange = {'max','min','mean','median','none'};

    % bgTextColor = [251 248 230]/255;
    ColorMinorGrid = [0.8 0.8 0.8];

    ann_typ_idx = 1;
    plot_options_str_idx = 2;
    x1_idx = 4;
    x2_idx = 5;
    y1_idx = 6;
    y2_idx = 7;
    n_cant_idx = 7;

    PwaveColor =   [210 255 189]/255;
    QRScplxColor = [255 200 255]/255;
    TwaveColor =   [255 240 170]/255;
    
    DelineationNotEditedColor = [ 0.8 0.8 0.8];
    DelineationEditedColor = [ 1 0 0 ];
    
    %% parse input

    p = inputParser;   % Create instance of inputParser class.
    p.addParamValue('ECG_header', [], @(x)(isstruct(x)) );
    p.addParamValue('QRS_locations', [], @(x)( isnumeric(x) && all(x > 0) ) );
    p.addParamValue('WinSize', [], @(x)( all(isnumeric(x)) ) );
    p.addParamValue('ECG_annotations', [], @(x)( isnumeric(x) || iscell(x) ) );
    p.addParamValue('ECG_delineation', [], @(x)(isstruct(x)) );
    p.addParamValue('MaxECGrange', [], @(x)(isnumeric(x) || ischar(x) && any(strcmpi(cModesMaxYrange, x)) ) );
    p.addParamValue('RowsCols', [], @(x)( all(isnumeric(x)) && all(x > 0) ) );
    p.addParamValue('FigureHdl', [], @(x)( ishandle(x) ) );
    p.addParamValue('AllowEdition', false, @(x)(islogical(x)) );

    try
        p.parse( varargin{:} );
    catch MyError
        fprintf(2, 'Incorrect argval/argvalue, use:\n' );
        fprintf(2, 'plot_ecg_mosaic( ECG, ''arg_name'', ''arg_value'');\n\n' );
        fprintf(2, 'Valid arg_name are:\n' );
        fprintf(2, '+ ECG_header\n' );
        fprintf(2, '+ QRS_locations\n' );
        fprintf(2, '+ ECG_annotations\n' );
        fprintf(2, '+ ECG_delineation\n' );
        fprintf(2, '+ WinSize\n' );
        fprintf(2, '+ MaxECGrange\n' );
        fprintf(2, '+ RowsCols\n' );
        fprintf(2, '+ FigureHdl\n' );
        fprintf(2, '\nOr run ''doc plot_ecg_mosaic'' for details.\n\n' );
        rethrow(MyError);    
    end

    QRS_locations = p.Results.QRS_locations;
    winSize = p.Results.WinSize;
    ECG_header = p.Results.ECG_header;
    annotations = p.Results.ECG_annotations;
    positions = p.Results.ECG_delineation;
    all_yranges = p.Results.MaxECGrange;
    rowscols = p.Results.RowsCols;
    fig_hdl = p.Results.FigureHdl;
    bAllowEdition = p.Results.AllowEdition;

    delete(p)

    if( isempty(fig_hdl) )
        fig_hdl = gcf;
    else
        figure(fig_hdl)
    end
    clf('reset')

    if( iscell(ECG) )

        [lECG, cant_sig ] = size(ECG{1});
        cant_realizations = length(ECG);

        ECG = cell2mat(reshape(ECG, 1, 1, cant_realizations));

    elseif( isnumeric(ECG) )

        if( ~isempty(QRS_locations) )
            ECG = pack_signal(ECG, QRS_locations, winSize);    
        end

        [lECG, cant_sig, cant_realizations] = size(ECG);

    end

    if( isempty(all_yranges) )
        bYrange_given = false;
        all_yranges = nan(cant_sig,2);
        rangeMode = 'none';
    else
        if(ischar(all_yranges))
            rangeMode = lower(all_yranges);
            all_yranges = nan(cant_sig,2);
            bYrange_given = false;
        else
            bYrange_given = true;
        end
    end

    % plot arrows when range was not given.
    bArrow = ~bYrange_given;

    if( isempty(ECG_header) )
        ECG_header.freq = nan;
        ECG_header.nsamp = lECG;
        ECG_header.nsig = cant_sig;
        ECG_header.desc = num2str(colvec(1:cant_sig));
        ECG_header.units = repmat('#', ECG_header.nsig , 1);
        ECG_header.gain = repmat(1, ECG_header.nsig , 1);
        ECG_header.adczero = repmat(0, ECG_header.nsig , 1);

        bGridScale = false;
    else
        if( ECG_header.nsig ~= cant_sig )
            ECG_header.nsig = cant_sig;
            ECG_header.units = repmat(ECG_header.units(1,:), ECG_header.nsig , 1);
            ECG_header.gain = repmat(ECG_header.gain(1), ECG_header.nsig , 1);
            ECG_header.adczero = repmat(ECG_header.adczero(1), ECG_header.nsig , 1);
        end
        bGridScale = true;

        [ECG ECG_header] = ADC2units(ECG, ECG_header, 'uv');
    end

    if( isempty(positions) )
        all_annotations = annotations;
    else

        if( isempty(annotations) )
            all_annotations = update_annotations_from_positions(positions);
        else
            all_annotations = cat(3, annotations, update_annotations_from_positions(positions) );
        end
        
    end

    if( isempty(all_annotations) )

        all_annotations = cell(cant_sig,1);
        anns_provided = false;

    else
        if(iscell(all_annotations))

            anns_provided = true;

            [lannotations, cant_params, cant_type_anns] = size(all_annotations);

            if(cant_params ~= n_cant_idx )
                error( 'Invalid number of params: \n\nCol1: Ann point location\nCol2: Ann line start\nCol3: Ann line end\nCol4: plot LineSpec or Prop Name/Val\nCol5: Ann legend\nCol6: Annotation type\n' );
            end

            if(lannotations == 1)
                %same ann for all signals
                all_annotations = repmat(all_annotations, cant_sig, 1);
            elseif( lannotations > 1 && lannotations < cant_sig )
                error([ 'Invalid number of annotations: ' num2str(lannotations) ' for ' num2str(cant_sig) ' signals.'  ])
            end

        else
            error('Invalid annotation format.')
        end
    end

    if( isempty(rowscols) || prod(rowscols) < cant_sig )

        prev_units = get(fig_hdl,'units');
        prev_size = get(fig_hdl,'outerposition');

        set(fig_hdl,'Visible', 'off');
        set(fig_hdl,'units','normalized');
        set(fig_hdl,'outerposition',[0 0.03 1 0.97]);
        set(fig_hdl,'units','pixels');
        screen_size = get(fig_hdl,'outerposition');

        % restore originals
        set(fig_hdl,'outerposition',prev_size);
        set(fig_hdl,'units',prev_units);
        set(fig_hdl,'Visible', 'on');

        target_aspect = screen_size(3)/screen_size(4);
        aux_val = [1:cant_sig; ceil(cant_sig./(1:cant_sig)) ] ;
        bExact = (cant_sig - prod(aux_val)) == 0;
        if( any(bExact(2:end-1)) )
            % non-trivial solution
            aux_val = aux_val(:, bExact);
        end
        aux_aspects = aux_val(1,:)./aux_val(2,:);
        [~, aux_idx] = sort(abs(aux_aspects - target_aspect));
        n_rows = aux_val(1,aux_idx(1));
        n_cols = aux_val(2,aux_idx(1));
    else
        n_rows = rowscols(1);
        n_cols = rowscols(2);
    end

    %%

    db_status = dbstatus();
    bRestoreErrorStatus = false;
    if( length(db_status) > 1 && strcmpi(db_status(end).cond, 'caught error')  )
        dbclear if caught error
        bRestoreErrorStatus = true;
    end

    ECG_hdl = [];
    axes_hdl = nan(cant_sig,1);
    anns_hdl = [];

    if( isnan(ECG_header.freq) )
        xTimeGridSpacing = round(lECG/5); % mseg
        xTimeGrid = 0:xTimeGridSpacing:lECG;
    else
        xTimeGridSpacing = 50; % mseg
        xTimeGrid = 0:(xTimeGridSpacing/1e3*ECG_header.freq):lECG;
    end


    plotXmin = 1;
    plotXmax = lECG;
    plotXrange = plotXmax - plotXmin;
    xTextOffset = 0.005*plotXrange;

    if(anns_provided)
        lineWidth = 1.3;
    else
        lineWidth = 1;
    end

    v_space = 0.02/(n_rows+1);
    h_space = 0.02/(n_cols+2);
    axe_width = 0.98/n_rows;
    axe_height = 0.98/n_cols;
    axe_x = h_space;
    axe_y = 1-axe_height-v_space;

    for ii = 1:cant_sig

        txthdl = [];

        axes_hdl(ii) = axes();

        % AUX ECG plot

        % ECG plot
        ECG_hdl = [ ECG_hdl; plot( squeeze(ECG(:,ii,:)), 'LineWidth', lineWidth )];

        if(~bYrange_given)
            all_yranges(ii,:) = rowvec(get(axes_hdl(ii), 'ylim'));
        end

        set(axes_hdl(ii), 'xlim', [1 lECG]);

        prev_units = get(axes_hdl(ii),'units');
        set(axes_hdl(ii),'units','pixels');
        axis_position = get(axes_hdl(ii),'Position');
        set(axes_hdl(ii),'units',prev_units);

        set(axes_hdl(ii), 'Box', 'off' );
        set(axes_hdl(ii), 'Xtick', [] );
        set(axes_hdl(ii), 'Ytick', [] );
        set(axes_hdl(ii), 'Xcolor', [1 1 1] );
        set(axes_hdl(ii), 'Ycolor', [1 1 1] );

        [x_loc, y_loc] = ind2sub([n_rows n_cols], ii);
        axe_x = h_space * x_loc + axe_width * (x_loc-1);
        axe_y = 1-( v_space * y_loc + axe_height * y_loc);

    %     aux_val = get(axes_hdl(ii), 'Position' );
    %     set(axes_hdl(ii), 'Position', [ aux_val(1:2) - (k_magnify-1)/2*aux_val(3:4) k_magnify*aux_val(3:4) ] );
        set(axes_hdl(ii), 'Position', [ axe_x axe_y axe_width axe_height ] );

        
    end

    max_yrange = [min(all_yranges(:,1)) max(all_yranges(:,2)) ];
    % min_range = min(all_yranges(:,2) - all_yranges(:,1));

    for ii = 1:cant_sig

        set(fig_hdl, 'CurrentAxes', axes_hdl(ii));
        hold(axes_hdl(ii), 'on');

        if(bGridScale) 
            % grid plot

            yTimeGridSpacing = (all_yranges(ii,2) - all_yranges(ii,1)) / 5; % uV

            yTimeGrid = max_yrange(1):yTimeGridSpacing:max_yrange(2);

            lyTimeGrid = length(yTimeGrid);
            if( lyTimeGrid > 10 )
                yTimeGrid = yTimeGrid( yTimeGrid >= all_yranges(ii,1) & yTimeGrid <= all_yranges(ii,2) );
            end

            aux_hdl = plot( repmat([1; lECG], 1, length(yTimeGrid)), repmat(yTimeGrid,2,1), 'Color', ColorMinorGrid, 'LineStyle', ':' );
            aux_hdl = [aux_hdl; plot( repmat(xTimeGrid,2,1), repmat(colvec(max_yrange), 1, length(xTimeGrid) ), 'Color', ColorMinorGrid, 'LineStyle', ':' ) ];

            set(aux_hdl,'ButtonDownFcn',@WindowButtonDown);            

        end

        % zero line
        aux_hdl = plot( [ xTimeGrid(1) xTimeGrid(end)], [0 0], '-k' );

        set(aux_hdl,'ButtonDownFcn',@WindowButtonDown);            
        
        hold(axes_hdl(ii), 'off');

    end

    if(~bYrange_given)
        switch(rangeMode)
            case 'max'
                all_yranges = repmat([ min(all_yranges(:,1)) max(all_yranges(:,2)) ], cant_sig, 1   );
            case 'min'
                all_yranges = repmat([ max(all_yranges(:,1)) min(all_yranges(:,2)) ], cant_sig, 1   );
            case 'mean'
                all_yranges = repmat([ mean(all_yranges(:,1)) mean(all_yranges(:,2)) ], cant_sig, 1   );
            case 'median'
                all_yranges = repmat([ median(all_yranges(:,1)) median(all_yranges(:,2)) ], cant_sig, 1   );
        end
    end


    ECG_lgd = [];
    
    for ii = 1:cant_sig
    %     aux_val = get(axes_hdl(ii), 'ylim');
    %     aux_val2 = (all_yranges - abs(diff(aux_val)))/2;
    %     set(axes_hdl(ii), 'ylim', aux_val + [-aux_val2 aux_val2] );
        set(axes_hdl(ii), 'ylim', all_yranges(ii,:) );

        set(fig_hdl, 'CurrentAxes', axes_hdl(ii));
        y_lims = get(axes_hdl(ii), 'ylim');
        x_lims = get(axes_hdl(ii), 'xlim');

        if( ii <= size(ECG_header.desc,1) )
            ECG_lgd = [ECG_lgd; text(x_lims(1) + 0.01*abs(diff(x_lims)), y_lims(2) - 0.05*abs(diff(y_lims)), ECG_header.desc(ii,:), 'EdgeColor', [0 0 0] )];
        else
            ECG_lgd = [ECG_lgd; text(x_lims(1) + 0.01*abs(diff(x_lims)), y_lims(2) - 0.05*abs(diff(y_lims)), ['Sig ' num2str(ii)], 'EdgeColor', [0 0 0] )];
        end

        plotXmin = x_lims(1);
        plotXmax = x_lims(2);
        plotXrange = plotXmax - plotXmin;
        plotYmin = y_lims(1);
        plotYmax = y_lims(2);
        plotYrange = plotYmax - plotYmin;
        yTextOffset = 0.01*plotYrange;

        yTimeGridSpacing = (all_yranges(ii,2) - all_yranges(ii,1)) / 5; % uV
        yTimeGrid = max_yrange(1):yTimeGridSpacing:max_yrange(2);

        aux_idx = find(yTimeGrid > plotYmin & yTimeGrid < plotYmax);

        width_legend = 0.18*plotXrange;
        height_legend = 0.1*plotYrange;
        left_legend = plotXmax - width_legend - 5*xTextOffset;
        bottom_legend = plotYmax - height_legend - 2*yTextOffset;

        % scale X
    %         plot( [-xTimeGridSpacing/2 -xTimeGridSpacing/2 xTimeGridSpacing/2 xTimeGridSpacing/2] + 5/7*lECG, plotYmin + 0.05*plotYrange + [yTextOffset 0 0 yTextOffset], 'Color', [0 0 0], 'LineWidth', 1 );
        if( bArrow )
            ECG_lgd = [ECG_lgd; arrow( [ xTimeGrid(2) plotYmin + 0.05*plotYrange], [ xTimeGrid(3) plotYmin + 0.05*plotYrange], 2, 0.5, [0 0 0], axes_hdl(ii) )];
        end

        if( isnan(ECG_header.freq) )
            aux_str =  [ num2str(xTimeGridSpacing) ' #' ];
        else
            aux_str =  [ num2str(xTimeGridSpacing) ' ms' ];
        end

        ECG_lgd = [ECG_lgd; text( mean(xTimeGrid(2:3)),  plotYmin + 0.05*plotYrange + 5*yTextOffset, aux_str, 'FontSize', 8, 'HorizontalAlignment', 'center') ];

        % scale Y
    %         plot( [ -xTextOffset 0 0 -xTextOffset] + 0.95*lECG, plotYmin + 5/7*plotYrange + [-yTimeGridSpacing/2 -yTimeGridSpacing/2 yTimeGridSpacing/2 yTimeGridSpacing/2], 'Color', [0 0 0], 'LineWidth', 1 );
        if( bArrow )
            ECG_lgd = [ECG_lgd; arrow( [0.95*lECG - xTextOffset; yTimeGrid(aux_idx(2))], [0.95*lECG - xTextOffset; yTimeGrid(aux_idx(3))], 2, 0.5, [0 0 0], axes_hdl(ii) )];
        end

        ECG_lgd = [ECG_lgd; text( 0.95*lECG - 10*xTextOffset, mean(yTimeGrid(aux_idx(2:3))), [ num2str(yTimeGridSpacing) ' ' ECG_header.units(ii,:) ], 'FontSize', 8, 'HorizontalAlignment', 'center', 'Rotation', 90)];

        txthdl = [];

        if( ~isempty(all_annotations{ii}) )

            update_annotations( ii )

        end

        if ~isempty(positions)
            ECG_lgd = [ ECG_lgd; patch([left_legend left_legend [left_legend left_legend]+width_legend left_legend ], [bottom_legend [bottom_legend bottom_legend]+height_legend bottom_legend bottom_legend], [1 1 1], 'EdgeColor', [0 0 0])];
            ECG_lgd = [ ECG_lgd; text( left_legend + width_legend/4, bottom_legend + height_legend/2, 'P', 'FontSize', 8, 'HorizontalAlignment', 'center', 'BackgroundColor', PwaveColor) ];
            ECG_lgd = [ ECG_lgd; text( left_legend + width_legend/2, bottom_legend + height_legend/2, 'QRS', 'FontSize', 8, 'HorizontalAlignment', 'center', 'BackgroundColor', QRScplxColor ) ];
            ECG_lgd = [ ECG_lgd; text( left_legend + width_legend*3/4, bottom_legend + height_legend/2, 'T', 'FontSize', 8, 'HorizontalAlignment', 'center', 'BackgroundColor', TwaveColor ) ];
        end
        
    end

    if(bRestoreErrorStatus)
        dbstop if caught error
    end
    
    my_timer = timer('TimerFcn',@timer_fcn, 'StopFcn', @timer_stop_fcn , 'StartDelay', 10);

    set(fig_hdl, ...
        'WindowScrollWheelFcn', {@ScrollWheel}, ...
        'WindowKeyPressFcn', {@KeyPress}, ...
        'WindowButtonUpFcn', {@WindowButtonUp} ...
        );            

    set(axes_hdl,...
        'ButtonDownFcn', {@WindowButtonDown} ...
        );            


    set(ECG_hdl,...
        'ButtonDownFcn', {@WindowButtonDown} ...
        );            

    if( isempty(positions) )
        for ii = 1:lWaveAnnotations
            manual_edited_position.(cWaveAnnotations{ii}) = nan;
        end
        manual_edited_position = repmat(manual_edited_position, ECG_header.nsig,1);
    else
        manual_edited_position = positions;
    end

    PrevStateWindowButtonMotionFcn = nan;
    bIsDragAllowed = false;
    lead_move_idx = [];
    axes_move = [];
    move_hdl = [];
    sample_move = [];
    bFixedWave = false;
    bPreserveFix = false;
    manual_edited_position_matrix = positions2matrix(manual_edited_position, cWaveAnnotations);
    bDelineationEdited = false;
    
    if( nargout == 0 )
        ECG_hdl = [];
    end
    
%% Internal functions    

    function timer_stop_fcn(obj,event_obj)
         
        if(bPreserveFix)
            % never stop when editing
            start(my_timer)
        end
        
    end

    function timer_fcn(obj,event_obj)
        
        delete(findobj('Tag', 'title_efimero' ));

        if(~bPreserveFix)
            % allow edition of the closer wave
            bFixedWave = false;
        end
        
    end

    function this_annotations = update_annotations_from_positions( this_positions )
        
        [min_vals_pw_x, min_vals_pw_y, max_vals_pw_x, max_vals_pw_y]     = CalcMinMax( this_positions, {'Pon' 'P' 'Poff'} );
        [min_vals_qrs_x, min_vals_qrs_y, max_vals_qrs_x, max_vals_qrs_y] = CalcMinMax( this_positions, {'QRSon' 'qrs' 'QRSoff'} );
        [min_vals_tw_x, min_vals_tw_y, max_vals_tw_x, max_vals_tw_y]     = CalcMinMax( this_positions, {'Ton' 'T' 'Toff'} );

        Pwave = cell(cant_sig,7);
        Pwave_max = cell(cant_sig,7);
        QRScplx = cell(cant_sig,7);
        Twave = cell(cant_sig,7);
        QRScplx_point = cell(cant_sig,7);
        Qwave = cell(cant_sig,7);
        Rwave = cell(cant_sig,7);
        Swave = cell(cant_sig,7);
        Twave_max = cell(cant_sig,7);

        Pwave(:,1) = {'patch'};
        Pwave(:,2) = { [ { 'FaceColor' }; ...   
                         { PwaveColor  } ]'}; 
        Pwave(:,[x1_idx y1_idx x2_idx y2_idx]) = [ num2cell(min_vals_pw_x), num2cell(min_vals_pw_y), num2cell(max_vals_pw_x), num2cell(max_vals_pw_y)];

        k_point_color = 0.9;
        
        Pwave_max(:,1) = {'point'};
        Pwave_max(:,2) = { [ { 'String' 'Marker' 'MarkerSize' 'Color' }; ...
                         {     'P'      'd'       5            k_point_color*PwaveColor    } ]'}; 
        Pwave_max(:,4) = {this_positions(:).P};
        
        QRScplx(:,1) = {'patch'};
        QRScplx(:,2) = { [ { 'FaceColor' }; ...   
                         { QRScplxColor  } ]'}; 
        QRScplx(:,[x1_idx y1_idx x2_idx y2_idx]) = [ num2cell(min_vals_qrs_x), num2cell(min_vals_qrs_y), num2cell(max_vals_qrs_x), num2cell(max_vals_qrs_y)];

%         QRScplx_point(:,1) = {'point'};
%         QRScplx_point(:,2) = { [ { 'String' 'Marker' 'MarkerSize' 'Color' }; ...
%                          {          'qrs'      'd'       5            k_point_color*QRScplxColor    } ]'}; 
%         QRScplx_point(:,4) = {this_positions(:).qrs};
        QRScplx_point = [];
        
        Qwave(:,1) = {'point'};
        Qwave(:,2) = { [ { 'String' 'Marker' 'MarkerSize' 'Color' }; ...
                         { 'Q'      'd'       5            k_point_color*QRScplxColor    } ]'}; 
        Qwave(:,4) = {this_positions(:).Q};
        
        Rwave(:,1) = {'point'};
        Rwave(:,2) = { [ { 'String' 'Marker' 'MarkerSize' 'Color' }; ...
                         { 'R'      'd'       5            k_point_color*QRScplxColor    } ]'}; 
        Rwave(:,4) = {this_positions(:).R};
        
        Swave(:,1) = {'point'};
        Swave(:,2) = { [ { 'String' 'Marker' 'MarkerSize' 'Color' }; ...
                         { 'S'      'd'       5            k_point_color*QRScplxColor    } ]'}; 
        Swave(:,4) = {this_positions(:).S};
        
        Twave(:,1) = {'patch'};
        Twave(:,2) = { [ { 'FaceColor' }; ...   
                         { TwaveColor  } ]'}; 
        Twave(:,[x1_idx y1_idx x2_idx y2_idx]) = [ num2cell(min_vals_tw_x), num2cell(min_vals_tw_y), num2cell(max_vals_tw_x), num2cell(max_vals_tw_y)];

        Twave_max(:,1) = {'point'};
        Twave_max(:,2) = { [ { 'String' 'Marker' 'MarkerSize' 'Color' }; ...
                         {     'T'      'd'       5            k_point_color*TwaveColor    } ]'}; 
        Twave_max(:,4) = {this_positions(:).T};
        
        this_annotations = cat(3, Pwave_max, QRScplx_point, Qwave, Rwave, Swave, Twave_max, Pwave, QRScplx, Twave);
        
    end
    
    function update_annotations( lead_idx )
        
        x_lims = get(axes_hdl(lead_idx), 'xlim');
        x_range = abs(diff(x_lims));
        y_lims = all_yranges(lead_idx,:);
        y_range = abs(diff(y_lims));

        for jj = 1:cant_type_anns

            switch( all_annotations{lead_idx,ann_typ_idx,jj} )

                case 'point'

                    if( ~isnan(all_annotations{lead_idx,x1_idx,jj}) )

                        aux_hdl = text_arrow(all_annotations{lead_idx,x1_idx,jj}, ECG(all_annotations{lead_idx,x1_idx,jj}, lead_idx), '', all_annotations{lead_idx,plot_options_str_idx,jj}, axes_hdl(lead_idx) );

                        set(aux_hdl, 'Tag', num2str(lead_idx) );
                        
                        txthdl = [ txthdl aux_hdl ];

                        uistack(aux_hdl,'bottom');

                        if(bAllowEdition)
                            set(aux_hdl,...
                                'ButtonDownFcn', {@WindowButtonDown} ...
                                );            
                        end

%                         aux_hdl = eval( [ 'plot( all_annotations{' num2str(lead_idx) ',' num2str(x1_idx) ',' num2str(jj) '}, ECG(all_annotations{' num2str(lead_idx) ',' num2str(x1_idx) ',' num2str(jj) '},' num2str(lead_idx) ',1)  ' all_annotations{lead_idx,plot_options_str_idx,jj} ')']);
%                         anns_hdl = [ anns_hdl; colvec(aux_hdl)];
%                         txthdl = [ txthdl; text( all_annotations{lead_idx,x1_idx,jj} + 0.03*x_range, ECG(all_annotations{lead_idx,x1_idx,jj}, lead_idx, 1), all_annotations{lead_idx,text_idx,jj}, 'BackgroundColor', bgTextColor )];
                    end

                case 'line'

                    if( isempty(all_annotations{lead_idx,x1_idx,jj}) )
                        all_annotations{lead_idx,x1_idx,jj} = repmat(x_lims(1),1,length(all_annotations{lead_idx,y1_idx,jj}));
                    end
                    if( isempty(all_annotations{lead_idx,x2_idx,jj}) )
                        all_annotations{lead_idx,x2_idx,jj} = repmat(x_lims(2),1,length(all_annotations{lead_idx,y1_idx,jj}));
                    end

                    if(cant_realizations > 1)
                        if( isempty(all_annotations{lead_idx,y1_idx,jj}) )
                            all_annotations{lead_idx,y1_idx,jj} = min(squeeze(ECG(all_annotations{lead_idx,x1_idx,jj},lead_idx,:)),[],2);
                        end
                        if( isempty(all_annotations{lead_idx,y2_idx,jj}) )
                            all_annotations{lead_idx,y2_idx,jj} = max(squeeze(ECG(all_annotations{lead_idx,x1_idx,jj},lead_idx,:)),[],2);
                        end
                    else
                        if( isempty(all_annotations{lead_idx,y1_idx,jj}) )
                            all_annotations{lead_idx,y1_idx,jj} = repmat(min(ECG(:,lead_idx)),1,length(all_annotations{lead_idx,x1_idx,jj}));
                        end
                        if( isempty(all_annotations{lead_idx,y2_idx,jj}) )
                            all_annotations{lead_idx,y2_idx,jj} = repmat(max(ECG(:,lead_idx)),1,length(all_annotations{lead_idx,x1_idx,jj}));
                        end
                    end

                    if( ~isempty(cell2mat(all_annotations(lead_idx,4:7,jj)))   )

                        aux_hdl = text_line(cell2mat(all_annotations(lead_idx,[x1_idx x2_idx],jj)), cell2mat(all_annotations(lead_idx,[y1_idx y2_idx],jj)), '', all_annotations{lead_idx,plot_options_str_idx,jj}, axes_hdl(lead_idx) );

                        set(aux_hdl, 'Tag', num2str(lead_idx) );
                        
                        txthdl = [ txthdl aux_hdl ];

                        for kk = aux_hdl
                            uistack(kk,'bottom');
                            if(bAllowEdition)
                                set(kk,...
                                    'ButtonDownFcn', {@WindowButtonDown} ...
                                    );            
                            end
                        end                        
%                         aux_hdl = eval( [ 'plot( cell2mat(all_annotations(' num2str(lead_idx) ', [' num2str(x1_idx) ' ' num2str(x2_idx) '],' num2str(jj) ')), cell2mat(all_annotations(' num2str(lead_idx) ', [' num2str(y1_idx) ' ' num2str(y2_idx) '],' num2str(jj) '))'  all_annotations{lead_idx,plot_options_str_idx,jj} ')']);
%                         anns_hdl = [ anns_hdl; colvec(aux_hdl)];
%                         if( ~strcmpi('', all_annotations{lead_idx,text_idx,jj}) )
%                             txthdl = [ txthdl; text(all_annotations{lead_idx,x1_idx,jj} + 0.03*x_range, all_annotations{lead_idx,y1_idx,jj} + 0.03*y_range, all_annotations{lead_idx,text_idx,jj}, 'VerticalAlignment', 'bottom','BackgroundColor', bgTextColor )];
% %                             set(txthdl(end), 'Rotation', round(180/pi*atan2( (all_annotations{lead_idx,y2_idx,jj} - all_annotations{lead_idx,y1_idx,jj}) * axis_position(4) / y_range , (all_annotations{lead_idx,x2_idx,jj} - all_annotations{lead_idx,x1_idx,jj}) * axis_position(3) / x_range )) );
% %                             set(txthdl(end), 'Rotation', 90 );
%                         end                    
                    end

                case 'patch'

                    if( isempty(all_annotations{lead_idx,x1_idx,jj}) )
                        all_annotations{lead_idx,x1_idx,jj} = repmat(x_lims(1),1,length(all_annotations{lead_idx,y1_idx,jj}));
                    end
                    if( isempty(all_annotations{lead_idx,x2_idx,jj}) )
                        all_annotations{lead_idx,x2_idx,jj} = repmat(x_lims(2),1,length(all_annotations{lead_idx,y1_idx,jj}));
                    end

                    if(cant_realizations > 1)
                        if( isempty(all_annotations{lead_idx,y1_idx,jj}) )
                            all_annotations{lead_idx,y1_idx,jj} = min(squeeze(ECG(all_annotations{lead_idx,x1_idx,jj},lead_idx,:)),[],2);
                        end
                        if( isempty(all_annotations{lead_idx,y2_idx,jj}) )
                            all_annotations{lead_idx,y2_idx,jj} = max(squeeze(ECG(all_annotations{lead_idx,x1_idx,jj},lead_idx,:)),[],2);
                        end
                    else
                        if( isempty(all_annotations{lead_idx,y1_idx,jj}) )
                            all_annotations{lead_idx,y1_idx,jj} = repmat(min(ECG(:,lead_idx)),1,length(all_annotations{lead_idx,x1_idx,jj}));
                        end
                        if( isempty(all_annotations{lead_idx,y2_idx,jj}) )
                            all_annotations{lead_idx,y2_idx,jj} = repmat(max(ECG(:,lead_idx)),1,length(all_annotations{lead_idx,x1_idx,jj}));
                        end
                    end                    

                    if( ~isempty(cell2mat(all_annotations(lead_idx,4:7,jj)))   )

                        aux_hdl = patch( [all_annotations{lead_idx,x1_idx,jj} all_annotations{lead_idx,x1_idx,jj} all_annotations{lead_idx,x2_idx,jj} all_annotations{lead_idx,x2_idx,jj} ], [ all_annotations{lead_idx,y1_idx,jj} all_annotations{lead_idx,y2_idx,jj} all_annotations{lead_idx,y2_idx,jj} all_annotations{lead_idx,y1_idx,jj} ], [1 1 1], 'EdgeColor', 'none');
                        aux_prop_vals = all_annotations{lead_idx,plot_options_str_idx,jj};
                        set( aux_hdl, aux_prop_vals(:,1)', aux_prop_vals(:,2)' );

                        set(aux_hdl, 'Tag', num2str(lead_idx) );
                        
                        if(bAllowEdition)
                            set(aux_hdl,...
                                'ButtonDownFcn', {@WindowButtonDown} ...
                                );            
                            
                        end
                        
                        txthdl = [ txthdl aux_hdl ];

                        uistack(aux_hdl,'bottom');

                    end                    

            end

        end
        
        
    end
    
    function ScrollWheel(obj,event_obj)

        if( ~bAllowEdition )
            return
        end
        
        if (event_obj.VerticalScrollCount > 0)
            waveMode_idx = waveMode_idx + 1;
            if(waveMode_idx > lWaveAnnotations)
                waveMode_idx = 1;
            end
        else
            waveMode_idx = waveMode_idx - 1;
            if(waveMode_idx < 1)
                waveMode_idx = lWaveAnnotations;
            end
        end
        
        bFixedWave = true;
        
        update_title_efimero( cWaveAnnotations{waveMode_idx}, 10 );
        
    end

    function KeyPress(obj,event_obj)

        if(strcmp(event_obj.Key,'e'))
            bAllowEdition = ~bAllowEdition;
            if(bAllowEdition)
                update_title_efimero( 'Edit mode ON', 2 )
                set(ECG_lgd, 'visible', 'off')
            else
                update_title_efimero( 'Edit mode OFF', 2 )
                set(ECG_lgd, 'visible', 'on')
            end
            
        elseif (strcmp(event_obj.Key,'s'))
            
            if( bDelineationEdited )
                % store manual_edited_position
                assignin('caller', 'manual_edited_position', manual_edited_position)
                fprintf(1, '[%s] Edited positions stored in caller workspace.\n', datestr(now) );
                set( fig_hdl, 'Color', DelineationNotEditedColor )                
            else
                fprintf(2, '[%s] No changes made, edit delineation before saving.\n', datestr(now) );
            end
            
        end

    end

    function [min_vals_x, min_vals_y, max_vals_x, max_vals_y] = CalcMinMax( this_annotation, field_names )

        if( isfield(this_annotation, field_names{1} ) )
            aux_on = cell2mat({this_annotation(:).(field_names{1})});
        else
            aux_on = [];
        end

        if( isfield(this_annotation, field_names{2} ) )
            aux_peak = cell2mat({this_annotation(:).(field_names{2})});
        else
            aux_peak = [];
        end

        if( isfield(this_annotation, field_names{3} ) )
            aux_off = cell2mat({this_annotation(:).(field_names{3})});
        else
            aux_off = [];
        end

        start_sample = 1;

        % wave start
        bOn = ~isnan(aux_on) & aux_on >= start_sample & aux_on <= start_sample+lECG;
        % wave end
        bOff = ~isnan(aux_off) & aux_off >= start_sample & aux_off <= start_sample+lECG;

        % wave peak
        bPeak = ~isnan(aux_peak) & aux_peak >= start_sample & aux_peak <= start_sample+lECG;

        % wave conection between start-end
        bOnOff = bOn & bOff;
        aux_complete_idx = find(bOnOff);

        aux_complete_idxx = num2cell(nan(ECG_header.nsig,1));
        % normal sampled version
        aux_complete_idxx(aux_complete_idx) = arrayfun( @(a)( max(1, aux_on(a)):min(lECG,aux_off(a)) ),aux_complete_idx, 'UniformOutput', false);
        if( any(~bOnOff & bOn & bPeak) )
            % on-peak
            aux_complete_idx = find( ~bOnOff & bOn & bPeak);
            aux_complete_idxx(aux_complete_idx) = arrayfun( @(a)( aux_on(a):aux_peak(a) ),aux_complete_idx, 'UniformOutput', false);
        end

        if( any(~bOnOff & bOff & bPeak) )
            % peak-off
            aux_complete_idx = find(~bOnOff & bPeak & bOff);
            aux_complete_idxx(aux_complete_idx) = arrayfun( @(a)( aux_peak(a):aux_off(a) ),aux_complete_idx, 'UniformOutput', false);
        end         

        max_vals_x = cellfun( @(a,lead)( a(end) ), aux_complete_idxx);
        max_vals_y = cellfun( @(a,lead)( max(dummyfun(ECG,a, lead))), aux_complete_idxx, num2cell((1:ECG_header.nsig)'));
        min_vals_x = cellfun( @(a,lead)( a(1) ), aux_complete_idxx);
        min_vals_y = cellfun( @(a,lead)( min(dummyfun(ECG,a, lead))), aux_complete_idxx, num2cell((1:ECG_header.nsig)'));

    end

    function fun_val = dummyfun(a,b,c)
        
        fun_val = nan;
        if( ~any(isnan(b)) && ~any(isnan(c)) )
            fun_val = a(b,c);
        end
        
    end

    function [lead_idx, obj ] = get_lead_idx(obj)
       
        lead_idx = find(obj == axes_hdl);
        
        if( isempty(lead_idx) )
            %click on the signal
            obj = get(obj, 'Parent');            
            lead_idx = find(obj == axes_hdl);
        end
        
    end

    function WindowButtonMotion(obj,event_obj)
        
        if (bIsDragAllowed)
        
            point = get(axes_move,'CurrentPoint');

            sample_move = round(point(1));

            sample_move = min( lECG, max(1, round(sample_move) ) );

            if( ~isempty(move_hdl) )
                delete(move_hdl)
            end

            move_hdl = plot(axes_move, sample_move, ECG(sample_move, lead_move_idx), '+r', 'MarkerSize', 8 );

        else
            set(fig_hdl, 'WindowButtonMotionFcn', PrevStateWindowButtonMotionFcn );
        end
    end

    function WindowButtonUp(obj,event_obj)
        
        if (bIsDragAllowed)

            bIsDragAllowed = false;
            bPreserveFix = false;
            
            set(fig_hdl, 'WindowButtonMotionFcn', PrevStateWindowButtonMotionFcn );
            
            hold(axes_move, 'off' )
            
            if( ~isempty(move_hdl) )
                delete(move_hdl)
                move_hdl = [];
            end
            
            manual_edited_position(lead_move_idx).((cWaveAnnotations{waveMode_idx})) = sample_move;
            
            aux_mat = manual_edited_position_matrix{lead_move_idx};
            aux_mat(waveMode_idx) = sample_move;
            manual_edited_position_matrix{lead_move_idx} = aux_mat;
            
            if( isempty(annotations) )
                all_annotations = update_annotations_from_positions(manual_edited_position);
            else
                all_annotations = cat(3, annotations, update_annotations_from_positions(manual_edited_position) );
            end
            
            delete(findobj('Tag', num2str(lead_move_idx) ))
            
            update_annotations( lead_move_idx )
            
            lead_move_idx = [];
            bDelineationEdited = true;
            set( fig_hdl, 'Color', DelineationEditedColor )                
            
        end
        
    end

    function WindowButtonDown(obj,event_obj)

        if( ~bAllowEdition )
            return
        end
        
        [lead_idx, obj] = get_lead_idx(obj);
        
        if (strcmp(get(fig_hdl,'SelectionType'),'alt'))
            % Delete annotation
            
            point = get(obj,'CurrentPoint');

            sample = round(point(1));

            if( ~bFixedWave )
                [~, waveMode_idx] = sort(abs(manual_edited_position_matrix{lead_idx} - sample));
                waveMode_idx = waveMode_idx(1);
            end
            
            manual_edited_position(lead_idx).((cWaveAnnotations{waveMode_idx})) = nan;
            
            aux_mat = manual_edited_position_matrix{lead_idx};
            aux_mat(waveMode_idx) = nan;
            manual_edited_position_matrix{lead_idx} = aux_mat;
            
            if( isempty(annotations) )
                all_annotations = update_annotations_from_positions(manual_edited_position);
            else
                all_annotations = cat(3, annotations, update_annotations_from_positions(manual_edited_position) );
            end
            
            delete(findobj('Tag', num2str(lead_idx) ))
            
            update_annotations( lead_idx )            

        elseif( strcmp(get(fig_hdl,'SelectionType'),'extend') )

            point = get(obj,'CurrentPoint');

            sample = round(point(1));

            if( bFixedWave )
                bPreserveFix = true;
            else
                [~, waveMode_idx] = sort(abs(manual_edited_position_matrix{lead_idx} - sample));
                
                if( sample < manual_edited_position(lead_idx).((cWaveAnnotations{waveMode_idx(1)})) )
                    waveMode_idx = max(1, waveMode_idx(1) - 1);
                elseif( sample == manual_edited_position(lead_idx).((cWaveAnnotations{waveMode_idx(1)})) )
                    waveMode_idx = waveMode_idx(1);
                else
                    waveMode_idx = min( lWaveAnnotations, waveMode_idx(1) + 1);
                end
                
            end            

            manual_edited_position(lead_idx).((cWaveAnnotations{waveMode_idx})) = sample;
            
            aux_mat = manual_edited_position_matrix{lead_idx};
            aux_mat(waveMode_idx) = sample;
            manual_edited_position_matrix{lead_idx} = aux_mat;
            
            if( isempty(annotations) )
                all_annotations = update_annotations_from_positions(manual_edited_position);
            else
                all_annotations = cat(3, annotations, update_annotations_from_positions(manual_edited_position) );
            end
            
            delete(findobj('Tag', num2str(lead_idx) ))
            
            update_annotations( lead_idx )
            
            bDelineationEdited = true;
            set( fig_hdl, 'Color', DelineationEditedColor )     
            
            
        elseif (strcmp(get(fig_hdl,'SelectionType'),'normal'))

            if (~bIsDragAllowed)

                axes_move = obj;
                lead_move_idx = lead_idx;
                bIsDragAllowed = true;

                point = get(axes_move,'CurrentPoint');

                sample = round(point(1));

                if( bFixedWave )
                    bPreserveFix = true;
                else
                    [~, waveMode_idx] = sort(abs(manual_edited_position_matrix{lead_move_idx} - sample));
                    waveMode_idx = waveMode_idx(1);
                end            
                
                sample_move = sample;
                
                PrevStateWindowButtonMotionFcn = get(fig_hdl, 'WindowButtonMotionFcn');
                set(fig_hdl, 'WindowButtonMotionFcn', @WindowButtonMotion);

                hold(axes_move, 'on' )
                
            end
                
        end

    end

    function update_title_efimero( strTitle, delay )
       
        delete(findobj('Tag', 'title_efimero' ))
        
        for hdl_idx = rowvec(axes_hdl) 
        
            set(fig_hdl, 'CurrentAxes', hdl_idx);

            y_lims = get(hdl_idx, 'ylim');
            x_lims = get(hdl_idx, 'xlim');

            plotXmin = x_lims(1);
            plotXmax = x_lims(2);
            plotXrange = plotXmax - plotXmin;
            plotYmin = y_lims(1);
            plotYmax = y_lims(2);
            plotYrange = plotYmax - plotYmin;
            yTextOffset = 0.01*plotYrange;

            width_legend = 0.2*plotXrange;
            height_legend = 0.1*plotYrange;
            left_legend = plotXmax - width_legend - 5*xTextOffset;
            bottom_legend = plotYmin;

            aux_hdl = text( left_legend + width_legend/2, bottom_legend + height_legend/2, strTitle, 'FontSize', 8, 'HorizontalAlignment', 'center', 'BackgroundColor', 'r' );

            set(aux_hdl, 'Tag', 'title_efimero')
            
        end
        
        if( ~isinf(delay) && strcmpi(my_timer.Running, 'off') )
            my_timer.StartDelay = delay;
            start(my_timer)
        end
        
    end


end