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

if do_usage $*; then
    usage "$script (stage file)+"
fi

if (( $# >= 1 )); then
    set_array  INPUT_FILES    $argv[1,-1]
else
    fatal "Wrong number of arguments."
fi

for (( i = 1; $i <= ${#INPUT_FILES}; ++i )); do

    set_variable  INPUT_FILE          ${INPUT_FILES[$i]}
    set_variable  GRAPH_TXT           ${TMPDIR:-/tmp}/graph\[$i\].txt
    set_variable  GRAPH_ROOT          ${TMPDIR:-/tmp}/graph\[$i\].root

    eval awk \'\/Processing stage\/,\/Elapsed\/ \{ print \$0 \}\'  $INPUT_FILE  \
        |  head -n -1                                                           \
        |  tail -n +2                                                           >  ${TMPDIR:-/tmp}/graph.txt

    if (( `cat ${TMPDIR:-/tmp}/graph.txt | wc -l` == 0 )); then
        cat $INPUT_FILE                                                         >  ${TMPDIR:-/tmp}/graph.txt
    fi
        
    cat ${TMPDIR:-/tmp}/graph.txt                                               \
        |  grep -v "warning"                                                    \
        |  awk '{ printf "%5d %8.3f\n", NR, $(NF) }'                            >  $GRAPH_TXT
    
    JGraph \
        -f  $GRAPH_TXT          \
        -o  $GRAPH_ROOT
done

GRAPH_ROOT=(${TMPDIR:-/tmp}/graph\[{1..${#INPUT_FILES}}\].root)

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

rm -f ${TMPDIR:-/tmp}/graph*.{txt,root}