Jpp
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
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 BACKGROUND ${SIRENE_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  JPlot1D \
40  -f "${INPUT_FILES//\.root/.root:job}" \
41  -\> "PDG code" \
42  -\^ "npe" \
43  -y "1e-2 1e5" -Y \
44  -S 0.8 -O P \
45  -T "Average number of photo-electrons per event by a track_in" \
46  -L TR \
47  -o $WORKDIR/hits_PDG.gif $BACKGROUND
48 fi
49 
50 if (( 1 )); then
51 
52  JPlot1D \
53  -f "${INPUT_FILES//\.root/.root:trks$}" \
54  -\> "PDG code" \
55  -\^ "number of tracks" \
56  -y "1e-3 1e1" -Y \
57  -S 0.8 -O P \
58  -T "Average occurrence of a track_in per event" \
59  -L TR \
60  -o $WORKDIR/trks_PDG.gif $BACKGROUND
61 fi
62 
63 if (( 1 )); then
64 
65  JPlot1D \
66  -f "${INPUT_FILES//\.root/.root:hits_per_E_in}" \
67  -y "0.5 1e8" -Y \
68  -XX \
69  -\> "E_{#nu} [GeV]" \
70  -\^ "npe" \
71  -S 0.8 \
72  -L TL \
73  -T "Average number of photo-elecrons per event inside the instrumented volume" \
74  -o $WORKDIR/hits_E.gif $BACKGROUND
75 fi
76 
77 if (( 1 )); then
78 
79  JPlot1D \
80  -f "${INPUT_FILES//\.root/.root:pmt}" \
81  -x "0.5 1e5" -X \
82  -y "1e-6 1e2" -Y \
83  -\> "npe" \
84  -\^ "number of hits [a.u.]" \
85  -L TR \
86  -T "number of photo-electrons per PMT" \
87  -o $WORKDIR/npe.gif $BACKGROUND
88 fi
89 
90 if (( 1 )); then
91 
92  set_variable TH2 nuExD
93  set_variable XMIN 2.0
94  set_variable XMAX 8.0
95  set_variable DX 0.5
96  set_variable ZMIN 1.0e-10
97  set_variable ZMAX 1.0e0
98 
99  for (( XL = $XMIN; $XL < $XMAX; XL += $DX )); do
100 
101  (( XR = $XL + $DX ))
102 
103  LABEL=`printf '10^{%3.1f} GeV < E_{#nu} < 10^{%3.1f} GeV' $XL $XR`
104 
105  JPlot1D \
106  -f "${INPUT_FILES//\.root/.root:${TH2}}" \
107  -P Y \
108  -x "$XL $XR" \
109  -y "0 750" \
110  -z "$ZMIN $ZMAX" -Z \
111  -\> "D [m]" \
112  -\^ "npe [a.u.]" \
113  -T "$LABEL" \
114  -L TR \
115  -o $WORKDIR/${TH2}_${XL}.gif -B
116  done
117 
118  montage \
119  -tile 4x3 \
120  -geometry +0+0 \
121  $WORKDIR/${TH2}_*.gif \ $WORKDIR/${TH2}.gif >& /dev/null
122 
123  rm -f $WORKDIR/${TH2}_*.gif
124 
125  if [[ "$BACKGROUND" != "-B" ]]; then
126  display $WORKDIR/${TH2}.gif
127  fi
128 fi
129 
130 if (( 1 )); then
131 
132  set_variable TH2 nuExc
133  set_variable XMIN 2.0
134  set_variable XMAX 8.0
135  set_variable DX 0.5
136  set_variable ZMIN 1.0e-7
137  set_variable ZMAX 1.0e-4
138 
139  for (( XL = $XMIN; $XL < $XMAX; XL += $DX )); do
140 
141  (( XR = $XL + $DX ))
142 
143  LABEL=`printf '10^{%3.1f} GeV < E_{#nu} < 10^{%3.1f} GeV' $XL $XR`
144 
145  JPlot1D \
146  -f "${INPUT_FILES//\.root/.root:${TH2}}" \
147  -P Y \
148  -x "$XL $XR" \
149  -z "$ZMIN $ZMAX" -Z \
150  -\> "cos(#theta_{0})" \
151  -\^ "npe [a.u.]" \
152  -T "$LABEL" \
153  -L TL \
154  -o $WORKDIR/${TH2}_${XL}.gif -B
155  done
156 
157  montage \
158  -tile 4x3 \
159  -geometry +0+0 \
160  $WORKDIR/${TH2}_*.gif \ $WORKDIR/${TH2}.gif >& /dev/null
161 
162  rm -f $WORKDIR/${TH2}_*.gif
163 
164  if [[ "$BACKGROUND" != "-B" ]]; then
165  display $WORKDIR/${TH2}.gif
166  fi
167 fi
168 
169 if (( 1 )); then
170 
171  set_variable TH2 nuDxc
172  set_variable XMIN 0.0
173  set_variable XMAX 240.0
174  set_variable DX 20
175  set_variable ZMIN 1.0e-1
176  set_variable ZMAX 1.0e+2
177 
178  typeset -Z3 CD
179 
180  for (( XL = $XMIN; $XL < $XMAX; XL += $DX )); do
181 
182  (( XR = $XL + $DX ))
183 
184  LABEL=`printf '%3.0f m < D < %3.0f m' $XL $XR`
185 
186  CD=`printf '%.0f' ${XL}`
187 
188  JPlot1D \
189  -f "${INPUT_FILES//\.root/.root:${TH2}}" \
190  -P Y \
191  -x "$XL $XR" \
192  -z "$ZMIN $ZMAX" -Z \
193  -\> "cos(#theta_{0})" \
194  -\^ "npe [a.u.]" \
195  -T "$LABEL" \
196  -L TL \
197  -o $WORKDIR/${TH2}_${CD}.gif -B
198  done
199 
200  montage \
201  -tile 4x3 \
202  -geometry +0+0 \
203  $WORKDIR/${TH2}_*.gif \ $WORKDIR/${TH2}.gif >& /dev/null
204 
205  rm -f $WORKDIR/${TH2}_*.gif
206 
207  if [[ "$BACKGROUND" != "-B" ]]; then
208  display $WORKDIR/${TH2}.gif
209  fi
210 fi
211 
212 if (( 1 )); then
213 
214  set_variable TH2 nuDxU
215  set_variable XMIN 0.0
216  set_variable XMAX 240.0
217  set_variable DX 20
218  set_variable ZMIN 1.0e-1
219  set_variable ZMAX 2.0e+1
220 
221  typeset -Z3 CD
222 
223  for (( XL = $XMIN; $XL < $XMAX; XL += $DX )); do
224 
225  (( XR = $XL + $DX ))
226 
227  LABEL=`printf '%3.0f m < D < %3.0f m' $XL $XR`
228 
229  CD=`printf '%.0f' ${XL}`
230 
231  JPlot1D \
232  -f "${INPUT_FILES//\.root/.root:${TH2}}" \
233  -P Y \
234  -x "$XL $XR" \
235  -z "$ZMIN $ZMAX" -Z \
236  -\> "cos(#theta_{1})" \
237  -\^ "npe [a.u.]" \
238  -T "$LABEL" \
239  -L TR \
240  -o $WORKDIR/${TH2}_${CD}.gif -B
241  done
242 
243  montage \
244  -tile 4x3 \
245  -geometry +0+0 \
246  $WORKDIR/${TH2}_*.gif \ $WORKDIR/${TH2}.gif >& /dev/null
247 
248  rm -f $WORKDIR/${TH2}_*.gif
249 
250  if [[ "$BACKGROUND" != "-B" ]]; then
251  display $WORKDIR/${TH2}.gif
252  fi
253 fi
254 
255 if (( 1 )); then
256 
257  set_variable TH2 nuDxT
258  set_variable XMIN 0.0
259  set_variable XMAX 240.0
260  set_variable DX 20
261  set_variable ZMIN 1.0e-3
262  set_variable ZMAX 2.0e+1
263 
264  typeset -Z3 CD
265 
266  for (( XL = $XMIN; $XL < $XMAX; XL += $DX )); do
267 
268  (( XR = $XL + $DX ))
269 
270  LABEL=`printf '%3.0f m < D < %3.0f m' $XL $XR`
271 
272  CD=`printf '%.0f' ${XL}`
273 
274  JPlot1D \
275  -f "${INPUT_FILES//\.root/.root:${TH2}}" \
276  -P Y \
277  -x "$XL $XR" \
278  -y "-50 +500" \
279  -z "$ZMIN $ZMAX" -Z \
280  -\> "#Deltat [ns]" \
281  -\^ "npe [a.u.]" \
282  -T "$LABEL" \
283  -L TR \
284  -o $WORKDIR/${TH2}_${CD}.gif -B
285  done
286 
287  montage \
288  -tile 4x3 \
289  -geometry +0+0 \
290  $WORKDIR/${TH2}_*.gif \ $WORKDIR/${TH2}.gif >& /dev/null
291 
292  rm -f $WORKDIR/${TH2}_*.gif
293 
294  if [[ "$BACKGROUND" != "-B" ]]; then
295  display $WORKDIR/${TH2}.gif
296  fi
297 
298 fi
299 
300 if (( 1 )); then
301 
302  set_variable TH2 muExR
303  set_variable XMIN 2.0
304  set_variable XMAX 8.0
305  set_variable DX 0.5
306  set_variable ZMIN 1.0e-5
307  set_variable ZMAX 5.0e0
308 
309  for (( XL = $XMIN; $XL < $XMAX; XL += $DX )); do
310 
311  (( XR = $XL + $DX ))
312 
313  LABEL=`printf '10^{%3.1f} GeV < E_{#mu} < 10^{%3.1f} GeV' $XL $XR`
314 
315  JPlot1D \
316  -f "${INPUT_FILES//\.root/.root:${TH2}}" \
317  -P Y \
318  -x "$XL $XR" \
319  -y "0 250" \
320  -z "$ZMIN $ZMAX" -Z \
321  -\> "R [m]" \
322  -\^ "npe [a.u.]" \
323  -T "$LABEL" \
324  -L TR \
325  -o $WORKDIR/${TH2}_${XL}.gif -B
326  done
327 
328  montage \
329  -tile 4x3 \
330  -geometry +0+0 \
331  $WORKDIR/${TH2}_*.gif \ $WORKDIR/${TH2}.gif >& /dev/null
332 
333  rm -f $WORKDIR/${TH2}_*.gif
334 
335  if [[ "$BACKGROUND" != "-B" ]]; then
336  display $WORKDIR/${TH2}.gif
337  fi
338 fi
339 
340 if (( 1 )); then
341 
342  set_variable TH2 muRxU
343  set_variable XMIN 0.0
344  set_variable XMAX 240.0
345  set_variable DX 20
346  set_variable ZMIN 1.0e-1
347  set_variable ZMAX 2.0e+1
348 
349  typeset -Z3 CD
350 
351  for (( XL = $XMIN; $XL < $XMAX; XL += $DX )); do
352 
353  (( XR = $XL + $DX ))
354 
355  LABEL=`printf '%3.0f m < R < %3.0f m' $XL $XR`
356 
357  CD=`printf '%.0f' ${XL}`
358 
359  JPlot1D \
360  -f "${INPUT_FILES//\.root/.root:${TH2}}" \
361  -P Y \
362  -x "$XL $XR" \
363  -z "$ZMIN $ZMAX" -Z \
364  -\> "cos(#theta_{1})" \
365  -\^ "npe [a.u.]" \
366  -T "$LABEL" \
367  -L TR \
368  -o $WORKDIR/${TH2}_${CD}.gif -B
369  done
370 
371  montage \
372  -tile 4x3 \
373  -geometry +0+0 \
374  $WORKDIR/${TH2}_*.gif \ $WORKDIR/${TH2}.gif >& /dev/null
375 
376  rm -f $WORKDIR/${TH2}_*.gif
377 
378  if [[ "$BACKGROUND" != "-B" ]]; then
379  display $WORKDIR/${TH2}.gif
380  fi
381 fi
382 
383 if (( 1 )); then
384 
385  set_variable TH2 muRxT
386  set_variable XMIN 0.0
387  set_variable XMAX 240.0
388  set_variable DX 20
389  set_variable ZMIN 1.0e-3
390  set_variable ZMAX 2.0e+1
391 
392  typeset -Z3 CD
393 
394  for (( XL = $XMIN; $XL < $XMAX; XL += $DX )); do
395 
396  (( XR = $XL + $DX ))
397 
398  LABEL=`printf '%3.0f m < R < %3.0f m' $XL $XR`
399 
400  CD=`printf '%.0f' ${XL}`
401 
402  JPlot1D \
403  -f "${INPUT_FILES//\.root/.root:${TH2}}" \
404  -P Y \
405  -x "$XL $XR" \
406  -y "-50 +500" \
407  -z "$ZMIN $ZMAX" -Z \
408  -\> "#Deltat [ns]" \
409  -\^ "npe [a.u.]" \
410  -T "$LABEL" \
411  -L TR \
412  -o $WORKDIR/${TH2}_${CD}.gif -B
413  done
414 
415  montage \
416  -tile 4x3 \
417  -geometry +0+0 \
418  $WORKDIR/${TH2}_*.gif \ $WORKDIR/${TH2}.gif >& /dev/null
419 
420  rm -f $WORKDIR/${TH2}_*.gif
421 
422  if [[ "$BACKGROUND" != "-B" ]]; then
423  display $WORKDIR/${TH2}.gif
424  fi
425 
426 fi
427 
428 
JSYSTEM::set_variable
void set_variable(const std::string &name, const std::string &value)
Set environment variable.
Definition: JSystemToolkit.hh:379
exit
exit
Definition: JPizza.sh:36
f
do set_array DAQHEADER JPrintDAQHeader f
Definition: JTuneHV.sh:74
Y
then fatal Wrong number of arguments fi set_variable STRING $argv[1] set_variable DETECTORXY_TXT $WORKDIR $DETECTORXY_TXT tail read X Y CHI2 RMS printf optimum n $X $Y $CHI2 $RMS awk v Y
Definition: detector-XY:fit1d.sh:32
echo
then echo
Definition: JDOMDAQDriver.sh:56
T
do set_variable OUTPUT_DIRECTORY $WORKDIR T
Definition: JCalibrateHeight.sh:61
SIRENE_DEBUG
then echo Variable JPP_DIR undefined exit fi source $JPP_DIR setenv sh $JPP_DIR zmodload zsh mathfunc set_variable DEBUG SIRENE_DEBUG
Definition: plot-Domino.sh:27
o
do $DIR JTransitTime o
Definition: JTransitTime.sh:44
usage
* usage
Definition: JLigier-local.sh:33
INPUT_FILES
then fatal Not enough arguments fi set_variable DETECTOR $argv[1] INPUT_FILES
Definition: JCalibrateStrings.sh:21
XMIN
then fatal Wrong number of arguments fi set_variable STRING $argv[1] set_variable DETECTOR_TXT $WORKDIR $DETECTOR_TXT tail read X Y Z RMS let XMIN
Definition: detector-XY:fit.sh:39
$argv
*set_variable DETECTOR $argv[1] set_variable INPUT_FILE $argv[2,-2] set_variable OUTPUT_FILE $argv[-1]
Definition: JDomino.sh:39
S
do set_variable SIGMA_NS set_variable OUTLIERS set_variable OUTPUT_FILE matrix[${ALPHA_DEG}\deg\] root $JPP JMatrixNZ a $DETECTOR f $INPUT_FILE o $OUTPUT_FILE S
Definition: JMatrixNZ.sh:56
$WORKDIR
then display $WORKDIR
Definition: plot-Domino.sh:127
version
version
Definition: JCalibratePMT.sh:7
WORKDIR
then usage $script< detector file >< inputfile > fi case set_variable WORKDIR
Definition: JLegolas.sh:28
null
esac $JPP_BIN JLogger sh $LOGGER until pgrep JGetMessage</dev/null > dev null
Definition: examples/JLogger/JLogger.sh:33
P
then $DIR JPlotNPE PDG P
Definition: JPlotNPE-PDG.sh:60
done
esac done
Definition: JAddHDE.sh:21
DEBUG
#define DEBUG(A)
Message macros.
Definition: JMessage.hh:62
Generated by   doxygen 1.8.5