JRootfitToGauss3D.sh

Go to the documentation of this file.

#!/bin/zsh
version=1.0
script=${0##*/}

source $JPP_DIR/setenv.sh $JPP_DIR


set_variable   DEBUG    1
set_variable   WORKDIR  `pwd`

if do_usage $*; then
    usage "$script [working directory]"
fi

if (( $# == 1 )); then
    set_variable  WORKDIR   $1
fi

set_variable  X  "-4.0 +4.0"
set_variable  Y  "-4.0 +4.0"
set_variable  Z  "-4.0 +4.0"
set_variable  N  4000000

if (( $N < 100000 )); then
    OPTION=count
else
    OPTION=value
fi

$JPP_DIR/examples/JROOT/JRootfitToGauss3D \
    -o $WORKDIR/fit.root                   \
    -x "$X"                                \
    -y "$Y"                                \
    -z "$Z"                                \
    -n "$N"                                \
    -O $OPTION                             \
    -d $DEBUG --!

set_variable FORMULA "[6] * "                     \
    "exp(-0.5 * (x-[0])*(x-[0]) / ([1]*[1])) * "  \
    "exp(-0.5 * (y-[2])*(y-[2]) / ([3]*[3])) * "  \
    "exp(-0.5 * (z-[4])*(z-[4]) / ([5]*[5])) / "  \
    "(2*pi*sqrt(2*pi)*[1]*[3]*[5])  +  [7]"

JFit3D \
    -f $WORKDIR/fit.root:h3                \
    -f $WORKDIR/fit.root:h3                \
    -o /dev/null                           \
    -x "$X"                                \
    -y "$Y"                                \
    -z "$Z"                                \
    -F "$FORMULA"                          \
    -@ "p0 = GetMean(1)"                   \
    -@ "p1 = GetStdDev(1) * 0.33"          \
    -@ "p2 = GetMean(2)"                   \
    -@ "p3 = GetStdDev(2) * 0.33"          \
    -@ "p4 = GetMean(3)"                   \
    -@ "p5 = GetStdDev(3) * 0.33"          \
    -@ "p6 = GetSumOfWeights() / (sqrt(2*pi)*2*pi)"   \
    -@ "p7 = GetMinimum"                   \
    -O NLQ                                 \
    -d $DEBUG