8 # ------------------------------------------------------------------------------------------
10 # Utility script to plot JDomino results.
12 # ------------------------------------------------------------------------------------------
14 if [ -z $JPP_DIR ]; then
15 echo "Variable JPP_DIR undefined."
19 source $JPP_DIR/setenv.sh $JPP_DIR
23 set_variable DEBUG ${DEBUG:-2} #
24 set_variable WORKDIR ${WORKDIR:-./} #
25 set_variable BACKGROUND ${BACKGROUND:-} # -B
27 if ( do_usage $* ); then
28 usage "$script (input file)+"
31 if (( $# == 0 )); then
32 fatal "No input files."
35 INPUT_FILES=($argv[*])
39 HISTOGRAMS=(`echo ${INPUT_FILES} | eval sed 's/\.root/\.root:job/g'`)
42 -f "${HISTOGRAMS[*]}" \
47 -T "Average number of photo-electrons per event by a track_in" \
49 -o $WORKDIR/hits_PDG.gif \
55 HISTOGRAMS=(`echo ${INPUT_FILES} | eval sed 's/\.root/\.root:trks\\$/g'`)
58 -f "${HISTOGRAMS[*]}" \
60 -\^ "number of tracks" \
63 -T "Average occurrence of a track_in per event" \
65 -o $WORKDIR/trks_PDG.gif \
71 HISTOGRAMS=(`echo ${INPUT_FILES} | eval sed 's/\.root/\.root:hits_per_E_in/g'`)
74 -f "${HISTOGRAMS[*]}" \
81 -T "Average number of photo-elecrons per event inside the instrumented volume" \
82 -o $WORKDIR/hits_E.gif \
88 HISTOGRAMS=(`echo ${INPUT_FILES} | eval sed 's/\.root/\.root:pmt/g'`)
91 -f "${HISTOGRAMS[*]}" \
95 -\^ "number of hits [a.u.]" \
97 -T "number of photo-electrons per PMT" \
104 set_variable TH2 nuExD
105 set_variable XMIN 2.0
106 set_variable XMAX 8.0
108 set_variable ZMIN 1.0e-10
109 set_variable ZMAX 1.0e0
111 HISTOGRAMS=(`echo ${INPUT_FILES} | eval sed 's/\.root/\.root:${TH2}/g'`)
113 for (( XL = $XMIN; $XL < $XMAX; XL += $DX )); do
117 LABEL=`printf '10^{%3.1f} GeV < E_{#nu} < 10^{%3.1f} GeV' $XL $XR`
120 -f "${HISTOGRAMS[*]}" \
124 -z "$ZMIN $ZMAX" -Z \
129 -o $WORKDIR/${TH2}_${XL}.gif \
136 $WORKDIR/${TH2}_*.gif \
139 rm -f $WORKDIR/${TH2}_*.gif
141 if [[ "$BACKGROUND" != "-B" ]]; then
142 display $WORKDIR/${TH2}.gif
148 set_variable TH2 nuExc
149 set_variable XMIN 2.0
150 set_variable XMAX 8.0
152 set_variable ZMIN 1.0e-7
153 set_variable ZMAX 1.0e-4
155 HISTOGRAMS=(`echo ${INPUT_FILES} | eval sed 's/\.root/\.root:${TH2}/g'`)
157 for (( XL = $XMIN; $XL < $XMAX; XL += $DX )); do
161 LABEL=`printf '10^{%3.1f} GeV < E_{#nu} < 10^{%3.1f} GeV' $XL $XR`
164 -f "${HISTOGRAMS[*]}" \
167 -z "$ZMIN $ZMAX" -Z \
168 -> "cos(#theta_{0})" \
172 -o $WORKDIR/${TH2}_${XL}.gif \
179 $WORKDIR/${TH2}_*.gif \
182 rm -f $WORKDIR/${TH2}_*.gif
184 if [[ "$BACKGROUND" != "-B" ]]; then
185 display $WORKDIR/${TH2}.gif
191 set_variable TH2 nuDxc
192 set_variable XMIN 0.0
193 set_variable XMAX 240.0
195 set_variable ZMIN 1.0e-1
196 set_variable ZMAX 1.0e+2
198 HISTOGRAMS=(`echo ${INPUT_FILES} | eval sed 's/\.root/\.root:${TH2}/g'`)
202 for (( XL = $XMIN; $XL < $XMAX; XL += $DX )); do
206 LABEL=`printf '%3.0f m < D < %3.0f m' $XL $XR`
208 CD=`printf '%.0f' ${XL}`
211 -f "${HISTOGRAMS[*]}" \
214 -z "$ZMIN $ZMAX" -Z \
215 -> "cos(#theta_{0})" \
219 -o $WORKDIR/${TH2}_${CD}.gif \
226 $WORKDIR/${TH2}_*.gif \
229 rm -f $WORKDIR/${TH2}_*.gif
231 if [[ "$BACKGROUND" != "-B" ]]; then
232 display $WORKDIR/${TH2}.gif
238 set_variable TH2 nuDxU
239 set_variable XMIN 0.0
240 set_variable XMAX 240.0
242 set_variable ZMIN 1.0e-1
243 set_variable ZMAX 2.0e+1
245 HISTOGRAMS=(`echo ${INPUT_FILES} | eval sed 's/\.root/\.root:${TH2}/g'`)
249 for (( XL = $XMIN; $XL < $XMAX; XL += $DX )); do
253 LABEL=`printf '%3.0f m < D < %3.0f m' $XL $XR`
255 CD=`printf '%.0f' ${XL}`
258 -f "${HISTOGRAMS[*]}" \
261 -z "$ZMIN $ZMAX" -Z \
262 -> "cos(#theta_{1})" \
266 -o $WORKDIR/${TH2}_${CD}.gif \
273 $WORKDIR/${TH2}_*.gif \
276 rm -f $WORKDIR/${TH2}_*.gif
278 if [[ "$BACKGROUND" != "-B" ]]; then
279 display $WORKDIR/${TH2}.gif
285 set_variable TH2 nuDxT
286 set_variable XMIN 0.0
287 set_variable XMAX 240.0
289 set_variable ZMIN 1.0e-3
290 set_variable ZMAX 2.0e+1
292 HISTOGRAMS=(`echo ${INPUT_FILES} | eval sed 's/\.root/\.root:${TH2}/g'`)
296 for (( XL = $XMIN; $XL < $XMAX; XL += $DX )); do
300 LABEL=`printf '%3.0f m < D < %3.0f m' $XL $XR`
302 CD=`printf '%.0f' ${XL}`
305 -f "${HISTOGRAMS[*]}" \
309 -z "$ZMIN $ZMAX" -Z \
314 -o $WORKDIR/${TH2}_${CD}.gif \
321 $WORKDIR/${TH2}_*.gif \
324 rm -f $WORKDIR/${TH2}_*.gif
326 if [[ "$BACKGROUND" != "-B" ]]; then
327 display $WORKDIR/${TH2}.gif
334 set_variable TH2 muExR
335 set_variable XMIN 2.0
336 set_variable XMAX 8.0
338 set_variable ZMIN 1.0e-5
339 set_variable ZMAX 5.0e0
341 HISTOGRAMS=(`echo ${INPUT_FILES} | eval sed 's/\.root/\.root:${TH2}/g'`)
343 for (( XL = $XMIN; $XL < $XMAX; XL += $DX )); do
347 LABEL=`printf '10^{%3.1f} GeV < E_{#mu} < 10^{%3.1f} GeV' $XL $XR`
350 -f "${HISTOGRAMS[*]}" \
354 -z "$ZMIN $ZMAX" -Z \
359 -o $WORKDIR/${TH2}_${XL}.gif \
366 $WORKDIR/${TH2}_*.gif \
369 rm -f $WORKDIR/${TH2}_*.gif
371 if [[ "$BACKGROUND" != "-B" ]]; then
372 display $WORKDIR/${TH2}.gif
378 set_variable TH2 muRxU
379 set_variable XMIN 0.0
380 set_variable XMAX 240.0
382 set_variable ZMIN 1.0e-1
383 set_variable ZMAX 2.0e+1
385 HISTOGRAMS=(`echo ${INPUT_FILES} | eval sed 's/\.root/\.root:${TH2}/g'`)
389 for (( XL = $XMIN; $XL < $XMAX; XL += $DX )); do
393 LABEL=`printf '%3.0f m < R < %3.0f m' $XL $XR`
395 CD=`printf '%.0f' ${XL}`
398 -f "${HISTOGRAMS[*]}" \
401 -z "$ZMIN $ZMAX" -Z \
402 -> "cos(#theta_{1})" \
406 -o $WORKDIR/${TH2}_${CD}.gif \
413 $WORKDIR/${TH2}_*.gif \
416 rm -f $WORKDIR/${TH2}_*.gif
418 if [[ "$BACKGROUND" != "-B" ]]; then
419 display $WORKDIR/${TH2}.gif
425 set_variable TH2 muRxT
426 set_variable XMIN 0.0
427 set_variable XMAX 240.0
429 set_variable ZMIN 1.0e-3
430 set_variable ZMAX 2.0e+1
432 HISTOGRAMS=(`echo ${INPUT_FILES} | eval sed 's/\.root/\.root:${TH2}/g'`)
436 for (( XL = $XMIN; $XL < $XMAX; XL += $DX )); do
440 LABEL=`printf '%3.0f m < R < %3.0f m' $XL $XR`
442 CD=`printf '%.0f' ${XL}`
445 -f "${HISTOGRAMS[*]}" \
449 -z "$ZMIN $ZMAX" -Z \
454 -o $WORKDIR/${TH2}_${CD}.gif \
461 $WORKDIR/${TH2}_*.gif \
464 rm -f $WORKDIR/${TH2}_*.gif
466 if [[ "$BACKGROUND" != "-B" ]]; then
467 display $WORKDIR/${TH2}.gif