pre-calibration_D0ORCA011.sh

Go to the documentation of this file.

#!/usr/bin/env zsh
script=${0##*/}
 
# --------------------------------------------------------------------------------------------
#
# Auxiliary script to derive the geometry of the D0ORCA011 (123) detector from other detectors.
#
# --------------------------------------------------------------------------------------------
 
if [ -z $JPP_DIR ]; then
    echo "Variable JPP_DIR undefined."
    exit
fi
 
source $JPP_DIR/setenv.sh $JPP_DIR >& /dev/null
 
set_variable  DETECTOR_ID            123
 
if do_usage $*; then
    usage "$script <source directory A> <source directory B>"                           \
          "\nThe source directory A corresponds to a pre-calibrated D0ORCA010 detector."  \
          "\nThe source directory B corresponds to a pre-calibrated D1ORCA011 detector."
fi
 
case $# in
    2)  set_variable DIR_DETECTOR_A       ${1:a}
        set_variable DIR_DETECTOR_B       ${2:a};;
    *)  fatal "Wrong number of arguments."
esac
 
set_variable  DETECTOR_A      $DIR_DETECTOR_A/detector.datx
set_variable  TRIPOD_A        $DIR_DETECTOR_A/tripod.txt
set_variable  DETECTOR_B      $DIR_DETECTOR_B/detector.datx
set_variable  TRIPOD_B        $DIR_DETECTOR_B/tripod.txt
 
if [[ ! -f $DETECTOR_A ]]; then
    fatal "Missing pre-calibrated D0ORCA010 input detector file"
fi
 
if [[ ! -f $DETECTOR_B ]]; then
    fatal "Missing pre-calibrated D1ORCA011 input detector file"
fi
 
if [[ ! -f $TRIPOD_B ]]; then
    fatal "Missing pre-calibrated D1ORCA011 input tripod file"
fi
 
# D1ORCA011 detector striped of all the new strings wrt D0ORCA011 and ID of D0ORCA011 assigned
JEditDetector \
    -a $DETECTOR_B                  \
    -@ "id=$DETECTOR_ID"            \
    -k "1 2 3 4 9 10 12 18 28"      \
    -o detector.datx 
 
# D0ORCA010 detector striped of all strings but the ones in D0ORCA011
JEditDetector \
    -a $DETECTOR_A                     \
    -k "1 2 3 9 10 11 12 19"           \
    -o detector_tmp.datx 
 
cat >tripod_tmp.txt<<EOF
EOF
 
# D0ORCA010 detector aligned with D1ORCA011 at the base level.
JAlignDetector \
    -A                                 \
    -O                                 \
    -a detector_tmp.datx               \
    -T tripod_tmp.txt                  \
    -b detector.datx                   \
    -r 0-0
 
# D0ORCA010 detector striped of all the strings but ORCA.0011
JEditDetector \
    -a detector_tmp.datx               \
    -k "11 19"                         \
    -o detector_tmp.datx 
 
# The two detector merged
JMergeDetector \
    -a detector.datx                   \
    -a detector_tmp.datx               \
    -o detector.datx 
 
JEditTripod \
    -f $TRIPOD_B                       \
    -r "6"                             \
    -o tripod.txt
 
rm -f tripod_tmp.txt
rm -f detector_tmp.datx
 
JMergeDisable \
    -f $DIR_DETECTOR_A/disable.txt     \
    -f $DIR_DETECTOR_B/disable.txt     \
    -o disable.txt
 
JEditTransmitter                       \
    -f $DIR_DETECTOR_B/transmitter.txt \
    -k "1 2 3 4 9 10 12 18 28"         \
    -o transmitter.txt
 
JEditTransmitter                       \
    -f $DIR_DETECTOR_A/transmitter.txt \
    -k "11 19"                         \
    -o transmitter_tmp.txt
 
JMergeTransmitter          \
    -f transmitter.txt     \
    -f transmitter_tmp.txt \
    -o transmitter.txt
 
rm -f transmitter_tmp.txt
 
JEditHydrophone                        \
    -f $DIR_DETECTOR_B/hydrophone.txt  \
    -k "1 2 3 4 9 10 12 18 28"         \
    -o hydrophone.txt
 
JEditHydrophone                        \
    -f $DIR_DETECTOR_A/hydrophone.txt  \
    -k "11 19"                         \
    -o hydrophone_tmp.txt
 
JMergeHydrophone          \
    -f hydrophone.txt     \
    -f hydrophone_tmp.txt \
    -o hydrophone.txt
 
rm -f hydrophone_tmp.txt
 
cat>waveform.txt<<EOF
# waveform identifier; emitter identifier
 16   3
-17   3
 23   5
-24   5
EOF