plot-perth.sh

Go to the documentation of this file.

#!/usr/bin/env 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   RECONSTRUCTION_PRECISION    0.1
set_variable:  FORMAT      GRAPHICS_FORMAT             gif
set_variable+  BATCH       GRAPHICS_BATCH              -B
 
if do_usage $*; then
    usage "$script <input file>"\
          "\nThe input file corresponds to the log file of JPerth."
fi
 
if (( $# != 1 )); then
    fatal "Wrong number of arguments."
fi
 
set_variable  INPUT_FILE  $1
 
grep "\‍(chi2\[.\]\|gradient\‍)" $INPUT_FILE  >  $WORKDIR/perth.txt
 
OUTPUT_FILES=()
 
awk '/chi2\[.\]/ { printf "%5d %9.5f\n", NR, -$3 }'                       $WORKDIR/perth.txt  >  $WORKDIR/g1.txt
awk '/gradient/  { printf "%5d %9.5f\n", NR,  $2 }'                       $WORKDIR/perth.txt  >  $WORKDIR/g2.txt
 
JGraph \
        -f $WORKDIR/g1.txt                       \
        -o $WORKDIR/g1.root                      \
        -T $TITLE
 
awk 'BEGIN { XMAX = -1.0e10; XMIN = +1.0e10 } END { printf "%9.5f %9.5f", XMIN, XMAX} { if ($1 < XMIN) { XMIN = $1 }; if ($1 > XMAX) { XMAX = $1 }}'  $WORKDIR/g1.txt  |  read XMIN XMAX
awk 'BEGIN { YMAX = -1.0e10; YMIN = +1.0e10 } END { printf "%9.5f %9.5f", YMIN, YMAX} { if ($2 < YMIN) { YMIN = $2 }; if ($2 > YMAX) { YMAX = $2 }}'  $WORKDIR/g1.txt  |  read YMIN YMAX
 
let "XMIN -= 1.0"
let "XMAX += 1.0"
 
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"
 
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/perth.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"                               \
        -\^ "likelihood"                         \
        -x "$XMIN $XMAX"                         \
        -y "$YMIN $YMAX"                         \
        -N "X 505"                               \
        -G Y                                     \
        -S 0.5                                   \
        -T "${${INPUT_FILE:t}:r}"                \
        -o  likelihood.$FORMAT                   $BATCH
 
JGraph \
        -f $WORKDIR/g2.txt                       \
        -o $WORKDIR/g2.root                      \
        -T $TITLE
 
awk 'BEGIN { YMAX = -1.0e10; YMIN = +1.0e10 } END { printf "%9.5f %9.5f", YMIN, YMAX} { if ($2 < YMIN) { YMIN = $2 }; if ($2 > YMAX) { YMAX = $2 }}'  $WORKDIR/g2.txt  |  read YMIN YMAX
 
let "Y    = $YMAX - $YMIN"
let "YMAX = floor(($YMAX + 0.1 * ($Y)) / $PRECISION + 1.0) * $PRECISION"
let "YMIN = 0.0"
 
JPlot1D \
        -w 1200x600                              \
        -f $WORKDIR/g2.root:\.\*                 \
        -> "step"                               \
        -\^ "|gradient|"                         \
        -x "$XMIN $XMAX"                         \
        -y "$YMIN $YMAX"                         \
        -N "X 505"                               \
        -G Y                                     \
        -S 0.5                                   \
        -T "${${INPUT_FILE:t}:r}"                \
        -o  gradient.$FORMAT                     $BATCH
 
rm -f $WORKDIR/perth.txt $WORKDIR/{g1,g2,L1}.{txt,root}