Jpp  19.1.0-rc.1
the software that should make you happy
JNarrabri.sh
Go to the documentation of this file.
1 #!/bin/zsh
2 script=${0##*/}
3 
4 if [ -z $JPP_DIR ]; then
5  echo "Variable JPP_DIR undefined."
6  exit
7 fi
8 
9 source $JPP_DIR/setenv.sh $JPP_DIR >& /dev/null
10 
11 set_variable: DEBUG DYNAMICS_DEBUG 2
12 set_variable: WORKDIR DYNAMICS_WORKDIR ./
13 set_variable: TMAX_S DYNAMICS_TMAX_S 1000.0
14 set_variable: YMIN DYNAMICS_YMIN 0.00
15 set_variable: YMAX DYNAMICS_YMAX 0.10
16 set_variable DIR $JPP_DIR/examples/JDynamics/
17 set_variable: TIMESTAMP GRAPHICS_TIMESTAMP utc
18 set_variable: FORMAT GRAPHICS_FORMAT gif
19 set_variable+ BATCH GRAPHICS_BATCH -B
20 
21 if do_usage $*; then
22  usage "$script <detector file> (input file)+"\
23  "\nInput files correspond to the output of JKatoomba[.sh]."
24 fi
25 
26 if (( $# < 2 )); then
27  fatal "Wrong number of arguments."
28 fi
29 
30 set_variable DETECTOR $argv[1]
31 set_variable INPUT_FILE $argv[2,-1]
32 
33 eval `JPrintDetector -a $DETECTOR -O SUMMARY`
34 
35 getMechanics.sh $DETECTOR
36 
37 $DIR/JNarrabri \
38  -a $DETECTOR \
39  -f "$INPUT_FILE[*]" \
40  -o $WORKDIR/narrabri.root \
41  -T $TMAX_S \
42  -d $DEBUG --!
43 
44 
45 typeset -Z 4 STRING
46 
47 for STRING in $STRINGS[*]; do
48 
49  JPrintResult \
50  -f "$WORKDIR/narrabri.root:^\[${STRING}\]\.tilt" \
51  -F "GetStdDev(1)" \
52  -F "GetStdDev(2)" \
53  -F "GetEntries - Integral" \
54  -F "GetEntries" | read -A RESULT
55 
56  printf "%04d %4.2f %4.2f %5.0f / %8.0f\n" $STRING $RESULT[*]
57 
58 done
59 
60 echo "Creating graphics output in directory $WORKDIR"
61 
62 JPlot1D \
63  -f $WORKDIR/narrabri.root:time \
64  -> "#Deltat [s]" \
65  -\^ "number of events [a.u.]" \
66  -XX \
67  -L TL \
68  -T "" \
69  -o $WORKDIR/time.$FORMAT $BATCH
70 
71 JPlot1D \
72  -w 1200x600 \
73  -f $WORKDIR/narrabri.root:G\.\*orientation \
74  -\^ "orientation [rad]" \
75  -y "-3.20 +3.20" \
76  -N "X 505" \
77  -t "$TIMESTAMP" \
78  -L BR \
79  -T "" \
80  -o $WORKDIR/orientation.$FORMAT $BATCH
81 
82 JPlot1D \
83  -w 1200x600 \
84  -f $WORKDIR/narrabri.root:G\.\*amplitude \
85  -\^ "amplitude [rad]" \
86  -N "X 505" \
87  -t "$TIMESTAMP" \
88  -L TR \
89  -T "" \
90  -o $WORKDIR/amplitude.$FORMAT $BATCH
91 
92 JPlot1D \
93  -w 1200x600 \
94  -f $WORKDIR/narrabri.root:H\.\*orientation \
95  -\^ "orientation [rad]" \
96  -y "-3.20 +3.20" \
97  -N "X 505" \
98  -t "$TIMESTAMP" \
99  -L BR \
100  -T "" \
101  -o $WORKDIR/orientation-L.$FORMAT $BATCH
102 
103 JPlot1D \
104  -w 1200x600 \
105  -f $WORKDIR/narrabri.root:H\.\*amplitude \
106  -\^ "amplitude [rad]" \
107  -N "X 505" \
108  -t "$TIMESTAMP" \
109  -L TR \
110  -T "" \
111  -o $WORKDIR/amplitude-L.$FORMAT $BATCH
112 
113 for STRING in $STRINGS[*]; do
114 
115  JPlot1D \
116  -w 1200x600 \
117  -f "$WORKDIR/narrabri.root:G\[${STRING}\]\.orientation" \
118  -f "$WORKDIR/narrabri.root:H\[${STRING}\]\.orientation" \
119  -\^ "orientation [rad]" \
120  -y "-3.20 +3.20" \
121  -N "X 505" \
122  -t "$TIMESTAMP" \
123  -T "[${STRING}]" \
124  -o $WORKDIR/orientation_${STRING}.$FORMAT $BATCH
125 
126  JPlot1D \
127  -w 1200x600 \
128  -f "$WORKDIR/narrabri.root:G\[${STRING}\]\.amplitude" \
129  -f "$WORKDIR/narrabri.root:H\[${STRING}\]\.amplitude" \
130  -\^ "amplitude [rad]" \
131  -y "$YMIN $YMAX" \
132  -N "X 505" \
133  -t "$TIMESTAMP" \
134  -T "[${STRING}]" \
135  -o $WORKDIR/amplitude_${STRING}.$FORMAT $BATCH
136 
137  JPlot2D \
138  -f "$WORKDIR/narrabri.root:^\[${STRING}\]\.tilt" \
139  -> "dTx/dt [mrad/10 min]" \
140  -< "dTy/dt [mrad/10 min]" \
141  -O COLZ \
142  -N "X 505" \
143  -N "Y 505" \
144  -s 101010 \
145  -T "[${STRING}]" \
146  -o $WORKDIR/tilt_${STRING}.$FORMAT $BATCH
147 done