Jpp  master_rocky
the software that should make you happy
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 SIRENE_DEBUG 2
24 set_variable: WORKDIR SIRENE_WORKDIR ./
25 set_variable: FORMAT GRAPHICS_FORMAT gif
26 set_variable+ BATCH GRAPHICS_BATCH -B
27 
28 if do_usage $*; then
29  usage "$script (input file)+"
30 fi
31 
32 if (( $# == 0 )); then
33  fatal "No input files."
34 fi
35 
36 INPUT_FILES=($argv[*])
37 
38 if (( 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
49 fi
50 
51 if (( 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
62 fi
63 
64 if (( 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
76 fi
77 
78 if (( 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
89 fi
90 
91 if (( 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
130 fi
131 
132 if (( 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
170 fi
171 
172 if (( 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
214 fi
215 
216 if (( 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
258 fi
259 
260 if (( 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 
304 fi
305 
306 if (( 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
345 fi
346 
347 if (( 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
389 fi
390 
391 if (( 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 
435 fi
436