tripod-3Z: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   FORMULA     "[0] + [1] * log(1.0 + 0.5*(x-[2])*(x-[2])/([3]*[3]))"  
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 [<tripod file>]"
fi

if (( $# > 1 )); then
    fatal "Wrong number of arguments."
fi

source JAcousticsToolkit.sh

CHECK_EXIT_CODE

set_variable  TRIPOD3Z_TXT        $WORKDIR/tripod3z.txt

if [[ ! -f $TRIPOD3Z_TXT ]] then
    fatal "No detector data file $TRIPOD3Z_TXT."
fi

grep -v '^#' $TRIPOD3Z_TXT | sort -gr -k 4,4 | tail -1 | read X Y Z CHI2 RMS

printf "optimum %7.3f %7.3f %7.3f %7.3f %7.3f\n" $X $Y $Z $CHI2 $RMS

awk                 '$1  ~ /^#/ { for (i = 2; i <= NF; ++i ) { printf " %s", $i }}' $TRIPOD3Z_TXT  |  read A B C

awk -v Y=$Y -v Z=$Z '$1 !~ /^#/ { if ( $2 == Y && $3 == Z )  { printf "%7.3f %7.3f %7.3f\n", $1, $4, $5 } }' $TRIPOD3Z_TXT  >!  $WORKDIR/tripod3z_${A}.txt
awk -v X=$X -v Z=$Z '$1 !~ /^#/ { if ( $1 == X && $3 == Z )  { printf "%7.3f %7.3f %7.3f\n", $2, $4, $5 } }' $TRIPOD3Z_TXT  >!  $WORKDIR/tripod3z_${B}.txt
awk -v X=$X -v Y=$Y '$1 !~ /^#/ { if ( $1 == X && $2 == Y )  { printf "%7.3f %7.3f %7.3f\n", $3, $4, $5 } }' $TRIPOD3Z_TXT  >!  $WORKDIR/tripod3z_${C}.txt

for TRIPOD in $A $B $C; do

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

    JGraph \
        -f $TRIPOD_TXT                            \
        -o $TRIPOD_ROOT

    sort -gr -k 2,2 $WORKDIR/tripod3z_${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 (( $# == 1 )); then
        JEditTripod -f $argv[1] -T "$TRIPOD add 0.0 0.0 $X" -o $argv[1]
    fi

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

    rm -f $TRIPOD_TXT $TRIPOD_ROOT $G1_TXT $G1_ROOT

done