function [FQRSout,QTout] = physionet2013(tm,ecg,paramStruct)
% template algorithm for Physionet/CinC competition 2013. This function can
% be used for events 1 and 2. Participants are free to modify any
% components of the code. However the function prototype must stay the same
% [FQRSout,QTout] = physionet2013(tm,ecg) where the inputs and outputs are specified
% below.
%
% inputs
%   tm :         Nx1 vector of time in milliseconds
%   ecg:         4x60000 (4 channels and 1min of signal at 1000Hz) matrix of
%                abdominal ecg channels.
%   paramStruct: param of the algorithm to use. if not inputed then the
%                default parameters are used. 
%                ***FOR SET-C of the challenge LEAVE EMPTY.***
% 
% output
%   FQRSout:    FQRS markers in seconds. Each marker indicates the position of one
%               of the FQRS detected by the algorithm.
%   QTout:      1x1 estimated fetal QTout duration
%
%
%
% Safe Foetus Monitoring Toolbox, version 1.0, Sept 2013
% Released under the GNU General Public License
%
% Copyright (C) 2013 Joachim Behar
% Oxford university, Intelligent Patient Monitoring Group - Oxford 2013
% joachim.behar@eng.ox.ac.uk
%
% Last updated : 02-08-2013
%
% 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 2 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.    

if nargin<2; error('physionet2013: wrong number of input arguments'); end;
if nargin<3; paramStruct=[]; end;
addpath(genpath('subfunctions/'));

try
    if isempty(paramStruct)
        % == for FQRS extraction
        paramStruct.method = 'FUSE';              % TS/TS-PCA/TS-EKF/DEFLATION/FUSE
        paramStruct.methodComplexity = 'COMPLEX'; % COMPLEX/SIMPLE
        paramStruct.smoothRR = 1;                 % postprocessing
        paramStruct.cinc_match = 1;               % twick to better match CinC scores
        paramStruct.debug = 0;
        paramStruct.fs = 1000;
        paramStruct.TSnbCycles = 20;              % nb of cycles to build template 
        paramStruct.EKFgain = 1;
        paramStruct.fhighNbCoeff = 5;
        paramStruct.fbasNbCoeff = 3;
        paramStruct.NotchCheck = 1;
        paramStruct.fbas = 10;
        paramStruct.fhigh = 99;
        paramStruct.STDcut = 29;                  % rubbish to fit to challenge scores (was 17)
        paramStruct.NbPC = 2;                     % nb of principal components in PCA
        paramStruct.detThres = 0.5;               % QRS detector threshold for FQRS
        paramStruct.qrsdet_wind = 15;             % qrs detection window: LENGTH_SEG/qrsdet_wind NEEDS TO BE AN INTEGER!!
        paramStruct.extract_qt = 0;
    end
    
    [FQRSout,~,outExtra] = physionet2013extract(tm,ecg,paramStruct);
catch ME
    for enb=1:length(ME.stack); disp(ME.stack(enb)); end;
    FQRSout=0.1:0.42:60;
    QTout=250;
end


try
    % == for QT extraction
    % the main reasons for doing the QT separately are:
    % (1) the baseline wander cutt-off must be lower than 10Hz otherwise
    % the T-wave is completely destroyed
    % (2) using BSS methods is not a good idea as part of the T-wave could
    % end in a secondary channel so we stay we time based methods.
    qt_paramStruct.method = 'TS-PCA';
    qt_paramStruct.methodComplexity = 'SIMPLE';  % MUST BE: SIMPLE (time based)
    qt_paramStruct.debug = 0; 
    qt_paramStruct.fs = 1000;
    qt_paramStruct.TSnbCycles = 20;
    qt_paramStruct.EKFgain = 1;
    qt_paramStruct.fhighNbCoeff = 5;
    qt_paramStruct.fbasNbCoeff = 3;
    qt_paramStruct.NotchCheck = 1;
    qt_paramStruct.fbas = 2; % a lower cut-off not to kill the T-wave
    qt_paramStruct.fhigh = 99;
    qt_paramStruct.NbPC = 2;
    qt_paramStruct.extract_qt = 1;    
    qt_paramStruct.rMQRS = outExtra.rMQRS; % use the MQRS returned by the previous step
    qt_paramStruct.rFQRS = outExtra.rFQRS; % use the FQRS returned by the previous step
    
    [~,QTout,~] = physionet2013extract(tm,ecg,qt_paramStruct);

catch ME
    for enb=1:length(ME.stack); disp(ME.stack(enb)); end;
    FQRSout=0.1:0.42:60;
    QTout=250;
end

end