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