hydrophone-phi:fit.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

set_variable  DEBUG               ${ACOUSTICS_DEBUG:-2}
set_variable  WORKDIR             ${ACOUSTICS_WORKDIR:-./}
set_variable  BACKGROUND          ${ACOUSTICS_BACKGROUND:+-B}
set_variable  FORMULA             "[0] + [1]*(sqrt(1.0 + 0.5*(x-[2])*(x-[2])/([3]*[3])) - 1.0)"

if ( do_usage $* ); then
    usage "$script <hydrophone file> [<string identifier>]"
fi

if (( $# > 2 )); then
    fatal "Wrong number of arguments."
fi

set_variable  HYDROPHONE_TXT      $1
set_variable  BASE_NAME           ${${HYDROPHONE_TXT%%\.*}##*/}
set_variable  HYDROPHONE_ROOT     $WORKDIR/${BASE_NAME}.root

source JAcoustics.sh   --

if [[ ! -f $HYDROPHONE_TXT ]]; then
    fatal "No sound hydrophone file $HYDROPHONE_TXT."
fi

JGraph \
    -f $HYDROPHONE_TXT                    \
    -o $HYDROPHONE_ROOT

sort -gr -k 2,2 $HYDROPHONE_TXT | tail -1 | read X Y

let XMIN="$X - 0.25"
let XMAX="$X + 0.25"

JFit \
    -f ${HYDROPHONE_ROOT}:\.\*      \
    -o $WORKDIR/fit.root            \
    -F "$FORMULA"                   \
    -x "$XMIN $XMAX"                \
    -@"p0 =   $Y"                   \
    -@"p1 =   1.0"                  \
    -@"p2 =   $X"                   \
    -@"p3 =   0.4"                  \
    -d $DEBUG

set_array P2 `JPrintFit -f $WORKDIR/fit.root:\.\* -@p2`

printf "result %10.5f +/- %10.5f\n" $P2[1] $P2[2]

if (( $# == 2 )); then

    set_variable  STRING              $2
    set_variable  PHI                 $P2[1]

    JEditHydrophone -f $WORKDIR/hydrophone.txt -S "$STRING rot $PHI"

else

    set_variable  STRING              ${BASE_NAME##*_}
fi

mv $WORKDIR/fit.root $HYDROPHONE_ROOT

set_variable  HYDROPHONE_GIF      $WORKDIR/${BASE_NAME}.gif

typeset -Z 4 STRING

JPlot1D \
    -f $HYDROPHONE_ROOT:\.\*            \
    -\> "#phi [rad]"                    \
    -\^ "RMS [ms]"                      \
    -T "$STRING"                        \
    -o $HYDROPHONE_GIF                  $BACKGROUND

rm -f  $HYDROPHONE_ROOT