pre-calibration_D0ARCA030.sh

Go to the documentation of this file.

#!/usr/bin/env zsh
script=${0##*/}
 
# --------------------------------------------------------------------------------------------
#
# Steering script for the global-fit-of-global-fits procedure of 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            232
set_array     RUNS                   22296-22310
set_variable: MUL                    ACOUSTICS_MUL          -1.21e-2
set_variable: THREADS                ACOUSTICS_THREADS       32
set_variable  SINGULARITY            -
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 D0ARCA028 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
 
set_variable  DETECTOR               $DIR/detector.datx
set_variable  TRIPOD                 $DIR/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
# Data from D0ARCA030_tripods.csv
9    +587646.419   +4017109.956   -3456.980
10   +587763.722   +4017253.398   -3453.894
13   +587510.740   +4016869.160   -3451.700
18   +587780.899   +4016532.953   -3454.403
19   +587955.229   +4016862.546   -3461.771
20   +587948.374   +4017173.149   -3461.153
21   +587241.776   +4016259.510   -3443.920
22   +587737.450   +4016020.395   -3446.611
EOF
 
eval `JPrintDetector -a $DETECTOR_INITIAL -O SUMMARY`;  FINAL_STRINGS=($STRINGS[*])
eval `JPrintDetector -a $DETECTOR         -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            -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  $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
 
print "FIXED STRINGS:" $FIXED_STRINGS
print "FINAL STRINGS:" $FINAL_STRINGS
print "FIXED TRIPODS:" $FIXED_TRIPODS
print "FINAL TRIPODS:" $FINAL_TRIPODS
 
 
cat>waveform.txt<<EOF
# waveform identifier; emitter identifier
 28   9
-29   9
 25  10
-26  10
 29  13
-30  13
 34  14
-35  14
 33  16
-34  16
 23  18
-24  18
 39  19
-40  19
 45  20
-46  20
 47  21
-48  21
 41  22
-42  22
EOF
 
cat>acoustics_trigger_parameters.txt<<EOF
Q             = 0.0;
TMax_s        = 0.035;
numberOfHits  = 420;
EOF
 
cat>disable.txt<<EOF
# disable transmissions
# $9136f45a-c2f3-450e-b450-ba8ea4497983$
# GIT 18.3.0-306-g8309c53d6-D
# ROOT 6.30/02
# application JEditDisable
# command /data1/dejong/km3net/GIT/Jpp/out//Linux/bin//JEditDisable -f disable.log -o disable.log -a detector.datx -q -d 2 --!
# namespace KM3NET
# system Linux verster.lorentz.leidenuniv.nl 6.6.8-200.fc39.x86_64 #1 SMP PREEMPT_DYNAMIC Thu Dec 21 04:01:49 UTC 2023 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 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 816919071
9 817595373
9 817605768
9 817608140
9 817802210
10 806649098
10 806657372
10 808430036
10 808437870
10 808468365
10 808474602
10 808964876
10 808974966
10 808977563
10 808984610
10 808984662
10 817606460
10 817608135
10 817616429
10 817802231
10 819852785
12 808454978
12 808474243
12 808491794
12 808982550
12 808995481
12 808998966
13 806649098
13 808430036
13 808430571
13 808447081
13 808474243
13 808474602
13 808480499
13 808952499
13 808964925
13 808974966
13 808981369
13 808984662
13 817802231
EOF
 
cat>transmitter.txt<<EOF
# Data from D0ARCA028_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
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
  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.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
 
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    0.5
stage  1B  0                2  250.0e-6  10.0   0  0.0005 0.2
stage  1C  ${SINGULARITY}2  0  250.0e-6  10.0  20  0.3       
stage  2A  ${SINGULARITY}2  0  100.0e-6  10.0  20  0.2    0.2
stage  2C  ${SINGULARITY}2  0  100.0e-6  10.0  20  0.2       
 
initialise
 
stage  2a  ${SINGULARITY}2  0  100.0e-6  10.0  20  0.2    0.2
stage  2c  ${SINGULARITY}2  0  100.0e-6  10.0  20  0.2       
stage  3a  ${SINGULARITY}2  0   50.0e-6  10.0  20  0.2    0.2
stage  3c  ${SINGULARITY}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