mechanics:plot.sh

Go to the documentation of this file.

#!/bin/zsh
script=${0##*/}
 
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
set_variable   ZRANGE      0.0
 
if do_usage $*; then
    usage "$script <string identifier> [z-range]"
fi
 
case $# in
    2) let           "ZRANGE = $argv[2]";&
    1) set_variable  STRING  $argv[1];;
    *) fatal "Wrong number of arguments."
esac
 
set_variable  STRING              $argv[1]
 
set_variable  MECHANICS_TXT       $WORKDIR/mechanicsab_${STRING}.txt
set_variable  MECHANICS_ROOT      $WORKDIR/mechanicsab_${STRING}.root
 
if [[ ! -f $MECHANICS_TXT ]] then
    fatal "No mechanics data file $MECHANICS_TXT."
fi
 
JGraph2D \
    -f $MECHANICS_TXT                            \
    -o $MECHANICS_ROOT
 
JPrintMechanics -M $WORKDIR/mechanics.txt -S $STRING | read DUMMY A B
 
awk -v A=$A -v B=$B 'BEGIN { DA=1e6; DB = 1e6; ZMIN = 1e6; } END { print ZMIN } \
    { \
    if (($1 - A >= -DA && $1 - A <= +DA) && ($2 - B >= -DB && $2 - B <= +DB)) { \
    ZMIN = $3;
    if ($1 > A) { DA = $1 - A} else { DA = A - $1};  \
    if ($2 > B) { DB = $2 - B} else { DB = B - $2};  \
    }}' $MECHANICS_TXT  |  read ZMIN
 
printf "current %4d %7.5f %7.3f %7.3f\n"  $STRING  $A  $B  $ZMIN
 
echo $A $B > ${TMPDIR:-/tmp}/graph_1.txt
 
JGraph \
    -f ${TMPDIR:-/tmp}/graph_1.txt               \
    -o ${TMPDIR:-/tmp}/graph_1.root
 
typeset -Z 4 STRING
 
sort -n  -k 3,3 $MECHANICS_TXT | tail -1 | read A B ZMAX RMS
sort -nr -k 3,3 $MECHANICS_TXT | tail -1 | read A B ZMIN RMS
 
printf "optimal %4d %7.5f %7.3f %7.3f\n"  $STRING  $A  $B  $ZMIN
 
echo $A $B > ${TMPDIR:-/tmp}/graph_2.txt
 
JGraph \
    -f ${TMPDIR:-/tmp}/graph_2.txt               \
    -o ${TMPDIR:-/tmp}/graph_2.root
 
if (( $ZRANGE > 0.0 )); then
    let "ZMAX = $ZMIN + $ZRANGE"
fi
 
JPlot2D \
    -f ${MECHANICS_ROOT}:\.\*                    \
    -f ${TMPDIR:-/tmp}/graph_1.root:\.\*         \
    -f ${TMPDIR:-/tmp}/graph_2.root:\.\*         \
    -> "a [m^{-1}]"                             \
    -< "b [m]"                                  \
    -\^ "#chi^{2}/NDF"                           \
    -N "X 505"                                   \
    -N "Y 510"                                   \
    -z "$ZMIN $ZMAX"                             \
    -O COLZ                                      \
    -G XY                                        \
    -S 1.5                                       \
    -T "[$STRING]"                               \
    -o $WORKDIR/mechanicsab_${STRING}.$FORMAT    $BATCH
 
rm -f  $MECHANICS_ROOT
rm -f ${TMPDIR:-/tmp}/graph_{1,2}.{txt,root}