back-propagation_D0ARCA028.sh

Go to the documentation of this file.

#!/bin/zsh
script=${0##*/}
 
# --------------------------------------------------------------------------------------------
#
# Steering script for the back-propagation procedure of the D0ARCA028 (160) detector from the
# D0ARCA030 (232) 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            160
set_array     RUNS                   18826-18833
set_variable: THREADS                ACOUSTICS_THREADS       8
set_variable: CALIBRATION_ARCHIVE    CALIBRATION_ARCHIVE     
set_variable  WORKDIR                `pwd`
set_variable  HOMEDIR                $WORKDIR
 
 
if do_usage $*; then
    usage "$script"\
          "\nRequires variable CALLIBRATION_ARCHIVE to be set for the user."
fi
 
source JAcousticsToolkit.sh
 
set_variable  DETECTOR               detector.datx
set_variable  TRIPOD_INITIAL         tripod_initial.txt
set_variable  TRIPOD_ABSOLUTE        tripod_absolute.txt
set_variable  TRIPOD                 tripod.txt
expand_array  RUNS
 
# Run in HOMEDIR
mkdir -p $HOMEDIR
cd       $HOMEDIR
 
 
# retrieve D0ARCA030 detector and tripod from calibration archive
getFile.sh $CALIBRATION_ARCHIVE 232 1 detector A_2.1.0 $WORKDIR/D0ARCA030_A_2.1.0
getFile.sh $CALIBRATION_ARCHIVE 232 1 tripod $WORKDIR/tripod_D0ARCA030
 
cp -p $WORKDIR/D0ARCA030_A_2.1.0.datx $WORKDIR/detector.datx
 
# define initial tripod for D0ARCA028 (to be optimised in the back-propagation procedure)
cat>$TRIPOD_INITIAL<<EOF
# Data from D0ARCA028_tripods.csv
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
13   +587510.740   +4016869.160   -3451.700
EOF
 
# Create TRIPOD for D0ARCA028 to be optimised in the back-propagation procedure.
typeset -A TRIPODS
 
JEditTripod   \
    -f $WORKDIR/tripod_D0ARCA030.txt   \
    -r "18 19 20 21 22"                \
    -o $TRIPOD_ABSOLUTE                \
 
get_tripods  $TRIPOD_INITIAL  TRIPODS;  FINAL_TRIPODS=(${(k)TRIPODS[*]})
get_tripods  $TRIPOD_ABSOLUTE 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_ABSOLUTE   -o tripod.txt
JEditTripod   -f tripod.txt  -k "${FINAL_TRIPODS}" -o tripod.txt
 
# Edit detector file to remove string from D0ARCA030 not present in D0ARCA028
JEditDetector \
    -a $DETECTOR           \
    -r "6 7"               \
    -@ "id=$DETECTOR_ID"   \
    -o $DETECTOR
 
# Create auxiliary files for D0ARCA028 back-propagation procedure.
cat>waveform.txt<<EOF
# waveform identifier; emitter identifier
 14  7
-15  7
 12  8
-13  8
 28  9
-29  9
 25  10
-26  10
 27  11
-28  11
 29  13
-30  13
 26  15
-27  15
 33  16
-34  16
 22  17
-23  17
EOF
 
cat>acoustics_trigger_parameters.txt<<EOF
Q             = 0.0;
TMax_s        = 0.020;
numberOfHits  = 400;
EOF
 
cat>disable.txt<<EOF
# disable transmissions
# -85c0-42be-8923-2b56d90c4ebf$
# GIT test-rotations-old-309-g7eec518ad-D
# ROOT 6.32.02
# application JEditDisable
# command /pbs/home/m/martin/Jpp_withVP/out//Linux/bin//JEditDisable -f disable.log -o disable.log -a detector.datx -q -d 2 --!
# namespace KM3NET
# system Linux ccwslurm2038 5.14.0-570.42.2.el9_6.x86_64 #1 SMP PREEMPT_DYNAMIC Mon Sep 8 12:15:13 EDT 2025 x86_64
2 808964925
3 808430571
3 808474602
3 808964925
3 808971936
3 808978646
3 808981369
4 808474602
4 808964925
7 806476519
7 806649098
7 808430571
7 808480499
7 808974928
7 808974966
7 808975010
7 808976351
7 808985071
7 816919166
7 816943531
7 817595446
7 817606460
7 819737799
8 806649098
8 808430571
8 808469556
8 808491794
8 808944706
8 808964757
8 808964845
8 808972687
8 808974966
8 808980070
8 808984610
8 817565754
8 817581876
8 817595428
8 817595516
8 817597135
8 817597212
8 817597230
8 817605775
8 817605786
8 817612559
8 817801419
9 806649098
9 806657372
9 808432849
9 808454978
9 808474602
9 808957378
9 808958576
9 808981206
9 808981369
9 808987098
9 809003775
9 816919071
9 817305565
9 817595373
9 817605768
9 817608140
9 817802210
10 806649098
10 806657372
10 808430036
10 808437870
10 808468365
10 808474602
10 808956926
10 808964876
10 808974966
10 808977563
10 808979220
10 808984610
10 808984662
10 817305565
10 817569715
10 817584003
10 817606460
10 817608135
10 817616429
10 817802231
10 819737799
10 819852785
12 808454978
12 808474243
12 808491794
12 808982550
12 808995481
12 808998966
13 806649098
13 808430036
13 808430571
13 808447081
13 808454978
13 808474243
13 808474602
13 808480499
13 808952499
13 808956926
13 808964925
13 808974966
13 808981369
13 808984662
13 817305565
13 817595519
13 817802231
EOF
 
cat>transmitter.txt<<EOF
# Data from D0ARCA028_baseAcousticBeacon.csv
15   32  0  -0.38  -0.67   0.48
16    8  0   0.44   0.63   0.48
17   18  0   0.12   0.76   0.48
EOF
 
cat>hydrophone.txt<<EOF
# Data from D0ARCA028_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
  12  0  -0.12   0.76   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.47   0.62   0.62
   5  0  -0.76  -0.13   0.62
   1  0  -0.20  -0.75   0.62
   3  0   0.12  -0.76   0.62
   4  0   0.62  -0.46   0.62
   8  0  -0.52  -0.57   0.62
  14  0  -0.20  -0.75   0.62
  15  0  -0.75  -0.20   0.62
  17  0   0.08   0.77   0.62
  18  0  -0.22  -0.74   0.62
  31  0   0.71   0.31   0.62
EOF
 
# Build acoustic events for D0ARCA028 pre-calibration runs
JAcousticsEventBuilder.sh  $DETECTOR  $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
 
# Do not fix string 9 and 31 as they were dead for D0ARCA030 but are active for D0ARCA028
fix  string  1 3 4 5 6 7 8 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 30 32
fix  tripod  13 9 10
 
stage  1A   -3  2  250.0e-6  10.0  20  0.5    0.5     
stage  2A   -3  2  100.0e-6  10.0  20  0.2    0.2
stage  3A   -3  2   50.0e-6  10.0  20  0.2    0.2
stage  4A    3  2   50.0e-6  10.0  20  0.2    0.2
 
EOF
 
JSydney.sh  $DETECTOR  $INPUT_FILES[*]  script.txt
 
# Apply creep correction to D0ARCA028 detector file
JKatoomba.sh  $DETECTOR  $INPUT_FILES[*]  katoomba.root
 
JCreep \
    -a $DETECTOR            \
    -f katoomba.root        \
    -A                      
 
if [ ! $HOMEDIR -ef $WORKDIR ]; then
    cp -vru * $HOMEDIR
fi