input-tables.sh

Go to the documentation of this file.

#!/usr/bin/env zsh
script=${0##*/}
 
# --------------------------------------------------------------------------------------------
#
# Auxiliary script to convert csv files from input-tables to input files.
#
# --------------------------------------------------------------------------------------------
 
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  DEBUG          2
set_variable  WORKDIR        ${TMPDIR:-/tmp}
 
if do_usage $*; then
    usage "$script <directory> <detector identifier>"\
          "\nThe directory corresponds to the GIT repository \"calibration/input-tables\""
fi
 
if (( $# != 2 )); then
    fatal "Wrong number of arguments."
fi
 
set_variable  DIR           $1
set_variable  DETOID        $2
 
if [[ ! -d $DIR/data/${DETOID} ]]; then
    fatal "No directory $DIR/data/${DETOID}"
fi
 
#
# Create tripod file.
#
# \param  1    input file
# \param  2    output file
#
function tripod()
{
    TRIPOD_CSV=$1
    TRIPOD_TXT=$2
 
    if [[ -f $TRIPOD_CSV ]]; then
 
        printf "# Data from %s\n" ${TRIPOD_CSV:t}                >  $TRIPOD_TXT
        egrep -v "^#" $TRIPOD_CSV  |  tail -n +2  |  tr ',' ' '  >> $TRIPOD_TXT
    fi
}
 
#
# Create hydrophone file.
#
# \param  1    input file
# \param  2    output file
#
function hydrophone()
{
    HYDROPHONE_CSV=$1
    HYDROPHONE_TXT=$2
 
    if [[ -f $HYDROPHONE_CSV ]]; then
 
        egrep -v "^#" $HYDROPHONE_CSV  |  tail -n +2  |  tr ',' ' '  >  $WORKDIR/hydrophones.txt
 
        printf "# Data from %s\n" ${HYDROPHONE_CSV:t}                >  $HYDROPHONE_TXT
 
        while read -A BUFFER; do
            
            set_variable STRING     $BUFFER[1]
            set_variable FLOOR      $BUFFER[2]
            set_variable X          $BUFFER[3]
            set_variable Y          $BUFFER[4]
            set_variable ANGLE_DEG  $BUFFER[5]
            set_variable Z          $BUFFER[6]
            
            let "X0 = -1.0 * $Y"
            let "Y0 = -1.0 * $X"
 
            let "CT = cos($ANGLE_DEG * $PI / 180.0)"
            let "ST = sin($ANGLE_DEG * $PI / 180.0)"
 
            let "X1 = $ST * $X0  +  $CT * $Y0"
            let "Y1 = $CT * $X0  -  $ST * $Y0"
            
            printf "%4d %2d %6.2f %6.2f %6.2f\n" $STRING $FLOOR $X1 $Y1 $Z  >>  $HYDROPHONE_TXT
 
        done < $WORKDIR/hydrophones.txt
 
        rm -f  $WORKDIR/hydrophones.txt
    fi
}
 
#
# Create transmitter file.
#
# \param  1    input file
# \param  2    output file
#
function transmitter()
{
    TRANSMITTER_CSV=$1
    TRANSMITTER_TXT=$2
 
    if [[ -f $TRANSMITTER_CSV ]]; then
 
        egrep -v "^#" $TRANSMITTER_CSV  |  tail -n +2  |  tr ',' ' '  >  $WORKDIR/transmitters.txt
 
        printf "# Data from %s\n" ${TRANSMITTER_CSV:t}                >  $TRANSMITTER_TXT
 
        while read -A BUFFER; do
 
            set_variable ID         $BUFFER[1]
            set_variable STRING     $BUFFER[2]
            set_variable FLOOR      $BUFFER[3]
            set_variable X          $BUFFER[4]
            set_variable Y          $BUFFER[5]
            set_variable ANGLE_DEG  $BUFFER[6]
            set_variable Z          $BUFFER[7]
 
            let "X0 = -1.0 * $Y"
            let "Y0 = -1.0 * $X"
 
            let "CT = cos($ANGLE_DEG * $PI / 180.0)"
            let "ST = sin($ANGLE_DEG * $PI / 180.0)"
 
            let "X1 = $ST * $X0  +  $CT * $Y0"
            let "Y1 = $CT * $X0  -  $ST * $Y0"
            
            printf "%2d %4d %2d %6.2f %6.2f %6.2f\n" $ID $STRING $FLOOR $X1 $Y1 $Z  >>  $TRANSMITTER_TXT
 
        done < $WORKDIR/transmitters.txt
 
        rm -f  $WORKDIR/transmitters.txt
    fi
}
 
#
# Create waveform file.
#
# \param  1    input file
# \param  2    output file
#
function waveform()
{
    WAVEFORM_CSV=$1
    WAVEFORM_TXT=$2
 
    if [[ -f $WAVEFORM_CSV ]]; then
 
        egrep -v "^#" $WAVEFORM_CSV  |  tail -n +2  |  tr ',' ' '  >  $WORKDIR/waveforms.txt
 
        printf "# Data from %s\n" ${WAVEFORM_CSV:t}                >  $WAVEFORM_TXT
 
        while read -A BUFFER; do
 
            set_variable WAVEFORM_ID  $BUFFER[1]
            set_variable EMITTER_ID   $BUFFER[2]
            
            printf "%3d %3d\n" $WAVEFORM_ID $EMITTER_ID  >>  $WAVEFORM_TXT
 
        done < $WORKDIR/waveforms.txt
 
        rm -f  $WORKDIR/waveforms.txt
    fi
}
 
#
# Create mechanics file.
#
# \param  1    detector identifier
# \param  2    output file
#
function mechanics()
{
    MECHANICS_TXT=$2
 
    if   [[ $1 =~ ORC ]]; then
 
        cat>$MECHANICS_TXT<<EOF
  -1   0.00311     85.966
   9   0.00087    267.054       
EOF
    elif [[ $1 =~ ARC ]]; then
 
        cat>$MECHANICS_TXT<<EOF
  -1   0.00094    294.291
EOF
    else
        fatal "Invalid detector identifier $1."
    fi
}
 
#
# Create sound velocity file.
#
# \param  1    detector identifier
# \param  2    output file
#
function sound_velocity()
{
    SOUND_VELOCITY_TXT=$2
 
    if   [[ $1 =~ ORC ]]; then
 
        cat>$SOUND_VELOCITY_TXT<<EOF
1538.00 -17e-3 -2000.00
EOF
    elif [[ $1 =~ ARC ]]; then
 
        cat>$SOUND_VELOCITY_TXT<<EOF
1541.00 -17e-3 -2000.00
EOF
    else
        fatal "Invalid detector identifier $1."
    fi
}
 
#
# Create disable file.
#
# \param  1    detector identifier
# \param  2    output file
#
function disable()
{
    DISABLE_TXT=$2
 
    cat>$DISABLE_TXT<<EOF
#list of bad transmissions
EOF
}
 
#
# Create trigger parameters file.
#
# \param  1    detector identifier
# \param  2    output file
#
function trigger_parameters()
{
    TRIGGER_PARAMETERS_TXT=$2
 
    cat>$TRIGGER_PARAMETERS_TXT<<EOF
Q             = 0.0;
TMax_s        = 0.020;
numberOfHits  = 300;
EOF
}
 
#
# Create fit parameters file.
#
# \param  1    detector identifier
# \param  2    output file
#
function fit_parameters()
{
    FIT_PARAMETERS_TXT=$2
 
    cat>$FIT_PARAMETERS_TXT<<EOF
Nmin         = 3;
Tmax_s       = 600.0;
chi2perNDF   = 100.0;
mestimator   = 1;
option       = 3;
sigma_s      = 50.0e-6;
stdev        = 10.0e+5;
EOF
}
 
tripod               $DIR/data/${DETOID}/${DETOID}_tripods.csv       tripod.txt
hydrophone           $DIR/data/${DETOID}/${DETOID}_hydrophones.csv   hydrophone.txt
transmitter          $DIR/data/${DETOID}/${DETOID}_transmitters.csv  transmitter.txt
waveform             $DIR/data/${DETOID}/${DETOID}_waveforms.csv     waveform.txt
mechanics            ${DETOID}                                       mechanics.txt 
sound_velocity       ${DETOID}                                       sound_velocity.txt 
disable              ${DETOID}                                       disable.txt 
trigger_parameters   ${DETOID}                                       acoustics_trigger_parameters.txt
fit_parameters       ${DETOID}                                       acoustics_fit_parameters.txt