pre-calibration_D0ORCA038.sh

Go to the documentation of this file.

#!/bin/zsh
script=${0##*/}
 
# --------------------------------------------------------------------------------------------
#
# Steering script for the global-fit-of-global-fits procedure of the D0ORCA038 (291) 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            291
set_array     RUNS                   29215-29222
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 D0ORCA033 detector, tripod from calibration archive
# MAKE SURE TO USE THE D0ORCA033_fixed BRANCH (NOT YET IN THE CALIBRATION ARCHIVE MAIN BRANCH)
getFile.sh $CALIBRATION_ARCHIVE 274 1 detector $WORKDIR/D0ORCA033
getFile.sh $CALIBRATION_ARCHIVE 274 1 tripod $WORKDIR/tripod_D0ORCA033
 
# retrieve D0ORCA033 mechanics, fit parameters, sound velocity
getFile.sh $JPP_DIR/inputs/ 274 1 mechanics $WORKDIR/mechanics
getFile.sh $JPP_DIR/inputs/ 274 1 acoustics_fit_parameters $WORKDIR/acoustics_fit_parameters
getFile.sh $JPP_DIR/inputs/ 274 1 sound_velocity $WORKDIR/sound_velocity
 
# Get a initial detector for D0ORCA038 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 D0ORCA038 input_tables
11 257208.80 4743611.20 -2440.19
12 257000.30 4743616.80 -2441.30
15 257222.80 4743333.20 -2438.49
16 256850.30 4743351.70 -2440.00
EOF
 
eval `JPrintDetector -a $DETECTOR_INITIAL -O SUMMARY`;  FINAL_STRINGS=($STRINGS[*])
eval `JPrintDetector -a D0ORCA033.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 D0ORCA033.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_D0ORCA033.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_D0ORCA033.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
 27  16
-28  16
EOF
 
# Changed to 500 based on investigation average number of hits per event in D0ORCA038 (around 550-600)
cat>acoustics_trigger_parameters.txt<<EOF
# Trigger parameters D0ORCA038
Q             = 0.0;
TMax_s        = 0.020;
numberOfHits  = 500;
EOF
 
# Updated disable file for D0ORCA038
cat>disable.txt<<EOF
# 78c0dfdb-6bd7-4d56-8d09-0774a88615c4$
# GIT 20.0.0-195-g190c9e876
# ROOT 6.32.02
# application JEditDisable
# command /pbs/throng/km3net/software/RHEL9/jpp/master/out//Linux/bin//JEditDisable -f new_disable.log -o new_disable.log -a detector.datx -q -d 2 --!
# namespace KM3NET
# system Linux cca013 5.14.0-611.49.1.el9_7.x86_64 #1 SMP PREEMPT_DYNAMIC Tue Apr 7 14:39:16 EDT 2026 x86_64
11 809161495
11 809175478
11 817317552
11 817331139
11 819737869
11 819740066
12 809161495
12 809521500
12 813709200
12 817320256
12 817355794
15 809076040
15 809161495
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 D0ORCA038 input_tables
   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
  72  0  -0.25   1.01   0.55
  73  0   0.46   0.93   0.55
  74  0   0.38   0.97   0.55
  48  0  -1.03  -0.16   0.55
  37  0  -1.01  -0.25   0.55
  26  0  -1.00   0.27   0.55
EOF
 
expand_array  RUNS
 
set_variable:  QUALITY_FACTOR  ACOUSTICS_QUALITYFACTOR  2.0e3
 
BUFFER=(`JRuns -D $DETECTOR_ID -F RUN`)
let "NUMBER_OF_ERRORS = 0"
for RUN in ${RUNS:*BUFFER}; do
    
    set_variable  INPUT_FILE   $WORKDIR/KM3NeT_${(l:8::0::0:)DETECTOR_ID}_${(l:8::0::0:)RUN}_toa.root
    set_variable  OUTPUT_FILE  $WORKDIR/KM3NeT_${(l:8::0::0:)DETECTOR_ID}_${(l:8::0::0:)RUN}_event.root
 
    if [[ ! -f $INPUT_FILE ]]; then
 
            JToA \
                -@ "detid  = `getDetector -D $DETECTOR_ID`"  \
                -@ "minrun = $RUN"                           \
                -@ "maxrun = $RUN"                           \
                -@ "QUALITYFACTOR >= $QUALITY_FACTOR"        \
                -o $INPUT_FILE                               \
                -d $DEBUG --!
 
        if (( $? != 0 )); then
            let "NUMBER_OF_ERRORS = $NUMBER_OF_ERRORS + 1"
        fi
    fi
 
    if [[ -f $INPUT_FILE ]] && [[ ! -f $OUTPUT_FILE ]]; then
 
            timer_start
 
            JAcousticsEventBuilder \
                -a $DETECTOR                                 \
                -f $INPUT_FILE                               \
                -T $WORKDIR/tripod.txt                       \
                -Y $WORKDIR/transmitter.txt                  \
            -V $WORKDIR/sound_velocity.txt               \
                -W $WORKDIR/waveform.txt                     \
                -@ $WORKDIR/acoustics_trigger_parameters.txt \
                -H $WORKDIR/hydrophone.txt                   \
                -o $OUTPUT_FILE                              \
                -d $DEBUG --!
 
        if (( $? != 0 )); then
            let "NUMBER_OF_ERRORS = $NUMBER_OF_ERRORS + 1"
        fi
 
            timer_stop
            timer_print
    fi
done
 
if (( $NUMBER_OF_ERRORS != 0 )); then
    error "Number of errors $NUMBER_OF_ERRORS"
fi
 
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 --!
 
 
# copy detector and tripod to align with initial ORCA033 and ROV measurements
cp -p $DETECTOR $WORKDIR/aligned_detector.datx
cp -p $TRIPOD $WORKDIR/aligned_tripod.txt
 
# Align with D0ORCA033
JAlignDetector \
    -a $WORKDIR/aligned_detector.datx \
    -b $WORKDIR/D0ORCA033.datx \
    -P "x y z phi tx ty" \
    -O \
    -T $WORKDIR/aligned_tripod.txt \
    -A 
 
#Align with ROV
JAlignDetector \
    -a $WORKDIR/aligned_detector.datx \
    -b $WORKDIR/$DETECTOR_INITIAL \
    -P "x y z" \
    -O \
    -T $WORKDIR/aligned_tripod.txt \
    -A 
 
# copy detector to apply z-adjustment
cp -p $WORKDIR/aligned_detector.datx $WORKDIR/aligned_detector_z_adjusted.datx
 
JAdjustZpos \
    -a $WORKDIR/aligned_detector_z_adjusted.datx \
    -f $WORKDIR/KM3NeT_00000291_*event.root \
    -o $WORKDIR/aligned_detector_z_adjusted.datx
 
if [ ! $HOMEDIR -ef $WORKDIR ]; then
    cp -vru * $HOMEDIR
fi