Jpp 19.3.0-rc.1
the software that should make you happy
Loading...
Searching...
No Matches
JRate.sh
Go to the documentation of this file.
1#!/usr/bin/env zsh
2#
3# \author mdejong
4#
5version=1.0
6script=${0##*/}
7
8# ------------------------------------------------------------------------------------------
9#
10# Utility script to plot event rate of astrophysical source.
11#
12# ------------------------------------------------------------------------------------------
13
14
15if [ -z $JPP_DIR ]; then
16 echo "Variable JPP_DIR undefined."
17 exit
18fi
19
20
21source $JPP_DIR/setenv.sh $JPP_DIR
22
23
24set_variable DEBUG 2
25set_variable WORKDIR ./
26set_variable INPUT_FILE $JPP_DATA/trigger_efficiency+background.root
27set_variable OUTPUT_FILE rate.root
28set_variable NUMBER_OF_BLOCKS 2
29set_variable: FORMAT GRAPHICS_FORMAT gif
30set_variable+ BATCH GRAPHICS_BATCH -B
31
32if do_usage $*; then
33 usage "$script [input file [output file]]"
34fi
35
36case $# in
37 2) set_variable OUTPUT_FILE $2;&
38 1) set_variable INPUT_FILE $1;;
39esac
40
41
42if (( 1 )); then
43
44 print_variable INPUT_FILE OUTPUT_FILE
45 check_input_file $INPUT_FILE
46
47 $JPP_DIR/examples/JAstronomy/JRate \
48 -f $INPUT_FILE \
49 -o $OUTPUT_FILE \
50 -N $NUMBER_OF_BLOCKS \
51 -d $DEBUG --!
52
53fi
54
55
56if (( 1 )); then
57
58 JPlot1D \
59 -f ${OUTPUT_FILE}:ct \
60 -> "cos(#theta)" \
61 -\^ "dP/d#Omega" \
62 -L TR
63
64fi
65
66
67JPlot1D \
68 -f ${OUTPUT_FILE}:h0 \
69 -f ${OUTPUT_FILE}:h1 \
70 -y "0.01 4" \
71 -XXY \
72 -> "E [GeV]" \
73 -\^ "Rate [events/year]" \
74 -L TR \
75 -o rate.$FORMAT $BATCH