JPseudoExperiment3D.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>"
fi
 
if (( $# != 1 )); then
    fatal "Wrong number of arguments."
fi
 
let "NUMBER_OF_TESTS  =  $argv[1]"
 
# signal PDF
 
let "X0 = 0.0"
let "XS = 0.3"
let "Y0 = 0.0"
let "YS = 0.3"
let "Z0 = 0.0"
let "ZS = 0.3"
 
set_variable FS  "exp(-0.5 * (x-$X0)*(x-$X0) / ($XS*$XS)) * exp(-0.5 * (y-$Y0)*(y-$Y0) / ($YS*$YS)) * exp(-0.5 * (y-$Z0)*(y-$Z0) / ($ZS*$ZS)) / 1560.654"
 
X_AXIS=(25 -1.0 +1.0)
Y_AXIS=(25 -1.0 +1.0)
Z_AXIS=(25 -1.0 +1.0)
 
$JPP_DIR/examples/JGizmo/JHistogram3D \
    -o signal.root             \
    -F "$FS"                   \
    -T "hs"                    \
    -x "$X_AXIS[*]"            \
    -y "$Y_AXIS[*]"            \
    -z "$Z_AXIS[*]"            \
    -d $DEBUG
 
# background PDF
 
let "NB = 3.0"
 
set_variable FB  "$NB / ($X_AXIS[1] * $Y_AXIS[1] * $Z_AXIS[1])"
 
$JPP_DIR/examples/JGizmo/JHistogram3D \
    -o background.root             \
    -F "$FB"                       \
    -T "hb"                        \
    -x "$X_AXIS[*]"                \
    -y "$Y_AXIS[*]"                \
    -z "$Z_AXIS[*]"                \
    -d $DEBUG
 
# pseudo experiments
 
NS=(0 2 4 6)
 
for N in $NS[*]; do
 
    $JPP_DIR/examples/JAstronomy/JPseudoExperiment \
        -o benchmark\[$N\].root    \
        -s signal.root:hs          \
        -b background.root:hb      \
        -n $NUMBER_OF_TESTS        \
        -N $N                      \
        -x "100 -10.0 +50.0"       \
        -d $DEBUG
done
 
JPlot1D \
    -fbenchmark\[${^NS}\].root:hl  \
    -Y                             \
    -> "likelihood ratio"         \
    -\^ "number of events"         \
    -T ""                          \
    -L TR                          \
    -o likelihood.$FORMAT          $BATCH
 
JPlot1D \
    -fbenchmark\[${^NS}\].root:hn  \
    -Y                             \
    -> "#DeltaN_{s}"              \
    -\^ "number of events"         \
    -T ""                          \
    -L TR                          \
    -o ns.$FORMAT                  $BATCH