5 source $JPP_DIR/setenv.sh $JPP_DIR
9 set_variable WORKDIR `pwd`
12 usage "$script [working directory]"
15 if (( $# == 1 )); then
16 set_variable WORKDIR $1
19 set_variable X "-4.0 +4.0"
20 set_variable Y "-4.0 +4.0"
21 set_variable Z "-4.0 +4.0"
22 set_variable N 4000000
24 if (( $N < 100000 )); then
30 $JPP_DIR/examples/JROOT/JRootfitToGauss3D \
31 -o $WORKDIR/fit.root \
39 set_variable FORMULA "[6] * " \
40 "exp(-0.5 * (x-[0])*(x-[0]) / ([1]*[1])) * " \
41 "exp(-0.5 * (y-[2])*(y-[2]) / ([3]*[3])) * " \
42 "exp(-0.5 * (z-[4])*(z-[4]) / ([5]*[5])) / " \
43 "(2*pi*sqrt(2*pi)*[1]*[3]*[5]) + [7]"
46 -f $WORKDIR/fit.root:h3 \
47 -f $WORKDIR/fit.root:h3 \
53 -@ "p0 = GetMean(1)" \
54 -@ "p1 = GetStdDev(1) * 0.33" \
55 -@ "p2 = GetMean(2)" \
56 -@ "p3 = GetStdDev(2) * 0.33" \
57 -@ "p4 = GetMean(3)" \
58 -@ "p5 = GetStdDev(3) * 0.33" \
59 -@ "p6 = GetSumOfWeights() / (sqrt(2*pi)*2*pi)" \
60 -@ "p7 = GetMinimum" \