8# ------------------------------------------------------------------------------------------
10# Utility script to plot JDomino results.
12# ------------------------------------------------------------------------------------------
14if [ -z $JPP_DIR ]; then
15 echo "Variable JPP_DIR undefined."
19source $JPP_DIR/setenv.sh $JPP_DIR
23set_variable: DEBUG SIRENE_DEBUG 2
24set_variable: WORKDIR SIRENE_WORKDIR ./
25set_variable: FORMAT GRAPHICS_FORMAT gif
26set_variable+ BATCH GRAPHICS_BATCH -B
29 usage "$script (input file)+"
33 fatal "No input files."
41 -f "${INPUT_FILES//\.root/.root:job}" \
46 -T "Average number of photo-electrons per event by a track_in" \
48 -o $WORKDIR/hits_PDG.$FORMAT $BATCH
54 -f "${INPUT_FILES//\.root/.root:trks$}" \
56 -\^ "number of tracks" \
59 -T "Average occurrence of a track_in per event" \
61 -o $WORKDIR/trks_PDG.$FORMAT $BATCH
67 -f "${INPUT_FILES//\.root/.root:hits_per_E_in}" \
74 -T "Average number of photo-elecrons per event inside the instrumented volume" \
75 -o $WORKDIR/hits_E.$FORMAT $BATCH
81 -f "${INPUT_FILES//\.root/.root:pmt}" \
85 -\^ "number of hits [a.u.]" \
87 -T "number of photo-electrons per PMT" \
88 -o $WORKDIR/npe.$FORMAT $BATCH
93 set_variable TH2 nuExD
97 set_variable ZMIN 1.0e-10
98 set_variable ZMAX 1.0e0
100 for (( XL = $XMIN; $XL < $XMAX; XL += $DX )); do
104 LABEL=`printf '10^{%3.1f} GeV < E_{#nu} < 10^{%3.1f} GeV' $XL $XR`
107 -f "${INPUT_FILES//\.root/.root:${TH2}}" \
111 -z "$ZMIN $ZMAX" -Z \
116 -o $WORKDIR/${TH2}_${XL}.$FORMAT $BATCH
122 $WORKDIR/${TH2}_*.$FORMAT \
123 $WORKDIR/${TH2}.$FORMAT >& /dev/null
125 rm -f $WORKDIR/${TH2}_*.$FORMAT
127 if [[ "$BATCH" != "-B" ]]; then
128 display $WORKDIR/${TH2}.$FORMAT
134 set_variable TH2 nuExc
135 set_variable XMIN 2.0
136 set_variable XMAX 8.0
138 set_variable ZMIN 1.0e-7
139 set_variable ZMAX 1.0e-4
141 for (( XL = $XMIN; $XL < $XMAX; XL += $DX )); do
145 LABEL=`printf '10^{%3.1f} GeV < E_{#nu} < 10^{%3.1f} GeV' $XL $XR`
148 -f "${INPUT_FILES//\.root/.root:${TH2}}" \
151 -z "$ZMIN $ZMAX" -Z \
152 -> "cos(#theta_{0})" \
156 -o $WORKDIR/${TH2}_${XL}.$FORMAT $BATCH
162 $WORKDIR/${TH2}_*.$FORMAT \
163 $WORKDIR/${TH2}.$FORMAT >& /dev/null
165 rm -f $WORKDIR/${TH2}_*.$FORMAT
167 if [[ "$BATCH" != "-B" ]]; then
168 display $WORKDIR/${TH2}.$FORMAT
174 set_variable TH2 nuDxc
175 set_variable XMIN 0.0
176 set_variable XMAX 240.0
178 set_variable ZMIN 1.0e-1
179 set_variable ZMAX 1.0e+2
183 for (( XL = $XMIN; $XL < $XMAX; XL += $DX )); do
187 LABEL=`printf '%3.0f m < D < %3.0f m' $XL $XR`
189 CD=`printf '%.0f' ${XL}`
192 -f "${INPUT_FILES//\.root/.root:${TH2}}" \
195 -z "$ZMIN $ZMAX" -Z \
196 -> "cos(#theta_{0})" \
200 -o $WORKDIR/${TH2}_${CD}.$FORMAT $BATCH
206 $WORKDIR/${TH2}_*.$FORMAT \
207 $WORKDIR/${TH2}.$FORMAT >& /dev/null
209 rm -f $WORKDIR/${TH2}_*.$FORMAT
211 if [[ "$BATCH" != "-B" ]]; then
212 display $WORKDIR/${TH2}.$FORMAT
218 set_variable TH2 nuDxU
219 set_variable XMIN 0.0
220 set_variable XMAX 240.0
222 set_variable ZMIN 1.0e-1
223 set_variable ZMAX 2.0e+1
227 for (( XL = $XMIN; $XL < $XMAX; XL += $DX )); do
231 LABEL=`printf '%3.0f m < D < %3.0f m' $XL $XR`
233 CD=`printf '%.0f' ${XL}`
236 -f "${INPUT_FILES//\.root/.root:${TH2}}" \
239 -z "$ZMIN $ZMAX" -Z \
240 -> "cos(#theta_{1})" \
244 -o $WORKDIR/${TH2}_${CD}.$FORMAT $BATCH
250 $WORKDIR/${TH2}_*.$FORMAT \
251 $WORKDIR/${TH2}.$FORMAT >& /dev/null
253 rm -f $WORKDIR/${TH2}_*.$FORMAT
255 if [[ "$BATCH" != "-B" ]]; then
256 display $WORKDIR/${TH2}.$FORMAT
262 set_variable TH2 nuDxT
263 set_variable XMIN 0.0
264 set_variable XMAX 240.0
266 set_variable ZMIN 1.0e-3
267 set_variable ZMAX 2.0e+1
271 for (( XL = $XMIN; $XL < $XMAX; XL += $DX )); do
275 LABEL=`printf '%3.0f m < D < %3.0f m' $XL $XR`
277 CD=`printf '%.0f' ${XL}`
280 -f "${INPUT_FILES//\.root/.root:${TH2}}" \
284 -z "$ZMIN $ZMAX" -Z \
289 -o $WORKDIR/${TH2}_${CD}.$FORMAT $BATCH
295 $WORKDIR/${TH2}_*.$FORMAT \
296 $WORKDIR/${TH2}.$FORMAT >& /dev/null
298 rm -f $WORKDIR/${TH2}_*.$FORMAT
300 if [[ "$BATCH" != "-B" ]]; then
301 display $WORKDIR/${TH2}.$FORMAT
308 set_variable TH2 muExR
309 set_variable XMIN 2.0
310 set_variable XMAX 8.0
312 set_variable ZMIN 1.0e-5
313 set_variable ZMAX 5.0e0
315 for (( XL = $XMIN; $XL < $XMAX; XL += $DX )); do
319 LABEL=`printf '10^{%3.1f} GeV < E_{#mu} < 10^{%3.1f} GeV' $XL $XR`
322 -f "${INPUT_FILES//\.root/.root:${TH2}}" \
326 -z "$ZMIN $ZMAX" -Z \
331 -o $WORKDIR/${TH2}_${XL}.$FORMAT $BATCH
337 $WORKDIR/${TH2}_*.$FORMAT \
338 $WORKDIR/${TH2}.$FORMAT >& /dev/null
340 rm -f $WORKDIR/${TH2}_*.$FORMAT
342 if [[ "$BATCH" != "-B" ]]; then
343 display $WORKDIR/${TH2}.$FORMAT
349 set_variable TH2 muRxU
350 set_variable XMIN 0.0
351 set_variable XMAX 240.0
353 set_variable ZMIN 1.0e-1
354 set_variable ZMAX 2.0e+1
358 for (( XL = $XMIN; $XL < $XMAX; XL += $DX )); do
362 LABEL=`printf '%3.0f m < R < %3.0f m' $XL $XR`
364 CD=`printf '%.0f' ${XL}`
367 -f "${INPUT_FILES//\.root/.root:${TH2}}" \
370 -z "$ZMIN $ZMAX" -Z \
371 -> "cos(#theta_{1})" \
375 -o $WORKDIR/${TH2}_${CD}.$FORMAT $BATCH
381 $WORKDIR/${TH2}_*.$FORMAT \
382 $WORKDIR/${TH2}.$FORMAT >& /dev/null
384 rm -f $WORKDIR/${TH2}_*.$FORMAT
386 if [[ "$BATCH" != "-B" ]]; then
387 display $WORKDIR/${TH2}.$FORMAT
393 set_variable TH2 muRxT
394 set_variable XMIN 0.0
395 set_variable XMAX 240.0
397 set_variable ZMIN 1.0e-3
398 set_variable ZMAX 2.0e+1
402 for (( XL = $XMIN; $XL < $XMAX; XL += $DX )); do
406 LABEL=`printf '%3.0f m < R < %3.0f m' $XL $XR`
408 CD=`printf '%.0f' ${XL}`
411 -f "${INPUT_FILES//\.root/.root:${TH2}}" \
415 -z "$ZMIN $ZMAX" -Z \
420 -o $WORKDIR/${TH2}_${CD}.$FORMAT $BATCH
426 $WORKDIR/${TH2}_*.$FORMAT \
427 $WORKDIR/${TH2}.$FORMAT >& /dev/null
429 rm -f $WORKDIR/${TH2}_*.$FORMAT
431 if [[ "$BATCH" != "-B" ]]; then
432 display $WORKDIR/${TH2}.$FORMAT