JRootfitToGauss.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        100
set_variable  X        "-3.5 +3.5"

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

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

set_variable FORMULA "[0] * exp(-0.5 * (x-[1])*(x-[1]) / ([2]*[2])) + [3]"

JFit \
    -f $WORKDIR/fit.root:h1                \
    -f $WORKDIR/fit.root:h1                \
    -o /dev/null                           \
    -x "$X"                                \
    -F "$FORMULA"                          \
    -@ "p0 = GetMaximum()"                 \
    -@ "p1 = GetMean()"                    \
    -@ "p2 = GetStdDev() * 0.66"           \
    -@ "p3 = GetMinimum  + 0.10"           \
    -O NLQ                                 \

JPlot1D \
    -f $WORKDIR/fit.root:h1                \
    -T ""                                  \
    -o f1.$FORMAT                          $BATCH