JMorphology2D.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
 
set_variable   FORMULA           "[0] * exp(-0.5 * (x-[1])*(x-[1]) / ([2]*[2])) * exp(-0.5 * (y-[1])*(y-[1]) / ([3]*[3]))"
 
UTC=(2025-01-01T00:00:00Z 2035-01-01T00:00:00Z)
 
MORPH_X_DEGS=(0.0 0.5)
MORPH_Y_DEGS=(0.0 0.5)
RESOLUTION_DEG=0.3
 
for MORPH_X_DEG in $MORPH_X_DEGS[*]; do
    for MORPH_Y_DEG in $MORPH_X_DEGS[*]; do
    
        $DIR/JMorphology \
            -o source.\[$MORPH_X_DEG.$MORPH_Y_DEG\].root                  \
            -U "$UTC[*]"                                                  \
            -n $NUMBER_OF_EVENTS                                          \
            -r "Formula exp(-0.5*x*x/($RESOLUTION_DEG*$RESOLUTION_DEG))"  \
            -M "Gauss2D $LOCATION_DEG $MORPH_X_DEG $MORPH_Y_DEG"          \
            -d $DEBUG
 
        JFit2D \
            -f source.\[$MORPH_X_DEG.$MORPH_Y_DEG\].root:h2          \
            -o fit.root                                              \
            -F "$FORMULA"                                            \
            -@ "p0 = GetMaximum"                                     \
            -@ "p1 = GetMean"                                        \
            -@ "p2 = 0.7*GetRMS(1)"                                  \
            -@ "p3 = 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 -@ p3))
 
        printf "morphology (%6.3f,%6.3f) [deg] resolution %6.3f [deg] =>\n"  $MORPH_X_DEG  $MORPH_Y_DEG  $RESOLUTION_DEG
        printf "(%6.3f,%6.3f) [deg]\n"  $((sqrt($MORPH_X_DEG*$MORPH_X_DEG + 0.5*$RESOLUTION_DEG*$RESOLUTION_DEG)))  $((sqrt($MORPH_Y_DEG*$MORPH_Y_DEG + 0.5*$RESOLUTION_DEG*$RESOLUTION_DEG)))
        printf "fit X-sigma = %6.3f [deg]\n" $SIGMA_X[1]
        printf "fit Y-sigma = %6.3f [deg]\n" $SIGMA_Y[1]
 
        rm -f fit.root >& /dev/null
    done
done
 
for MORPH_X_DEG in $MORPH_X_DEGS[*]; do
    for MORPH_Y_DEG in $MORPH_X_DEGS[*]; do
 
        JPlot2D \
            -f source.\[$MORPH_X_DEG.$MORPH_Y_DEG\].root:h2          \
            -> "#DeltaRA [deg]"                                     \
            -< "#DeltaDEC [deg]"                                    \
            -O COLZ                                                  \
            -N "X 505"                                               \
            -N "Y 505"                                               \
            -s 11000                                                 \
            -T "source ($MORPH_X_DEG,$MORPH_Y_DEG) [deg]"            \
            -o source.\[$MORPH_X_DEG.$MORPH_Y_DEG\].$FORMAT          $BATCH
    done
done