examples/JAcoustics/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
 
zmodload zsh/mathfunc
 
set_variable   DIR                $JPP_DIR/examples/JAcoustics/
set_variable   WORKDIR            ./
set_variable:  DEBUG              ACOUSTICS_DEBUG              2
set_variable+  HYDROPHONE         ACOUSTICS_HYDROPHONE         -H hydrophone.txt
set_variable+  TRANSMITTER        ACOUSTICS_TRANSMITTER        -Y transmitter.txt
set_variable:  SIGMA_S            ACOUSTICS_SIGMA_S            50.0E-6
set_variable:  NUMBER_OF_EVENTS   ACOUSTICS_NUMBER_OF_EVENTS   1000
set_variable:  FORMAT             GRAPHICS_FORMAT              gif
set_variable+  BATCH              GRAPHICS_BATCH               -B
 
if do_usage $*; then
    usage "$script <detector file> <tripod file>"
fi
 
if (( $# != 2 )); then
    fatal "Wrong number of arguments."
fi
 
set_variable  DETECTOR            $argv[1]
set_variable  TRIPOD              $argv[2]
 
if [[ ! -f $DETECTOR ]]; then
    fatal "No file $DETECTOR."
fi
 
if [[ ! -f $TRIPOD ]]; then
    fatal "No file $TRIPOD."
fi
 
eval `JPrintDetector -a $DETECTOR -O IDENTIFIER`
eval `JPrintDetector -a $DETECTOR -O SUMMARY`
 
if (( 1 )); then
 
    JDrawDetector2D \
        -a $DETECTOR                              \
        -T $TRIPOD                                \
        -L TL                                     \
        -o footprint.$FORMAT                      $BATCH
fi
 
if ( ! reuse_file $WORKDIR/katoomba.root ); then
 
    $DIR/JKatoomba \
        -a $DETECTOR                              \
        -o $WORKDIR/katoomba.root                 \
        -T $TRIPOD                                \
        ${TRANSMITTER}                            \
        -n $NUMBER_OF_EVENTS                      \
        -@ acoustics_fit_parameters.txt           \
        -@ "sigma_s    = $SIGMA_S;"               \
        -M mechanics.txt                          \
        ${HYDROPHONE}                             \
        -S 12345678                               \
        -d $DEBUG --!
fi
 
if (( 1 )); then
 
    rm -f $WORKDIR/strings.txt
 
    for STRING in ${STRINGS[*]}; do
 
        RESULT=(`JPrintResult -f "$WORKDIR/katoomba.root:string\[$STRING\]" -F "GetRMS(1)" -F "GetRMS(2)"`)
 
        printf "%4d  %9.5f  %9.5f [mrad]\n"  $STRING  $RESULT[*]
 
        printf "%9.5f\n"  $((sqrt($RESULT[1]*$RESULT[2] + $RESULT[1]*$RESULT[2]))) >> $WORKDIR/strings.txt
    done
 
    $JPP_DIR/examples/JGizmo/JHistogram1D \
        -x "50 0 +2.5"                                            \
        -f $WORKDIR/strings.txt                                   \
        -o $WORKDIR/strings.root                                  
 
    JPlot1D \
        -f $WORKDIR/strings.root:\.\*                             \
        -> "#sigma_{T} [mrad]"                                   \
        -\^ "number of strings"                                   \
        -T ""                                                     \
        -o strings.$FORMAT                                        $BATCH
    
    rm -f $WORKDIR/strings.{txt,root}
fi
 
if (( 1 )); then
 
    JPlot1D \
        -f "$WORKDIR/katoomba.root:cpu"                           \
        -> "cpu [ms]"                                            \
        -\^ "number of events [a.u.]"                             \
        -T "$FIT"                                                 \
        -o cpu.$FORMAT                                            $BATCH
 
    JPlot1D \
        -f "$WORKDIR/katoomba.root:chi2"                          \
        -x "0 2"                                                  \
        -> "chi2/NDF"                                            \
        -\^ "number of events [a.u.]"                             \
        -s 110                                                    \
        -T "$FIT"                                                 \
        -o chi2.$FORMAT                                           $BATCH
fi
 
if (( 1 )); then
 
    source JAcousticsToolkit.sh
 
    typeset -A EMITTERS
 
    get_tripods+       tripod.txt       EMITTERS
    get_transmitters+  transmitter.txt  EMITTERS
 
    for EMITTER in % ${(@k)EMITTERS}; do
 
        JPlot1D \
            -f "$WORKDIR/katoomba.root:emitter\[${EMITTER}\]"         \
            -> "#Deltat [s]"                                         \
            -\^ "number of events [a.u.]"                             \
            -N "X 505"                                                \
            -s 111010                                                 \
            -T "[$EMITTER]"                                           \
            -o TX${EMITTER}.$FORMAT                                   $BATCH
    done
fi
 
if (( 1 )); then
 
    JEllipse \
        -o $WORKDIR/ellipse.root                                      \
        -p "0.0 0.0 1.0"                                              \
        -@ "line.color = 2"                                           \
        -@ "line.width = 2"
 
    for STRING in % $STRINGS[*]; do
 
        JPlot2D \
            -f "$WORKDIR/katoomba.root:string\[${STRING}\]"           \
            -f  $WORKDIR/ellipse.root:\.\*                            \
            -O COLZ                                                   \
            -> "#DeltaTx [mrad]"                                     \
            -< "#DeltaTy [mrad]"                                     \
            -N "X 505"                                                \
            -N "Y 505"                                                \
            -p 57                                                     \
            -s 101010                                                 \
            -T "[$STRING]"                                            \
            -o RX${STRING}.$FORMAT                                    $BATCH
    done
fi