Jpp 19.3.0-rc.2
the software that should make you happy
Loading...
Searching...
No Matches
plot-katoomba.sh
Go to the documentation of this file.
1#!/usr/bin/env zsh
2script=${0##*/}
3
4if [ -z $JPP_DIR ]; then
5 echo "Variable JPP_DIR undefined."
6 exit
7fi
8
9source $JPP_DIR/setenv.sh $JPP_DIR >& /dev/null
10
11set_variable: DEBUG ACOUSTICS_DEBUG 2
12set_variable: WORKDIR ACOUSTICS_WORKDIR ./
13set_variable: FORMAT GRAPHICS_FORMAT gif
14set_variable+ BATCH GRAPHICS_BATCH -B
15
16if do_usage $*; then
17 usage "$script <detector file> (input file)+"
18fi
19
20if (( $# < 2 )); then
21 fatal "Wrong number of arguments."
22fi
23
24set_variable DETECTOR $argv[1]
25set_array INPUT_FILES $argv[2,-1]
26
27eval `JPrintDetector -a $DETECTOR -O IDENTIFIER`
28eval `JPrintDetector -a $DETECTOR -O SUMMARY`
29
30JAcoustics.sh $DETECTOR_ID
31
32source JAcousticsToolkit.sh
33
34CHECK_EXIT_CODE
35
36
37JPlot1D \
38 -f${^INPUT_FILES}:chi2 \
39 -x "0 10.0" \
40 -> "#chi^{2}/NDF" \
41 -\^ "number of events [a.u.]" \
42 -T "" \
43 -o chi2.$FORMAT $BATCH
44
45JPlot1D \
46 -f${^INPUT_FILES}:hn \
47 -XX \
48 -> "number of iterations" \
49 -\^ "number of events [a.u.]" \
50 -T "" \
51 -o step.$FORMAT $BATCH
52
53
54typeset -A EMITTERS
55
56get_tripods+ $WORKDIR/tripod.txt EMITTERS
57get_transmitters+ $WORKDIR/transmitter.txt EMITTERS
58
59for EMITTER in ${(@k)EMITTERS}; do
60
61 for (( i = 1; $i <= ${#INPUT_FILES}; ++i )); do
62
63 JOpera2D \
64 -f "${INPUT_FILES[$i]}:ha\[${EMITTER}\]" \
65 -f "${INPUT_FILES[$i]}:hb\[${EMITTER}\]" \
66 -u Subtract \
67 -o $WORKDIR/subtract.root >& /dev/null
68
69 JOpera2D \
70 -f "subtract.root:Subtract" \
71 -f "${INPUT_FILES[$i]}:ha\[${EMITTER}\]" \
72 -u Divide \
73 -o $WORKDIR/ratio.root >& /dev/null
74
75 JPlot2D \
76 -f $WORKDIR/ratio.root:\.\* \
77 -> "string" \
78 -< "floor" \
79 -\^ "inefficiency" \
80 -z "1.0e-4 1.05" -Z \
81 -O COLZ \
82 -T "$EMITTER" \
83 -o ratio\[${EMITTER}\]_${i}.$FORMAT $BATCH
84
85 rm -f $WORKDIR/subtract.root
86 rm -f $WORKDIR/ratio.root
87
88 done
89done