JKatoomba.sh

Go to the documentation of this file.

#!/bin/zsh
script=${0##*/}
 
if [ -z $JPP_DIR ]; then
    echo "Variable JPP_DIR undefined."
    exit
fi

source $JPP_DIR/setenv.sh $JPP_DIR >& /dev/null

set_variable  DEBUG               2
set_variable  DIR                 /user/mjg/Jpp/JAcoustics.examples/
set_variable  WORKDIR             /tmp/
set_variable  BACKGROUND          #-B

if ( do_usage $* ); then
    usage "$script"
fi

function clean()
{
    rm -f $WORKDIR/detector.detx
    rm -f $WORKDIR/tripod.txt
    rm -f $WORKDIR/tripod.txt
    rm -f $WORKDIR/cpu.gif
    rm -f $WORKDIR/chi2.gif
    rm -f $WORKDIR/TX*.gif
    rm -f $WORKDIR/RX*.gif
}

if [[ ! -f $WORKDIR/detector.detx ]]; then

    set_variable numberOfLines          2
    set_variable numberOfFloors        18
    set_variable floorDistance_m       36
    set_variable heightOfFirstFloor_m  72
    set_variable detectorType           1

    JDetector \
        -D 1000                                  \
        -V V3                                    \
        -@ "detectorType    = $detectorType"     \
        -@ "numberOfLines   = $numberOfLines"    \
        -@ "numberOfFloors  = $numberOfFloors"   \
        -@ "floorDistance_m = $floorDistance_m"  \
        -@ "lineDistance_m  = $lineDistance_m"   \
        -z $heightOfFirstFloor_m                 \
        -f footprint.txt                         \
        -o $WORKDIR/detector.detx -d 1 --!

fi

if [[ ! -f $WORKDIR/tripod.txt ]]; then

    cat>$WORKDIR/tripod.txt<<EOF
# identifier UTM position (x,y,z)
#
  12    587888 4017079 -3450.0
  14    587436 4016992 -3450.0
  16    587843 4016727 -3450.0
EOF
fi

if (( 1 )); then

    JDrawDetector2D \
        -a $WORKDIR/detector.detx                 \
        -T $WORKDIR/tripod.txt                    \
        -L                                        \
        -o $WORKDIR/footprint.gif       $BACKGROUND
fi

if ( ! reuse_file $WORKDIR/katoomba.root ); then

    $DIR/JKatoomba \
        -a $WORKDIR/detector.detx                 \
        -o $WORKDIR/katoomba.root                 \
        -T $WORKDIR/tripod.txt                    \
        -n 1000                                   \
        -N 10                                     \
        -M 0                                      \
        -S 12345678                               \
        -d $DEBUG --!
fi

if (( 1 )); then

    JPlot1D \
        -f "$WORKDIR/katoomba.root:cpu"                           \
        -XX                                                       \
        -\> "cpu [#mus]"                                          \
        -\^ "number of events [a.u.]"                             \
        -T "$FIT"                                                 \
        -o $WORKDIR/cpu.gif                             $BACKGROUND

    JPlot1D \
        -f "$WORKDIR/katoomba.root:chi2"                          \
        -\> "chi2/NDF"                                            \
        -\^ "number of events [a.u.]"                             \
        -s 110                                                    \
        -T "$FIT"                                                 \
        -o $WORKDIR/chi2.gif                            $BACKGROUND
fi

if (( 1 )); then

    for EMITTER in 12 14 16; do

        JPlot1D \
            -f "$WORKDIR/katoomba.root:emitter\[${EMITTER}\]"         \
            -\> "#Deltat [s]"                                         \
            -\^ "number of events [a.u.]"                             \
            -N "X 505"                                                \
            -s 111010                                                 \
            -T "?"                                                    \
            -o $WORKDIR/TX${EMITTER}.gif                    $BACKGROUND
    done
fi

if (( 1 )); then

    for STRING in 1 2; do

        JPlot2D \
            -f "$WORKDIR/katoomba.root:string\[${STRING}\]"           \
            -O COLZ                                                   \
            -\> "Tx"                                                  \
            -\< "Ty"                                                  \
            -N "X 505"                                                \
            -N "Y 505"                                                \
            -s 101010                                                 \
            -T "?"                                                    \
            -o $WORKDIR/RX${STRING}.gif                     $BACKGROUND
    done
fi