Jpp  master_rocky-40-g5f0272dcd
the software that should make you happy
JRootfitToGauss2D.sh
Go to the documentation of this file.
1 #!/bin/zsh
2 version=1.0
3 script=${0##*/}
4 
5 source $JPP_DIR/setenv.sh $JPP_DIR
6 
7 
8 set_variable DEBUG 1
9 set_variable WORKDIR `pwd`
10 set_variable: FORMAT GRAPHICS_FORMAT gif
11 set_variable+ BATCH GRAPHICS_BATCH -B
12 
13 if do_usage $*; then
14  usage "$script [working directory]"
15 fi
16 
17 if (( $# == 1 )); then
18  set_variable WORKDIR $1
19 fi
20 
21 set_variable N 1000
22 set_variable X "-3.5 +3.5"
23 set_variable Y "-3.5 +3.5"
24 
25 $JPP_DIR/examples/JROOT/JRootfitToGauss2D \
26  -o $WORKDIR/fit.root \
27  -x "$X" \
28  -y "$Y" \
29  -n "$N" \
30  -w \
31  -d 0 >& /dev/null
32 
33 $JPP_DIR/examples/JROOT/JRootfitToGauss2D \
34  -f $WORKDIR/fit.root \
35  -o /dev/null \
36  -x "$X" \
37  -y "$Y" \
38  -d $DEBUG
39 
40 set_variable FORMULA "[4] * exp(-0.5 * (x-[0])*(x-[0]) / ([1]*[1])) * exp(-0.5 * (y-[2])*(y-[2]) / ([3]*[3])) / (2*pi*[1]*[3]) + [5]"
41 
42 JFit2D \
43  -f $WORKDIR/fit.root:h2 \
44  -f $WORKDIR/fit.root:h2 \
45  -o /dev/null \
46  -x "$X" \
47  -y "$Y" \
48  -F "$FORMULA" \
49  -@ "p0 = GetMean(1)" \
50  -@ "p1 = GetStdDev(1) * 0.50" \
51  -@ "p2 = GetMean(2)" \
52  -@ "p3 = GetStdDev(2) * 0.50" \
53  -@ "p4 = GetSumOfWeights() / (2*pi)" \
54  -@ "p5 = GetMinimum" \
55  -O NLQ \
56  -d $DEBUG
57 
58 JPlot2D \
59  -f $WORKDIR/fit.root:h2 \
60  -x "$X" \
61  -y "$Y" \
62  -O COLZ \
63  -T "" \
64  -o f2.$FORMAT $BATCH