hydrophone-phi: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+  ALIGN       ACOUSTICS_ALIGN        Y     
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 <string identifier> [<hydrophone file>]"
fi
 
if (( $# == 0 || $# > 2 )); then
    fatal "Wrong number of arguments."
fi
 
set_variable  STRING              $argv[1]
 
set_variable  HYDROPHONE_TXT      $WORKDIR/hydrophonephi_${STRING}.txt
set_variable  HYDROPHONE_ROOT     $WORKDIR/hydrophonephi_${STRING}.root
 
if [[ ! -f $HYDROPHONE_TXT ]]; then
    fatal "No hydrophone data file $HYDROPHONE_TXT."
fi
 
sort -gr -k 2,2 $HYDROPHONE_TXT | tail -1 | read X Y
 
if [[ "$ALIGN" == "Y" ]]; then
   
    echo "Apply alignment."
 
    awk -v X=$X 'BEGIN { PI = atan2(0.0,-1.0) } { \
        x=($1 - X);  \
        if (x < -PI) { x += 2*PI };
        if (x > +PI) { x -= 2*PI };
        printf "%7.3f %12.5f\n", x, $2 \
        } ' $HYDROPHONE_TXT > ${TMPDIR:-/tmp}/hydrophone.txt
 
    X1=0.0
else
    
    cp $HYDROPHONE_TXT ${TMPDIR:-/tmp}/hydrophone.txt
 
    X1=$X
    X=0.0
fi
 
JGraph \
    -f ${TMPDIR:-/tmp}/hydrophone.txt          \
    -o $HYDROPHONE_ROOT
 
let "XMIN = $X1 - 0.35"
let "XMAX = $X1 + 0.35"
 
JFit \
    -f ${HYDROPHONE_ROOT}:\.\*      \
    -o $WORKDIR/fit.root            \
    -F "$FORMULA"                   \
    -x "$XMIN $XMAX"                \
    -@"p0 =   $Y"                   \
    -@"p1 =   1.0"                  \
    -@"p2 =   $X1"                  \
    -@"p3 =   0.4"                  \
    -d $DEBUG
 
set_array P2 `JPrintFit -f $WORKDIR/fit.root:\.\* -@p2`
 
P2[1]=$(($P2[1] + $X))
 
printf "result %10.5f +/- %10.5f [rad]\n" $P2[1] $P2[2]
printf "result %10.3f +/- %10.3f [deg]\n" $(($P2[1] * 180.0 / 3.1415)) $(($P2[2] * 180.0 / 3.1415))
 
if (( $# == 2 )); then
 
    set_variable  HYDROPHONE     $2
    set_variable  PHI            $P2[1]
 
    JEditHydrophone -f $HYDROPHONE -S "$STRING rot $PHI" -o $HYDROPHONE
fi
 
mv $WORKDIR/fit.root $HYDROPHONE_ROOT
 
typeset -Z 4 STRING
 
JPlot1D \
    -f ${HYDROPHONE_ROOT}:\.\*                     \
    -> "#phi [rad]"                               \
    -\^ "RMS [ms]"                                 \
    -T "[$STRING]"                                 \
    -o $WORKDIR/hydrophonephi_${STRING}.$FORMAT    $BATCH
 
rm -f  $HYDROPHONE_ROOT
rm -f  ${TMPDIR:-/tmp}/hydrophone.txt