Jpp  15.0.0
the software that should make you happy
 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: 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 \ $WORKDIR/${TH2}.$FORMAT >& /dev/null
123 
124  rm -f $WORKDIR/${TH2}_*.$FORMAT
125 
126  if [[ "$BATCH" != "-B" ]]; then
127  display $WORKDIR/${TH2}.$FORMAT
128  fi
129 fi
130 
131 if (( 1 )); then
132 
133  set_variable TH2 nuExc
134  set_variable XMIN 2.0
135  set_variable XMAX 8.0
136  set_variable DX 0.5
137  set_variable ZMIN 1.0e-7
138  set_variable ZMAX 1.0e-4
139 
140  for (( XL = $XMIN; $XL < $XMAX; XL += $DX )); do
141 
142  (( XR = $XL + $DX ))
143 
144  LABEL=`printf '10^{%3.1f} GeV < E_{#nu} < 10^{%3.1f} GeV' $XL $XR`
145 
146  JPlot1D \
147  -f "${INPUT_FILES//\.root/.root:${TH2}}" \
148  -P Y \
149  -x "$XL $XR" \
150  -z "$ZMIN $ZMAX" -Z \
151  -\> "cos(#theta_{0})" \
152  -\^ "npe [a.u.]" \
153  -T "$LABEL" \
154  -L TL \
155  -o $WORKDIR/${TH2}_${XL}.$FORMAT $BATCH
156  done
157 
158  montage \
159  -tile 4x3 \
160  -geometry +0+0 \
161  $WORKDIR/${TH2}_*.$FORMAT \ $WORKDIR/${TH2}.$FORMAT >& /dev/null
162 
163  rm -f $WORKDIR/${TH2}_*.$FORMAT
164 
165  if [[ "$BATCH" != "-B" ]]; then
166  display $WORKDIR/${TH2}.$FORMAT
167  fi
168 fi
169 
170 if (( 1 )); then
171 
172  set_variable TH2 nuDxc
173  set_variable XMIN 0.0
174  set_variable XMAX 240.0
175  set_variable DX 20
176  set_variable ZMIN 1.0e-1
177  set_variable ZMAX 1.0e+2
178 
179  typeset -Z3 CD
180 
181  for (( XL = $XMIN; $XL < $XMAX; XL += $DX )); do
182 
183  (( XR = $XL + $DX ))
184 
185  LABEL=`printf '%3.0f m < D < %3.0f m' $XL $XR`
186 
187  CD=`printf '%.0f' ${XL}`
188 
189  JPlot1D \
190  -f "${INPUT_FILES//\.root/.root:${TH2}}" \
191  -P Y \
192  -x "$XL $XR" \
193  -z "$ZMIN $ZMAX" -Z \
194  -\> "cos(#theta_{0})" \
195  -\^ "npe [a.u.]" \
196  -T "$LABEL" \
197  -L TL \
198  -o $WORKDIR/${TH2}_${CD}.$FORMAT $BATCH
199  done
200 
201  montage \
202  -tile 4x3 \
203  -geometry +0+0 \
204  $WORKDIR/${TH2}_*.$FORMAT \ $WORKDIR/${TH2}.$FORMAT >& /dev/null
205 
206  rm -f $WORKDIR/${TH2}_*.$FORMAT
207 
208  if [[ "$BATCH" != "-B" ]]; then
209  display $WORKDIR/${TH2}.$FORMAT
210  fi
211 fi
212 
213 if (( 1 )); then
214 
215  set_variable TH2 nuDxU
216  set_variable XMIN 0.0
217  set_variable XMAX 240.0
218  set_variable DX 20
219  set_variable ZMIN 1.0e-1
220  set_variable ZMAX 2.0e+1
221 
222  typeset -Z3 CD
223 
224  for (( XL = $XMIN; $XL < $XMAX; XL += $DX )); do
225 
226  (( XR = $XL + $DX ))
227 
228  LABEL=`printf '%3.0f m < D < %3.0f m' $XL $XR`
229 
230  CD=`printf '%.0f' ${XL}`
231 
232  JPlot1D \
233  -f "${INPUT_FILES//\.root/.root:${TH2}}" \
234  -P Y \
235  -x "$XL $XR" \
236  -z "$ZMIN $ZMAX" -Z \
237  -\> "cos(#theta_{1})" \
238  -\^ "npe [a.u.]" \
239  -T "$LABEL" \
240  -L TR \
241  -o $WORKDIR/${TH2}_${CD}.$FORMAT $BATCH
242  done
243 
244  montage \
245  -tile 4x3 \
246  -geometry +0+0 \
247  $WORKDIR/${TH2}_*.$FORMAT \ $WORKDIR/${TH2}.$FORMAT >& /dev/null
248 
249  rm -f $WORKDIR/${TH2}_*.$FORMAT
250 
251  if [[ "$BATCH" != "-B" ]]; then
252  display $WORKDIR/${TH2}.$FORMAT
253  fi
254 fi
255 
256 if (( 1 )); then
257 
258  set_variable TH2 nuDxT
259  set_variable XMIN 0.0
260  set_variable XMAX 240.0
261  set_variable DX 20
262  set_variable ZMIN 1.0e-3
263  set_variable ZMAX 2.0e+1
264 
265  typeset -Z3 CD
266 
267  for (( XL = $XMIN; $XL < $XMAX; XL += $DX )); do
268 
269  (( XR = $XL + $DX ))
270 
271  LABEL=`printf '%3.0f m < D < %3.0f m' $XL $XR`
272 
273  CD=`printf '%.0f' ${XL}`
274 
275  JPlot1D \
276  -f "${INPUT_FILES//\.root/.root:${TH2}}" \
277  -P Y \
278  -x "$XL $XR" \
279  -y "-50 +500" \
280  -z "$ZMIN $ZMAX" -Z \
281  -\> "#Deltat [ns]" \
282  -\^ "npe [a.u.]" \
283  -T "$LABEL" \
284  -L TR \
285  -o $WORKDIR/${TH2}_${CD}.$FORMAT $BATCH
286  done
287 
288  montage \
289  -tile 4x3 \
290  -geometry +0+0 \
291  $WORKDIR/${TH2}_*.$FORMAT \ $WORKDIR/${TH2}.$FORMAT >& /dev/null
292 
293  rm -f $WORKDIR/${TH2}_*.$FORMAT
294 
295  if [[ "$BATCH" != "-B" ]]; then
296  display $WORKDIR/${TH2}.$FORMAT
297  fi
298 
299 fi
300 
301 if (( 1 )); then
302 
303  set_variable TH2 muExR
304  set_variable XMIN 2.0
305  set_variable XMAX 8.0
306  set_variable DX 0.5
307  set_variable ZMIN 1.0e-5
308  set_variable ZMAX 5.0e0
309 
310  for (( XL = $XMIN; $XL < $XMAX; XL += $DX )); do
311 
312  (( XR = $XL + $DX ))
313 
314  LABEL=`printf '10^{%3.1f} GeV < E_{#mu} < 10^{%3.1f} GeV' $XL $XR`
315 
316  JPlot1D \
317  -f "${INPUT_FILES//\.root/.root:${TH2}}" \
318  -P Y \
319  -x "$XL $XR" \
320  -y "0 250" \
321  -z "$ZMIN $ZMAX" -Z \
322  -\> "R [m]" \
323  -\^ "npe [a.u.]" \
324  -T "$LABEL" \
325  -L TR \
326  -o $WORKDIR/${TH2}_${XL}.$FORMAT $BATCH
327  done
328 
329  montage \
330  -tile 4x3 \
331  -geometry +0+0 \
332  $WORKDIR/${TH2}_*.$FORMAT \ $WORKDIR/${TH2}.$FORMAT >& /dev/null
333 
334  rm -f $WORKDIR/${TH2}_*.$FORMAT
335 
336  if [[ "$BATCH" != "-B" ]]; then
337  display $WORKDIR/${TH2}.$FORMAT
338  fi
339 fi
340 
341 if (( 1 )); then
342 
343  set_variable TH2 muRxU
344  set_variable XMIN 0.0
345  set_variable XMAX 240.0
346  set_variable DX 20
347  set_variable ZMIN 1.0e-1
348  set_variable ZMAX 2.0e+1
349 
350  typeset -Z3 CD
351 
352  for (( XL = $XMIN; $XL < $XMAX; XL += $DX )); do
353 
354  (( XR = $XL + $DX ))
355 
356  LABEL=`printf '%3.0f m < R < %3.0f m' $XL $XR`
357 
358  CD=`printf '%.0f' ${XL}`
359 
360  JPlot1D \
361  -f "${INPUT_FILES//\.root/.root:${TH2}}" \
362  -P Y \
363  -x "$XL $XR" \
364  -z "$ZMIN $ZMAX" -Z \
365  -\> "cos(#theta_{1})" \
366  -\^ "npe [a.u.]" \
367  -T "$LABEL" \
368  -L TR \
369  -o $WORKDIR/${TH2}_${CD}.$FORMAT $BATCH
370  done
371 
372  montage \
373  -tile 4x3 \
374  -geometry +0+0 \
375  $WORKDIR/${TH2}_*.$FORMAT \ $WORKDIR/${TH2}.$FORMAT >& /dev/null
376 
377  rm -f $WORKDIR/${TH2}_*.$FORMAT
378 
379  if [[ "$BATCH" != "-B" ]]; then
380  display $WORKDIR/${TH2}.$FORMAT
381  fi
382 fi
383 
384 if (( 1 )); then
385 
386  set_variable TH2 muRxT
387  set_variable XMIN 0.0
388  set_variable XMAX 240.0
389  set_variable DX 20
390  set_variable ZMIN 1.0e-3
391  set_variable ZMAX 2.0e+1
392 
393  typeset -Z3 CD
394 
395  for (( XL = $XMIN; $XL < $XMAX; XL += $DX )); do
396 
397  (( XR = $XL + $DX ))
398 
399  LABEL=`printf '%3.0f m < R < %3.0f m' $XL $XR`
400 
401  CD=`printf '%.0f' ${XL}`
402 
403  JPlot1D \
404  -f "${INPUT_FILES//\.root/.root:${TH2}}" \
405  -P Y \
406  -x "$XL $XR" \
407  -y "-50 +500" \
408  -z "$ZMIN $ZMAX" -Z \
409  -\> "#Deltat [ns]" \
410  -\^ "npe [a.u.]" \
411  -T "$LABEL" \
412  -L TR \
413  -o $WORKDIR/${TH2}_${CD}.$FORMAT $BATCH
414  done
415 
416  montage \
417  -tile 4x3 \
418  -geometry +0+0 \
419  $WORKDIR/${TH2}_*.$FORMAT \ $WORKDIR/${TH2}.$FORMAT >& /dev/null
420 
421  rm -f $WORKDIR/${TH2}_*.$FORMAT
422 
423  if [[ "$BATCH" != "-B" ]]; then
424  display $WORKDIR/${TH2}.$FORMAT
425  fi
426 
427 fi
428 
429 
$argv
then usage $script< detector file >< detectorfile > nIf the range of floors is the first detector file is aligned to the second before the comparison fi case set_variable RANGE $argv[3]
Definition: JCompareDetector.sh:31
f
then JMuonPostfit f
Definition: JMuonPostfit.sh:118
WORKDIR
then usage $script< detector file >< inputfile > nUtility script to create PDF and CDF of transition time distribution fi case set_variable WORKDIR
Definition: JLegolas.sh:31
exit
exit
Definition: JPizza.sh:36
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:33
echo
then echo
Definition: JDOMDAQDriver.sh:56
set_variable
then echo Variable JPP_DIR undefined exit fi source $JPP_DIR setenv sh $JPP_DIR &dev null set_variable
Definition: JAcoustics.sh:20
T
do set_variable OUTPUT_DIRECTORY $WORKDIR T
Definition: JCalibrateHeight.sh:61
X
then break fi done getCenter read X Y Z let X
Definition: acoustics-fit.sh:387
o
do $DIR JTransitTime o
Definition: JTransitTime.sh:44
usage
* usage
Definition: JLigier-local.sh:34
INPUT_FILES
set_array INPUT_FILES
Definition: JCalibrateStrings.sh:26
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:58
DEBUG
then usage $script< input_file >< detector_file > fi set_variable OUTPUT_DIR set_variable SELECTOR JDAQTimesliceL1 set_variable DEBUG case set_variable DEBUG
Definition: JCalibratePMT.sh:21
$WORKDIR
then display $WORKDIR
Definition: plot-Domino.sh:128
Z
do set_variable MODULE getModule a $WORKDIR detector_a datx L $STRING JEditDetector a $WORKDIR detector_a datx M $MODULE setz o $WORKDIR detector_a datx JEditDetector a $WORKDIR detector_b datx M $MODULE setz o $WORKDIR detector_b datx done echo Output stored at $WORKDIR detector_a datx and $WORKDIR tripod_a txt JDrawDetector2D a $WORKDIR detector_a datx a $WORKDIR detector_b datx L BL o detector $FORMAT $BATCH JDrawDetector2D T $WORKDIR tripod_a txt T $WORKDIR tripod_b txt L BL o tripod $FORMAT $BATCH JCompareDetector a $WORKDIR detector_a datx b $WORKDIR detector_b datx o $WORKDIR abc root &dev null for KEY in X Y Z
Definition: JAlignDetector.sh:81
XMIN
then let XMIN
Definition: detector-XY:fit.sh:67
done
set_variable NUMBER_OF_ITERATIONS set_variable EPSILON cat acoustics_fit_parameters txt<< EOF $CONFIGURATION[*]Nmin=3;sigma_s=100.0e-6;stdev=10.0;mestimator=0;fixStrings=0;EOF for STRING in $STRINGS[*];do#fit stretching and(z) position of given string set_variable DETECTOR_TMP/tmp/detector_A.datx JEditDetector-a $DETECTOR-o $DETECTOR_TMP-r $STRING JEditDetector-a $DETECTOR-o $DETECTOR-k $STRING for MUL in 0.005 0.001;do DX_M=0.2 for((N=0;$N< $NUMBER_OF_ITERATIONS;++N));do CHI2[3]=$CHI2[1] fitPositionOfString $STRING Z $DX_M fitStretchingOfString $STRING $MUL if(($CHI2[3]-$CHI2[1]< $EPSILON));then break fi done if(($N >=$NUMBER_OF_ITERATIONS));then printf"warning: reached maximum number of iterations %d - converenge %7.3f\n"$N $(($CHI2[3]-$CHI2[1])) fi done JMergeDetector-a $DETECTOR-a $DETECTOR_TMP-o $DETECTOR rm-f $DETECTOR_TMP JConvertDetectorFormat-a $DETECTOR-o $DETECTOR-r-d 0 > &dev null done
Definition: acoustics-fit.sh:427
version
version
Definition: JCalibratePMT.sh:7
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:62
Generated by   doxygen 1.8.5