#! /bin/sh

F=0
T=24
I=4
P=1
PLT_OPT=""
STT_OPT=""
STPLOT_TYPE="uncorrected"

while [ "$1" != "" ]
do
OPT="$1"
case $OPT in
    -c) STT_OPT="-c"; STPLOT_TYPE="baseline-corrected" ;;
    -r) shift; R=$1 ;;
    -f) shift; F=$1 ;;
    -t) shift; T=$1 ;;
    -i) shift; I=$1 ;;
    -p) shift; P=$1 ;;
    -T) shift; PLT_OPT="-T $1" ;;
esac
shift
done

if [ "x$R" = "x" ]
then
	echo $*
	echo R = x"$R"x
	echo usage: stplot -r RECORD [ -f FROM -t TO -i TICK-INTERVAL -T lw ]
	exit
fi

sttrend -r $R -a 16a -f $F:0:0 -t $T:0:0 $STT_OPT >stt.$R
if [ \! -s ste.$R ]
then
    stevents -r $R | sort -n >ste.$R
fi
plt $PLT_OPT stt.$R 0 2 -W .1 .75 .95 .95 -g "grid sub none" \
    -x"" -y"ST0 (uV, $STPLOT_TYPE)" \
    -t"Record $R                                                                                                                                                - $P -" \
    -xa $F $T 1 " " 1 0 -ya -500 500 100 - 2 $F
plt $PLT_OPT ste.$R -pl 0 1 0 2 -W .1 .625 .95 .725 -seym \
    -x"" -y"Events" -t"" \
    -xa $F $T 1 " " 1 0 -ya -50 50 5 " " 3 $F
plt $PLT_OPT stt.$R 0 3 -W .1 .4 .95 .6 -g "grid sub none" -se \
    -x"" -y"ST1 (uV, $STPLOT_TYPE)" -t"" \
    -xa $F $T 1 " " 1 0 -ya -500 500 100 - 2 $F
plt $PLT_OPT stt.$R 0 1 -W .1 .05 .95 .3 -g "grid sub none" -se \
    -x"Elapsed time (hours)" -y"HR (bpm)" -t"" \
    -xa $F $T 1 - $I 0 -ya 0 200 50 - 2 $F
rm stt.$R
