acoustics-fit-toolkit.sh File Reference

Go to the source code of this file.

Functions
then usage source $script nThis script contains auxiliary functions for the pre calibration of specific detectors nIt should be sourced by the detector specific scripts fi source JAcousticsToolkit sh set_variable WORKDIR pwd set_variable TMPDIR $WORKDIR typeset A TRIPODS typeset a CHI2 set_variable PRECISION set_variable EPSILON set_variable RADIUS_M function	initialise ()
function	clean ()
function	fixStrings ()
function	fixTripods ()
function	getCenter ()
function	setCenter ()
function	getChi2 ()
function	fitPositionOfString ()
function	fitStretchingOfString ()
function	fitPositionOfTripod ()
function	fitPositionOfModule ()
function	fitPositionOfAnchor ()
function	fitRotationOfAnchor ()
function	stage_0 ()
function	stage_1A ()
function	stage_1B ()
function	stage_2 ()
function	stage_3 ()
function	stage_D ()

Variables
	script

Function Documentation

then usage source $script nThis script contains auxiliary functions for the pre calibration of specific detectors nIt should be sourced by the detector specific scripts fi source JAcousticsToolkit sh set_variable WORKDIR pwd set_variable TMPDIR $WORKDIR typeset A TRIPODS typeset a CHI2 set_variable PRECISION set_variable EPSILON set_variable RADIUS_M function initialise

(

)

Definition at line 66 of file acoustics-fit-toolkit.sh.

{
    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  TRIPODS

    TRANSMITTERS=(`get_strings_with_transmitter $WORKDIR/transmitter.txt`)

    mkdir -p  $TMPDIR
}

function clean

(

)

Definition at line 86 of file acoustics-fit-toolkit.sh.

{
    rm -rf $TMPDIR >& /dev/null
    
    rm -f  {acoustics_fit_parameters,acoustics_trigger_parameters}.txt

    JAcoustics.sh  $DETECTOR_ID    
}

function fixStrings

(

)

Definition at line 101 of file acoustics-fit-toolkit.sh.

{
    typeset -a BUFFER

    BUFFER=(`echo $*`)

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

function fixTripods

(

)

Definition at line 119 of file acoustics-fit-toolkit.sh.

{
    typeset -a BUFFER

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

function getCenter

(

)

Definition at line 138 of file acoustics-fit-toolkit.sh.

{
    typeset -A __TRIPODS__

    get_tripods $TRIPOD __TRIPODS__

    let "X = 0.0"
    let "Y = 0.0"
    let "Z = 0.0"

    if (( ${#__TRIPODS__} >= 1 )); then
        
        for _X _Y _Z in `echo ${(@v)__TRIPODS__}`; do
            let "X = $X + $_X"
            let "Y = $Y + $_Y"
            let "Z = $Z + $_Z"
        done

        let "X = $X / ${#__TRIPODS__}"
        let "Y = $Y / ${#__TRIPODS__}"
        let "Z = $Z / ${#__TRIPODS__}"
    fi

    printf "%15.5f %15.5f %15.5f"  $X  $Y  $Z   
}

function setCenter

(

)

Definition at line 172 of file acoustics-fit-toolkit.sh.

{
    getCenter | read X Y Z

    let "X = $X - $1"
    let "Y = $Y - $2"
    let "Z = $Z - $3"

    printf "Fix center %15.5f %15.5f %15.5f\n"  $X  $Y  $Z
    
    JEditTripod   -f $TRIPOD   -T "-1 sub $X $Y $Z" -o $TRIPOD   -d 0  >&  /dev/null
    JEditDetector -a $DETECTOR -S "-1 sub $X $Y $Z" -o $DETECTOR -d 0  >&  /dev/null
}

function getChi2

(

)

Definition at line 192 of file acoustics-fit-toolkit.sh.

{
    for (( m = 0; $m != 3; ++m )); do

        rm -f $TMPDIR/katoomba.{root,log} >& /dev/null

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

        ERROR=$?

        if (( $ERROR == 0 )); then
        
            let "VALUE = $(JPrintResult -f $TMPDIR/katoomba.root:chi2 -F GetMean)" >>& $TMPDIR/katoomba.log

            ERROR=$?

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

                eval $1=$VALUE

                return
            else
                printf "warning: exit code let %d\n" $ERROR
            fi
        else
            printf "warning: exit code JKatoomba %d\n" $ERROR
        fi

        # backup error status

        set_variable  UDIR  `mktemp -d $WORKDIR/XXXXXX`
        
        cp -p $TMPDIR/katoomba.log      $UDIR
        cp -p $TMPDIR/katoomba.root     $UDIR
        cp -p $DETECTOR                 $UDIR
        cp -p $TRIPOD                   $UDIR
        cp -p $WORKDIR/transmitter.txt  $UDIR

        printf "warning: backup %s\n" $UDIR
    done

    printf "warning: no valid chi2 -> exit\n"
    exit 1
}

function fitPositionOfString

(

)

Definition at line 260 of file acoustics-fit-toolkit.sh.

{
    typeset -A BUFFER

    BUFFER=(X 0.0 Y 0.0 Z 0.0)
    
    BUFFER[$2]=$3

    if [[ "$2" == "Z" ]]; then
        OPTION=-s                  # optical modules only
    else
        OPTION=-S                  # all modules
    fi

    for (( ; $N != $NUMBER_OF_ITERATIONS; ++N )); do

        cp $DETECTOR $TMPDIR/detector.${DETECTOR:e}
        
        JEditDetector -a $DETECTOR $OPTION "$1 add ${(v)BUFFER}" -o $DETECTOR -q -d 0 >& /dev/null

        getChi2 CHI2\[2\]

        if (( $CHI2[2] >= $CHI2[1] - $PRECISION )); then

            mv $TMPDIR/detector.${DETECTOR:e} $DETECTOR

            break
        fi

        CHI2[1]=$CHI2[2]
    done

    for (( ; $N != $NUMBER_OF_ITERATIONS; ++N )); do
        
        cp $DETECTOR $TMPDIR/detector.${DETECTOR:e}

        JEditDetector -a $DETECTOR $OPTION "$1 sub ${(v)BUFFER}" -o $DETECTOR -q -d 0 >& /dev/null

        getChi2 CHI2\[2\]

        if (( $CHI2[2] >= $CHI2[1] - $PRECISION )); then

            mv $TMPDIR/detector.${DETECTOR:e} $DETECTOR

            break
        fi

        CHI2[1]=$CHI2[2]
    done

    printf "string     %04d  %s %6.4f        %6d %8.4f\n" $1 $2 $3 $N $CHI2[1]
}

function fitStretchingOfString

(

)

Definition at line 320 of file acoustics-fit-toolkit.sh.

{
    let "Z = $CAN_ZMAX_M * $2"

    for (( ; $N != $NUMBER_OF_ITERATIONS; ++N )); do

        cp $DETECTOR $TMPDIR/detector.${DETECTOR:e}

        JEditDetector -a $DETECTOR -s "$1 mul $2" -s "$1 sub 0.0 0.0 $Z" -o $DETECTOR -q -d 0 >& /dev/null

        getChi2 CHI2\[2\]

        if (( $CHI2[2] >= $CHI2[1] - $PRECISION )); then

            mv $TMPDIR/detector.${DETECTOR:e} $DETECTOR
            
            break
        fi

        CHI2[1]=$CHI2[2]
    done

    for (( ; $N != $NUMBER_OF_ITERATIONS; ++N )); do
        
        cp $DETECTOR $TMPDIR/detector.${DETECTOR:e}

        JEditDetector -a $DETECTOR -s "$1 add 0.0 0.0 $Z" -s "$1 div $2" -o $DETECTOR -q -d 0 >& /dev/null

        getChi2 CHI2\[2\]

        if (( $CHI2[2] >= $CHI2[1] - $PRECISION )); then

            mv $TMPDIR/detector.${DETECTOR:e} $DETECTOR

            break
        fi

        CHI2[1]=$CHI2[2]
    done

    printf "string     %04d  M %6.4f        %6d %8.4f\n" $1 $2 $N $CHI2[1]
}

function fitPositionOfTripod

(

)

Definition at line 371 of file acoustics-fit-toolkit.sh.

{
    typeset -A BUFFER

    BUFFER=(X 0.0 Y 0.0 Z 0.0)

    BUFFER[$2]=$3
    
    for (( ; $N != $NUMBER_OF_ITERATIONS; ++N )); do

        cp $TRIPOD $TMPDIR/tripod.txt
        
        JEditTripod -f $TRIPOD -T "$1 add ${(v)BUFFER}" -o $TRIPOD -q -d 0 >& /dev/null

        getChi2 CHI2\[2\]

        if (( $CHI2[2] >= $CHI2[1] - $PRECISION )); then

            mv $TMPDIR/tripod.txt $TRIPOD

            break
        fi

        CHI2[1]=$CHI2[2]
    done

    for (( ; $N != $NUMBER_OF_ITERATIONS; ++N )); do
        
        cp $TRIPOD $TMPDIR/tripod.txt

        JEditTripod -f $TRIPOD -T "$1 sub ${(v)BUFFER}" -o $TRIPOD -q -d 0 >& /dev/null

        getChi2 CHI2\[2\]

        if (( $CHI2[2] >= $CHI2[1] - $PRECISION )); then

            mv $TMPDIR/tripod.txt $TRIPOD

            break
        fi

        CHI2[1]=$CHI2[2]
    done

    printf "tripod       %2d  %s %6.4f        %6d %8.4f\n" $1 $2 $3 $N $CHI2[1]
}

function fitPositionOfModule

(

)

Definition at line 426 of file acoustics-fit-toolkit.sh.

{
    typeset -A BUFFER

    BUFFER=(X 0.0 Y 0.0 Z $3)

    set_variable  NUMBER_OF_STEPS   10       # limit number of steps
    set_variable  MODULE            `getModule -a $DETECTOR -L "$1 $2"`

    for (( i=0 ; $i != $NUMBER_OF_STEPS && $N != $NUMBER_OF_ITERATIONS; ++i, ++N )); do
        
        cp $DETECTOR $TMPDIR/detector.${DETECTOR:e}

        JEditDetector -a $DETECTOR -M "$MODULE add ${(v)BUFFER}" -o $DETECTOR -q -d 0 >& /dev/null

        getChi2 CHI2\[2\]

        if (( $CHI2[2] >= $CHI2[1] - $PRECISION )); then

            mv $TMPDIR/detector.${DETECTOR:e} $DETECTOR
            
            break
        fi

        CHI2[1]=$CHI2[2]
    done

    for (( i=0 ; $i != $NUMBER_OF_STEPS && $N != $NUMBER_OF_ITERATIONS; ++i, ++N )); do
        
        cp $DETECTOR $TMPDIR/detector.${DETECTOR:e}

        JEditDetector -a $DETECTOR -M "$MODULE sub ${(v)BUFFER}" -o $DETECTOR -q -d 0 >& /dev/null

        getChi2 CHI2\[2\]

        if (( $CHI2[2] >= $CHI2[1] - $PRECISION )); then

            mv $TMPDIR/detector.${DETECTOR:e} $DETECTOR

            break
        fi

        CHI2[1]=$CHI2[2]
    done

    printf "module     %04d.%02d %6.4f        %6d %8.4f\n" $1 $2 $3 $N $CHI2[1]
}

function fitPositionOfAnchor

(

)

Definition at line 486 of file acoustics-fit-toolkit.sh.

{
    typeset -A BUFFER

    BUFFER=(X 0.0 Y 0.0 Z 0.0)

    BUFFER[$2]=$3

    if [[ "$2" == "Z" ]]; then

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

        ADD="-M $MODULE add ${(v)BUFFER}"                  # base module only
        SUB="-M $MODULE sub ${(v)BUFFER}"                  # base module only
    else
        ADD="-S $1      add ${(v)BUFFER}"                  # all modules
        SUB="-S $1      sub ${(v)BUFFER}"                  # all modules
    fi

    for (( ; $N != $NUMBER_OF_ITERATIONS; ++N )); do

        cp $DETECTOR $TMPDIR/detector.${DETECTOR:e}

        JEditDetector -a $DETECTOR $ADD -o $DETECTOR -q -d 0 >& /dev/null

        getChi2 CHI2\[2\]

        if (( $CHI2[2] >= $CHI2[1] - $PRECISION )); then

            mv $TMPDIR/detector.${DETECTOR:e} $DETECTOR

            break
        fi

        CHI2[1]=$CHI2[2]
    done

    for (( ; $N != $NUMBER_OF_ITERATIONS; ++N )); do

        cp $DETECTOR $TMPDIR/detector.${DETECTOR:e}

        JEditDetector -a $DETECTOR $SUB -o $DETECTOR -q -d 0 >& /dev/null

        getChi2 CHI2\[2\]

        if (( $CHI2[2] >= $CHI2[1] - $PRECISION )); then

            mv $TMPDIR/detector.${DETECTOR:e} $DETECTOR

            break
        fi

        CHI2[1]=$CHI2[2]
    done

    printf "anchor     %04d  %s %6.4f        %6d %8.4f\n" $1 $2 $3 $N $CHI2[1]
}

function fitRotationOfAnchor

(

)

Definition at line 554 of file acoustics-fit-toolkit.sh.

{
    let "ROT = $2"

    for (( ; $N != $NUMBER_OF_ITERATIONS; ++N )); do

        cp $WORKDIR/hydrophone.txt  $TMPDIR/hydrophone.txt
        cp $WORKDIR/transmitter.txt $TMPDIR/transmitter.txt
        
        JEditHydrophone  -f $WORKDIR/hydrophone.txt  -S "$1 rot $ROT" -o $WORKDIR/hydrophone.txt  -q -d 0 >& /dev/null
        JEditTransmitter -f $WORKDIR/transmitter.txt -S "$1 rot $ROT" -o $WORKDIR/transmitter.txt -q -d 0 >& /dev/null

        getChi2 CHI2\[2\]

        if (( $CHI2[2] >= $CHI2[1] - $PRECISION )); then

            mv $TMPDIR/hydrophone.txt  $WORKDIR/hydrophone.txt
            mv $TMPDIR/transmitter.txt $WORKDIR/transmitter.txt

            break
        fi

        CHI2[1]=$CHI2[2]
    done

    let "ROT = -1.0 * $ROT"

    for (( ; $N != $NUMBER_OF_ITERATIONS; ++N )); do

        cp $WORKDIR/hydrophone.txt  $TMPDIR/hydrophone.txt
        cp $WORKDIR/transmitter.txt $TMPDIR/transmitter.txt
        
        JEditHydrophone  -f $WORKDIR/hydrophone.txt  -S "$1 rot $ROT" -o $WORKDIR/hydrophone.txt  -q -d 0 >& /dev/null
        JEditTransmitter -f $WORKDIR/transmitter.txt -S "$1 rot $ROT" -o $WORKDIR/transmitter.txt -q -d 0 >& /dev/null

        getChi2 CHI2\[2\]

        if (( $CHI2[2] >= $CHI2[1] - $PRECISION )); then

            mv $TMPDIR/hydrophone.txt  $WORKDIR/hydrophone.txt
            mv $TMPDIR/transmitter.txt $WORKDIR/transmitter.txt

            break
        fi

        CHI2[1]=$CHI2[2]
    done

    printf "string     %04d  R %6.4f        %6d %8.4f\n" $1 $2 $N $CHI2[1]
}

function stage_0

(

)

Definition at line 611 of file acoustics-fit-toolkit.sh.

{
    cat>$TMPDIR/acoustics_fit_parameters.txt<<EOF
`egrep Tmax_s $WORKDIR/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
    
    getChi2 CHI2\[1\]

    set_variable  NUMBER_OF_ITERATIONS       1000

    for DX_M in 0.5 0.2 0.1; do              # determine (x,y,z) positions of tripods

        let "A_RAD = $DX_M / $RADIUS_M"      # use maximal horizontal distance between T-bar and emitter/hydrophone

        for (( N = 0; $N < $NUMBER_OF_ITERATIONS; ++N )); do

            CHI2[3]=$CHI2[1]
            
            for ID in ${(k)TRIPODS}; do
                fitPositionOfTripod $ID X $DX_M
                fitPositionOfTripod $ID Y $DX_M
                fitPositionOfTripod $ID Z $DX_M
            done

            for STRING in $TRANSMITTERS[*]; do  
                fitPositionOfAnchor $STRING X $DX_M
                fitPositionOfAnchor $STRING Y $DX_M
                fitPositionOfAnchor $STRING Z $DX_M
                fitRotationOfAnchor $STRING $A_RAD
            done

            if (( $CHI2[3] - $CHI2[1] < $EPSILON )); then
                break
            fi
        done

        if (( $N >= $NUMBER_OF_ITERATIONS )); then
            printf "warning: reached maximum number of iterations %d - convergence %7.3f\n" $N $(($CHI2[3] - $CHI2[1]))
        fi
    done

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

function stage_1A

(

)

Definition at line 674 of file acoustics-fit-toolkit.sh.

{
    cat>$TMPDIR/acoustics_fit_parameters.txt<<EOF
`egrep Tmax_s $WORKDIR/acoustics_fit_parameters.txt`
Nmin        =   3;
sigma_s     = 250.0e-6;
stdev       =  10.0;
mestimator  =   2;
option      =   0;
EOF

    getChi2 CHI2\[1\]

    set_variable  NUMBER_OF_ITERATIONS       1000

    for DX_M in 0.5 0.2; do                  # determine (x,y,z) positions of strings and tripods

        let "A_RAD = $DX_M / $RADIUS_M"      # use maximal horizontal distance between T-bar and emitter/hydrophone

        for (( N = 0; $N < $NUMBER_OF_ITERATIONS; ++N )); do

            CHI2[3]=$CHI2[1]
            
            for STRING in $STRINGS[*]; do       
                fitPositionOfString $STRING X $DX_M
                fitPositionOfString $STRING Y $DX_M
                fitPositionOfString $STRING Z $DX_M
            done

            for ID in ${(k)TRIPODS}; do
                fitPositionOfTripod $ID X $DX_M
                fitPositionOfTripod $ID Y $DX_M
                fitPositionOfTripod $ID Z $DX_M
            done

            for STRING in $TRANSMITTERS[*]; do  
                fitPositionOfAnchor $STRING X $DX_M
                fitPositionOfAnchor $STRING Y $DX_M
                fitPositionOfAnchor $STRING Z $DX_M
                fitRotationOfAnchor $STRING $A_RAD
            done

            if (( $CHI2[3] - $CHI2[1] < $EPSILON )); then
                break
            fi
        done

        if (( $N >= $NUMBER_OF_ITERATIONS )); then
            printf "warning: reached maximum number of iterations %d - convergence %7.3f\n" $N $(($CHI2[3] - $CHI2[1]))
        fi
    done

    getChi2 CHI2\[1\]

    set_variable  NUMBER_OF_ITERATIONS       2000

    for DX_M in 0.1; do                      # determine (z) positions of modules
        
        for (( N = 0; $N < $NUMBER_OF_ITERATIONS; ++N )); do
            
            CHI2[3]=$CHI2[1]

            for STRING in $STRINGS[*]; do
                for (( FLOOR = 1; $FLOOR <= 18; FLOOR += 1 )); do
                    fitPositionOfModule $STRING $FLOOR $DX_M
                done
            done

            if (( $CHI2[3] - $CHI2[1] < $EPSILON )); then
                break
            fi
        done

        if (( $N >= $NUMBER_OF_ITERATIONS )); then
            printf "warning: reached maximum number of iterations %d - convergence %7.3f\n" $N $(($CHI2[3] - $CHI2[1]))
        fi
    done
}

function stage_1B

(

)

Definition at line 759 of file acoustics-fit-toolkit.sh.

{
    cat>$TMPDIR/acoustics_fit_parameters.txt<<EOF
`egrep Tmax_s $WORKDIR/acoustics_fit_parameters.txt`
Nmin        =   3;
sigma_s     = 100.0e-6;
stdev       =  10.0;
mestimator  =   0;
option      =   1;
EOF
 
    for STRING in $STRINGS[*]; do            # determine stretching and (z) position of given string          
        
        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

        getChi2 CHI2\[1\]

        set_variable  NUMBER_OF_ITERATIONS       1000

        for MUL in 0.005 0.001; do

            DX_M=0.2

            for (( N = 0; $N < $NUMBER_OF_ITERATIONS; ++N )); do

                CHI2[3]=$CHI2[1]

                fitPositionOfString   $STRING Z  $DX_M
                fitStretchingOfString $STRING $MUL

                if (( $CHI2[3] - $CHI2[1] < $EPSILON )); then
                    break
                fi
            done

            if (( $N >= $NUMBER_OF_ITERATIONS )); then
                printf "warning: reached maximum number of iterations %d - convergence %7.3f\n" $N $(($CHI2[3] - $CHI2[1]))
            fi
        done

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

function stage_2

(

)

Definition at line 810 of file acoustics-fit-toolkit.sh.

{
    cat>$TMPDIR/acoustics_fit_parameters.txt<<EOF
`egrep Tmax_s $WORKDIR/acoustics_fit_parameters.txt`
Nmin        =   3;
sigma_s     = 100.0e-6;
stdev       =  10.0;
mestimator  =   0;
option      =   1;
EOF

    getChi2 CHI2\[1\]

    set_variable  NUMBER_OF_ITERATIONS       1000

    for DX_M in 0.2; do                      # determine (x,y,z) positions of strings and tripods

        let "A_RAD = $DX_M / $RADIUS_M"      # use maximal horizontal distance between T-bar and emitter/hydrophone
        
        for (( N = 0; $N < $NUMBER_OF_ITERATIONS; ++N )); do

            CHI2[3]=$CHI2[1]
            
            for STRING in $STRINGS[*]; do       
                fitPositionOfString $STRING X $DX_M
                fitPositionOfString $STRING Y $DX_M
                fitPositionOfString $STRING Z $DX_M
            done

            for ID in ${(k)TRIPODS}; do
                fitPositionOfTripod $ID X $DX_M
                fitPositionOfTripod $ID Y $DX_M
                fitPositionOfTripod $ID Z $DX_M
            done

            for STRING in $TRANSMITTERS[*]; do  
                fitPositionOfAnchor $STRING X $DX_M
                fitPositionOfAnchor $STRING Y $DX_M
                fitPositionOfAnchor $STRING Z $DX_M
                fitRotationOfAnchor $STRING $A_RAD
            done

            if (( $CHI2[3] - $CHI2[1] < $EPSILON )); then
                break
            fi
        done

        if (( $N >= $NUMBER_OF_ITERATIONS )); then
            warning "reached maximum number of iterations $N - convergence $(($CHI2[3] - $CHI2[1]))"
        fi
    done

    getChi2 CHI2\[1\]

    set_variable  NUMBER_OF_ITERATIONS       2000

    for DX_M in 0.1; do                      # determine (z) positions of modules and (x,y) positions of strings

        for (( N = 0; $N < $NUMBER_OF_ITERATIONS; ++N )); do

            CHI2[3]=$CHI2[1]
            
            for STRING in $STRINGS[*]; do       

                for (( FLOOR = 1; $FLOOR <= 18; FLOOR += 1 )); do
                    fitPositionOfModule $STRING $FLOOR $DX_M
                done

                fitPositionOfString $STRING X $DX_M
                fitPositionOfString $STRING Y $DX_M
            done

            if (( $CHI2[3] - $CHI2[1] < $EPSILON )); then
                break
            fi
        done

        if (( $N >= $NUMBER_OF_ITERATIONS )); then
            warning "reached maximum number of iterations $N - convergence $(($CHI2[3] - $CHI2[1]))"
        fi
    done
}

function stage_3

(

)

Definition at line 898 of file acoustics-fit-toolkit.sh.

{
    cat>$TMPDIR/acoustics_fit_parameters.txt<<EOF
`egrep Tmax_s $WORKDIR/acoustics_fit_parameters.txt`
Nmin        =   3;
sigma_s     =  50.0e-6;
stdev       =  10.0;
mestimator  =   0;
option      =   1;
EOF

    getChi2 CHI2\[1\]

    set_variable  NUMBER_OF_ITERATIONS       1000

    for DX_M in 0.10 0.05; do                # determine (x,y,z) positions of strings and tripods

        let "A_RAD = $DX_M / $RADIUS_M"      # use maximal horizontal distance between T-bar and emitter/hydrophone
        
        for (( N = 0; $N < $NUMBER_OF_ITERATIONS; ++N )); do

            CHI2[3]=$CHI2[1]
            
            for STRING in $STRINGS[*]; do       
                fitPositionOfString $STRING X $DX_M
                fitPositionOfString $STRING Y $DX_M
                fitPositionOfString $STRING Z $DX_M
            done

            for ID in ${(k)TRIPODS}; do
                fitPositionOfTripod $ID X $DX_M
                fitPositionOfTripod $ID Y $DX_M
                fitPositionOfTripod $ID Z $DX_M
            done

            for STRING in $TRANSMITTERS[*]; do  
                fitPositionOfAnchor $STRING X $DX_M
                fitPositionOfAnchor $STRING Y $DX_M
                fitPositionOfAnchor $STRING Z $DX_M
                fitRotationOfAnchor $STRING $A_RAD
            done

            if (( $CHI2[3] - $CHI2[1] < $EPSILON )); then
                break
            fi
        done

        if (( $N >= $NUMBER_OF_ITERATIONS )); then
            warning "reached maximum number of iterations $N - convergence $(($CHI2[3] - $CHI2[1]))"
        fi
    done

    getChi2 CHI2\[1\]

    set_variable  NUMBER_OF_ITERATIONS       2000

    for DX_M in 0.05; do                     # determine (z) positions of modules and (x,y) positions of strings

        for (( N = 0; $N < $NUMBER_OF_ITERATIONS; ++N )); do

            CHI2[3]=$CHI2[1]
            
            for STRING in $STRINGS[*]; do       

                for (( FLOOR = 1; $FLOOR <= 18; FLOOR += 1 )); do
                    fitPositionOfModule $STRING $FLOOR $DX_M
                done

                fitPositionOfString $STRING X $DX_M
                fitPositionOfString $STRING Y $DX_M
            done

            if (( $CHI2[3] - $CHI2[1] < $EPSILON )); then
                break
            fi
        done

        if (( $N >= $NUMBER_OF_ITERATIONS )); then
            warning "reached maximum number of iterations $N - convergence $(($CHI2[3] - $CHI2[1]))"
        fi
    done
}

function stage_D

(

)

Definition at line 988 of file acoustics-fit-toolkit.sh.

{
    cat>$TMPDIR/acoustics_fit_parameters.txt<<EOF
`egrep Tmax_s $WORKDIR/acoustics_fit_parameters.txt`
Nmin        =   3;
sigma_s     =  50.0e-6;
stdev       =  10.0;
mestimator  =   0;
option      =   1;
EOF

    getChi2 CHI2\[1\]

    set_variable  HYDROPHONE                 
    set_variable  NUMBER_OF_ITERATIONS       2000

    for DX_M in 0.10; do                     # determine rotation of string and z-position of anchors

        let "A_RAD = $DX_M / $RADIUS_M"      # use maximal horizontal distance between T-bar and emitter/hydrophone

        for (( N = 0; $N < $NUMBER_OF_ITERATIONS; ++N )); do

            CHI2[3]=$CHI2[1]

            for STRING in $STRINGS[*]; do
                fitRotationOfAnchor $STRING $A_RAD
                fitPositionOfAnchor $STRING Z $DX_M
            done

            if (( $CHI2[3] - $CHI2[1] < $EPSILON )); then
                break
            fi
        done

        if (( $N >= $NUMBER_OF_ITERATIONS )); then
            warning "reached maximum number of iterations $N - convergence $(($CHI2[3] - $CHI2[1]))"
        fi
    done

    getChi2 CHI2\[1\]

    set_variable  NUMBER_OF_ITERATIONS       2000

    for DX_M in 0.05; do                       # determine position of strings and tripods

        for (( N = 0; $N < $NUMBER_OF_ITERATIONS; ++N )); do

            CHI2[3]=$CHI2[1]

            for STRING in $STRINGS[*]; do
                fitPositionOfString $STRING X $DX_M
                fitPositionOfString $STRING Y $DX_M
                fitPositionOfString $STRING Z $DX_M
                fitPositionOfAnchor $STRING Z $DX_M
            done

            for ID in ${(k)TRIPODS}; do
                fitPositionOfTripod $ID X $DX_M
                fitPositionOfTripod $ID Y $DX_M
                fitPositionOfTripod $ID Z $DX_M
            done

            if (( $CHI2[3] - $CHI2[1] < $EPSILON )); then
                break
            fi
        done

        if (( $N >= $NUMBER_OF_ITERATIONS )); then
            warning "reached maximum number of iterations $N - convergence $(($CHI2[3] - $CHI2[1]))"
        fi
    done

    getChi2 CHI2\[1\]

    set_variable  NUMBER_OF_ITERATIONS       1000

    for DX_M in 0.05; do                     # determine (z) positions of modules and (x,y) positions of strings

        for (( N = 0; $N < $NUMBER_OF_ITERATIONS; ++N )); do

            CHI2[3]=$CHI2[1]

            for STRING in $STRINGS[*]; do

                for (( FLOOR = 1; $FLOOR <= 18; FLOOR += 1 )); do
                    fitPositionOfModule $STRING $FLOOR $DX_M
                done

                fitPositionOfString $STRING X $DX_M
                fitPositionOfString $STRING Y $DX_M
            done

            if (( $CHI2[3] - $CHI2[1] < $EPSILON )); then
                break
            fi
        done

        if (( $N >= $NUMBER_OF_ITERATIONS )); then
            warning "reached maximum number of iterations $N - convergence $(($CHI2[3] - $CHI2[1]))"
        fi
    done
}

Variable Documentation

script

Definition at line 2 of file acoustics-fit-toolkit.sh.