JMorphologyHn.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   1000000
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)
 
RESOLUTION_DEG=0.2
 
set_variable FORMULA      "[0] * exp(-0.5 * (x-[1])*(x-[1]) / ([3]*[3])) * exp(-0.5 * (y-[2])*(y-[2]) / ([4]*[4]))"
 
$JPP_DIR/examples/JGizmo/JHistogram2D \
    -o ha.root                                               \
    -T "user"                                                \
    -F "$FORMULA"                                            \
    -@ "p0 = 1.0; p1 = +0.1; p2 = +0.5; p3 = 0.3; p4 = 0.3"  \
    -x "300 -1.5 +1.5"                                       \
    -y "300 -1.5 +1.5"                                       \
    -d $DEBUG
 
$JPP_DIR/examples/JGizmo/JHistogram2D \
    -o hb.root                                               \
    -T "user"                                                \
    -F "$FORMULA"                                            \
    -@ "p0 = 3.0; p1 = -0.5; p2 = -0.2; p3 = 0.2; p4 = 0.2"  \
    -x "300 -1.5 +1.5"                                       \
    -y "300 -1.5 +1.5"                                       \
    -d $DEBUG
 
$JPP_DIR/examples/JGizmo/JHistogram2D \
    -o hc.root                                               \
    -T "user"                                                \
    -F "$FORMULA"                                            \
    -@ "p0 = 1.5; p1 = +0.7; p2 = -0.3; p3 = 0.6; p4 = 0.2"  \
    -x "300 -1.5 +1.5"                                       \
    -y "300 -1.5 +1.5"                                       \
    -d $DEBUG
 
JOpera2D \
    -fh{a,b,c}.root:user                                     \
    -o histogram.root                                        \
    -u Add                                                   \
    -O user                                                  
    
$DIR/JMorphology \
    -o source.histogram.root                                 \
    -U "$UTC[*]"                                             \
    -n $NUMBER_OF_EVENTS                                     \
    -r "Gauss $RESOLUTION_DEG"                               \
    -M "Histogram $LOCATION_DEG histogram.root user"         \
    -x "300 -3.5 +3.5"                                       \
    -y "300 -3.5 +3.5"                                       \
    -d $DEBUG --!
 
JPlot2D \
    -f source.histogram.root:h2                              \
    -> "#DeltaRA [deg]"                                     \
    -< "#DeltaDEC [deg]"                                    \
    -G XY                                                    \
    -O COLZ                                                  \
    -N "X 505"                                               \
    -N "Y 505"                                               \
    -T ""                                                    \
    -p 57                                                    \
    -o source.histogram.$FORMAT                              $BATCH
 
rm -f h{a,b,c}.root   >& /dev/null
rm -f histogram.root  >& /dev/null