Jpp  19.1.0-rc.1
the software that should make you happy
pre-calibration_D0ARCA028.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 D0ARCA028 (160) 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 160
18 set_array RUNS 18826-18833
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 D0ARCA021 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_D0ARCA021.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 expand_array RUNS
53 
54 
55 # Run in HOMEDIR
56 
57 mkdir -p $HOMEDIR
58 cd $HOMEDIR
59 
60 JDetectorDB \
61  -D $DETECTOR_ID \
62  -r $RUNS[1] \
63  -V "" \
64  -WW \
65  -o $DETECTOR_INITIAL
66 
67 JEditDetector \
68  -a $DETECTOR \
69  -r "14 15" \
70  -o detector_tmp.datx --!
71 
72 eval `JPrintDetector -a $DETECTOR_INITIAL -O SUMMARY`; FINAL_STRINGS=($STRINGS[*])
73 eval `JPrintDetector -a detector_tmp.datx -O SUMMARY`; FIXED_STRINGS=($STRINGS[*])
74 
75 # Apply modifications.
76 
77 JEditDetector \
78  -a $DETECTOR_INITIAL \
79  -s "-1 mul $MUL" \
80  -o detector.datx --!
81 
82 eval `JPrintDetector -a detector.datx -O SUMMARY`
83 
84 for STRING in ${STRINGS[*]}; do
85 
86  set_variable MODULE `getModule -a detector.datx -L "$STRING 0"`
87 
88  JEditDetector \
89  -a detector.datx \
90  -M "$MODULE setz -11.7" \
91  -o detector.datx --!
92 done
93 
94 JEditDetector -a detector.datx -r "$FIXED_STRINGS[*]" -o detector.datx --!
95 JMergeDetector -a detector.datx -a detector_tmp.datx -o detector.datx --!
96 
97 JEditDetector \
98  -a detector.datx \
99  -M "817603901 lower $((+60.0 * $PI / 180.0 ))" \
100  -o detector.datx --!
101 
102 JEditTripod \
103  -f $TRIPOD \
104  -r "12" \
105  -o tripod.txt
106 
107 cat>acoustics_trigger_parameters.txt<<EOF
108 Q = 0.0;
109 TMax_s = 0.020;
110 numberOfHits = 400;
111 EOF
112 
113 cat>disable.txt<<EOF
114 # disable transmissions
115 # $9136f45a-c2f3-450e-b450-ba8ea4497983$
116 # GIT 18.3.0-306-g8309c53d6-D
117 # ROOT 6.30/02
118 # application JEditDisable
119 # command /data1/dejong/km3net/GIT/Jpp/out//Linux/bin//JEditDisable -f disable.log -o disable.log -a detector.datx -q -d 2 --!
120 # namespace KM3NET
121 # 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
122 2 808964925
123 3 808430571
124 3 808474602
125 3 808964925
126 3 808971936
127 3 808978646
128 3 808981369
129 4 808474602
130 4 808964925
131 7 806476519
132 7 806649098
133 7 808430571
134 7 808480499
135 7 808974928
136 7 808974966
137 7 808975010
138 7 808976351
139 7 808985071
140 7 816919166
141 7 816943531
142 7 817595446
143 7 817606460
144 7 819737799
145 8 806649098
146 8 808430571
147 8 808469556
148 8 808491794
149 8 808944706
150 8 808964757
151 8 808964845
152 8 808972687
153 8 808974966
154 8 808980070
155 8 808984610
156 8 817565754
157 8 817581876
158 8 817595516
159 8 817597135
160 8 817597212
161 8 817597230
162 8 817605775
163 8 817605786
164 8 817612559
165 8 817801419
166 9 806649098
167 9 806657372
168 9 808432849
169 9 808454978
170 9 808474602
171 9 808957378
172 9 808958576
173 9 808981206
174 9 808981369
175 9 808987098
176 9 816919071
177 9 817595373
178 9 817605768
179 9 817608140
180 9 817802210
181 10 806649098
182 10 806657372
183 10 808430036
184 10 808437870
185 10 808468365
186 10 808474602
187 10 808964876
188 10 808974966
189 10 808977563
190 10 808984610
191 10 808984662
192 10 817606460
193 10 817608135
194 10 817616429
195 10 817802231
196 10 819852785
197 12 808454978
198 12 808474243
199 12 808491794
200 12 808982550
201 12 808995481
202 12 808998966
203 13 806649098
204 13 808430036
205 13 808430571
206 13 808447081
207 13 808474243
208 13 808474602
209 13 808480499
210 13 808952499
211 13 808964925
212 13 808974966
213 13 808981369
214 13 808984662
215 13 817802231
216 EOF
217 
218 cat>transmitter.txt<<EOF
219 # Data from D0ARCA028_baseAcousticBeacon.csv
220 11 16 0 0.03 -0.77 0.48
221 14 5 0 0.74 0.23 0.48
222 15 32 0 -0.38 -0.67 0.48
223 16 8 0 0.44 0.63 0.48
224 17 18 0 0.12 0.76 0.48
225 EOF
226 
227 cat>hydrophone.txt<<EOF
228 # Data from D0ARCA028_hydrophones.csv
229  9 -1 -0.57 -0.42 0.57
230  10 0 -0.67 0.39 0.62
231  11 0 -0.67 0.39 0.62
232  13 0 -0.68 -0.37 0.62
233  14 0 -0.20 -0.75 0.62
234  12 0 -0.12 0.76 0.62
235  15 0 -0.72 -0.27 0.62
236  19 0 -0.16 -0.76 0.62
237  16 0 0.07 0.77 0.62
238  20 0 0.09 -0.77 0.62
239  23 0 -0.59 -0.49 0.62
240  24 0 -0.68 -0.37 0.62
241  21 0 -0.46 0.62 0.62
242  22 0 -0.26 0.73 0.62
243  25 0 -0.74 -0.22 0.62
244  30 0 -0.08 -0.77 0.62
245  26 0 -0.06 0.77 0.62
246  27 0 -0.21 0.74 0.62
247  28 0 0.12 0.76 0.62
248  32 0 0.47 0.62 0.62
249  5 0 -0.76 -0.13 0.62
250  1 0 -0.20 -0.75 0.62
251  3 0 0.12 -0.76 0.62
252  4 0 0.62 -0.46 0.62
253  8 0 -0.52 -0.57 0.62
254  14 0 -0.20 -0.75 0.62
255  15 0 -0.75 -0.20 0.62
256  17 0 0.08 0.77 0.62
257  18 0 -0.22 -0.74 0.62
258  31 0 0.71 0.31 0.62
259 EOF
260 
261 JAcousticsEventBuilder.sh detector.datx $RUNS[*]
262 
263 INPUT_FILES=(`ls KM3NeT_${(l:8::0::0:)DETECTOR_ID}_0*${^RUNS}_event.root`)
264 
265 
266 # Run in WORKDIR
267 
268 cd $WORKDIR
269 
270 if [ ! $HOMEDIR -ef $WORKDIR ]; then
271  cp -p $HOMEDIR/$ACOUSTICS_DETECTOR $WORKDIR
272  cp -p $HOMEDIR/${^ACOUSTICS_KEYS}.txt $WORKDIR
273  cp -p $HOMEDIR/${^ACOUSTICS_AUXS}.txt $WORKDIR
274  cp -p $HOMEDIR/${^INPUT_FILES} $WORKDIR
275 fi
276 
277 cat>script.txt<<EOF
278 
279 fix string $FIXED_STRINGS[*]
280 fix tripod $FIXED_TRIPODS[*]
281 
282 stage 1A 0 2 250.0e-6 10.0 20 0.5 0.5
283 stage 1B 0 2 250.0e-6 10.0 0 0.0005 0.2
284 stage 1C ${SINGULARITY}2 0 250.0e-6 10.0 20 0.3
285 stage 2A ${SINGULARITY}2 0 100.0e-6 10.0 20 0.2 0.2
286 stage 2C ${SINGULARITY}2 0 100.0e-6 10.0 20 0.2
287 
288 initialise
289 
290 stage 2a ${SINGULARITY}2 0 100.0e-6 10.0 20 0.2 0.2
291 stage 2c ${SINGULARITY}2 0 100.0e-6 10.0 20 0.2
292 stage 3a ${SINGULARITY}2 0 50.0e-6 10.0 20 0.2 0.2
293 stage 3c ${SINGULARITY}2 0 50.0e-6 10.0 20 0.2
294 EOF
295 
296 JSydney.sh detector.datx $INPUT_FILES[*] script.txt
297 
298 if [ ! $HOMEDIR -ef $WORKDIR ]; then
299  cp -vru * $HOMEDIR
300 fi
Generated by doxygen 1.9.1