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   ALIGN       ${ACOUSTICS_ALIGN}     
set_variable   FORMULA     "[0] + [1]*(sqrt(1.0 + 0.5*(x-[2])*(x-[2])/([3]*[3])) - 1.0)"  

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" != "" ]]; then
   
    awk -v X=$X 'BEGIN { PI = atan2(0.0,-1.0) } { \
        x=($1 - X);  \
        if (x < -PI) { x += 2PI };
        if (x > +PI) { x -= 2*PI };
        printf "%7.3f %12.5f\n", x, $2 \
        } ' $HYDROPHONE_TXT > /tmp/hydrophone.txt

    X1=0.0
else
    
    cp $HYDROPHONE_TXT /tmp/hydrophone.txt

    X1=$X
    X=0.0
fi

JGraph \
    -f /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\n" $P2[1] $P2[2]

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

    set_variable  HYDROPHONE_FILE     $2
    set_variable  PHI                 $P2[1]

    JEditHydrophone -f $HYDROPHONE_FILE -S "$STRING rot $PHI"
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}.gif   $BACKGROUND

rm -f  $HYDROPHONE_ROOT
rm -f  /tmp/hydrophone.txt