JSeaBird.sh

Go to the documentation of this file.

#!/usr/bin/env zsh
script=${0##*/}
 
if [ -z $JPP_DIR ]; then
    echo "Variable JPP_DIR undefined."
    exit
fi
 
source $JPP_DIR/setenv.sh $JPP_DIR #>& /dev/null
 
set_variable   DEBUG           3
set_variable   DIR             $JPP_DIR/examples/JAcoustics/
set_variable:  FORMAT          GRAPHICS_FORMAT        gif
set_variable+  BATCH           GRAPHICS_BATCH         -B
 
set_variable   SOUND_VELOCITY  1541.00 -17e-3 -2000.00
 
let "v0 = 1541.00"
let "b  = -17e-3"
let "z0 = -2000.00"
let "za = -2500.00"
let "zb = -3450.00"
 
if do_usage $*; then
    usage "$script <input file>"\
          "\nSee https://git.km3net.de/calibration/input_tables/data/sound/ for example files."
fi
 
case $# in
    1)  set_variable INPUT_FILE      $argv[1];;
    *)  fatal "Wrong number of arguments.";;
esac
 
 
$DIR/JSeaBird \
    -f $INPUT_FILE                       \
    -o seabird.root                      \
    -V "$SOUND_VELOCITY"                 \
    -d $DEBUG --!
 
$DIR/JSeaBird \
    -f $INPUT_FILE                       \
    -o seabird+.root                     \
    -V "$SOUND_VELOCITY"                 \
    -O                                   \
    -d $DEBUG --!
 
if (( 1 )); then
 
    JPlot1D \
        -f seabird.root:hz                   \
        -f seabird+.root:hz                  \
        -> "#Deltaz [m]"                    \
        -\^ "number of entries"              \
        -T ""                                \
        -o z.$FORMAT                         $BATCH
fi
if (( 1 )); then
 
    JPlot1D \
        -f seabird.root:v1                   \
        -f seabird.root:v2                   \
        -> "#Deltav [m/s]"                  \
        -\^ "number of entries"              \
        -T ""                                \
        -L TL                                \
        -o v.$FORMAT                         $BATCH
fi
if (( 1 )); then
 
    JLine \
        -o line.root                         \
        -@ "color = 2"                       \
        -p "$((-1 * $za)) $(($v0 + $b*($za - $z0))) $((-1 * $zb)) $(($v0 + $b*($zb - $z0)))" 
 
    JPlot1D \
        -@ "titleOffset  = 1.2"              \
        -w 600x600                           \
        -x "2500 3450"                       \
        -y "1540 1570"                       \
        -f seabird.root:ga                   \
        -f line.root:\.\*                    \
        -> "depth [m]"                      \
        -\^ "v [m/s]"                        \
        -N "X 505"                           \
        -N "Y 505"                           \
        -S 0.6                               \
        -T ""                                \
        -o ga.$FORMAT                        $BATCH
fi
if (( 1 )); then
 
    JPlot1D \
        -@ "titleOffset  = 0.6"              \
        -@ "bottomMargin = 0.2"              \
        -@ "fontSize     = 4.5"              \
        -@ "tickLength   = 4.5"              \
        -w 600x300                           \
        -f seabird.root:gc                   \
        -f seabird.root:gd                   \
        -y "-1.0 +1.0"                       \
        -> "depth [m]"                      \
        -\^ "#Deltav [m/s]"                  \
        -N "X 505"                           \
        -N "Y 505"                           \
        -S 0.4                               \
        -T ""                                \
        -o gc.$FORMAT                        $BATCH
fi
if (( 1 )); then
 
    JFit \
        -f seabird.root:ga                   \
        -o fit.root                          \
        -F "[0] + [1]*(x - 2000.0)"          \
        -x "2500 3450"                       \
        -@ "p0 = 1540.0"                     \
        -@ "p1 = 17e-3"                      
 
    JPlot1D \
        -@ "titleOffset  = 1.5"              \
        -f fit.root:\.\*                     \
        -x "2500 3450"                       \
        -y "1540 1570"                       \
        -> "depth [m]"                      \
        -\^ "v [m/s]"                        \
        -N "X 505"                           \
        -N "Y 505"                           \
        -S 0.6                               \
        -T ""                                \
        -o fit.$FORMAT                       $BATCH
fi