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  BACKGROUND          ${ACOUSTICS_BACKGROUND:+-B}
set_variable  FORMULA             "[0] + [1] * log(1.0 + 0.5*(x-[2])*(x-[2])/([3]*[3]))"

if ( do_usage $* ); then
    usage "$script [<tripod file>]"
fi

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

source JAcoustics.sh --

CHECK_EXIT_CODE

set_variable  TRIPOD3Z_TXT        $WORKDIR/tripod3z.txt

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

sort -gr -k 4,4 $TRIPOD3Z_TXT | 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

    JGraph \
        -f $TRIPOD_TXT                            \
        -o $TRIPOD_ROOT

    sort -gr -k 2,2 $WORKDIR/tripod3z_${TRIPOD}.txt | tail -1 | read X Y EY

    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`

    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 $P2[1] `JPrintFit -f $WORKDIR/fit.root:\.\* -x "$P2[1]"`

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


    mv $WORKDIR/fit.root $TRIPOD_ROOT

    sort -gr -k 1,1 $WORKDIR/tripod3z_${TRIPOD}.txt | tail -1 | read XMIN Y EY
    sort -g  -k 1,1 $WORKDIR/tripod3z_${TRIPOD}.txt | tail -1 | read XMAX Y EY

    let XMIN="$XMIN - 0.05"
    let XMAX="$XMAX + 0.05"

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

    rm -f $WORKDIR/tripod3z_${TRIPOD}.txt
    rm -f $WORKDIR/tripod3z_${TRIPOD}.root

done