/* rrlist.c		Joe Mietus		Apr 1 2011 */

/*

rrlist.c		Joe Mietus		Dec 19 2001
_______________________________________________________________________________
rrlist: List R-R intervals
Copyright (C) 2002 Joe Mietus

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.

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.

You may contact the author by e-mail (joe at physionet dot org) or postal mail
(MIT Room E25-505A, Cambridge, MA 02139 USA).  For updates to this software,
please visit PhysioNet (http://www.physionet.org/).
_______________________________________________________________________________

List R-R intervals

Usage: rrlist annotator tape [options]

options are :

    [-f start] : begin at time 'start'
    [-t end] : end at time 'end'
    [-l length] : output for duration 'length'
    [-h] : output time in hours in first column
    [-m] : output time in minutes in first column
    [-s] : output time in sec in first column
    [-a annotation] : list only intervals between consecutive annotations

*/

#include <stdio.h>
#include <wfdb/wfdb.h>
#include <wfdb/ecgmap.h>


main(argc, argv)	
int argc;
char *argv[];
{
    int i, annotation;
    int hflag, mflag, sflag, aflag;
    long start, end, length;
    double sps;
    struct WFDB_anninfo ai[1];
    struct WFDB_ann annot[2];

    if (argc < 3) {
        usage(argv[0]);
	exit(1);
    }

    hflag = mflag = sflag = aflag = 0;
    start = end = length = 0L;

    sps = sampfreq(argv[2]);
    ai[0].name = argv[1];
    ai[0].stat = WFDB_READ;
    if (annopen(argv[2], ai, 1) < 0)
       exit(2);

    for (i=2; ++i < argc && *argv[i] == '-'; ) {
        switch(argv[i][1]) {
	    case 'f': start = strtim(argv[++i]);
	              break;
	    case 't': end = strtim(argv[++i]);
	              break;
	    case 'l': length = strtim(argv[++i]);
	              break;
	    case 'h': hflag = 1;
	              break;
	    case 'm': mflag = 1;
	              break;
	    case 's': sflag = 1;
	              break;
	    case 'a': annotation = strann(argv[++i]);
	              aflag = 1;
	              break;
	    default:  usage(argv[0]);
	              exit(1);
        }
    }
    if (end == 0L && length != 0L)
        end = start + length;

    if (iannsettime(start) < 0)
        exit(2);

    start = 0L;

    while (getann(0, &annot[0]) >= 0) {
        if (isqrs(annot[0].anntyp))
        if ((!aflag && isqrs(annot[0].anntyp)) ||
            (aflag && annot[0].anntyp == annotation))
            break;
    }

    while (getann(0, &annot[1]) >= 0 && 
           (annot[1].time <= end || end == 0L)) {
	if ((!aflag && isqrs(annot[1].anntyp)) ||
            (aflag && annot[0].anntyp == annotation && 
             annot[1].anntyp == annotation)) {

	    if (hflag)
		printf("%.8f ", (annot[1].time-start)/(sps*3600));
	    else if (mflag)
		printf("%.6f ", (annot[1].time-start)/(sps*60));
	    else if (sflag)
		printf("%.3f ", (annot[1].time-start)/sps);
	    printf("%.3f %s\n", (annot[1].time - annot[0].time)/sps,
		                annstr(annot[1].anntyp));
	}
        if (isqrs(annot[1].anntyp)) {
	    annot[0] = annot[1];
	}
    }

    exit(0);
}


usage(prog)
char *prog;
{
    fprintf(stderr, "Usage : %s annotator tape [options]\n\n", prog);
    fprintf(stderr, " options are :\n");
    fprintf(stderr, " [-f start] : begin at time 'start'\n");
    fprintf(stderr, " [-t end] : end at time 'end'\n");
    fprintf(stderr, " [-l length] : output for duration 'length'\n");
    fprintf(stderr, " [-h] : output time in hours in first column\n");
    fprintf(stderr, " [-m] : output time in minutes in first column\n");
    fprintf(stderr, " [-s] : output time in seconds in first column\n");
    fprintf(stderr, " [-a annotation] : ");
    fprintf(stderr, " list only intervals between consecutive annotations\n");
}