JGandalfFitToGauss.sh

Go to the documentation of this file.

#!/usr/bin/env 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  ${TMPDIR:-/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