JPseudoExperimentUpperLimit.sh

Go to the documentation of this file.

#!/bin/env zsh
#
# \author mdejong
#
version=1.0
script=${0##*/}
 
source $JPP_DIR/setenv.sh $JPP_DIR >& /dev/null
 
set_variable   DEBUG       2
set_variable:  FORMAT      GRAPHICS_FORMAT        gif
set_variable+  BATCH       GRAPHICS_BATCH         -B
 
if do_usage $*; then
    usage "$script <number of tests> <signal> <background>"
fi
 
if (( $# != 3 )); then
    fatal "Wrong number of arguments."
fi
 
let "NUMBER_OF_TESTS  =  $argv[1]"
let "NS               =  $argv[2] * 1.0"
let "NB               =  $argv[3] * 1.0"
 
if (( 1 )); then
 
    X_AXIS=(1 -0.5 +0.5)
 
    set_variable FS  "$NS / $X_AXIS[1]"
    set_variable FB  "$NB / $X_AXIS[1]"
 
else
 
    X_AXIS=(50 -1.0 +1.0)
 
    let "X0 = 0.0"
    let "XS = 0.3"
 
    set_variable FS  "$NS * (exp(-0.5 * (x-$X0)*(x-$X0) / ($XS*$XS)) / (sqrt(2.0*$PI) * $XS)) * ($X_AXIS[3] - $X_AXIS[2]) / $X_AXIS[1]"
    set_variable FB  "$NB / $X_AXIS[1]"
fi
 
$JPP_DIR/examples/JGizmo/JHistogram1D \
    -o signal.root                 \
    -F "$FS"                       \
    -T "hs"                        \
    -x "$X_AXIS[*]"                \
    -d 0
 
$JPP_DIR/examples/JGizmo/JHistogram1D \
    -o background.root             \
    -F "$FB"                       \
    -T "hb"                        \
    -x "$X_AXIS[*]"                \
    -d 0
 
HS=signal.root:hs
HB=background.root:hb
 
timer_start
 
$JPP_DIR/examples/JAstronomy/JPseudoExperimentUpperLimit \
        -o benchmark.root          \
        -E "$HS $HB"               \
        -n $NUMBER_OF_TESTS        \
        -x "500 -1.0e1 +1.0e1"     \
        -d $DEBUG
 
timer_stop
timer_print
 
JPlot1D \
    -f benchmark.root:h0           \
    -f benchmark.root:h1           \
    -> "signal"                   \
    -T ""                          \
    -L TR                          \
    -o likelihood.$FORMAT          $BATCH
 
JPlot1D \
    -f $HB                         \
    -f $HS                         \
    -T ""                          \
    -o benchmark.$FORMAT           $BATCH