#!/bin/csh -f
# $Id: run,v 2.0 94/07/17 18:33:45 dfk Exp $

set File = $1
set Option = $2
set TICSperMSEC = 25000.
#set verbose
#set echo

onintr cleanup

zcat $File | srt-extract $Option | sort +1n +2n -T . \
    | analyze temp > /dev/null
dm "x1/$TICSperMSEC" < temp | cdfhere > cdf_real
dm "x2/$TICSperMSEC" < temp | cdfhere > cdf_model

dm x1 < cdf_real > sort_real
dm x1 < cdf_model > sort_model

/u/dfk/bin/plot -lines -t "$File cdfs" cdf_real cdf_model
/u/dfk/bin/plot -lines -x 0:50 -t "$File cdfs" cdf_real cdf_model

abut sort_real sort_model | dm "(x2 - x1)^2" | stats mean | dm "sqrt(x1)" > demerit
stats mean < sort_real > mean
echo -n mean:
cat mean
echo -n demerit figure:
cat demerit
echo -n demerit percentage:
abut demerit mean | dm "if x2 = 0 then 'error' else (x1/x2 * 100.)"

abut sort_real sort_model | running-demerit > run_demerit

/u/dfk/bin/plot -lines -t "$File running demerit percentage" run_demerit

set name=$File:t
set name=$name:r
set name=$name:r

# rm -f temp
mv demerit     Data/$name.demerit
mv mean        Data/$name.mean
mv cdf_model   Data/$name.model.cdf
mv cdf_real    Data/$name.real.cdf
mv run_demerit Data/$name.run_demerit

exit 0

cleanup:
rm -f temp demerit mean cdf_model cdf_real 
rm -f sort_model sort_real run_demerit
exit 1
