JNuisance.sh

Go to the documentation of this file.

#!/bin/env zsh
script=${0##*/}
 
if [ -z $JPP_DIR ]; then
    echo "Variable JPP_DIR undefined."
    exit
fi
 
source $JPP_DIR/setenv.sh $JPP_DIR
 
set_variable   DIR                $JPP_DIR/examples/JAstronomy
set_variable   DEBUG              2
set_variable   NUMBER_OF_EVENTS   100000
set_variable:  FORMAT             GRAPHICS_FORMAT       gif
set_variable+  BATCH              GRAPHICS_BATCH        -B
 
if do_usage $*; then
    usage "$script"
fi
 
set_variable   FORMULA           "[0] * exp(-0.5 * (x-[1])*(x-[1]) / ([2]*[2]))"
set_variable   SIGMA             0.3
 
$DIR/JNuisance \
        -o nuisance.root           \
        -n $NUMBER_OF_EVENTS       \
        -u "Gauss $SIGMA"          \
        -d $DEBUG
 
JFit \
        -f nuisance.root:h0        \
    -o fit.root                \
    -F "$FORMULA"              \
    -@ "p0 = GetMaximum"       \
    -@ "p1 = 1.0"              \
    -@ "p2 = 0.7*GetRMS"       \
        -O L                       \
    -d $DEBUG                  >& /dev/null
 
printf "nuisance %6.3f => %6.3f +/- %6.3f\n"  $SIGMA   $(JPrintFit -f fit.root:h0 -@ p2)
 
rm -f nuisance.root fit.root