JAAnet.sh

Go to the documentation of this file.

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

# ------------------------------------------------------------------------------------------
#
#                         Utility script to test JAAnet.
#
# ------------------------------------------------------------------------------------------


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


source $JPP_DIR/setenv.sh $JPP_DIR


set_variable  DEBUG          2
set_variable  WORKDIR        ./
set_variable  INPUT_FILE     $WORKDIR/aanet.root
set_variable  OUTPUT_FILE    $WORKDIR/histogram.root

if ( do_usage $* ); then
    usage "$script [input file [output file]]"
fi

case $# in
    2) OUTPUT_FILE=$2;&
    1) INPUT_FILE=$1;;
esac

print_variable     INPUT_FILE OUTPUT_FILE
check_input_file   $INPUT_FILE


if (( 1 )) ; then

    JAAnet \
        -f $INPUT_FILE               \
        -o $OUTPUT_FILE              \
        -d $DEBUG

fi


if (( 1 )); then

    JPlot1D \
        -f ${OUTPUT_FILE}:hx                   \
        -\> "angle [deg]"                      \
        -\^ "number of events [a.u]"           \
        -XX                                    \
    -o angle.gif

fi

if (( 1 )); then

    JPlot1D \
        -f ${OUTPUT_FILE}:hd                   \
        -\> "distance [m]"                     \
        -\^ "number of events [a.u]"           \
    -o distance.gif

fi

if (( 1 )); then

    JPlot1D \
        -f ${OUTPUT_FILE}:ht                   \
        -\> "#Delta t [ns]"                    \
        -\^ "number of events [a.u]"        \
    -o time.gif

fi

if (( 1 )); then

    JFit \
        -f ${OUTPUT_FILE}:he                            \
        -F "[0]*exp(-0.5*(x-[1])*(x-[1])/([2]*[2]))"    \
        -@ "p0 = GetMaximum"                            \
        -@ "p1 = 0"                                     \
        -@ "p2 = 0.3"                                   \
        -x "-0.8 +0.8"

    JPlot1D \
        -f fit.root:he                                  \
        -\> "log10(E_{fit}/E_{true})"                   \
        -\^ "number of events [a.u]"                    \
    -o eratio.gif

fi