tripod-Z: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 <tripod identifier> [<tripod file>]"
fi

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

set_variable  TRIPOD              $argv[1]

source JAcoustics.sh --

CHECK_EXIT_CODE

set_variable  TRIPOD_TXT          $WORKDIR/tripodz_${TRIPOD}.txt
set_variable  TRIPOD_ROOT         $WORKDIR/tripodz_${TRIPOD}.root
set_variable  G1_TXT              $WORKDIR/g1.txt
set_variable  G1_ROOT             $WORKDIR/g1.root

if [[ ! -f $TRIPOD_TXT ]] then
    fatal "No tripod data file  $TRIPOD_TXT."
fi

JGraph \
    -f $TRIPOD_TXT                            \
    -o $TRIPOD_ROOT

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

echo $X $Y > $G1_TXT

JGraph \
    -f $G1_TXT                                \
    -o $G1_ROOT

awk '
    BEGIN {
    XMIN = +1e3
    XMAX = -1e3
    }
    END { print XMIN, XMAX }
    {
    if ( $1 < XMIN ) { XMIN = $1 }
    if ( $1 > XMAX ) { XMAX = $1 }
    }
    ' $TRIPOD_TXT | read XMIN XMAX

if (( $X > $XMIN && $X < $XMAX )); then

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

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

    set_variable X  $P2[1]

    printf "nominal %2d %7.3f %7.3f\n" $TRIPOD 0.0 `JPrintFit -f $WORKDIR/fit.root:\.\* -x "0.0"`
    printf "optimal %2d %7.3f %7.3f\n" $TRIPOD $X  `JPrintFit -f $WORKDIR/fit.root:\.\* -x "$X"`

    mv $WORKDIR/fit.root $TRIPOD_ROOT

else

    echo   "Optimum value at limit; -> no fit."
    printf "nominal %2d %7.3f %7.3f\n" $TRIPOD 0.0 `awk -v X=0.0 '{ if ($1 == X) { print $2 }}' $TRIPOD_TXT`
    printf "optimal %2d %7.3f %7.3f\n" $TRIPOD $X  $Y

fi
    
if (( $# == 2 )); then
    JEditTripod -f $argv[2] -T "$TRIPOD add 0.0 0.0 $X"
fi

JPlot1D \
    -f $TRIPOD_ROOT:\.\*                      \
    -f ${G1_ROOT}:\.\*                        \
    -\> "#Deltaz [m]"                         \
    -\^ "#chi^{2}/NDF"                        \
    -T "[$TRIPOD]"                            \
    -o $WORKDIR/tripodz_${TRIPOD}.gif         $BACKGROUND

rm -f  $TRIPOD_ROOT $G1_TXT $G1_ROOT