Jpp
JGizmo/JHistogram2D.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 JHistogram2D, JFit2D and JPlot2D.
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 /tmp/
26 
27 if ( do_usage $* ); then
28  usage "$script [working directory]"
29 fi
30 
31 case $# in
32  1) WORKDIR=$1;;
33 esac
34 
35 
36 if (( 1 )); then
37 
38  set_variable FORMULA "[0] * exp(-0.5 * (x-[1])*(x-[1]) / ([2]*[2])) * exp(-0.5 * (y-[1])*(y-[1]) / ([2]*[2]))"
39  set_variable OUTPUT_FILE histogram.root
40 
41  JHistogram2D \
42  -o $WORKDIR/$OUTPUT_FILE \
43  -F "$FORMULA" \
44  -="p0 = 1.0; p1 = 0.0; p2 = 2.0" \
45  -x "100 -10 +10" \
46  -y "100 -10 +10" \
47  -n 100000
48 
49 
50  JFit2D \
51  -f $WORKDIR/${OUTPUT_FILE}:h0 \
52  -o $WORKDIR/fit.root \
53  -F "$FORMULA" \
54  -@"p0 = GetMaximum" \
55  -@"p1 = GetMean" \
56  -@"p2 = GetRMS" \
57  -d $DEBUG
58 
59 
60  JPlot2D \
61  -f $WORKDIR/fit.root:h0 \
62  -O LEGO
63 
64 fi