JMorphology.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)
 
MORPHOLOGY_DEGS=(0.0 0.5)
RESOLUTION_DEGS=(0.0 0.1 0.3)
 
for MORPHOLOGY_DEG in $MORPHOLOGY_DEGS[*]; do
    for RESOLUTION_DEG in $RESOLUTION_DEGS[*]; do
    
        $DIR/JMorphology \
            -o source.$MORPHOLOGY_DEG.\[$RESOLUTION_DEG\].root       \
            -U "$UTC[*]"                                             \
            -n $NUMBER_OF_EVENTS                                     \
            -r "Gauss $RESOLUTION_DEG"                               \
            -M "Gauss $LOCATION_DEG $MORPHOLOGY_DEG"                 \
            -d $DEBUG
 
 
        JFit2D \
            -f source.$MORPHOLOGY_DEG.\[$RESOLUTION_DEG\].root:h2    \
            -o fit.root                                              \
            -F "$FORMULA"                                            \
            -@ "p0 = GetMaximum"                                     \
            -@ "p1 = 0.0"                                            \
            -@ "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 [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 fit.root >& /dev/null
    done
done
 
for MORPHOLOGY_DEG in $MORPHOLOGY_DEGS[*]; do
 
    JPlot1D \
        -fsource.$MORPHOLOGY_DEG.\[${^RESOLUTION_DEGS}\].root:ha   \
        -XX                                                        \
        -Y                                                         \
        -> "angle [deg]"                                          \
        -\^ "number of events [a.u]"                               \
        -L TR                                                      \
        -T "morphology $MORPHOLOGY_DEG [deg]"                      \
        -o angle.$MORPHOLOGY_DEG.$FORMAT                           $BATCH
done