Jpp test-rotations-old
the software that should make you happy
Loading...
Searching...
No Matches
JRootfitToGauss2D.sh
Go to the documentation of this file.
1#!/usr/bin/env zsh
2version=1.0
3script=${0##*/}
4
5source $JPP_DIR/setenv.sh $JPP_DIR
6
7
8set_variable DEBUG 1
9set_variable WORKDIR `pwd`
10set_variable: FORMAT GRAPHICS_FORMAT gif
11set_variable+ BATCH GRAPHICS_BATCH -B
12
13if do_usage $*; then
14 usage "$script [working directory]"
15fi
16
17if (( $# == 1 )); then
18 set_variable WORKDIR $1
19fi
20
21set_variable N 1000
22set_variable X "-3.5 +3.5"
23set_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
40set_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
42JFit2D \
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
58JPlot2D \
59 -f $WORKDIR/fit.root:h2 \
60 -x "$X" \
61 -y "$Y" \
62 -O COLZ \
63 -T "" \
64 -o f2.$FORMAT $BATCH