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

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

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


set_variable  GRAPH_TXT           /tmp/graph.txt
set_variable  GRAPH_ROOT          /tmp/graph.root

if [[ $STAGE == "A" ]] || (( $STAGE > 0 )); then

    eval awk \'\/Processing stage $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

else

    let STAGE="-1 * $STAGE"
    
    egrep -v "(Processing|Elapsed)"  $INPUT_FILE                                       \
        |  awk '{ printf "%5d %8.3f\n", NR, $(NF) }'                                   >  $GRAPH_TXT
fi

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  "$STAGE"                \
    -o stage-${STAGE}.$FORMAT   $BATCH