Jpp test-rotations-old
the software that should make you happy
Loading...
Searching...
No Matches
plot-Domino.sh
Go to the documentation of this file.
1#!/usr/bin/env zsh
2#
3# \author mdejong
4#
5version=1.0
6script=${0##*/}
7
8# ------------------------------------------------------------------------------------------
9#
10# Utility script to plot JDomino results.
11#
12# ------------------------------------------------------------------------------------------
13
14if [ -z $JPP_DIR ]; then
15 echo "Variable JPP_DIR undefined."
16 exit
17fi
18
19source $JPP_DIR/setenv.sh $JPP_DIR
20
21zmodload zsh/mathfunc
22
23set_variable: DEBUG SIRENE_DEBUG 2
24set_variable: WORKDIR SIRENE_WORKDIR ./
25set_variable: FORMAT GRAPHICS_FORMAT gif
26set_variable+ BATCH GRAPHICS_BATCH -B
27
28if do_usage $*; then
29 usage "$script (input file)+"
30fi
31
32if (( $# == 0 )); then
33 fatal "No input files."
34fi
35
36INPUT_FILES=($argv[*])
37
38if (( 1 )); then
39
40 JPlot1D \
41 -f "${INPUT_FILES//\.root/.root:job}" \
42 -> "PDG code" \
43 -\^ "npe" \
44 -y "1e-2 1e5" -Y \
45 -S 0.8 -O P \
46 -T "Average number of photo-electrons per event by a track_in" \
47 -L TR \
48 -o $WORKDIR/hits_PDG.$FORMAT $BATCH
49fi
50
51if (( 1 )); then
52
53 JPlot1D \
54 -f "${INPUT_FILES//\.root/.root:trks$}" \
55 -> "PDG code" \
56 -\^ "number of tracks" \
57 -y "1e-3 1e1" -Y \
58 -S 0.8 -O P \
59 -T "Average occurrence of a track_in per event" \
60 -L TR \
61 -o $WORKDIR/trks_PDG.$FORMAT $BATCH
62fi
63
64if (( 1 )); then
65
66 JPlot1D \
67 -f "${INPUT_FILES//\.root/.root:hits_per_E_in}" \
68 -y "0.5 1e8" -Y \
69 -XX \
70 -> "E_{#nu} [GeV]" \
71 -\^ "npe" \
72 -S 0.8 \
73 -L TL \
74 -T "Average number of photo-elecrons per event inside the instrumented volume" \
75 -o $WORKDIR/hits_E.$FORMAT $BATCH
76fi
77
78if (( 1 )); then
79
80 JPlot1D \
81 -f "${INPUT_FILES//\.root/.root:pmt}" \
82 -x "0.5 1e5" -X \
83 -y "1e-6 1e2" -Y \
84 -> "npe" \
85 -\^ "number of hits [a.u.]" \
86 -L TR \
87 -T "number of photo-electrons per PMT" \
88 -o $WORKDIR/npe.$FORMAT $BATCH
89fi
90
91if (( 1 )); then
92
93 set_variable TH2 nuExD
94 set_variable XMIN 2.0
95 set_variable XMAX 8.0
96 set_variable DX 0.5
97 set_variable ZMIN 1.0e-10
98 set_variable ZMAX 1.0e0
99
100 for (( XL = $XMIN; $XL < $XMAX; XL += $DX )); do
101
102 (( XR = $XL + $DX ))
103
104 LABEL=`printf '10^{%3.1f} GeV < E_{#nu} < 10^{%3.1f} GeV' $XL $XR`
105
106 JPlot1D \
107 -f "${INPUT_FILES//\.root/.root:${TH2}}" \
108 -P Y \
109 -x "$XL $XR" \
110 -y "0 750" \
111 -z "$ZMIN $ZMAX" -Z \
112 -> "D [m]" \
113 -\^ "npe [a.u.]" \
114 -T "$LABEL" \
115 -L TR \
116 -o $WORKDIR/${TH2}_${XL}.$FORMAT $BATCH
117 done
118
119 montage \
120 -tile 4x3 \
121 -geometry +0+0 \
122 $WORKDIR/${TH2}_*.$FORMAT \
123 $WORKDIR/${TH2}.$FORMAT >& /dev/null
124
125 rm -f $WORKDIR/${TH2}_*.$FORMAT
126
127 if [[ "$BATCH" != "-B" ]]; then
128 display $WORKDIR/${TH2}.$FORMAT
129 fi
130fi
131
132if (( 1 )); then
133
134 set_variable TH2 nuExc
135 set_variable XMIN 2.0
136 set_variable XMAX 8.0
137 set_variable DX 0.5
138 set_variable ZMIN 1.0e-7
139 set_variable ZMAX 1.0e-4
140
141 for (( XL = $XMIN; $XL < $XMAX; XL += $DX )); do
142
143 (( XR = $XL + $DX ))
144
145 LABEL=`printf '10^{%3.1f} GeV < E_{#nu} < 10^{%3.1f} GeV' $XL $XR`
146
147 JPlot1D \
148 -f "${INPUT_FILES//\.root/.root:${TH2}}" \
149 -P Y \
150 -x "$XL $XR" \
151 -z "$ZMIN $ZMAX" -Z \
152 -> "cos(#theta_{0})" \
153 -\^ "npe [a.u.]" \
154 -T "$LABEL" \
155 -L TL \
156 -o $WORKDIR/${TH2}_${XL}.$FORMAT $BATCH
157 done
158
159 montage \
160 -tile 4x3 \
161 -geometry +0+0 \
162 $WORKDIR/${TH2}_*.$FORMAT \
163 $WORKDIR/${TH2}.$FORMAT >& /dev/null
164
165 rm -f $WORKDIR/${TH2}_*.$FORMAT
166
167 if [[ "$BATCH" != "-B" ]]; then
168 display $WORKDIR/${TH2}.$FORMAT
169 fi
170fi
171
172if (( 1 )); then
173
174 set_variable TH2 nuDxc
175 set_variable XMIN 0.0
176 set_variable XMAX 240.0
177 set_variable DX 20
178 set_variable ZMIN 1.0e-1
179 set_variable ZMAX 1.0e+2
180
181 typeset -Z3 CD
182
183 for (( XL = $XMIN; $XL < $XMAX; XL += $DX )); do
184
185 (( XR = $XL + $DX ))
186
187 LABEL=`printf '%3.0f m < D < %3.0f m' $XL $XR`
188
189 CD=`printf '%.0f' ${XL}`
190
191 JPlot1D \
192 -f "${INPUT_FILES//\.root/.root:${TH2}}" \
193 -P Y \
194 -x "$XL $XR" \
195 -z "$ZMIN $ZMAX" -Z \
196 -> "cos(#theta_{0})" \
197 -\^ "npe [a.u.]" \
198 -T "$LABEL" \
199 -L TL \
200 -o $WORKDIR/${TH2}_${CD}.$FORMAT $BATCH
201 done
202
203 montage \
204 -tile 4x3 \
205 -geometry +0+0 \
206 $WORKDIR/${TH2}_*.$FORMAT \
207 $WORKDIR/${TH2}.$FORMAT >& /dev/null
208
209 rm -f $WORKDIR/${TH2}_*.$FORMAT
210
211 if [[ "$BATCH" != "-B" ]]; then
212 display $WORKDIR/${TH2}.$FORMAT
213 fi
214fi
215
216if (( 1 )); then
217
218 set_variable TH2 nuDxU
219 set_variable XMIN 0.0
220 set_variable XMAX 240.0
221 set_variable DX 20
222 set_variable ZMIN 1.0e-1
223 set_variable ZMAX 2.0e+1
224
225 typeset -Z3 CD
226
227 for (( XL = $XMIN; $XL < $XMAX; XL += $DX )); do
228
229 (( XR = $XL + $DX ))
230
231 LABEL=`printf '%3.0f m < D < %3.0f m' $XL $XR`
232
233 CD=`printf '%.0f' ${XL}`
234
235 JPlot1D \
236 -f "${INPUT_FILES//\.root/.root:${TH2}}" \
237 -P Y \
238 -x "$XL $XR" \
239 -z "$ZMIN $ZMAX" -Z \
240 -> "cos(#theta_{1})" \
241 -\^ "npe [a.u.]" \
242 -T "$LABEL" \
243 -L TR \
244 -o $WORKDIR/${TH2}_${CD}.$FORMAT $BATCH
245 done
246
247 montage \
248 -tile 4x3 \
249 -geometry +0+0 \
250 $WORKDIR/${TH2}_*.$FORMAT \
251 $WORKDIR/${TH2}.$FORMAT >& /dev/null
252
253 rm -f $WORKDIR/${TH2}_*.$FORMAT
254
255 if [[ "$BATCH" != "-B" ]]; then
256 display $WORKDIR/${TH2}.$FORMAT
257 fi
258fi
259
260if (( 1 )); then
261
262 set_variable TH2 nuDxT
263 set_variable XMIN 0.0
264 set_variable XMAX 240.0
265 set_variable DX 20
266 set_variable ZMIN 1.0e-3
267 set_variable ZMAX 2.0e+1
268
269 typeset -Z3 CD
270
271 for (( XL = $XMIN; $XL < $XMAX; XL += $DX )); do
272
273 (( XR = $XL + $DX ))
274
275 LABEL=`printf '%3.0f m < D < %3.0f m' $XL $XR`
276
277 CD=`printf '%.0f' ${XL}`
278
279 JPlot1D \
280 -f "${INPUT_FILES//\.root/.root:${TH2}}" \
281 -P Y \
282 -x "$XL $XR" \
283 -y "-50 +500" \
284 -z "$ZMIN $ZMAX" -Z \
285 -> "#Deltat [ns]" \
286 -\^ "npe [a.u.]" \
287 -T "$LABEL" \
288 -L TR \
289 -o $WORKDIR/${TH2}_${CD}.$FORMAT $BATCH
290 done
291
292 montage \
293 -tile 4x3 \
294 -geometry +0+0 \
295 $WORKDIR/${TH2}_*.$FORMAT \
296 $WORKDIR/${TH2}.$FORMAT >& /dev/null
297
298 rm -f $WORKDIR/${TH2}_*.$FORMAT
299
300 if [[ "$BATCH" != "-B" ]]; then
301 display $WORKDIR/${TH2}.$FORMAT
302 fi
303
304fi
305
306if (( 1 )); then
307
308 set_variable TH2 muExR
309 set_variable XMIN 2.0
310 set_variable XMAX 8.0
311 set_variable DX 0.5
312 set_variable ZMIN 1.0e-5
313 set_variable ZMAX 5.0e0
314
315 for (( XL = $XMIN; $XL < $XMAX; XL += $DX )); do
316
317 (( XR = $XL + $DX ))
318
319 LABEL=`printf '10^{%3.1f} GeV < E_{#mu} < 10^{%3.1f} GeV' $XL $XR`
320
321 JPlot1D \
322 -f "${INPUT_FILES//\.root/.root:${TH2}}" \
323 -P Y \
324 -x "$XL $XR" \
325 -y "0 250" \
326 -z "$ZMIN $ZMAX" -Z \
327 -> "R [m]" \
328 -\^ "npe [a.u.]" \
329 -T "$LABEL" \
330 -L TR \
331 -o $WORKDIR/${TH2}_${XL}.$FORMAT $BATCH
332 done
333
334 montage \
335 -tile 4x3 \
336 -geometry +0+0 \
337 $WORKDIR/${TH2}_*.$FORMAT \
338 $WORKDIR/${TH2}.$FORMAT >& /dev/null
339
340 rm -f $WORKDIR/${TH2}_*.$FORMAT
341
342 if [[ "$BATCH" != "-B" ]]; then
343 display $WORKDIR/${TH2}.$FORMAT
344 fi
345fi
346
347if (( 1 )); then
348
349 set_variable TH2 muRxU
350 set_variable XMIN 0.0
351 set_variable XMAX 240.0
352 set_variable DX 20
353 set_variable ZMIN 1.0e-1
354 set_variable ZMAX 2.0e+1
355
356 typeset -Z3 CD
357
358 for (( XL = $XMIN; $XL < $XMAX; XL += $DX )); do
359
360 (( XR = $XL + $DX ))
361
362 LABEL=`printf '%3.0f m < R < %3.0f m' $XL $XR`
363
364 CD=`printf '%.0f' ${XL}`
365
366 JPlot1D \
367 -f "${INPUT_FILES//\.root/.root:${TH2}}" \
368 -P Y \
369 -x "$XL $XR" \
370 -z "$ZMIN $ZMAX" -Z \
371 -> "cos(#theta_{1})" \
372 -\^ "npe [a.u.]" \
373 -T "$LABEL" \
374 -L TR \
375 -o $WORKDIR/${TH2}_${CD}.$FORMAT $BATCH
376 done
377
378 montage \
379 -tile 4x3 \
380 -geometry +0+0 \
381 $WORKDIR/${TH2}_*.$FORMAT \
382 $WORKDIR/${TH2}.$FORMAT >& /dev/null
383
384 rm -f $WORKDIR/${TH2}_*.$FORMAT
385
386 if [[ "$BATCH" != "-B" ]]; then
387 display $WORKDIR/${TH2}.$FORMAT
388 fi
389fi
390
391if (( 1 )); then
392
393 set_variable TH2 muRxT
394 set_variable XMIN 0.0
395 set_variable XMAX 240.0
396 set_variable DX 20
397 set_variable ZMIN 1.0e-3
398 set_variable ZMAX 2.0e+1
399
400 typeset -Z3 CD
401
402 for (( XL = $XMIN; $XL < $XMAX; XL += $DX )); do
403
404 (( XR = $XL + $DX ))
405
406 LABEL=`printf '%3.0f m < R < %3.0f m' $XL $XR`
407
408 CD=`printf '%.0f' ${XL}`
409
410 JPlot1D \
411 -f "${INPUT_FILES//\.root/.root:${TH2}}" \
412 -P Y \
413 -x "$XL $XR" \
414 -y "-50 +500" \
415 -z "$ZMIN $ZMAX" -Z \
416 -> "#Deltat [ns]" \
417 -\^ "npe [a.u.]" \
418 -T "$LABEL" \
419 -L TR \
420 -o $WORKDIR/${TH2}_${CD}.$FORMAT $BATCH
421 done
422
423 montage \
424 -tile 4x3 \
425 -geometry +0+0 \
426 $WORKDIR/${TH2}_*.$FORMAT \
427 $WORKDIR/${TH2}.$FORMAT >& /dev/null
428
429 rm -f $WORKDIR/${TH2}_*.$FORMAT
430
431 if [[ "$BATCH" != "-B" ]]; then
432 display $WORKDIR/${TH2}.$FORMAT
433 fi
434
435fi
436