sound-V:fit.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       ACOUSTICS_DEBUG        2
set_variable:  WORKDIR     ACOUSTICS_WORKDIR      ./
set_variable   FORMULA     "[0] + [1]*(sqrt(1.0 + 0.5*(x-[2])*(x-[2])/([3]*[3])) - 1.0)"  
set_variable:  FORMAT      GRAPHICS_FORMAT        gif
set_variable+  BATCH       GRAPHICS_BATCH         -B
 
if do_usage $*; then
    usage "$script"
fi
 
if (( $# != 0 )); then
    fatal "Wrong number of arguments."
fi
 
set_variable  VELOCITY_TXT        $WORKDIR/soundv.txt
set_variable  VELOCITY_ROOT       $WORKDIR/soundv.root
 
if [[ ! -f $VELOCITY_TXT ]]; then
    fatal "No sound velocity file $VELOCITY_TXT."
fi
 
JGraph \
    -f $VELOCITY_TXT                    \
    -o $VELOCITY_ROOT
 
sort -gr -k 2,2 $VELOCITY_TXT | tail -1 | read X Y RMS
 
JFit \
    -f ${VELOCITY_ROOT}:\.\*        \
    -o $WORKDIR/fit.root            \
    -F "$FORMULA"                   \
    -@"p0 =   $Y"                   \
    -@"p1 =   1.0"                  \
    -@"p2 =   $X"                   \
    -@"p3 =   1.0"                  \
    -d $DEBUG
 
set_array P2 `JPrintFit -f $WORKDIR/fit.root:\.\* -@p2`
 
printf "result %10.5f +/- %10.5f\n" $P2[1] $P2[2]
 
mv $WORKDIR/fit.root $VELOCITY_ROOT
 
JPlot1D \
    -f $VELOCITY_ROOT:\.\*              \
    -> "V [m/s]"                       \
    -\^ "#chi^{2}/NDF"                  \
    -T ""                               \
    -o $WORKDIR/velocity.$FORMAT        $BATCH
 
rm -f  $VELOCITY_ROOT