JStarTrek.sh

Go to the documentation of this file.

#!/usr/bin/env 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   OUTPUT_FILE       startrek.root                                     
set_variable   NUMBER_OF_BLOCKS  2                                             
set_variable:  FORMAT            GRAPHICS_FORMAT        gif
set_variable+  BATCH             GRAPHICS_BATCH         -B
 
if do_usage $*; then
    usage "$script <input file>"
fi
 
if (( $# != 1 )); then
    fatal "Wrong number of arguments."
fi
 
set_variable  INPUT_FILE   $1
 
 
if (( 1 )); then 
 
    print_variable INPUT_FILE OUTPUT_FILE
    check_input_file  $INPUT_FILE
    
    $JPP_DIR/examples/JAstronomy/JStarTrek \
            -f $INPUT_FILE          \
            -o $OUTPUT_FILE         \
            -N $NUMBER_OF_BLOCKS    \
        -o $OUTPUT_FILE         \
            -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.$FORMAT          $BATCH