JRate.sh

Go to the documentation of this file.

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

# ------------------------------------------------------------------------------------------
#
#                   Utility script to plot event rate of astrophysical source.
#
# ------------------------------------------------------------------------------------------


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        $JPP_DATA/trigger_efficiency+background.root
set_variable     OUTPUT_FILE       rate.root
set_variable     NUMBER_OF_BLOCKS  2


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

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


if (( 1 )); then 

    print_variable INPUT_FILE OUTPUT_FILE
    check_input_file  $INPUT_FILE
    
    $JPP_DIR/examples/JAstronomy/JRate \
        -f $INPUT_FILE          \
        -o $OUTPUT_FILE         \
        -N $NUMBER_OF_BLOCKS    \
        -d $DEBUG --!
    
fi


if (( 1 )); then

    JPlot1D \
        -f ${OUTPUT_FILE}:ct     \
        -\> "cos(#theta)"        \
        -\^ "dP/d#Omega"         \
        -L TR

fi    


JPlot1D \
    -f ${OUTPUT_FILE}:h0     \
    -f ${OUTPUT_FILE}:h1     \
    -y "0.01 4"              \
    -XXY                     \
    -\> "E [GeV]"            \
    -\^ "Rate [events/year]" \
    -L TR                    \
    -o rate.gif