upperlimit_mset.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 (( $# != 2 )); then
    fatal "Wrong number of arguments."
fi
 
set_variable  CONFIG      $argv[1]
let  "NUMBER_OF_TESTS  =  $argv[2]"
 
# pseudo experiments
 
if (( $NUMBER_OF_TESTS > 0 )); then
 
    timer_start
 
    $JPP_DIR/examples/JAstronomy/JGen2UpperLimit \
        -o benchmark.root        \
        -E $CONFIG               \
        -n $NUMBER_OF_TESTS      \
        -M $M_SIZE               \
        -R $SNR                  \
        -x "500 -1.0e2 +2.0e1"   \
        -S 0                     \
        -d $DEBUG                >& pseudo-experiment-limit.log
    timer_stop
    timer_print
 
    awk '/<mu>/ { print $(NF) }' pseudo-experiment-limit.log | read UPPER_LIMIT
 
    SIGMAREL=`echo "sqrt( (1./$NUMBER_OF_TESTS) + 1./($NUMBER_OF_TESTS*5.)*0.9/0.1 )" | bc -l` 
    SIGMA=`echo "$UPPER_LIMIT*$SIGMAREL" | bc -l`  
 
    printf "<mu> %1.2e +- %1.2e\n" $UPPER_LIMIT $SIGMA
fi