Jpp 20.0.0-rc.2
the software that should make you happy
Loading...
Searching...
No Matches
pre-calibration_DBORCA018.sh
Go to the documentation of this file.
1#!/bin/zsh
2script=${0##*/}
3
4# --------------------------------------------------------------------------------------------
5#
6# Steering script for the global-fit-of-global-fits procedure of the D0ORCA018 (148) detector.
7# It includes ORCA.0027 and ORCA.0028. For use during the run period for when they were disentangled
8#
9# --------------------------------------------------------------------------------------------
10
11if [ -z $JPP_DIR ]; then
12 echo "Variable JPP_DIR undefined."
13 exit
14fi
15
16source $JPP_DIR/setenv.sh $JPP_DIR >& /dev/null
17
18set_variable DETECTOR_ID 148
19set_array RUNS 19441-19450
20set_variable: THREADS ACOUSTICS_THREADS 5
21set_variable SINGULARITY -
22set_variable WORKDIR `pwd`
23set_variable HOMEDIR $WORKDIR
24set_variable ACOUSTICS_TRANSMITTER
25
26if do_usage $*; then
27 usage "$script <source directory> [user directory]"\
28 "\nThe source directory corresponds to a pre-calibrated D1ORCA015 detector."\
29 "\nThe user directory should be specified when this script is submitted in a batch queue."
30fi
31
32case $# in
33 2) set_variable HOMEDIR $2;&
34 1) set_variable DIR ${1:a};;
35 *) fatal "Wrong number of arguments."
36esac
37
38source JAcousticsToolkit.sh
39
40set_variable DETECTOR $DIR/detector.datx
41set_variable TRIPOD $DIR/tripod.txt
42set_variable DETECTOR_INITIAL detector_initial.datx
43set_variable TRIPOD_INITIAL tripod_initial.txt
44expand_array RUNS
45
46
47# Run in HOMEDIR
48
49mkdir -p $HOMEDIR
50cd $HOMEDIR
51
52JDetectorDB \
53 -D $DETECTOR_ID \
54 -r $RUNS[1] \
55 -WW \
56 -o $DETECTOR_INITIAL
57
58cat>$TRIPOD_INITIAL<<EOF
59# position of emitter 5 taken from https://elog.km3net.de/Operations+FR/8271 (RAB1v2)
60# position of emitter 6 taken from https://elog.km3net.de/Operations+FR/8731 (RAB2)
61# position of emitter 8 taken from https://elog.km3net.de/Operations+FR/9603 (RAB3)
625 +257040.800 +4743317.600 -2439.900
636 +256807.500 +4743636.500 -2442.500
648 +257118.100 +4743781.100 -2440.300
65EOF
66
67JEditDetector \
68 -a $DETECTOR \
69 -o detector_tmp.datx --!
70
71eval `JPrintDetector -a $DETECTOR_INITIAL -O SUMMARY`; FINAL_STRINGS=($STRINGS[*])
72eval `JPrintDetector -a detector_tmp.datx -O SUMMARY`; FIXED_STRINGS=($STRINGS[*])
73
74# Apply modifications.
75
76JEditDetector \
77 -a $DETECTOR_INITIAL \
78 -o detector.datx --!
79
80eval `JPrintDetector -a detector.datx -O SUMMARY`
81
82#for STRING in ${STRINGS[*]}; do
83#
84# set_variable MODULE `getModule -a detector.datx -L "$STRING 0"`
85#
86# JEditDetector \
87# -a detector.datx \
88# -M "$MODULE setz -2.9" \
89# -o detector.datx --!
90#done
91
92JEditDetector -a detector.datx -r "$FIXED_STRINGS[*]" -o detector.datx --!
93JMergeDetector -a detector.datx -a detector_tmp.datx -o detector.datx --!
94JEditDetector -a detector.datx -k "$FINAL_STRINGS[*]" -o detector.datx --!
95
96JEditDetector \
97 -a detector.datx \
98 -M "817315169 swap 14 18" \
99 -o detector.datx --!
100
101typeset -A TRIPODS
102
103get_tripods $TRIPOD_INITIAL TRIPODS; FINAL_TRIPODS=(${(k)TRIPODS[*]})
104get_tripods $TRIPOD TRIPODS; FIXED_TRIPODS=(${(k)TRIPODS[*]})
105
106cp -p $TRIPOD_INITIAL tripod.txt
107
108JEditTripod -f tripod.txt -r "${FIXED_TRIPODS}" -o tripod.txt
109JMergeTripod -f tripod.txt -f $TRIPOD -o tripod.txt
110JEditTripod -f tripod.txt -k "${FINAL_TRIPODS}" -o tripod.txt
111
112cat>waveform.txt<<EOF
113# waveform identifier; emitter identifier
114 23 5
115-24 5
116 28 6
117-29 6
118 16 8
119-17 8
120 31 9
121-32 9
122EOF
123
124cat>acoustics_trigger_parameters.txt<<EOF
125Q = 0.0;
126TMax_s = 0.015;
127numberOfHits = 210;
128EOF
129
130cat>disable.txt<<EOF
131# disable transmissions
132# $36807518-43bb-453f-aa3c-75a68d558d7f$
133# GIT master_rocky-301-g22ecddb76
134# ROOT 6.32.02
135# application JEditDisable
136# command /sps/km3net/users/rclark/DynamicCalib/Jpp/out//Linux/bin//JEditDisable -f disable.log -o disable.log -a detector.datx -q -d 2 --!
137# namespace KM3NET
138# system Linux ccwslurm0360 5.14.0-427.37.1.el9_4.x86_64 #1 SMP PREEMPT_DYNAMIC Fri Sep 13 12:41:50 EDT 2024 x86_64
1391 809521500
1402 808472260
1412 808966194
1422 808977319
1433 808966194
1443 813553576
1455 806487231
1465 808978694
1475 816935869
1486 817301284
1496 817329086
1508 806483369
1518 806487231
1528 808964883
1538 817338744
1548 817605247
155EOF
156
157cat>transmitter.txt<<EOF
158# Data from D0ORCA018_baseAcousticBeacon.csv
159 9 27 0 -0.74 -0.64 0.71
160EOF
161
162cat>hydrophone.txt<<EOF
163# Data from D0ORCA018_hydrophones.csv
164 1 -1 -0.93 0.46 0.55
165 2 0 0.02 1.04 0.55
166 3 0 0.46 0.93 0.55
167 9 0 -0.75 -0.72 0.55
168 10 -1 0.11 1.04 0.55
169 12 0 0.46 0.93 0.55
170 4 0 0.20 1.02 0.55
171 18 0 0.61 0.84 0.55
172 28 0 0.69 0.78 0.55
173 11 0 -0.07 1.04 0.55
174 19 0 0.46 0.93 0.55
175 20 0 0.61 0.84 0.55
176 29 0 0.46 0.93 0.55
177 30 0 0.46 0.93 0.55
178 21 0 0.11 1.04 0.55
179 17 0 0.46 0.93 0.55
180 27 0 1.04 0.07 0.55
181 31 0 0.75 0.72 0.55
182EOF
183
184JAcousticsEventBuilder.sh detector.datx $RUNS[*]
185
186INPUT_FILES=(`ls KM3NeT_${(l:8::0::0:)DETECTOR_ID}_0*${^RUNS}_event.root`)
187
188
189# Run in WORKDIR
190
191cd $WORKDIR
192
193if [ ! $HOMEDIR -ef $WORKDIR ]; then
194 cp -p $HOMEDIR/$ACOUSTICS_DETECTOR $WORKDIR
195 cp -p $HOMEDIR/${^ACOUSTICS_KEYS}.txt $WORKDIR
196 cp -p $HOMEDIR/${^ACOUSTICS_AUXS}.txt $WORKDIR
197 cp -p $HOMEDIR/${^INPUT_FILES} $WORKDIR
198fi
199
200cat>script.txt<<EOF
201
202fix string $FIXED_STRINGS[*]
203fix tripod $FIXED_TRIPODS[*]
204
205stage 1B ${SINGULARITY}1 0 100.0e-6 10.0 0 0.001 0.1
206stage 2A ${SINGULARITY}1 0 100.0e-6 10.0 20 0.2 0.5
207stage 2C ${SINGULARITY}1 0 100.0e-6 10.0 20 0.2
208
209initialise
210
211stage 2a ${SINGULARITY}1 0 100.0e-6 10.0 20 0.2 0.5
212stage 2c ${SINGULARITY}1 0 100.0e-6 10.0 20 0.2
213stage 3a ${SINGULARITY}1 0 50.0e-6 10.0 20 0.2 0.5
214stage 3c ${SINGULARITY}1 0 50.0e-6 10.0 20 0.2
215EOF
216
217JSydney.sh detector.datx $INPUT_FILES[*] script.txt
218
219if [ ! $HOMEDIR -ef $WORKDIR ]; then
220 cp -vru * $HOMEDIR
221fi