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