JRootfitToGauss2D.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`
set_variable:  FORMAT      GRAPHICS_FORMAT        gif
set_variable+  BATCH       GRAPHICS_BATCH         -B

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

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

set_variable  N        1000
set_variable  X        "-3.5 +3.5"
set_variable  Y        "-3.5 +3.5"

$JPP_DIR/examples/JROOT/JRootfitToGauss2D \
    -o $WORKDIR/fit.root                   \
    -x "$X"                                \
    -y "$Y"                                \
    -n "$N"                                \
    -w                                     \
    -d 0                                   >& /dev/null

$JPP_DIR/examples/JROOT/JRootfitToGauss2D \
    -f $WORKDIR/fit.root                   \
    -o /dev/null                           \
    -x "$X"                                \
    -y "$Y"                                \
    -d $DEBUG

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]"

JFit2D \
    -f $WORKDIR/fit.root:h2                \
    -f $WORKDIR/fit.root:h2                \
    -o /dev/null                           \
    -x "$X"                                \
    -y "$Y"                                \
    -F "$FORMULA"                          \
    -@ "p0 = GetMean(1)"                   \
    -@ "p1 = GetStdDev(1) * 0.50"          \
    -@ "p2 = GetMean(2)"                   \
    -@ "p3 = GetStdDev(2) * 0.50"          \
    -@ "p4 = GetSumOfWeights() / (2*pi)"   \
    -@ "p5 = GetMinimum"                   \
    -O NLQ                                 \
    -d $DEBUG

JPlot2D \
    -f $WORKDIR/fit.root:h2                \
    -x "$X"                                \
    -y "$Y"                                \
    -O COLZ                                \
    -T ""                                  \
    -o f2.$FORMAT                          $BATCH