hydrophone-Z: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   FORMULA     "[0] + [1]*((x-[2])*(x-[2])/([3]*[3]))"  
set_variable:  FORMAT      GRAPHICS_FORMAT        gif
set_variable+  BATCH       GRAPHICS_BATCH         -B
 
if do_usage $*; then
    usage "$script <string identifier> [<detector file>]"
fi
 
if (( $# == 0 || $# > 2 )); then
    fatal "Wrong number of arguments."
fi
 
set_variable  STRING              $argv[1]
 
set_variable  HYDROPHONE_TXT      $WORKDIR/hydrophonez_${STRING}.txt
set_variable  HYDROPHONE_ROOT     $WORKDIR/hydrophonez_${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
 
JGraph \
    -f $HYDROPHONE_TXT              \
    -o $HYDROPHONE_ROOT
 
let "XMIN = $X - 2.5"
let "XMAX = $X + 2.5"
 
JFit \
    -f ${HYDROPHONE_ROOT}:\.\*      \
    -o $WORKDIR/fit.root            \
    -F "$FORMULA"                   \
    -x "$XMIN $XMAX"                \
    -@"p0 =   $Y"                   \
    -@"p1 =   0.001"                \
    -@"p2 =   $X"                   \
    -@"p3 =   0.3"                  \
    -d $DEBUG
 
set_array P2 `JPrintFit -f $WORKDIR/fit.root:\.\* -@p2`
 
printf "result %7.3f +/- %7.3f [m]\n" $P2[1] $P2[2]
 
if (( $# == 2 )); then
 
    set_variable  DETECTOR            $2
    set_variable  MODULE              `getModule -a $DETECTOR -L "$STRING 0"`
    set_variable  Z                   $P2[1]
 
    JEditDetector -a $DETECTOR -o $DETECTOR -M "$MODULE add 0.0 0.0 $Z"
fi
 
mv $WORKDIR/fit.root $HYDROPHONE_ROOT
 
typeset -Z 4 STRING
 
JPlot1D \
    -f ${HYDROPHONE_ROOT}:\.\*                     \
    -> "#Deltaz [m]"                              \
    -\^ "RMS [ms]"                                 \
    -T "[$STRING]"                                 \
    -o $WORKDIR/hydrophonez_${STRING}.$FORMAT      $BATCH
 
rm -f  $HYDROPHONE_ROOT
rm -f  ${TMPDIR:-/tmp}/hydrophone.txt