pre-calibration_D0ARCA021.sh

Go to the documentation of this file.

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

# --------------------------------------------------------------------------------------------
#
# Steering script for the global-fit-of-global-fits procedure of the D0ARCA021 (133) detector.
#
# --------------------------------------------------------------------------------------------

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            133
set_array     RUNS                   13290-13297
set_variable: MUL                    ACOUSTICS_MUL    -1.21e-2
set_variable: THREADS                ACOUSTICS_THREADS    5
set_variable  WORKDIR                `pwd`
set_variable  HOMEDIR                $WORKDIR

if do_usage $*; then
    usage "$script <source directory> [user directory]"\
          "\nThe source directory corresponds to a pre-calibrated D0ARCA009 detector."\
          "\nThe user directory should be specified when this script is submitted in a batch queue."
fi

case $# in
    2)  set_variable HOMEDIR              $2;&
    1)  set_variable DIR                  ${1:a};;
    *)  fatal "Wrong number of arguments."
esac

source JAcousticsToolkit.sh


# post-calibration of input detector

mkdir -p $WORKDIR/post-calibration
pushd    $WORKDIR/post-calibration

post-calibration_D0ARCA009.sh  $DIR

popd

set_variable  DETECTOR               $WORKDIR/post-calibration/detector.datx
set_variable  TRIPOD                 $WORKDIR/post-calibration/tripod.txt
set_variable  DETECTOR_INITIAL       detector_initial.datx
set_variable  TRIPOD_INITIAL         tripod_initial.txt
expand_array  RUNS


# Run in HOMEDIR

mkdir -p $HOMEDIR
cd       $HOMEDIR

JDetectorDB \
    -D $DETECTOR_ID       \
    -r $RUNS[1]           \
    -V ""                 \
    -WW                   \
    -o $DETECTOR_INITIAL

cat>$TRIPOD_INITIAL<<EOF
# Very preliminary acoustic emitter positions from survey
7   +587198.628  +4016228.693  -3433.306
8   +588257.896  +4017034.442  -3457.240
9   +587646.419  +4017109.956  -3456.980
10  +587763.722  +4017253.398  -3453.894
12  +587600.000  +4017500.000  -3457.240
13  +587510.740  +4016869.160  -3451.700
EOF

JEditDetector \
    -a $DETECTOR           \
    -r "19"                \
    -o detector_tmp.datx --!

eval `JPrintDetector -a $DETECTOR_INITIAL -O SUMMARY`;  FINAL_STRINGS=($STRINGS[*])
eval `JPrintDetector -a detector_tmp.datx -O SUMMARY`;  FIXED_STRINGS=($STRINGS[*])

# Apply modifications.

JEditDetector \
    -a $DETECTOR_INITIAL     \
    -s "-1 mul $MUL"         \
    -o detector.datx --!

eval `JPrintDetector -a detector.datx -O SUMMARY`

for STRING in ${STRINGS[*]}; do

    set_variable  MODULE  `getModule -a detector.datx -L "$STRING 0"`

    JEditDetector \
        -a detector.datx         \
        -M "$MODULE setz -11.7"  \
        -o detector.datx --!
done

JEditDetector   -a detector.datx  -r "$FIXED_STRINGS[*]"  -o detector.datx --!
JMergeDetector  -a detector.datx  -a detector_tmp.datx    -o detector.datx --!
JEditDetector   -a detector.datx  -k "$FINAL_STRINGS[*]"  -o detector.datx --!

JEditDetector \
    -a detector.datx                 \
    -M "817802210 lower $((-60.0 * $PI / 180.0 ))" \
    -o detector.datx --!            

typeset -A TRIPODS

get_tripods  $TRIPOD_INITIAL  TRIPODS;  FINAL_TRIPODS=(${(k)TRIPODS[*]})
get_tripods  $TRIPOD          TRIPODS;  FIXED_TRIPODS=(${(k)TRIPODS[*]})

cp -p  $TRIPOD_INITIAL  tripod.txt

JEditTripod   -f tripod.txt  -r "${FIXED_TRIPODS}" -o tripod.txt
JMergeTripod  -f tripod.txt  -f $TRIPOD            -o tripod.txt
JEditTripod   -f tripod.txt  -k "${FINAL_TRIPODS}" -o tripod.txt

cat>acoustics_trigger_parameters.txt<<EOF
Q             = 0.0;
TMax_s        = 0.020;
numberOfHits  = 300;
EOF

cat>disable.txt<<EOF
# disable transmissions
# <emitter identifier> <module identifier>
-1 808981121
-1 808430572
7 808961261
EOF

cat>transmitter.txt<<EOF
# Data from D0ARCA021_baseAcousticBeacon.csv
11   16  0   0.03  -0.77   0.48
14    5  0   0.74   0.23   0.48
15   32  0  -0.38  -0.67   0.48
EOF

cat>hydrophone.txt<<EOF
# Data from D0ARCA021_hydrophones.csv
   9 -1  -0.57  -0.42   0.57
  10  0  -0.67   0.39   0.62
  11  0  -0.67   0.39   0.62
  13  0  -0.68  -0.37   0.62
  14  0  -0.20  -0.75   0.62
  12  0  -0.12   0.76   0.62
  15  0  -0.72  -0.27   0.62
  19  0  -0.16  -0.76   0.62
  16  0   0.07   0.77   0.62
  20  0   0.09  -0.77   0.62
  23  0  -0.59  -0.49   0.62
  24  0  -0.68  -0.37   0.62
  21  0  -0.46   0.62   0.62
  22  0  -0.26   0.73   0.62
  25  0  -0.74  -0.22   0.62
  30  0  -0.08  -0.77   0.62
  26  0  -0.06   0.77   0.62
  27  0  -0.21   0.74   0.62
  28  0   0.12   0.76   0.62
  32  0  -0.77   0.06   0.62
   5  0  -0.77   0.06   0.62
EOF

JAcousticsEventBuilder.sh  detector.datx  $RUNS[*]

INPUT_FILES=(`ls KM3NeT_${(l:8::0::0:)DETECTOR_ID}_0*${^RUNS}_event.root`)


# Run in WORKDIR

cd  $WORKDIR

if [ ! $HOMEDIR -ef $WORKDIR ]; then
    cp -p $HOMEDIR/$ACOUSTICS_DETECTOR     $WORKDIR
    cp -p $HOMEDIR/${^ACOUSTICS_KEYS}.txt  $WORKDIR
    cp -p $HOMEDIR/${^ACOUSTICS_AUXS}.txt  $WORKDIR
    cp -p $HOMEDIR/${^INPUT_FILES}         $WORKDIR
fi

cat>script.txt<<EOF

fix  string  $FIXED_STRINGS[*]
fix  tripod  $FIXED_TRIPODS[*]

stage   0  0  2  250.0e-6  10.0  20  0.5       
stage  1A  0  2  250.0e-6  10.0  20  0.5    1.0
stage  1B  0  2  250.0e-6  10.0   0  0.0005 0.2
stage  1C  2  0  250.0e-6  10.0  20  0.3       
stage  2A  2  0  100.0e-6  10.0  20  0.2    0.4
stage  2C  2  0  100.0e-6  10.0  20  0.2       

initialise

stage  2a  2  0  100.0e-6  10.0  20  0.2    0.4
stage  2c  2  0  100.0e-6  10.0  20  0.2       
stage  3a  2  0   50.0e-6  10.0  20  0.2    0.4
stage  3c  2  0   50.0e-6  10.0  20  0.2       
EOF

JSydney.sh  detector.datx  $INPUT_FILES[*]  script.txt

if [ ! $HOMEDIR -ef $WORKDIR ]; then
    cp -vru * $HOMEDIR
fi