JPizza.sh

Go to the documentation of this file.

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

# ------------------------------------------------------------------------------------------
#
#                         Utility script to run JPizza.
#
# ------------------------------------------------------------------------------------------


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      ${TMPDIR:-/tmp}/                
set_variable   OUTPUT_FILE  $WORKDIR/pizza.root  


if do_usage $*; then
    usage "$script (input file)+ <output file>"\
          "\nNote that if more than one input file is specified, all other arguments must be provided."
fi

case $# in
    2) set_variable OUTPUT_FILE  $2;&
    1) set_variable INPUT_FILE   $1;;
    0) exit;;
    *) set_variable INPUT_FILE   $argv[1,-2]
       set_variable OUTPUT_FILE  $argv[-1];;
esac


if (( 1 )); then

    JPizza \
        -f ${INPUT_FILE}                \
        -o ${OUTPUT_FILE}               \
        -d ${DEBUG} --!

    JPrintDomino \
        -f ${OUTPUT_FILE}               

fi


if (( 1 )); then 

    JPlot1D \
        -f ${OUTPUT_FILE}:hn              \
        -y "0.5 1.0e3" -Y                 \
        -\> "number of muons"             \
        -\^ "number of events [a.u.]"     \
        -T ""

fi


if (( 1 )); then 

    JPlot1D \
        -f ${OUTPUT_FILE}:he              \
        -y "0.5 1.0e3" -Y                 \
        -\> "E_{#mu}"                     \
        -\^ "number of events [a.u.]"     \
        -T ""

fi


if (( 1 )); then

    JPlot2D \
        -f ${OUTPUT_FILE}:h2              \
        -O COLZ                           \
        -\> "R^{2} [m^{2}]"               \
        -\< "z [m]"

fi