acoustics-fit-toolkit.sh

Go to the documentation of this file.

#!/bin/zsh
script=${0##*/}

# -----------------------------------------------------------------------------------------------------
#
# Toolkit for the determination of the fixed parameters in the system, a.k.a. "pre-calibration".
# The list of parameters includes (x,y,z) positions of strings and tripods, stretching of strings,
# z-positions of the optical modules, orientations of the strings and possibly z-positions of anchors.
#
# The value of a given parameter is optimised by minimising the average chi2/NDF
# of a series of global fits applied to a predefined set of acoustic events.
# For each parameter, there exists a designated function.
#
# There are also steering functions to consistently determine the optimal values for a set of parameters.
# In these, the current accuracy of the parameter values is taken into account.
# These are therefore referred to as "stages".
# The pre-calibration thus consists of a sequence of stages.
# For each detector, a complete pre-calibration procedure is implemented in a designated script.
#
# Mandatory input variables:
#
# DETECTOR              # detector file
# INPUT_FILES           # list of input files (output of JAcousticsEventBuilder[.sh])
#
# Optional input variables:
#
# HYDROPHONE            # if defined, use hydrophones in global fit; else do not use hydrophones
#
# In the following:
#
# STRINGS corresponds to list of strings to be fitted; this can be changed with function fixStrings.
# TRIPODS corresponds to list of tripods to be fitted; this can be changed with function fixTripods.
#
# TRANSMITTERS corresponds to list of strings which have an emitter on the anchor (maybe empty).
# For these strings, the orientations and the z-positions of the anchors will separately be fitted.
# This list will also be changed with function fixStrings.
# -----------------------------------------------------------------------------------------------------

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

if [[ ! $ZSH_EVAL_CONTEXT =~ :file$ ]]; then
    echo "Auxiliary script for the pre-calibrations."
    echo "This script should be sourced by the steering script for a specific detector."
    echo "Consult the script '$0' for more details." 
    exit
fi

source JAcousticsToolkit.sh

set_variable   TMPDIR     `pwd`/.$$                     # temporary directory
set_variable:  DEBUG      ACOUSTICS_DEBUG      3        # debug level
set_variable:  JOBS       ACOUSTICS_JOBS       0        # number of parallel threads
set_variable:  SLEEP_US   ACOUSTICS_SLEEP_US   100      # sleep time [us] of main thread when data queue exceeds limit

typeset -A     TRIPODS               # tripods
typeset -a     CHI2                  # chi2 values; index [1] corresponds to best value

let "NUMBER_OF_ITERATIONS = 100"     # maximal number of iterations
let "EPSILON              = 5.0E-4"  # maximal distance to minimum of chi2
let "RADIUS_M             = 1.0"     # maximal horizontal distance between T-bar and emitter/hydrophone
let "PRECISION            = 10"      # number of decimals of chi2 evaluation


# -----------------------------------------------------------------------------------------------------
# Initialise.
#
# Install input files in working directory and create temporary directory TMPDIR.
#
# -----------------------------------------------------------------------------------------------------
function initialise()
{
    eval `JPrintDetector -a $DETECTOR -O IDENTIFIER`
    eval `JPrintDetector -a $DETECTOR -O SUMMARY`
    eval `JPrintDetector -a $DETECTOR -O CAN`

    JAcoustics.sh  $DETECTOR_ID

    get_tripods  tripod.txt  TRIPODS

    TRANSMITTERS=(`get_strings_with_transmitter transmitter.txt`)

    if [[ "${HYDROPHONE+YES}" == "YES" ]]; then
        HYDROPHONES=(`get_strings_with_hydrophone hydrophone.txt`)
    else
        HYDROPHONES=()
    fi
        
    mkdir -p  $TMPDIR
}


# -----------------------------------------------------------------------------------------------------
# Backup files in given directory.
#
# \param  1       directory
# -----------------------------------------------------------------------------------------------------
function backup()
{
    mkdir -p $1

    cp -p  {$DETECTOR,tripod.txt,transmitter.txt,hydrophone.txt} $1

}


# -----------------------------------------------------------------------------------------------------
# Remove temporary directory and restore input files in working directory.
#
# -----------------------------------------------------------------------------------------------------
function clean()
{
    rm -rf $TMPDIR >& /dev/null
    
    rm -f  {acoustics_fit_parameters,acoustics_trigger_parameters}.txt

    JAcoustics.sh  $DETECTOR_ID    
}


# -----------------------------------------------------------------------------------------------------
# Fix strings.
#
# \param  1-N     identifiers
# -----------------------------------------------------------------------------------------------------
function fixStrings()
{
    typeset -a BUFFER

    BUFFER=(`echo $*`)

    if (( ${#BUFFER} > 0 )); then
        STRINGS=(${STRINGS:|BUFFER})
        TRANSMITTERS=(${TRANSMITTERS:|BUFFER})
        HYDROPHONES=(${HYDROPHONES:|BUFFER})
    fi
}


# -----------------------------------------------------------------------------------------------------
# Fix tripods.
#
# \param  1-N     identifiers
# -----------------------------------------------------------------------------------------------------
function fixTripods()
{
    typeset -a BUFFER

    BUFFER=(`echo $*`)
    
    if (( ${#BUFFER} > 0 )); then
        for ID in $BUFFER[*]; do
            unset "TRIPODS[${ID}]"
        done
    fi
}


# -----------------------------------------------------------------------------------------------------
# Evaluate current chi2.
#
# \param  1       variable containing chi2 value on return
# -----------------------------------------------------------------------------------------------------
function getChi2()
{
    rm -f $TMPDIR/katoomba.root >& /dev/null

    if (( $JOBS == 0 )); then

        JKatoomba \
            -a $DETECTOR                              \
            -f "$INPUT_FILES[*]"                      \
            -o $TMPDIR/katoomba.root                  \
            -T tripod.txt                             \
            -Y transmitter.txt                        \
            -V sound_velocity.txt                     \
            -M mechanics.txt                          \
            -@ $TMPDIR/acoustics_fit_parameters.txt   \
            -! disable.txt                            \
            ${HYDROPHONE+-H hydrophone.txt}           \
            -u                                        \
            -d 0 --!                                  >& $TMPDIR/job.log

    else

        JFremantle \
            -a $DETECTOR                              \
            -f "$INPUT_FILES[*]"                      \
            -o $TMPDIR/katoomba.root                  \
            -T tripod.txt                             \
            -Y transmitter.txt                        \
            -V sound_velocity.txt                     \
            -M mechanics.txt                          \
            -@ $TMPDIR/acoustics_fit_parameters.txt   \
            -! disable.txt                            \
            ${HYDROPHONE+-H hydrophone.txt}           \
            -u                                        \
            -N $JOBS                                  \
            -s $SLEEP_US                              \
            -d 0 --!                                  >& $TMPDIR/job.log
    fi
        
    if (( $? != 0 )); then
        fatal "non-zero exit code, see $TMPDIR/job.log"
    fi

    let "UV = $($JPP_DIR/examples/JAcoustics/JPrintChi2 -f $TMPDIR/katoomba.root -p $PRECISION)"

    if (( $? != 0 )); then
        fatal "non-zero exit code let."
    fi
    
    eval $1=$UV
}


# -----------------------------------------------------------------------------------------------------
# Stage 0.
# This stage can be used to determine the (x,y,z) positions of tripods when a number of strings is fixed.
#
# \param  1       step size
# \param  2       maximum number of extra steps
# -----------------------------------------------------------------------------------------------------
function stage_0()
{
    cat>$TMPDIR/acoustics_fit_parameters.txt<<EOF
`egrep Tmax_s acoustics_fit_parameters.txt`
Nmin        =   3;
sigma_s     = 100.0e-6;
stdev       =  10.0;
mestimator  =   0;
option      =   1;
EOF

    if (( ${#STRINGS} > 0 )); then
        JEditDetector  -a $DETECTOR  -o $TMPDIR/detector_0.datx  -k "$STRINGS[*]" -q -d 0 >& /dev/null
        JEditDetector  -a $DETECTOR  -o $DETECTOR                -r "$STRINGS[*]" -q -d 0 >& /dev/null
    fi

    PARAMETERS=()

    for ID in ${(k)TRIPODS}; do
        PARAMETERS['JEditTripod -f tripod.txt -T "$ID addx %" -o tripod.txt -q -d 0 >& /dev/null']=$1
        PARAMETERS['JEditTripod -f tripod.txt -T "$ID addy %" -o tripod.txt -q -d 0 >& /dev/null']=$1
        PARAMETERS['JEditTripod -f tripod.txt -T "$ID addz %" -o tripod.txt -q -d 0 >& /dev/null']=$1
    done

    gradient $2 Y

    printf "detector                         %6d %8.4f\n" $N $Y

    if (( ${#STRINGS} > 0 )); then
        JMergeDetector  -a $DETECTOR  -a $TMPDIR/detector_0.datx  -o $DETECTOR >& /dev/null
    fi
}


# -----------------------------------------------------------------------------------------------------
# Stage 1a.
# This stage can be used to roughly determine the positions of the strings and tripods.
# During this procedure, the strings are kept vertical in the global fits.
#
# \param  1       step size strings
# \param  2       step size tripods
# \param  3       maximum number of extra steps
# -----------------------------------------------------------------------------------------------------
function stage_1a()
{
    cat>$TMPDIR/acoustics_fit_parameters.txt<<EOF
`egrep Tmax_s acoustics_fit_parameters.txt`
Nmin        =   3;
sigma_s     = 250.0e-6;
stdev       =  10.0;
mestimator  =   2;
option      =   0;
EOF

    PARAMETERS=()

    for STRING in $STRINGS[*]; do
        PARAMETERS['JEditDetector -a $DETECTOR -S "$STRING addx %" -o $DETECTOR -q -d 0 >& /dev/null']=$1
        PARAMETERS['JEditDetector -a $DETECTOR -S "$STRING addy %" -o $DETECTOR -q -d 0 >& /dev/null']=$1
        PARAMETERS['JEditDetector -a $DETECTOR -s "$STRING addz %" -o $DETECTOR -q -d 0 >& /dev/null']=$1
    done

    for ID in ${(k)TRIPODS}; do
        PARAMETERS['JEditTripod -f tripod.txt -T "$ID addx %" -o tripod.txt -q -d 0 >& /dev/null']=$2
        PARAMETERS['JEditTripod -f tripod.txt -T "$ID addy %" -o tripod.txt -q -d 0 >& /dev/null']=$2
        PARAMETERS['JEditTripod -f tripod.txt -T "$ID addz %" -o tripod.txt -q -d 0 >& /dev/null']=$2
    done

    for STRING in ${HYDROPHONES:|TRANSMITTERS}; do
        PARAMETERS['JEditHydrophone -f hydrophone.txt -S "$STRING rot %" -o hydrophone.txt -q -d 0 >& /dev/null; JEditTransmitter -f transmitter.txt -S "$STRING rot %" -o transmitter.txt -q -d 0 >& /dev/null']=$(($1 / $RADIUS_M))
    done

    for STRING in $TRANSMITTERS[*]; do

        set_variable  MODULE  `getModule -a $DETECTOR -L "$STRING 0"`

        PARAMETERS['JEditDetector -a $DETECTOR -M "$MODULE addz %" -o $DETECTOR -q -d 0 >& /dev/null']=$1
        PARAMETERS['JEditHydrophone -f hydrophone.txt -S "$STRING rot %" -o hydrophone.txt -q -d 0 >& /dev/null; JEditTransmitter -f transmitter.txt -S "$STRING rot %" -o transmitter.txt -q -d 0 >& /dev/null']=$(($1 / $RADIUS_M))
    done

    gradient $3 Y

    printf "detector                         %6d %8.4f\n" $N $Y
}


# -----------------------------------------------------------------------------------------------------
# Stage 1b.
# This stage can be used to roughly determine the z-positions of the modules.
# During this procedure, the strings are kept vertical in the global fits.
#
# \param  1       step size
# \param  2       maximum number of extra steps
# -----------------------------------------------------------------------------------------------------
function stage_1b()
{
    cat>$TMPDIR/acoustics_fit_parameters.txt<<EOF
`egrep Tmax_s acoustics_fit_parameters.txt`
Nmin        =   3;
sigma_s     = 250.0e-6;
stdev       =  10.0;
mestimator  =   2;
option      =   0;
EOF

    PARAMETERS=()

    for STRING in $STRINGS[*]; do       

        for (( FLOOR = 1; $FLOOR <= 18; FLOOR += 1 )); do

            set_variable  MODULE  `getModule -a $DETECTOR -L "$STRING $FLOOR"`

            PARAMETERS['JEditDetector -a $DETECTOR -M "$MODULE addz %" -o $DETECTOR -q -d 0 >& /dev/null']=$1
        done
    done
        
    gradient $2 Y
        
    printf "detector                         %6d %8.4f\n" $N $Y
}


# -----------------------------------------------------------------------------------------------------
# Stage 1B.
# This stage can be used to simultaneously determine the stretching and z-positions of individual strings.
#
# \param  1       step size stretching
# \param  2       step size z-position
# \param  3       maximum number of extra steps
# -----------------------------------------------------------------------------------------------------
function stage_1B()
{
    cat>$TMPDIR/acoustics_fit_parameters.txt<<EOF
`egrep Tmax_s acoustics_fit_parameters.txt`
Nmin        =   3;
sigma_s     = 100.0e-6;
stdev       =  10.0;
mestimator  =   0;
option      =   1;
EOF

    for STRING in $STRINGS[*]; do
        
        JEditDetector  -a $DETECTOR  -o $TMPDIR/detector_1.datx  -r "$STRING" -q -d 0 >& /dev/null
        JEditDetector  -a $DETECTOR  -o $DETECTOR                -k "$STRING" -q -d 0 >& /dev/null

        PARAMETERS=()

        PARAMETERS['JEditDetector -a $DETECTOR -s "$STRING mul %"  -o $DETECTOR -q -d 0 >& /dev/null; JEditDetector -a $DETECTOR -s "$STRING subz \$\(\(% * $CAN_ZMAX_M\)\)"  -o $DETECTOR -q -d 0 >& /dev/null']=$1
        PARAMETERS['JEditDetector -a $DETECTOR -s "$STRING addz %" -o $DETECTOR -q -d 0 >& /dev/null']=$2

        gradient $3 Y

        printf "string     %04d                  %6d %8.4f\n" $STRING $N $Y

        JMergeDetector  -a $DETECTOR  -a $TMPDIR/detector_1.datx  -o $DETECTOR >& /dev/null
    done
}


# -----------------------------------------------------------------------------------------------------
# Stage 2a.
# This stage can be used to improve the determination the positions of the strings and tripods.
#
# \param  1       step size strings
# \param  2       step size tripods
# \param  3       maximum number of extra steps
# -----------------------------------------------------------------------------------------------------
function stage_2a()
{
    cat>$TMPDIR/acoustics_fit_parameters.txt<<EOF
`egrep Tmax_s acoustics_fit_parameters.txt`
Nmin        =   3;
sigma_s     = 100.0e-6;
stdev       =  10.0;
mestimator  =   0;
option      =   1;
EOF

    PARAMETERS=()

    for STRING in $STRINGS[*]; do
        PARAMETERS['JEditDetector -a $DETECTOR -S "$STRING addx %" -o $DETECTOR -q -d 0 >& /dev/null']=$1
        PARAMETERS['JEditDetector -a $DETECTOR -S "$STRING addy %" -o $DETECTOR -q -d 0 >& /dev/null']=$1
        PARAMETERS['JEditDetector -a $DETECTOR -s "$STRING addz %" -o $DETECTOR -q -d 0 >& /dev/null']=$1
    done

    for ID in ${(k)TRIPODS}; do
        PARAMETERS['JEditTripod -f tripod.txt -T "$ID addx %" -o tripod.txt -q -d 0 >& /dev/null']=$2
        PARAMETERS['JEditTripod -f tripod.txt -T "$ID addy %" -o tripod.txt -q -d 0 >& /dev/null']=$2
        PARAMETERS['JEditTripod -f tripod.txt -T "$ID addz %" -o tripod.txt -q -d 0 >& /dev/null']=$2
    done

    for STRING in ${HYDROPHONES:|TRANSMITTERS}; do
        PARAMETERS['JEditHydrophone -f hydrophone.txt -S "$STRING rot %" -o hydrophone.txt -q -d 0 >& /dev/null; JEditTransmitter -f transmitter.txt -S "$STRING rot %" -o transmitter.txt -q -d 0 >& /dev/null']=$(($1 / $RADIUS_M))
    done

    for STRING in $TRANSMITTERS[*]; do

        set_variable  MODULE  `getModule -a $DETECTOR -L "$STRING 0"`

        PARAMETERS['JEditDetector -a $DETECTOR -M "$MODULE addz %" -o $DETECTOR -q -d 0 >& /dev/null']=$1
        PARAMETERS['JEditHydrophone -f hydrophone.txt -S "$STRING rot %" -o hydrophone.txt -q -d 0 >& /dev/null; JEditTransmitter -f transmitter.txt -S "$STRING rot %" -o transmitter.txt -q -d 0 >& /dev/null']=$(($1 / $RADIUS_M))
    done

    gradient $3 Y

    printf "detector                         %6d %8.4f\n" $N $Y
}


# -----------------------------------------------------------------------------------------------------
# Stage 2b.
# This stage can be used to improve the determination the z-positions of the modules.
#
# \param  1       step size
# \param  2       maximum number of extra steps
# -----------------------------------------------------------------------------------------------------
function stage_2b()
{
    cat>$TMPDIR/acoustics_fit_parameters.txt<<EOF
`egrep Tmax_s acoustics_fit_parameters.txt`
Nmin        =   3;
sigma_s     = 100.0e-6;
stdev       =  10.0;
mestimator  =   0;
option      =   1;
EOF

    PARAMETERS=()

    for STRING in $STRINGS[*]; do       

        for (( FLOOR = 1; $FLOOR <= 18; FLOOR += 1 )); do

            set_variable  MODULE  `getModule -a $DETECTOR -L "$STRING $FLOOR"`

            PARAMETERS['JEditDetector -a $DETECTOR -M "$MODULE addz %" -o $DETECTOR -q -d 0 >& /dev/null']=$1
        done

        #PARAMETERS['JEditDetector -a $DETECTOR -S "$STRING addx %" -o $DETECTOR -q -d 0 >& /dev/null']=$1
        #PARAMETERS['JEditDetector -a $DETECTOR -S "$STRING addy %" -o $DETECTOR -q -d 0 >& /dev/null']=$1
    done
        
    gradient $2 Y
        
    printf "detector                         %6d %8.4f\n" $N $Y
}


# -----------------------------------------------------------------------------------------------------
# Stage 3a.
# This stage can be used to finalise the determination of the positions of the strings and tripods.
#
# \param  1       step size strings
# \param  2       step size tripods
# \param  3       maximum number of extra steps
# -----------------------------------------------------------------------------------------------------
function stage_3a()
{
   cat>$TMPDIR/acoustics_fit_parameters.txt<<EOF
`egrep Tmax_s acoustics_fit_parameters.txt`
Nmin        =   3;
sigma_s     =  50.0e-6;
stdev       =  10.0;
mestimator  =   0;
option      =   1;
EOF

    PARAMETERS=()

    for STRING in $STRINGS[*]; do
        PARAMETERS['JEditDetector -a $DETECTOR -S "$STRING addx %" -o $DETECTOR -q -d 0 >& /dev/null']=$1
        PARAMETERS['JEditDetector -a $DETECTOR -S "$STRING addy %" -o $DETECTOR -q -d 0 >& /dev/null']=$1
        PARAMETERS['JEditDetector -a $DETECTOR -s "$STRING addz %" -o $DETECTOR -q -d 0 >& /dev/null']=$1
    done

    for ID in ${(k)TRIPODS}; do
        PARAMETERS['JEditTripod -f tripod.txt -T "$ID addx %" -o tripod.txt -q -d 0 >& /dev/null']=$2
        PARAMETERS['JEditTripod -f tripod.txt -T "$ID addy %" -o tripod.txt -q -d 0 >& /dev/null']=$2
        PARAMETERS['JEditTripod -f tripod.txt -T "$ID addz %" -o tripod.txt -q -d 0 >& /dev/null']=$2
    done

    for STRING in ${HYDROPHONES:|TRANSMITTERS}; do
        PARAMETERS['JEditHydrophone -f hydrophone.txt -S "$STRING rot %" -o hydrophone.txt -q -d 0 >& /dev/null; JEditTransmitter -f transmitter.txt -S "$STRING rot %" -o transmitter.txt -q -d 0 >& /dev/null']=$(($1 / $RADIUS_M))
    done

    for STRING in $TRANSMITTERS[*]; do

        set_variable  MODULE  `getModule -a $DETECTOR -L "$STRING 0"`

        PARAMETERS['JEditDetector -a $DETECTOR -M "$MODULE addz %" -o $DETECTOR -q -d 0 >& /dev/null']=$1
        PARAMETERS['JEditHydrophone -f hydrophone.txt -S "$STRING rot %" -o hydrophone.txt -q -d 0 >& /dev/null; JEditTransmitter -f transmitter.txt -S "$STRING rot %" -o transmitter.txt -q -d 0 >& /dev/null']=$(($1 / $RADIUS_M))
    done

    gradient $3 Y

    printf "detector                         %6d %8.4f\n" $N $Y
}


# -----------------------------------------------------------------------------------------------------
# Stage 3b.
# This stage can be used to finalise the determination of the z-positions of the modules.
#
# \param  1       step size
# \param  2       maximum number of extra steps
# -----------------------------------------------------------------------------------------------------
function stage_3b()
{
    cat>$TMPDIR/acoustics_fit_parameters.txt<<EOF
`egrep Tmax_s acoustics_fit_parameters.txt`
Nmin        =   3;
sigma_s     =  50.0e-6;
stdev       =  10.0;
mestimator  =   0;
option      =   1;
EOF

    PARAMETERS=()

    for STRING in $STRINGS[*]; do       

        for (( FLOOR = 1; $FLOOR <= 18; FLOOR += 1 )); do

            set_variable  MODULE  `getModule -a $DETECTOR -L "$STRING $FLOOR"`

            PARAMETERS['JEditDetector -a $DETECTOR -M "$MODULE addz %" -o $DETECTOR -q -d 0 >& /dev/null']=$1
        done

        #PARAMETERS['JEditDetector -a $DETECTOR -S "$STRING addx %" -o $DETECTOR -q -d 0 >& /dev/null']=$1
        #PARAMETERS['JEditDetector -a $DETECTOR -S "$STRING addy %" -o $DETECTOR -q -d 0 >& /dev/null']=$1
    done

    gradient $2 Y

    printf "detector                         %6d %8.4f\n" $N $Y
}


# -----------------------------------------------------------------------------------------------------
# Stage D.
# This stage can be used to determine the rotations of the strings and z-positions of the anchors
# including the response of hydrophones.
# It includes a determination of the positions of the strings, tripods and modules.
#
# \param  1       step size strings
# \param  2       step size tripods
# \param  3       maximum number of extra steps
# -----------------------------------------------------------------------------------------------------
function stage_D()
{
    cat>$TMPDIR/acoustics_fit_parameters.txt<<EOF
`egrep Tmax_s acoustics_fit_parameters.txt`
Nmin        =   3;
sigma_s     =  50.0e-6;
stdev       =  10.0;
mestimator  =   0;
option      =   1;
EOF

    set_variable  HYDROPHONE                 

    PARAMETERS=()

    for STRING in $STRINGS[*]; do

        set_variable  MODULE  `getModule -a $DETECTOR -L "$STRING 0"`

        PARAMETERS['JEditDetector -a $DETECTOR -S "$STRING addx %" -o $DETECTOR -q -d 0 >& /dev/null']=$1
        PARAMETERS['JEditDetector -a $DETECTOR -S "$STRING addy %" -o $DETECTOR -q -d 0 >& /dev/null']=$1
        PARAMETERS['JEditDetector -a $DETECTOR -s "$STRING addz %" -o $DETECTOR -q -d 0 >& /dev/null']=$1
        PARAMETERS['JEditDetector -a $DETECTOR -M "$MODULE addz %" -o $DETECTOR -q -d 0 >& /dev/null']=$1

        PARAMETERS['JEditHydrophone -f hydrophone.txt -S "$STRING rot %" -o hydrophone.txt -q -d 0 >& /dev/null; JEditTransmitter -f transmitter.txt -S "$STRING rot %" -o transmitter.txt -q -d 0 >& /dev/null']=$(($1 / $RADIUS_M))
    done

    for ID in ${(k)TRIPODS}; do
        PARAMETERS['JEditTripod -f tripod.txt -T "$ID addx %" -o tripod.txt -q -d 0 >& /dev/null']=$2
        PARAMETERS['JEditTripod -f tripod.txt -T "$ID addy %" -o tripod.txt -q -d 0 >& /dev/null']=$2
        PARAMETERS['JEditTripod -f tripod.txt -T "$ID addz %" -o tripod.txt -q -d 0 >& /dev/null']=$2
    done

    for STRING in ${HYDROPHONES:|TRANSMITTERS}; do
        PARAMETERS['JEditHydrophone -f hydrophone.txt -S "$STRING rot %" -o hydrophone.txt -q -d 0 >& /dev/null; JEditTransmitter -f transmitter.txt -S "$STRING rot %" -o transmitter.txt -q -d 0 >& /dev/null']=$(($1 / $RADIUS_M))
    done

    for STRING in $TRANSMITTERS[*]; do

        set_variable  MODULE  `getModule -a $DETECTOR -L "$STRING 0"`

        PARAMETERS['JEditDetector -a $DETECTOR -M "$MODULE addz %" -o $DETECTOR -q -d 0 >& /dev/null']=$1
        PARAMETERS['JEditHydrophone -f hydrophone.txt -S "$STRING rot %" -o hydrophone.txt -q -d 0 >& /dev/null; JEditTransmitter -f transmitter.txt -S "$STRING rot %" -o transmitter.txt -q -d 0 >& /dev/null']=$(($1 / $RADIUS_M))
    done

    gradient $3 Y

    printf "detector                         %6d %8.4f\n" $N $Y
}