pre-calibration_D0ARCA006.sh

Go to the documentation of this file.

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

# --------------------------------------------------------------------------------------------
#
# Steering script for the global-fit-of-global-fits procedure of the D0ARCA006 (75) 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            75
set_array     RUNS                   9700 9777 9884 9984 10000 10113 10200 10263
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 [user directory]"\
          "\nThe user directory should be specified when this script is submitted in a batch queue."
fi

case $# in
    1)  set_variable HOMEDIR              $1;;
    0)  ;;
    *)  fatal "Wrong number of arguments."
esac

source JAcousticsToolkit.sh

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       \
        -@ "tcal = A03166911" \
        -@ "pcal = A02980895" \
        -@ "rcal = A02980908" \
        -V ""                 \
        -WW                   \
        -o $DETECTOR_INITIAL

cat>$TRIPOD_INITIAL<<EOF
# See https://git.km3net.de/calibration/input_tables for the origin(s) of the values
2  +587441.700 +4017014.700 -3459.600
3  +587848.700 +4016749.700 -3459.600
4  +587657.170 +4017113.940 -3459.600
5  +587906.500 +4017094.300 -3459.600
EOF

# 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                 \
    -M "806481218 swap  9 11"        \
    -o detector.datx --!            

cp -p  $TRIPOD_INITIAL  tripod.txt

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

cat>disable.txt<<EOF
# disable transmissions
# <emitter identifier> <module identifier>
4 808430036
3 808430571
4 808447091
2 808964925
3 808964925
3 808978646
3 808474602
4 808974966
4 808957378
3 808447094
3 808470373
2 808961261
3 808976233
5 808976233
3 808985803
EOF

cat>transmitter.txt<<EOF
# $script
EOF

cat>hydrophone.txt<<EOF
# Data from D0ARCA006_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
  18  0  -0.70  -0.34   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
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
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