hydrophone-XY:plot.sh

Go to the documentation of this file.

#!/bin/zsh
script=${0##*/}
 
if [ -z $JPP_DIR ]; then
    echo "Variable JPP_DIR undefined."
    exit
fi
 
source $JPP_DIR/setenv.sh $JPP_DIR >& /dev/null
 
zmodload zsh/mathfunc
 
set_variable:  DEBUG       ACOUSTICS_DEBUG        2
set_variable:  WORKDIR     ACOUSTICS_WORKDIR      ./
set_variable:  FORMAT      GRAPHICS_FORMAT        gif
set_variable+  BATCH       GRAPHICS_BATCH         -B
 
if do_usage $*; then
    usage "$script <string identifier>"
fi
 
if (( $# != 1 )); then
    fatal "Wrong number of arguments."
fi
 
set_variable  STRING              $argv[1]
 
set_variable  HYDROPHONE_TXT      $WORKDIR/hydrophonexy_${STRING}.txt
set_variable  HYDROPHONE_ROOT     $WORKDIR/hydrophonexy_${STRING}.root
 
if [[ ! -f $HYDROPHONE_TXT ]] then
    fatal "No hydrophone data file $HYDROPHONE_TXT."
fi
 
JGraph2D \
    -f $HYDROPHONE_TXT                        \
    -o $HYDROPHONE_ROOT
 
JPrintHydrophone -f $WORKDIR/hydrophone.txt -S $STRING | read DUMMY FLOOR X Y Z
 
JEllipse -o ${TMPDIR:-/tmp}/ellipse.root -p "0.0 0.0 $((sqrt($X*$X + $Y*$Y)))"
 
echo $X $Y > ${TMPDIR:-/tmp}/graph.txt
 
JGraph \
    -f ${TMPDIR:-/tmp}/graph.txt                          \
    -o ${TMPDIR:-/tmp}/graph.root
 
typeset -Z 4 STRING
 
JPlot2D \
    -f ${HYDROPHONE_ROOT}:\.\*                    \
    -f ${TMPDIR:-/tmp}/ellipse.root:\.\*                     \
    -f ${TMPDIR:-/tmp}/graph.root:\.\*                       \
    -> "#Deltax [m]"                             \
    -< "#Deltay [m]"                             \
    -\^ "RMS [ms]"                                \
    -N "X 505"                                    \
    -N "Y 505"                                    \
    -Z                                            \
    -O COLZ                                       \
    -T "[$STRING]"                                \
    -o $WORKDIR/hydrophonexy_${STRING}.$FORMAT    $BATCH
 
rm -f  $HYDROPHONE_ROOT
rm -f ${TMPDIR:-/tmp}/ellipse.root
rm -f ${TMPDIR:-/tmp}/graph.{txt,root}