JSoundVelocity.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             $JPP_DIR/examples/JAcoustics/  
set_variable   WORKDIR         ./                             
set_variable   QUANTILE        0.0                            
set_variable   TMAX_S          0.1                            
set_variable   NUMBER_OF_HITS  15                             

set_variable   VELOCITY_TXT    $WORKDIR/velocity.txt          
set_variable   VELOCITY_ROOT   $WORKDIR/velocity.root         

if do_usage $*; then
    usage "$script <detector> <run> [option]"\
          "\nwhere <option> can be clean."
fi

case $# in
    3)  set_variable OPTION       $argv[3];&
    2)  set_variable RUN          $argv[2];
        set_variable DETECTOR     $argv[1];;
    1)  set_variable OPTION       $argv[1];;
    *)  fatal "Wrong number of arguments.";;
esac

function clean()
{
    rm -f $WORKDIR/detector.detx
    rm -f $WORKDIR/tripod.txt
    rm -f $WORKDIR/base.txt
    rm -f $WORKDIR/toashort.root
    rm -f $WORKDIR/event.root
    rm -f $VELOCITY_TXT
    rm -f $VELOCITY_ROOT
    rm -f $WORKDIR/velocity.gif
}

case "$OPTION" in
    clean) clean;;
    "")    ;;
    *)     fatal "Wrong option $OPTION.";;
esac

JCookie.sh

if [[ -n "$DETECTOR" && -n "$RUN" ]]; then

# create input files

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

        JDetectorDB \
            -D $DETECTOR                               \
            -r $RUN                                    \
            -o $WORKDIR/detector.detx --!
    fi

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

        cat>$WORKDIR/tripod.txt<<EOF
# some comment
  12    587888 4017079 -3450.0
  14    587436 4016992 -3450.0
  16    587843 4016727 -3450.0
EOF
    fi

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

        $DIR/JBaseModuleWriter \
            -a $WORKDIR/detector.detx                 \
            -o $WORKDIR/base.txt
    fi

    if [[ ! -f $WORKDIR/toashort.root ]]; then

# create intemediate ROOT file for fast access

        JConvertDB \
            -q toashort                                \
            -@ "detid  = `getDetector -D $DETECTOR`"   \
            -@ "minrun = $RUN"                         \
            -@ "maxrun = $RUN"                         \
            -o $WORKDIR/toashort.root                  \
            -d $DEBUG --!
    fi

    if ( ! reuse_file $VELOCITY_TXT ); then

        for A in \
            1500.0 1510.0 1520.0                       \
            1530.0 1532.0 1534.0 1536.0 1538.0         \
            1540.0 1542.0 1544.0 1546.0 1548.0         \
            1550.0 1560.0 1570.0 1580.0; do

            rm -f event.root

            JAcousticsEventBuilder \
                -a $WORKDIR/detector.detx              \
                -f $WORKDIR/toashort.root              \
                -@ "Q            = $QUANTILE;"         \
                -@ "TMax_s       = $TMAX_S;"           \
                -@ "numberOfHits = $NUMBER_OF_HITS;"   \
                -T $WORKDIR/tripod.txt                 \
                -B $WORKDIR/base.txt                   \
                -o $WORKDIR/event.root                 \
                -V "$A"                                \
                -d $DEBUG --!

            RMS=`$DIR/JPrintRMS -f $WORKDIR/event.root -N 1 -d0`

            let RMS="$RMS*1000"

            echo "$A $RMS" >> $VELOCITY_TXT

        done
    fi

    JGraph \
        -f $VELOCITY_TXT                  \
        -o $VELOCITY_ROOT

    JPlot1D \
        -f $VELOCITY_ROOT:\.\*            \
        -x "1490 1590"                    \
        -\> "V [m/s]"                     \
        -\^ "RMS [ms]"                    \
        -T "$DETECTOR [$RUN,$RUN]"        \
        -o $WORKDIR/velocity.gif

fi