% -------------------------------------------------------------------------------------------------
% plotSgnMrkN.M: Plot signals (on multiple horizontal strips) with vertical markers.
%                --> With signal normalization
%                sgn    = signal matrix
%                pqrs   = vector of time indexes (markers)
%                freq   = sampling frequency
%                name   = figure title
%                nsecr  = number of seconds per row
%                nbands = number of strips per page
%
% Copyright (C) 2004 Maurizio Varanini, Clinical Physiology Institute, CNR, Pisa, Italy
%
% 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 "as is" and "as available" in the hope that it will be useful,
% but WITHOUT ANY WARRANTY of any kind; 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, write to the Free Software Foundation, Inc.,
% 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA.
%
% For any comment or bug report, please send e-mail to: maurizio.varanini@ifc.cnr.it
% -------------------------------------------------------------------------------------------------

function PlotSgnMrkN(sgn,pqrs,freq,name, nband, nsecr)
% convert column to row if necessary
[m,n] = size(sgn);  if n<m,   sgn = sgn.';   end
[m,n] = size(pqrs); if n<m,  pqrs = pqrs.';   end
yminc=min(sgn,[],2);
ymaxc=max(sgn,[],2);
nsgn=size(sgn,1);
nsampt=size(sgn,2);
if(nargin<3) freq=1; end
nsect=nsampt/freq;
if(nargin<4) name=''; end
if(nargin<5) nband=6; end
% if(nargin<6) nsecr=10; end
if(nargin<6) nsecr=nsect/nband; end
nsampr=nsecr*freq;   % number of samples per row
isampi=1;
isampf=0;
vmrkcol=['m', 'b', 'r', 'g', 'k'];

while isampf < nsampt
    figure
    for i=0:nband-1
        subplot(nband,1,i+1);
        isampf = isampi+nsampr-1;
        if(isampf > nsampt) isampf=nsampt; end;
%        yminc=min(sgn(:,isampi:isampf),[],2);
%        ymaxc=max(sgn(:,isampi:isampf),[],2);
        ymin=-1; ymax=1;
        yn=(ymax-ymin)./(ymaxc-yminc);
        yoffset=(1*(ymax-ymin)).*((nsgn:-1:1)-1)';
        clear ysgn;
        for is=1:nsgn
            ysgn(is,:)= yoffset(is) +ymin + yn(is)*(sgn(is,isampi:isampf)-yminc(is));
        end
        yming=ymin;
        ymaxg=ymax+max(yoffset);
        xmin=isampi/freq;
        xmax=(isampi+nsampr)/freq;
        plot([isampi:isampf]/freq ,ysgn);
        axis tight;
        if(i==0) ht=title(name); set(ht,'Interpreter','none'); end;
        XLIM([xmin, xmax]);
        YLIM([yming, ymaxg]);
        if(~isempty(pqrs))
            hold on;
            for im=1:size(pqrs,1);  % loop on marker type
                vmrk=pqrs(im,:);
                mrkcol=vmrkcol(mod(im,size(vmrkcol,2)));
                iqrsint=find((isampi < vmrk)&(isampf > vmrk));
                pqrsint=vmrk(iqrsint); 
                fxpqrs=[pqrsint; pqrsint]/freq;
                fypqrs = ones(size(fxpqrs)); fypqrs(1,:)=yming; fypqrs(2,:)=ymaxg;
                plot(fxpqrs,fypqrs, [mrkcol ':+']);
            end
        end
        if(isampf == nsampt) break; end;
        isampi = isampf+1;
    end
end