#!/usr/bin/csh -f

# creates the distance-plot of the two input files
# for all three combos (DE1, DE2, CDE)
# in LOG-LOG scales

if ( $#argv != 2 ) then
    echo "USAGE: $0 fname1 fname2 "
    exit
else
  while ( $1 =~ -* )
      echo "flag $1"
      shift
  end
  set fname1 = $1
  set fname2 = $2
  echo $fname1 $fname2 
endif


# set NAME = ${fname1}-${fname2}-out

set NAME = ${fname1}-${fname2}.pcp-ll
set FINAL = ${fname1}-${fname2}.pcp-ll
set FINAL1 = ${fname1}-${fname1}.pcp-ll
set FINAL2 = ${fname2}-${fname2}.pcp-ll

# set FINAL = ${fname1}-${fname2}.pcp-llnf
# set FINAL1 = ${fname1}-${fname1}.pcp-llnf
# set FINAL2 = ${fname2}-${fname2}.pcp-llnf

# set FINAL = ${fname1}-${fname2}-out.final
# set FINAL1 = ${fname1}-${fname1}-out.final
# set FINAL2 = ${fname2}-${fname2}-out.final
\rm -f $FINAL $FINAL1 $FINAL2

# npairs.pl $1 $2 | sort -n | awk '{print $0, NR }'  |\
#            flip12.pl | choptoplot | flip12.pl |\
#               awk -f dezero.awk | loglog.pl |\
#       killFlats.pl > $FINAL

perl5 npairs.pl $1 $2 | perl5 lbucketize.pl |\
    perl5 dezero.pl | perl5 loglog.pl > $FINAL
    # perl5 killFlats.pl 

perl5 npairs.pl $1 $1 | perl5 lbucketize.pl |\
    perl5 dezero.pl | perl5 loglog.pl > $FINAL1
    # perl5 killFlats.pl 

perl5 npairs.pl $2 $2 | perl5 lbucketize.pl |\
    perl5 dezero.pl | perl5 loglog.pl > $FINAL2
    # perl5 killFlats.pl 


echo " *** creating the plot "
# put everything into a eps plot
newmygnuplot -l "PCP self- and cross- plots" \
   -slope -o ${NAME}-all.eps -xlabel "log(dist)" \
   -ylabel "log(count-of-pairs)" \
   -s points \
   $FINAL \
   $FINAL1 $FINAL2

gv ${NAME}-all.eps
