plot-acoustics-fit.sh

Go to the documentation of this file.

#!/bin/zsh
script=${0##*/}

# --------------------------------------------------------------------------------------------
#
# Auxiliary script to plot chi2 evolution of acoustics-fit.sh
#
# --------------------------------------------------------------------------------------------

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:  FORMAT      GRAPHICS_FORMAT        gif
set_variable+  BATCH       GRAPHICS_BATCH         -B

if do_usage $*; then
    usage "$script <stage file> [YMIN YMAX]"
fi

case $# in
    3)  set_variable  YMAX          $argv[3];
        set_variable  YMIN          $argv[2];&
    1)  set_variable  INPUT_FILE    $argv[1];;
    *)  fatal "Wrong number of arguments."
esac

set_variable  GRAPH_TXT           ${TMPDIR:-/tmp}/graph.txt
set_variable  GRAPH_ROOT          ${TMPDIR:-/tmp}/graph.root

eval awk \'\/Processing stage\/,\/Elapsed\/ \{ print \$0 \}\'  $INPUT_FILE  \
    |  head -n -1                                                                  \
    |  tail -n +2                                                                  \
    |  grep -v "warning"                                                           \
    |  awk '{ printf "%5d %8.3f\n", NR, $(NF) }'                                   >  $GRAPH_TXT

JGraph \
    -f  $GRAPH_TXT          \
    -o  $GRAPH_ROOT

if [[ "$YMIN" != "" && "$YMAX" != "" ]]; then
    set_variable OPTION   -y "$YMIN $YMAX"
fi

JPlot1D \
    -f  $GRAPH_ROOT:\.\*                     \
    $OPTION                                  \
    -> "step"                               \
    -\^ "#chi^{2}/NDF"                       \
    -S  0.5                                  \
    -T  ""                                   \
    -o  ${${${INPUT_FILE}##*/}%%.*}.$FORMAT  $BATCH