JGandalfFitToGauss.sh

Go to the documentation of this file.

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

# ------------------------------------------------------------------------------------------
#
#                         Example script to test JGandalf.
#
# ------------------------------------------------------------------------------------------

if [ -z $JPP_DIR ]; then
    echo "Variable JPP_DIR undefined."
    exit
fi

source $JPP_DIR/setenv.sh $JPP_DIR

set_variable     DEBUG           1
set_variable     WORKDIR         /tmp/

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

case $# in
    1) set_variable WORKDIR $1;;
esac

set_variable  DIR                $JPP_DIR/examples/JFit/
set_variable  OUTPUT_FILE        $WORKDIR/gandalf.root
set_variable  NUMBER_OF_EVENTS   1000
set_variable  GAUSS              "0.0  1.0  1000.0 100.0"
set_variable  PRECISION          "0.05 0.05   20.0  20.0"

$DIR/JGandalfFitToGauss \
    -o $OUTPUT_FILE            \
    -n $NUMBER_OF_EVENTS       \
    -@ "$GAUSS"                \
    -e "$PRECISION"            \
    -d $DEBUG --!

typeset -A X_LABEL

X_LABEL[ha]=mean
X_LABEL[hb]=sigma
X_LABEL[hc]=signal
X_LABEL[hd]=background

for H1 in ha hb hc hd; do
    JPlot1D \
        -f ${OUTPUT_FILE}:${H1}         \
        -\> "$X_LABEL[${H1}]"           \
        -\^ "number of events [a.u]"    \
        -s 111111 -T ""
done