pre-calibration_D0ORCA007.sh

Go to the documentation of this file.

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

# --------------------------------------------------------------------------------------------
#
# Auxiliary script to derive the geometry of the D0ORCA007 (110) 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            110

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 D0ORCA007 and ID of D0ORCA007 assigned
JEditDetector \
    -a $DETECTOR_B                  \
    -@ "id=$DETECTOR_ID"            \
    -k "1 2 3 9 10 12"              \
    -o detector.datx 

# D0ORCA010 detector striped of all strings but the ones in D0ORCA007
JEditDetector \
    -a $DETECTOR_A                     \
    -k "1 2 3 9 10 11 12"              \
    -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"                            \
    -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

echo "# $script"                           >  disable.txt
egrep -v "^#" $DIR_DETECTOR_A/disable.txt  >> disable.txt
egrep -v "^#" $DIR_DETECTOR_B/disable.txt  >> disable.txt

JEditTransmitter \
    -f $DIR_DETECTOR_A/transmitter.txt \
    -k "1 2 3 9 10 11 12"              \
    -o transmitter.txt

JEditHydrophone \
    -f $DIR_DETECTOR_A/hydrophone.txt  \
    -k "1 2 3 9 10 11 12"              \
    -o hydrophone.txt