upper-limit.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:  SNR         ASTRONOMY_SNR          0.0
set_variable:  M_SIZE      ASTRONOMY_M_SIZE       10000
set_variable:  BOOST       ASTRONOMY_BOOST        1.0
set_variable:  DEBUG       ASTRONOMY_DEBUG        1
set_variable:  FORMAT      GRAPHICS_FORMAT        gif
set_variable+  BATCH       GRAPHICS_BATCH         -B
 
if do_usage $*; then
    usage "$script <signal histogram> <background histogram> <number of tests>"\
          "\nThe histograms correspond to \"<file name>:<histogram name>\"."
fi
 
if (( $# != 3 )); then
    fatal "Wrong number of arguments."
fi
 
set_variable  HS          $argv[1]
set_variable  HB          $argv[2]
let  "NUMBER_OF_TESTS  =  $argv[3]"
 
# pseudo experiments
 
if (( $NUMBER_OF_TESTS > 0 )); then
 
    timer_start
 
    $JPP_DIR/examples/JAstronomy/JPseudoExperimentUpperLimit \
        -o benchmark.root        \
        -E "$HS $HB"             \
        -n $NUMBER_OF_TESTS      \
        -B $BOOST                \
        -M $M_SIZE               \
        -R $SNR                  \
        -x "500 -1.0e2 +2.0e1"   \
        -S 0                     \
        -d $DEBUG                #>& pseudo-experiment-limit.log &
 
    timer_stop
    timer_print
fi