plot-acoustics-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

zmodload zsh/mathfunc

set_variable   DEBUG       2
set_variable   WORKDIR     ${TMPDIR:-/tmp}/
set_variable:  PRECISION   ACOUSTICS_PRECISION    0.1
set_variable:  FORMAT      GRAPHICS_FORMAT        gif
set_variable+  BATCH       GRAPHICS_BATCH         -B

if do_usage $*; then
    usage "$script <directory>"\
          "\nInput directory corresponds to location of stage files of a pre-calibration procedure."
fi

if (( $# == 1 )); then 
    set_variable  DIR  $1
else
    fatal "Wrong number of arguments."
fi

for INPUT_FILE in `ls -rt $DIR/stage-*.log`; do

    grep "\(chi2\[.\]\|gradient\)" $INPUT_FILE  >  $WORKDIR/stage.txt

    awk 'BEGIN { YMAX = 0.0; YMIN = 1.0e10 } END { printf "%9.5f %9.5f", YMIN, YMAX} /chi2/ { if ($3 < 1.0e10) { if ($3 < YMIN) { YMIN = $3}; if ($3 > YMAX) { YMAX = $3 }}}'  $WORKDIR/stage.txt  |  read YMIN YMAX

    let "Y    = $YMAX - $YMIN"
    let "YMAX = floor(($YMAX + 0.1 * ($Y)) / $PRECISION + 1.0) * $PRECISION"
    let "YMIN = floor(($YMIN - 0.1 * ($Y)) / $PRECISION - 0.0) * $PRECISION"

    if (( $YMIN < 0.0 )); then
        let "YMIN = 0.0"
    fi

    OUTPUT_FILES=()

    awk '/chi2\[.\]/ { printf "%5d %9.5f\n", NR, $3 }'                        $WORKDIR/stage.txt  >  $WORKDIR/g1.txt

    JGraph \
        -f $WORKDIR/g1.txt                       \
        -o $WORKDIR/g1.root                      \
        -T $TITLE

    OUTPUT_FILES+=$WORKDIR/g1.root

    awk \
        -v YMIN=$YMIN \
        -v YMAX=$YMAX \
        '/gradient/  { printf "%5d %9.5f %5d %9.5f\n", NR, YMIN, NR, YMAX }'  $WORKDIR/stage.txt  >  $WORKDIR/L1.txt  

    JLine \
        -p $WORKDIR/L1.txt                       \
        -o $WORKDIR/L1.root                      \
        -@ "width = 1"                           \
        -@ "style = 2"                           

    OUTPUT_FILES+=$WORKDIR/L1.root

    JPlot1D \
        -w 1200x600                              \
        -f${^OUTPUT_FILES}":.*"                  \
        -\> "step"                               \
        -\^ "#chi^{2}/NDF"                       \
        -y "$YMIN $YMAX"                         \
        -N "X 505"                               \
        -G Y                                     \
        -S 0.5                                   \
        -T "${${INPUT_FILE:t}:r}"                \
        -o  ${${INPUT_FILE:t}:r}.$FORMAT         $BATCH

    rm -f $WORKDIR/stage.txt $WORKDIR/g1.root $WORKDIR/L1.root
done