JRose.sh

Go to the documentation of this file.

#!/bin/zsh
#
#
# \author mdejong
#
version=1.0
script=${0##*/}

if [ -z $JPP_DIR ]; then
    echo "Variable JPP_DIR undefined."
    exit
fi

source $JPP_DIR/setenv.sh $JPP_DIR

set_variable:  DEBUG       COMPASS_DEBUG                2
set_variable:  WORKDIR     COMPASS_WORKDIR              ./
set_variable   DIR         $JPP_DIR/examples/JCompass/  
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 parameters."
fi

set_variable     INPUT_FILE              $argv[1]

$DIR/JRose \
    -f ${INPUT_FILE}                      \
    -o $WORKDIR/rose.root                 \
    -d $DEBUG

JPlot1D \
    -f "$WORKDIR/rose.root:HC.*"          \
    -\> "#chi^{2}/NDF"                    \
    -XX                                   \
    -T ""                                 \
    -L TR                                 \
    -o $WORKDIR/chi2.$FORMAT              $BATCH

JPlot1D \
    -f "$WORKDIR/rose.root:H0.*"          \
    -\> "residual [deg]"                  \
    -T "Q0"                               \
    -N "X 505"                            \
    -L TR                                 \
    -o $WORKDIR/Q0.$FORMAT                $BATCH

JPlot1D \
    -f "$WORKDIR/rose.root:H1.*"          \
    -\> "residual [deg/m]"                \
    -T "Q1"                               \
    -N "X 505"                            \
    -L TR                                 \
    -o $WORKDIR/Q1.$FORMAT                $BATCH

JPlot1D \
    -w 1200x600                           \
    -f "$WORKDIR/rose.root:.*\.Q0\.twist" \
    -\^ "twist [rad]"                     \
    -y "-3.2 +3.2"                        \
    -t "$TIMESTAMP"                       \
    -N "X 505"                            \
    -T "Q0"                               \
    -L TL                                 \
    -o $WORKDIR/twist.Q0.$FORMAT          $BATCH

JPlot1D \
    -w 1200x600                           \
    -f "$WORKDIR/rose.root:.*\.Q1\.twist" \
    -y "-0.005 +0.005"                    \
    -\^ "twist [rad/m]"                   \
    -t "$TIMESTAMP"                       \
    -N "X 505"                            \
    -T "Q1"                               \
    -L TL                                 \
    -o $WORKDIR/twist.Q1.$FORMAT          $BATCH

JPlot1D \
    -w 1200x600                           \
    -f "$WORKDIR/rose.root:.*\.Q0\.swing" \
    -\^ "swing [rad]"                     \
    -y "0.0 0.04"                         \
    -t "$TIMESTAMP"                       \
    -N "X 505"                            \
    -T "Q0"                               \
    -L TL                                 \
    -o $WORKDIR/swing.Q0.$FORMAT          $BATCH

JPlot1D \
    -w 1200x600                           \
    -f "$WORKDIR/rose.root:.*\.Q0\.atan2" \
    -\^ "atan2 [rad]"                     \
    -y "-3.2 +3.2"                        \
    -t "$TIMESTAMP"                       \
    -N "X 505"                            \
    -T "Q0"                               \
    -L TL                                 \
    -o $WORKDIR/atan2.Q0.$FORMAT          $BATCH