Jpp
JAAnet.sh
Go to the documentation of this file.
1 #!/bin/zsh
2 #
3 # \author mdejong
4 #
5 version=1.0
6 script=${0##*/}
7 
8 # ------------------------------------------------------------------------------------------
9 #
10 # Utility script to test JAAnet.
11 #
12 # ------------------------------------------------------------------------------------------
13 
14 
15 if [ -z $JPP_DIR ]; then
16  echo "Variable JPP_DIR undefined."
17  exit
18 fi
19 
20 
21 source $JPP_DIR/setenv.sh $JPP_DIR
22 
23 
24 set_variable DEBUG 2
25 set_variable WORKDIR ./
26 set_variable INPUT_FILE $WORKDIR/aanet.root
27 set_variable OUTPUT_FILE $WORKDIR/histogram.root
28 
29 if ( do_usage $* ); then
30  usage "$script [input file [output file]]"
31 fi
32 
33 case $# in
34  2) OUTPUT_FILE=$2;&
35  1) INPUT_FILE=$1;;
36 esac
37 
38 print_variable INPUT_FILE OUTPUT_FILE
39 check_input_file $INPUT_FILE
40 
41 
42 if (( 1 )) ; then
43 
44  JAAnet \
45  -f $INPUT_FILE \
46  -o $OUTPUT_FILE \
47  -d $DEBUG
48 
49 fi
50 
51 
52 if (( 1 )); then
53 
54  JPlot1D \
55  -f ${OUTPUT_FILE}:hx \
56  -> "angle [deg]" \
57  -\^ "number of events [a.u]" \
58  -XX \
59  -o angle.gif
60 
61 fi
62 
63 if (( 1 )); then
64 
65  JPlot1D \
66  -f ${OUTPUT_FILE}:hd \
67  -> "distance [m]" \
68  -\^ "number of events [a.u]" \
69  -o distance.gif
70 
71 fi
72 
73 if (( 1 )); then
74 
75  JPlot1D \
76  -f ${OUTPUT_FILE}:ht \
77  -> "#Delta t [ns]" \
78  -\^ "number of events [a.u]" \
79  -o time.gif
80 
81 fi
82 
83 if (( 1 )); then
84 
85  JFit \
86  -f ${OUTPUT_FILE}:he \
87  -F "[0]*exp(-0.5*(x-[1])*(x-[1])/([2]*[2]))" \
88  -@ "p0 = GetMaximum" \
89  -@ "p1 = 0" \
90  -@ "p2 = 0.3" \
91  -x "-0.8 +0.8"
92 
93  JPlot1D \
94  -f fit.root:he \
95  -> "log10(E_{fit}/E_{true})" \
96  -\^ "number of events [a.u]" \
97  -o eratio.gif
98 
99 fi
100 
Generated by   doxygen 1.8.16