Jpp  master_rocky-40-g5f0272dcd
the software that should make you happy
pre-calibration_D0ARCA021.sh
Go to the documentation of this file.
1 #!/bin/zsh
2 script=${0##*/}
3 
4 # --------------------------------------------------------------------------------------------
5 #
6 # Steering script for the global-fit-of-global-fits procedure of the D0ARCA021 (133) detector.
7 #
8 # --------------------------------------------------------------------------------------------
9 
10 if [ -z $JPP_DIR ]; then
11  echo "Variable JPP_DIR undefined."
12  exit
13 fi
14 
15 source $JPP_DIR/setenv.sh $JPP_DIR >& /dev/null
16 
17 set_variable DETECTOR_ID 133
18 set_array RUNS 13290-13297
19 set_variable: MUL ACOUSTICS_MUL -1.21e-2
20 set_variable: THREADS ACOUSTICS_THREADS 5
21 set_variable SINGULARITY -
22 set_variable WORKDIR `pwd`
23 set_variable HOMEDIR $WORKDIR
24 
25 if do_usage $*; then
26  usage "$script <source directory> [user directory]"\
27  "\nThe source directory corresponds to a pre-calibrated D0ARCA009 detector."\
28  "\nThe user directory should be specified when this script is submitted in a batch queue."
29 fi
30 
31 case $# in
32  2) set_variable HOMEDIR $2;&
33  1) set_variable DIR ${1:a};;
34  *) fatal "Wrong number of arguments."
35 esac
36 
37 source JAcousticsToolkit.sh
38 
39 
40 # post-calibration of input detector
41 
42 mkdir -p $WORKDIR/post-calibration
43 pushd $WORKDIR/post-calibration
44 
45 post-calibration_D0ARCA009.sh $DIR
46 
47 popd
48 
49 set_variable DETECTOR $WORKDIR/post-calibration/detector.datx
50 set_variable TRIPOD $WORKDIR/post-calibration/tripod.txt
51 set_variable DETECTOR_INITIAL detector_initial.datx
52 set_variable TRIPOD_INITIAL tripod_initial.txt
53 expand_array RUNS
54 
55 
56 # Run in HOMEDIR
57 
58 mkdir -p $HOMEDIR
59 cd $HOMEDIR
60 
61 JDetectorDB \
62  -D $DETECTOR_ID \
63  -r $RUNS[1] \
64  -V "" \
65  -WW \
66  -o $DETECTOR_INITIAL
67 
68 cat>$TRIPOD_INITIAL<<EOF
69 # Very preliminary acoustic emitter positions from survey
70 7 +587198.628 +4016228.693 -3433.306
71 8 +588257.896 +4017034.442 -3457.240
72 9 +587646.419 +4017109.956 -3456.980
73 10 +587763.722 +4017253.398 -3453.894
74 12 +587600.000 +4017500.000 -3457.240
75 13 +587510.740 +4016869.160 -3451.700
76 EOF
77 
78 JEditDetector \
79  -a $DETECTOR \
80  -r "19" \
81  -o detector_tmp.datx --!
82 
83 eval `JPrintDetector -a $DETECTOR_INITIAL -O SUMMARY`; FINAL_STRINGS=($STRINGS[*])
84 eval `JPrintDetector -a detector_tmp.datx -O SUMMARY`; FIXED_STRINGS=($STRINGS[*])
85 
86 # Apply modifications.
87 
88 JEditDetector \
89  -a $DETECTOR_INITIAL \
90  -s "-1 mul $MUL" \
91  -o detector.datx --!
92 
93 eval `JPrintDetector -a detector.datx -O SUMMARY`
94 
95 for STRING in ${STRINGS[*]}; do
96 
97  set_variable MODULE `getModule -a detector.datx -L "$STRING 0"`
98 
99  JEditDetector \
100  -a detector.datx \
101  -M "$MODULE setz -11.7" \
102  -o detector.datx --!
103 done
104 
105 JEditDetector -a detector.datx -r "$FIXED_STRINGS[*]" -o detector.datx --!
106 JMergeDetector -a detector.datx -a detector_tmp.datx -o detector.datx --!
107 JEditDetector -a detector.datx -k "$FINAL_STRINGS[*]" -o detector.datx --!
108 
109 JEditDetector \
110  -a detector.datx \
111  -M "817802210 lower $((-60.0 * $PI / 180.0 ))" \
112  -o detector.datx --!
113 
114 typeset -A TRIPODS
115 
116 get_tripods $TRIPOD_INITIAL TRIPODS; FINAL_TRIPODS=(${(k)TRIPODS[*]})
117 get_tripods $TRIPOD TRIPODS; FIXED_TRIPODS=(${(k)TRIPODS[*]})
118 
119 cp -p $TRIPOD_INITIAL tripod.txt
120 
121 JEditTripod -f tripod.txt -r "${FIXED_TRIPODS}" -o tripod.txt
122 JMergeTripod -f tripod.txt -f $TRIPOD -o tripod.txt
123 JEditTripod -f tripod.txt -k "${FINAL_TRIPODS}" -o tripod.txt
124 
125 cat>acoustics_trigger_parameters.txt<<EOF
126 Q = 0.0;
127 TMax_s = 0.020;
128 numberOfHits = 300;
129 EOF
130 
131 cat>disable.txt<<EOF
132 # disable transmissions
133 # $e807fd96-3473-48a0-a5d8-5cbe0cc1ebc4$
134 # GIT 18.3.0-305-g4e82b04f4-D
135 # ROOT 6.30/02
136 # application JEditDisable
137 # command /data1/dejong/km3net/GIT/Jpp/out//Linux/bin//JEditDisable -f disable.log -o disable.log -a detector.datx -q -d 2 --!
138 # namespace KM3NET
139 # 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
140 2 808470373
141 2 808961261
142 2 808964925
143 2 808976233
144 3 808430571
145 3 808447094
146 3 808451782
147 3 808470373
148 3 808474602
149 3 808961261
150 3 808964925
151 3 808971936
152 3 808976233
153 3 808978646
154 3 808981369
155 3 808985803
156 4 808474602
157 4 808961261
158 4 808964925
159 5 808447094
160 5 808451782
161 5 808976207
162 5 808976233
163 6 808447094
164 6 808964808
165 7 806476519
166 7 808961261
167 7 808974966
168 7 808976319
169 8 808430571
170 8 808956926
171 8 808964757
172 8 808974966
173 8 808980070
174 8 817612559
175 9 808454978
176 9 808474602
177 9 808958576
178 9 808981206
179 9 808981369
180 9 808987098
181 9 816919071
182 9 817595373
183 9 817802210
184 10 808430036
185 10 808437870
186 10 808956926
187 10 808964876
188 10 808971766
189 10 808977563
190 10 817616429
191 10 817802231
192 12 808454978
193 12 808474243
194 12 808491794
195 12 808982550
196 12 808995481
197 12 808998966
198 13 808430036
199 13 808430571
200 13 808447081
201 13 808451689
202 13 808470373
203 13 808474243
204 13 808474602
205 13 808952499
206 13 808964925
207 13 808974966
208 13 808976233
209 13 808976352
210 13 808981369
211 13 808984662
212 13 817802231
213 EOF
214 
215 cat>transmitter.txt<<EOF
216 # Data from D0ARCA021_baseAcousticBeacon.csv
217 11 16 0 0.03 -0.77 0.48
218 14 5 0 0.74 0.23 0.48
219 15 32 0 -0.38 -0.67 0.48
220 EOF
221 
222 cat>hydrophone.txt<<EOF
223 # Data from D0ARCA021_hydrophones.csv
224  9 -1 -0.57 -0.42 0.57
225  10 0 -0.67 0.39 0.62
226  11 0 -0.67 0.39 0.62
227  13 0 -0.68 -0.37 0.62
228  14 0 -0.20 -0.75 0.62
229  12 0 -0.12 0.76 0.62
230  15 0 -0.72 -0.27 0.62
231  19 0 -0.16 -0.76 0.62
232  16 0 0.07 0.77 0.62
233  20 0 0.09 -0.77 0.62
234  23 0 -0.59 -0.49 0.62
235  24 0 -0.68 -0.37 0.62
236  21 0 -0.46 0.62 0.62
237  22 0 -0.26 0.73 0.62
238  25 0 -0.74 -0.22 0.62
239  30 0 -0.08 -0.77 0.62
240  26 0 -0.06 0.77 0.62
241  27 0 -0.21 0.74 0.62
242  28 0 0.12 0.76 0.62
243  32 0 -0.77 0.06 0.62
244  5 0 -0.77 0.06 0.62
245 EOF
246 
247 JAcousticsEventBuilder.sh detector.datx $RUNS[*]
248 
249 INPUT_FILES=(`ls KM3NeT_${(l:8::0::0:)DETECTOR_ID}_0*${^RUNS}_event.root`)
250 
251 
252 # Run in WORKDIR
253 
254 cd $WORKDIR
255 
256 if [ ! $HOMEDIR -ef $WORKDIR ]; then
257  cp -p $HOMEDIR/$ACOUSTICS_DETECTOR $WORKDIR
258  cp -p $HOMEDIR/${^ACOUSTICS_KEYS}.txt $WORKDIR
259  cp -p $HOMEDIR/${^ACOUSTICS_AUXS}.txt $WORKDIR
260  cp -p $HOMEDIR/${^INPUT_FILES} $WORKDIR
261 fi
262 
263 cat>script.txt<<EOF
264 
265 fix string $FIXED_STRINGS[*]
266 fix tripod $FIXED_TRIPODS[*]
267 
268 stage 0 0 2 250.0e-6 10.0 20 0.5
269 stage 1A 0 2 250.0e-6 10.0 20 0.5 0.5
270 stage 1B 0 2 250.0e-6 10.0 0 0.0005 0.2
271 stage 1C ${SINGULARITY}2 0 250.0e-6 10.0 20 0.3
272 stage 2A ${SINGULARITY}2 0 100.0e-6 10.0 20 0.2 0.2
273 stage 2C ${SINGULARITY}2 0 100.0e-6 10.0 20 0.2
274 
275 initialise
276 
277 stage 2a ${SINGULARITY}2 0 100.0e-6 10.0 20 0.2 0.2
278 stage 2c ${SINGULARITY}2 0 100.0e-6 10.0 20 0.2
279 stage 3a ${SINGULARITY}2 0 50.0e-6 10.0 20 0.2 0.2
280 stage 3c ${SINGULARITY}2 0 50.0e-6 10.0 20 0.2
281 EOF
282 
283 JSydney.sh detector.datx $INPUT_FILES[*] script.txt
284 
285 if [ ! $HOMEDIR -ef $WORKDIR ]; then
286  cp -vru * $HOMEDIR
287 fi
Generated by doxygen 1.9.1