function handles = update_plots(handles)
%
% Author(s):    Jonathan Moeyersons       (Jonathan.Moeyersons@esat.kuleuven.be)
%               Sabine Van Huffel         (Sabine.Vanhuffel@esat.kuleuven.be)
%               Carolina Varon            (Carolina.Varon@esat.kuleuven.be)
%
% Version History:
% - 06/05/2019   JM      Initial version
%
% Copyright (c) 2019,  Jonathan Moeyersons, KULeuven-ESAT-STADIUS 
%
% This software is made available for non commercial research purposes only 
% under the GNU General Public License. However, notwithstanding any 
% provision of the GNU General Public License, this software may not be 
% used for commercial purposes without explicit written permission after 
% contacting jonathan.moeyersons@esat.kuleuven.be
%
% This program is free software: you can redistribute it and/or modify
% it under the terms of the GNU General Public License as published by
% the Free Software Foundation, either version 3 of the License, or
% (at your option) any later version.
% 
% This program is distributed in the hope that it will be useful,
% but WITHOUT ANY WARRANTY; without even the implied warranty of
% MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
% GNU General Public License for more details.
% 
% You should have received a copy of the GNU General Public License
% along with this program.  If not, see <https://www.gnu.org/licenses/>.
% Get the listbox value
listbox_value = get(handles.listbox_channels,'value')-1; % 0=All, 1=Channel 1,...

% Define the looping variable
if listbox_value == 0
    loop_var = 1:handles.data.channels;
else % Only the selection
    loop_var = listbox_value;
end

% Define sampling frequency
fs = handles.data.fs;

% Get the start time
start = max([1 round(fs*seconds(handles.data.start_analysis))]);

% Get the stop time
stop = min([size(handles.data.signal.filtered,1) round(fs*seconds(handles.data.start_analysis+handles.data.duration_analysis))]);

% Get the time
time = seconds((1:size(handles.data.signal.filtered,1))/fs); 
time.Format = 'hh:mm:ss';

% Clear the axes
cla(handles.axes_ECG)
cla(handles.axes_tachogram)

% Adjust axis
for ii = 1:handles.data.channels
    % Loaded signal
    handles.graph.signal.original.(handles.labels{ii}) = plot(handles.axes_ECG,time(start:stop),handles.data.signal.original(start:stop,ii),...
        'color',[handles.color(ii,:),0.2],...
        'visible','off');

    % Filtered signal
    handles.graph.signal.filtered.(handles.labels{ii}) = plot(handles.axes_ECG,time(start:stop),handles.data.signal.filtered(start:stop,ii),...
        'color',handles.color(ii,:),...
        'visible','off',...
        'tag',handles.labels{ii},...
        'UIContextMenu', handles.context_lead_options);

    if ~isempty(handles.data.R)
        % Define the R-peaks that remain
        temp = handles.data.R{ii} > time(start);
        temp = temp + (handles.data.R{ii} < time(stop));
        R_peak = handles.data.R{ii}(temp == 2);
        
        if ~isempty(R_peak)
            % Define the RR-intervals
            RR_int = seconds(diff(R_peak))*1000;
            
            % Adjust the tachogram based on the units selection
            contents = get(handles.popupmenu_units_tachogram,'string');

            switch contents{get(handles.popupmenu_units_tachogram,'Value')}
                case 'HR (bpm)'
                    % Change the data
                    RR_int = 60000./RR_int;
            end

            % Plot the R-peaks on the ECG-graph
            handles.graph.R.(handles.labels{ii}) = plot(handles.axes_ECG,...
                R_peak,...
                handles.data.signal.filtered(round(seconds(R_peak)*fs),ii),...
                'o',...
                'color',handles.color(ii,:),...
                'markerfacecolor',handles.color(ii,:),...
                'visible','off',...
                'tag',handles.labels{ii},...
                'UIContextMenu', handles.context_R_peak_options);

            % Create the tachogram
            handles.graph.RR.(handles.labels{ii}) = plot(handles.axes_tachogram,...
                R_peak(2:end),...
                RR_int,...
                '-o',...
                'color',handles.color(ii,:),...
                'markerfacecolor',handles.color(ii,:),...
                'visible','off',...
                'tag',handles.labels{ii});
        end
    end
end

% Make the signals visible
for ii = loop_var
    set(handles.graph.signal.filtered.(handles.labels{ii}),'visible','on')

    if handles.checkbox_show_original.Value
        set(handles.graph.signal.original.(handles.labels{ii}),'visible','on')
    end

    if ~isempty(handles.data.R)
        set(handles.graph.R.(handles.labels{ii}),'visible','on')
        set(handles.graph.RR.(handles.labels{ii}),'visible','on')
    end
end