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     BACKGROUND             ${COMPASS_BACKGROUND:+-B}
set_variable     DIR                    $JPP_DIR/examples/JCompass/

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.gif                  $BACKGROUND

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

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

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.gif              $BACKGROUND

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.gif              $BACKGROUND

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.gif              $BACKGROUND

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.gif              $BACKGROUND