pre-calibration_D0ORCA033.sh

Go to the documentation of this file.

#!/bin/zsh
script=${0##*/}
 
# --------------------------------------------------------------------------------------------
#
# Steering script for the global-fit-of-global-fits procedure of the D0ORCA033 (274) detector.
#
# --------------------------------------------------------------------------------------------
 
if [ -z $JPP_DIR ]; then
    echo "Variable JPP_DIR undefined."
    exit
fi
 
if [ -z $CALIBRATION_ARCHIVE ]; then
    echo "Variable CALIBRATION_ARCHIVE undefined."
    exit
fi
 
source $JPP_DIR/setenv.sh $JPP_DIR >& /dev/null
 
set_variable  DETECTOR_ID            274
set_array     RUNS                   28672-28679
set_variable: THREADS                ACOUSTICS_THREADS       32
set_variable: CALIBRATION_ARCHIVE    CALIBRATION_ARCHIVE
set_variable  SINGULARITY            -
set_variable:  WORKDIR               ACOUSTICS_WORKDIR     ./
set_variable  HOMEDIR                $WORKDIR
 
 
if do_usage $*; then
    usage "$script"\
          "\nRequires variable CALIBRATION_ARCHIVE to be set for the user."
fi
 
source JAcousticsToolkit.sh
 
set_variable  DETECTOR               $WORKDIR/detector.datx
set_variable  TRIPOD                 $WORKDIR/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
 
# retrieve D2ORCA024 detector and tripod from calibration archive
# MAKE SURE TO USE THE D1ORCA024 BRANCH (D2ORCA024 IS NOT YET IN THE CALIBRATION ARCHIVE MAIN BRANCH)
getFile.sh $CALIBRATION_ARCHIVE 244 1 detector $WORKDIR/D2ORCA024
getFile.sh $CALIBRATION_ARCHIVE 244 1 tripod $WORKDIR/tripod_D2ORCA024
 
# retrieve D2ORCA024 mechanics, fit parameters, sound velocity
getFile.sh $JPP_DIR/inputs/ 244 1 mechanics $WORKDIR/mechanics
getFile.sh $JPP_DIR/inputs/ 244 1 acoustics_fit_parameters $WORKDIR/acoustics_fit_parameters
getFile.sh $JPP_DIR/inputs/ 244 1 sound_velocity $WORKDIR/sound_velocity
 
# Get a initial detector for D0ORCA033 from the JDetectorDB, to be optimised in the global-fit-of-global-fits procedure.
JDetectorDB \
    -D $DETECTOR_ID       \
    -r $RUNS[1]           \
    -V ""                 \
    -WW                   \
    -o $DETECTOR_INITIAL
 
# Edited tripod 11 and 15 initial positions to adjust for 3.5 m vertical height above sea floor
cat>$TRIPOD_INITIAL<<EOF
# Data from D0ORCA033_tripods.csv
11 257208.80 4743611.20 -2440.20
12 257000.30 4743616.80 -2441.30
15 257222.80 4743333.20 -2438.50 
EOF
 
#JEditDetector \
#    -a $DETECTOR           \
#    -r "1 27 28"                 \
#    -o detector_tmp.datx --!
 
eval `JPrintDetector -a $DETECTOR_INITIAL -O SUMMARY`;  FINAL_STRINGS=($STRINGS[*])
eval `JPrintDetector -a D2ORCA024.datx -O SUMMARY`;  FIXED_STRINGS=($STRINGS[*])
 
# Apply modifications.
 
JEditDetector \
    -a $DETECTOR_INITIAL \
    -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 -2.9"   \
        -o detector.datx --!
done
 
JEditDetector   -a detector.datx  -r "$FIXED_STRINGS[*]"  -o detector.datx --!
JMergeDetector  -a detector.datx  -a D2ORCA024.datx       -o detector.datx --!
JEditDetector   -a detector.datx  -k "$FINAL_STRINGS[*]"  -o detector.datx --!
 
typeset -A TRIPODS
 
get_tripods  $TRIPOD_INITIAL              TRIPODS;  FINAL_TRIPODS=(${(k)TRIPODS[*]})
get_tripods  $WORKDIR/tripod_D2ORCA024.txt   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 $WORKDIR/tripod_D2ORCA024.txt    -o tripod.txt
JEditTripod   -f tripod.txt  -k "${FINAL_TRIPODS}" -o tripod.txt
 
 
cat>waveform.txt<<EOF
# Data from D2ORCA024_waveforms.csv
 35   9
-36   9
 34  10
-35  10
 28  11
-29  11
 48  12
-49  12
 49  14
-50  14
 16  15
-17  15
EOF
 
# Changed to 450 based on investigation average number of hits per event in D0ORCA033 (around 500)
cat>acoustics_trigger_parameters.txt<<EOF
# Trigger parameters D0ORCA033
Q             = 0.0;
TMax_s        = 0.020;
numberOfHits  = 450;
EOF
 
# Updated disable file for D0ORCA033
cat>disable.txt<<EOF
# 1e50f729-5c90-4c3f-973b-3068e31580a5$
# GIT 19.1.0
# ROOT 6.32.02
# application JEditDisable
# command /pbs/throng/km3net/software/RHEL9/jpp/v19.1.0/out//Linux/bin//JEditDisable -f new_disable.txt -o new_disable.txt -a detector.datx -q -d 2 --!
# namespace KM3NET
# system Linux cca013 5.14.0-611.27.1.el9_7.x86_64 #1 SMP PREEMPT_DYNAMIC Tue Jan 20 05:28:20 EST 2026 x86_64
11 808489117
11 813524183
11 813524234
11 813700098
11 817301284
11 817317552
11 817331139
11 817352720
11 819740066
12 808972598
12 809521500
12 813709200
12 817355794
15 808981510
15 811760695
15 813672377
15 817295057
15 817301191
EOF
 
cat>transmitter.txt<<EOF
# Data from D0ORCA033
 9   27  0  -0.74  -0.64   0.71
10   41  0  -0.40  -0.89   0.71
14   50  0  -0.26  -0.94   0.71
EOF
 
cat>hydrophone.txt<<EOF
# Data from D1ORCA024_hydrophones.csv combined with new D0ORCA033 values
   2  0   0.02   1.04   0.55
   3  0   0.46   0.93   0.55
   9  0  -0.75  -0.72   0.55
  10 -1   0.11   1.04   0.55
  12  0   0.46   0.93   0.55
   4  0   0.20   1.02   0.55
  18  0   0.61   0.84   0.55
  28  0   0.69   0.78   0.55
  11  0  -0.07   1.04   0.55
  19  0   0.46   0.93   0.55
  20  0   0.61   0.84   0.55
  29  0   0.46   0.93   0.55
  30  0   0.46   0.93   0.55
  21  0   0.11   1.04   0.55
  17  0   0.46   0.93   0.55
  27  0   1.04   0.07   0.55
  31  0   0.75   0.72   0.55
  39  0   0.46   0.93   0.55
  40  0   0.46   0.93   0.55
  41  0   0.91   0.50   0.55
  42  0   0.29   1.00   0.55
  38  0   0.89  -0.54   0.55
   1  0  -0.25   1.01   0.55
  49  0   0.32  -0.99   0.55
  50  0   0.83   0.62   0.55
  51  0   0.54   0.89   0.55
  52  0   0.82   0.64   0.55
  53  0   0.43   0.95   0.55
  60  0  -0.16   1.03   0.55
  61  0  -0.42   0.95   0.55
  62  0   1.01   0.25   0.55
  63  0   0.20   1.02   0.55
  64 -1  -0.07   1.04   0.55
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  ${SINGULARITY}1  0  100.0e-6  10.0  10  0.5       
stage  1B  ${SINGULARITY}1  0  100.0e-6  10.0   0  0.001  0.1
stage  2A  ${SINGULARITY}1  0  100.0e-6  10.0  20  0.2    0.5
stage  2C  ${SINGULARITY}1  0  100.0e-6  10.0  20  0.2       
 
initialise
 
stage  2a  ${SINGULARITY}1  0  100.0e-6  10.0  20  0.2    0.5
stage  2c  ${SINGULARITY}1  0  100.0e-6  10.0  20  0.2       
stage  3a  ${SINGULARITY}1  0   50.0e-6  10.0  20  0.2    0.5
stage  3c  ${SINGULARITY}1  0   50.0e-6  10.0  20  0.2       
EOF
 
JSydney \
    -a $DETECTOR                               \
    -f "$INPUT_FILES[*]"                       \
    -T tripod.txt                              \
    -Y transmitter.txt                         \
    -V sound_velocity.txt                      \
    -M mechanics.txt                           \
    -@ acoustics_fit_parameters.txt            \
    -! disable.txt                             \
    -N $THREADS                                \
    -s script.txt                              \
    -d $DEBUG --!
 
if [ ! $HOMEDIR -ef $WORKDIR ]; then
    cp -vru * $HOMEDIR
fi