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

if do_usage $*; then
    usage "$script <string identifier>"
fi

if (( $# != 1 )); then
    fatal "Wrong number of arguments."
fi

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

echo $A $B > /tmp/graph.txt

JGraph \
    -f /tmp/graph.txt                          \
    -o /tmp/graph.root

typeset -Z 4 STRING

set_variable  DZ                  1.50

sort -nr -k 3,3 $MECHANICS_TXT | tail -1 | read A B CHI2 RMS

let ZMIN="$CHI2"
let ZMAX="$CHI2 + $DZ"

JPlot2D \
    -f ${MECHANICS_ROOT}:\.\*                    \
    -f /tmp/graph.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 /tmp/graph.{txt,root}