JMorphology+background.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
 
zmodload zsh/mathfunc
 
set_variable   DIR                $JPP_DIR/examples/JAstronomy
set_variable   DEBUG              2
set_variable   LOCATION_DEG       -39.77 -101.72
set_variable   NUMBER_OF_EVENTS   35
set_variable   NUMBER_OF_SIGMAS   2.0
set_variable   BACKGROUND         50000                 # number of events per year
set_variable   OMEGA              2.0e-5                # solid angle
set_variable:  FORMAT             GRAPHICS_FORMAT       gif
set_variable+  BATCH              GRAPHICS_BATCH        -B
 
if do_usage $*; then
    usage "$script"
fi
 
UTC=(2025-01-01T00:00:00Z 2035-01-01T00:00:00Z)
 
MORPHOLOGY_DEG=0.0
RESOLUTION_DEG=0.4
 
let "R = (1.0 * $BACKGROUND) / (365.0 * 24.0 * 60.0 * 60.0) / (2.0 * $PI)"
 
XMIN=-3.5
XMAX=+3.5
YMIN=-3.5
YMAX=+3.5
 
$DIR/JMorphology \
    -@ "ct = -0.1 +1.0"                                        \
    -o source+background.root                                  \
    -U "$UTC[*]"                                               \
    -n $NUMBER_OF_EVENTS                                       \
    -r "Gauss $RESOLUTION_DEG"                                 \
    -x "200 $XMIN $XMAX"                                       \
    -y "200 $YMIN $YMAX"                                       \
    -M "Point $LOCATION_DEG"                                   \
    -B $R                                                      \
    -O $OMEGA                                                  \
    -d $DEBUG --!
 
JPlot2D \
    -f source+background.root:g1                               \
    -f source+background.root:g2                               \
    -> "#DeltaRA [deg]"                                       \
    -< "#DeltaDEC [deg]"                                      \
    -x "$XMIN $XMAX"                                           \
    -y "$YMIN $YMAX"                                           \
    -G XY                                                      \
    -O COLZ                                                    \
    -N "X 505"                                                 \
    -N "Y 505"                                                 \
    -S 0.4                                                     \
    -T "$NUMBER_OF_EVENTS"                                     \
    -o skymap+background.$FORMAT                               $BATCH
 
JPrintRange1D -f source+background.root:h1 | read XMIN YMIN XMAX YMAX
 
let "XS = cos($NUMBER_OF_SIGMAS * sqrt($RESOLUTION_DEG*$RESOLUTION_DEG + $MORPHOLOGY_DEG*$MORPHOLOGY_DEG) * $PI / 180.0)"
 
JFit \
    -f source+background.root:h1                               \
    -o fit.root                                                \
    -x "$XMIN $XS"                                             \
    -F "exp([0])"                                              \
    -@ "p0 = $((log($YMIN + 0.5)))"                            \
    -O L                                                       \
    -d $DEBUG 
 
JPrintFit -f fit.root:\.\* -@ p0 | read -A BUFFER
 
printf "Fitted background %7.3f +/- %7.3f\n" $((exp($BUFFER[1]))) $((exp($BUFFER[1]) * $BUFFER[2]))
 
JPlot1D \
    -f fit.root:h1                                             \
    -> "cos(#alpha)"                                          \
    -\^ "number of events"                                     \
    -N "X 505"                                                 \
    -T "$NUMBER_OF_EVENTS"                                     \
    -o source+background.$FORMAT                               $BATCH
 
rm -f fit.root >& /dev/null