JMorphologyH2.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       +45.0 -90.0
set_variable   NUMBER_OF_EVENTS   100000
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.5
RESOLUTION_DEG=0.2
 
set_variable F1      "[0] * exp(-0.5 * (x-[1])*(x-[1]) / ([2]*[2]))"
set_variable F2      "[0] * exp(-0.5 * (x-[1])*(x-[1]) / ([2]*[2])) * exp(-0.5 * (y-[3])*(y-[3]) / ([4]*[4]))"
 
$JPP_DIR/examples/JGizmo/JHistogram1D \
    -o resolution.root                                       \
    -T "user"                                                \
    -F "$F1"                                                 \
    -@ "p0 = 1.0; p1 = 0.0; p2 = $RESOLUTION_DEG"            \
    -x "300 0.0 +1.5"                                        \
    -d $DEBUG
 
$JPP_DIR/examples/JGizmo/JHistogram2D \
    -o histogram.root                                        \
    -T "user"                                                \
    -F "$F2"                                                 \
    -@ "p0 = 1.0; p1 = 0.0; p3 = 0.0"                        \
    -@ "p2 = $(($MORPHOLOGY_DEG / sqrt(2.0)))"               \
    -@ "p4 = $(($MORPHOLOGY_DEG / sqrt(2.0)))"               \
    -x "300 -1.5 +1.5"                                       \
    -y "300 -1.5 +1.5"                                       \
    -d $DEBUG
 
$DIR/JMorphology \
    -o source.histogram.root                                 \
    -U "$UTC[*]"                                             \
    -n $NUMBER_OF_EVENTS                                     \
    -r "Histogram resolution.root user"                      \
    -M "Histogram $LOCATION_DEG histogram.root user"         \
    -d $DEBUG --!
 
JFit2D \
    -f source.histogram.root:h2                              \
    -o fit.root                                              \
    -F "$F2"                                                 \
    -@ "p0 = GetMaximum"                                     \
    -@ "p1 = GetMean(1)"                                     \
    -@ "p2 = 0.7*GetRMS(1)"                                  \
    -@ "p3 = GetMean(2)"                                     \
    -@ "p4 = 0.7*GetRMS(2)"                                  \
    -O L                                                     \
    -d $DEBUG                                                >& /dev/null
 
SIGMA_X=($(JPrintFit -f fit.root:h2 -@ p2))
SIGMA_Y=($(JPrintFit -f fit.root:h2 -@ p4))
 
printf "morphology %6.3f [deg] resolution %6.3f [deg] => %6.3f [deg]\n"  $MORPHOLOGY_DEG  $RESOLUTION_DEG  $((sqrt($MORPHOLOGY_DEG*$MORPHOLOGY_DEG + $RESOLUTION_DEG*$RESOLUTION_DEG)/sqrt(2.0)))
printf "fit X-sigma = %6.3f [deg]\n" $SIGMA_X[1]
printf "fit Y-sigma = %6.3f [deg]\n" $SIGMA_Y[1]
 
rm -f {resolution,histogram,fit}.root >& /dev/null
 
JPlot2D \
    -f source.histogram.root:h2                              \
    -> "#DeltaRA [deg]"                                     \
    -< "#DeltaDEC [deg]"                                    \
    -G XY                                                    \
    -O COLZ                                                  \
    -N "X 505"                                               \
    -N "Y 505"                                               \
    -T ""                                                    \
    -o source.histogram.$FORMAT                              $BATCH