JMorphologyBn.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)
 
RESOLUTION_DEG=0.2
 
X1=+1.5
Y1=-1.0
S1=0.2
W1=0.2
 
X2=-1.0
Y2=+0.5
S2=0.3
W2=0.8
 
$DIR/JMorphology \
    -o source.binary.root                                      \
    -U "$UTC[*]"                                               \
    -n $NUMBER_OF_EVENTS                                       \
    -r "Gauss $RESOLUTION_DEG"                                 \
    -M "Binary $LOCATION_DEG $X1 $Y1 $W1 $S1 $X2 $Y2 $W2 $S2"  \
    -d $DEBUG --!
 
JPlot2D \
    -f source.binary.root:h2                                   \
    -> "#DeltaRA [deg]"                                       \
    -< "#DeltaDEC [deg]"                                      \
    -G XY                                                      \
    -O COLZ                                                    \
    -N "X 505"                                                 \
    -N "Y 505"                                                 \
    -T "source ($X1,$Y1) ($X2,$Y2) [deg]"                      \
    -o source.binary.$FORMAT                                   $BATCH