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  
set_variable:  FORMAT       GRAPHICS_FORMAT          gif
set_variable+  BATCH        GRAPHICS_BATCH           -B

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

case $# in
    2) set_variable OUTPUT_FILE $2;&
    1) set_variable 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.$FORMAT                       $BATCH

fi

if (( 1 )); then

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

fi

if (( 1 )); then

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

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.$FORMAT                               $BATCH

fi