plot-perth.sh

Go to the documentation of this file.

#!/bin/zsh
 
if [ -z $JPP_DIR ]; then
    echo "Variable JPP_DIR undefined."
    exit
fi
 
source $JPP_DIR/setenv.sh $JPP_DIR
 
zmodload zsh/mathfunc
 
set_variable   DEBUG       2
set_variable   WORKDIR     ${TMPDIR:-/tmp}/
set_variable:  PRECISION   ACOUSTICS_PRECISION    1.0
set_variable:  FORMAT      GRAPHICS_FORMAT        gif
set_variable+  BATCH       GRAPHICS_BATCH         -B
 
if do_usage $*; then
    usage "$script <input file>"
fi
 
if (( $# != 1 )); then
    fatal "Wrong number of arguments."
fi
 
set_variable INPUT_FILE      $argv[1]
 
grep "\‍(chi2\[.\]\|gradient\‍)" $INPUT_FILE  >  $WORKDIR/perth.txt
 
awk 'BEGIN { YMAX = -1.0e10; YMIN = +1.0e10 }   \
    END { printf "%9.5f %9.5f", YMIN, YMAX}     \
    /chi2/ { if ($3 < YMIN) { YMIN = $3}; if ($3 > YMAX) { YMAX = $3 }}'  $WORKDIR/perth.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"
 
awk '/chi2\[.\]/ { printf "%5d %9.5f\n", NR, $3 }'                        $WORKDIR/perth.txt  >  $WORKDIR/g1.txt
 
JGraph \
    -f $WORKDIR/g1.txt                       \
    -o $WORKDIR/g1.root                      \
    -T $TITLE
 
awk \
    -v YMIN=$YMIN \
    -v YMAX=$YMAX \
    '/gradient/  { printf "%5d %9.5f %5d %9.5f\n", NR, YMIN, NR, YMAX }'  $WORKDIR/perth.txt  >  $WORKDIR/L1.txt  
 
JLine \
    -p $WORKDIR/L1.txt                       \
    -o $WORKDIR/L1.root                      \
    -@ "width = 1"                           \
    -@ "style = 2"                           
 
JPlot1D \
    -w 1200x600                              \
    -f $WORKDIR/g1.root:\.\*                 \
    -f $WORKDIR/L1.root:\.\*                 \
    -> "step"                               \
    -\^ "#chi^{2}"                           \
    -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/g1.root $WORKDIR/L1.root  $WORKDIR/perth.txt