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  BACKGROUND          ${ACOUSTICS_BACKGROUND:+-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                  0.20

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 505"                                \
    -z "$ZMIN $ZMAX"                          \
    -O COLZ                                   \
    -S 1.5                                    \
    -T "[$STRING]"                            \
    -o $WORKDIR/mechanicsab_${STRING}.gif     $BACKGROUND

rm -f  $MECHANICS_ROOT
rm -f /tmp/graph.{txt,root}