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