Jpp
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
plot-canberra.sh
Go to the documentation of this file.
1 #!/bin/zsh
2 script=${0##*/}
3 
4 zmodload zsh/mathfunc
5 
6 if [ -z $JPP_DIR ]; then
7  echo "Variable JPP_DIR undefined."
8  exit
9 fi
10 
11 source $JPP_DIR/setenv.sh $JPP_DIR >& /dev/null
12 
14 set_variable WORKDIR ${ACOUSTICS_WORKDIR:-./}
15 set_variable BACKGROUND ${ACOUSTICS_BACKGROUND:+-B}
16 set_variable FORMULA "[0] * exp(-0.5 * (x-[1])*(x-[1]) / ([2]*[2]))"
17 set_variable NUMBER_OF_ENTRIES 100
18 set_variable TMIN_MS -1.0
19 set_variable TMAX_MS +1.0
20 
21 if ( do_usage $* ); then
22  usage "$script <detector file> <input file>"
23 fi
24 
25 if (( $# != 2 )); then
26  fatal "Wrong number of arguments."
27 fi
28 
31 
32 eval `JPrintDetector -a $DETECTOR -O IDENTIFIER`
33 eval `JPrintDetector -a $DETECTOR -O SUMMARY`
34 
35 source JAcoustics.sh $DETECTOR_ID
36 
37 
38 typeset -A TRIPODS
39 
40 get_tripods $WORKDIR/tripod.txt TRIPODS
41 
42 for EMITTER in ${(@k)TRIPODS}; do
43 
44  eval JCanberra \
45  -a $DETECTOR \
47  -o \"$WORKDIR/canberra\[${EMITTER}\].root\" \
48  -T `print_tripods` \
49  -V `print_velocity` \
50  -E $EMITTER \
51  -d $DEBUG --!
52 done
53 
54 
55 HISTOGRAMS=(`ls $WORKDIR/canberra\[*\].root`)
56 
57 JOpera1D \
58  -f "${HISTOGRAMS//\.root/.root:.*toa}" \
59  -u Add \
60  -o $WORKDIR/sum.root >& /dev/null
61 
62 set_variable X `JPrintResult -f $WORKDIR/sum.root:\.\* -F GetMean`
63 set_variable RMS `JPrintResult -f $WORKDIR/sum.root:\.\* -F GetRMS`
64 
65 let XMIN="$X - 5.0*$RMS"
66 let XMAX="$X + 5.0*$RMS"
67 
68 JFit \
69  -f $WORKDIR/sum.root:\.\* \
70  -F "$FORMULA" \
71  -@ "p0 = GetMaximum" \
72  -@ "p1 = $X" \
73  -@ "p2 = $RMS" \
74  -x "$XMIN $XMAX" \
75  -o $WORKDIR/fit.root \
76  -d $DEBUG
77 
78 let XMIN="$TMIN_MS * 1.0e-3"
79 let XMAX="$TMAX_MS * 1.0e-3"
80 
81 JPlot1D \
82  -f $WORKDIR/fit.root:\.\* \
83  -\> "#Deltat [s]" \
84  -x "$XMIN $XMAX" \
85  -N "X 505" \
86  -T "" \
87  -s 0 \
88  -o $WORKDIR/canberra.gif $BACKGROUND
89 
90 rm -f $WORKDIR/sum.root
91 rm -f $WORKDIR/fit.root
92 
93 attach getModule -a $DETECTOR
94 
95 typeset -Z 4 STRING
96 typeset -Z 2 FLOOR
97 
98 set_variable FIRST_FLOOR 1
99 
100 for STRING in $STRINGS[*]; do
101 
102  echo -n "Creating graphics for string $STRING.."
103 
104  HISTOGRAMS=()
105 
106  for EMITTER in ${(@k)TRIPODS}; do
107 
108  set_variable STRING_TXT $WORKDIR/string_${STRING}\[${EMITTER}\].txt
109  set_variable STRING_ROOT $WORKDIR/string_${STRING}\[${EMITTER}\].root
110 
111  rm -f $STRING_TXT
112 
113  for (( FLOOR = $FIRST_FLOOR; $FLOOR <= $LAST_FLOOR; FLOOR += 1 )); do
114 
115  sput $STRING $FLOOR
116  sget MODULE
117 
118  set_variable N `JPrintResult -f $WORKDIR/canberra\[$EMITTER\].root:${MODULE} -F GetEntries`
119 
120  if [[ -n "$N" ]] && (( $N >= $NUMBER_OF_ENTRIES )); then
121 
122  set_variable X `JPrintResult -f $WORKDIR/canberra\[$EMITTER\].root:${MODULE} -F GetMean`
123  set_variable RMS `JPrintResult -f $WORKDIR/canberra\[$EMITTER\].root:${MODULE} -F GetRMS`
124 
125  let XMIN="$X - 3.0*$RMS"
126  let XMAX="$X + 3.0*$RMS"
127 
128  JFit \
129  -f $WORKDIR/canberra\[$EMITTER\].root:${MODULE} \
130  -F "$FORMULA" \
131  -@ "p0 = GetMaximum" \
132  -@ "p1 = $X" \
133  -@ "p2 = $RMS" \
134  -x "$XMIN $XMAX" \
135  -O L \
136  -o $WORKDIR/fit.root \
137  -d $DEBUG
138 
139  if ( `JPrintFit -f $WORKDIR/fit.root:\.\* -@p0 >& /dev/null` ); then
140 
141  set_array P1 `JPrintFit -f $WORKDIR/fit.root:\.\* -@p1`
142  set_array P2 `JPrintFit -f $WORKDIR/fit.root:\.\* -@p2`
143 
144  ((P1[1] = $P1[1] * 1.0e3)) # [ms]
145  ((P2[1] = abs($P2[1]) * 1.0e3)) # [ms]
146 
147  printf "%2d %9.6f %9.6f\n" $FLOOR $P1[1] $P2[1] >> $STRING_TXT
148  fi
149 
150  rm -f $WORKDIR/fit.root
151  fi
152  done
153 
154  JGraph \
155  -f $STRING_TXT \
156  -o $STRING_ROOT
157 
158  HISTOGRAMS+="${STRING_ROOT}:.*"
159 
160  done
161 
162  JPlot1D \
163  -f "$HISTOGRAMS[*]" \
164  -w 800x400 \
165  -x "0.5 18.5" \
166  -y "$TMIN_MS $TMAX_MS" \
167  -\> "floor" \
168  -\^ "#Deltat [ms]" \
169  -N "X 505" \
170  -T "[$STRING]" \
171  -L TR \
172  -o $WORKDIR/string_${STRING}.gif $BACKGROUND
173 
174  echo " -> $WORKDIR/string_${STRING}.gif"
175 
176 done
177 
178 detach
function get_tripods()
Definition: JAcoustics.sh:100
data_type w[N+1][M+1]
Definition: JPolint.hh:741
TPaveText * p1
void set_variable(const std::string &name, const std::string &value)
Set environment variable.
do eval JCanberra a $DETECTOR f $INPUT_FILE o $WORKDIR canberra[${EMITTER}\] root T print_tripods V print_velocity E $EMITTER d $DEBUG!done HISTOGRAMS
clean eval JPrintDetector a $DETECTOR O IDENTIFIER eval JPrintDetector a $DETECTOR O SUMMARY set_variable STRING
set_variable INPUT_FILE
then fatal No sound hydrophone file $HYDROPHONE_TXT fi JGraph f $HYDROPHONE_TXT o $HYDROPHONE_ROOT sort gr k
do set_variable STRING_TXT $WORKDIR string_
exit
Definition: JPizza.sh:36
do set_array DAQHEADER JPrintDAQHeader f
Definition: JTuneHV.sh:79
then echo
then fatal Invalid string $STRING
do set_variable OUTPUT_DIRECTORY $WORKDIR T
do $DIR JTransitTime o
Definition: JTransitTime.sh:44
then awk F
* usage
then JFileTuna f $INPUT_FILE
Definition: file-Tuna.sh:66
then fatal Wrong number of arguments fi set_variable STRING $argv[1] set_variable DETECTOR_TXT $WORKDIR $DETECTOR_TXT tail read X Y Z RMS let XMIN
then echo Variable JPP_DIR undefined exit fi source $JPP_DIR setenv sh $JPP_DIR &dev null set_variable ACOUSTICS_DEBUG
Definition: JAcoustics.sh:17
p2
Definition: module-Z:fit.sh:48
*set_variable DETECTOR $argv[1] set_variable INPUT_FILE $argv[2,-2] set_variable OUTPUT_FILE $argv[-1]
Definition: JDomino.sh:39
then JMuonMCEvt f $INPUT_FILE o $INTERMEDIATE_FILE d
Definition: JMuonPath.sh:45
then $JPP_DIR software JCalibrate JCalibrateToT a
Definition: JTuneHV.sh:108
alias put_queue eval echo n
Definition: qlib.csh:19
then display $WORKDIR
Definition: plot-Domino.sh:127
set_variable DETECTOR
then fatal Wrong number of arguments fi set_variable DETECTOR $argv[1] set_variable INPUT_FILE $argv[2] eval JPrintDetector a $DETECTOR O IDENTIFIER eval JPrintDetector a $DETECTOR O SUMMARY source JAcoustics sh $DETECTOR_ID typeset A TRIPODS get_tripods $WORKDIR tripod txt TRIPODS for EMITTER in
Definition: JCanberra.sh:36
then usage $script[input file[working directory[option]]] nWhere option can be N
Definition: JMuonPostfit.sh:37
script
Definition: JAcoustics.sh:2
then usage $script< detector file >< inputfile > fi case set_variable WORKDIR
Definition: JLegolas.sh:28
source $JPP_DIR setenv csh $JPP_DIR eval JShellParser o a A
esac $JPP_BIN JLogger sh $LOGGER until pgrep JGetMessage</dev/null > dev null
const JModule & getModule(const JDetector &detector, const JModuleLocation &location)
find module with a given string and floor number
do eval JCanberra a $DETECTOR f $INPUT_FILE o $WORKDIR canberra[${EMITTER}\] root T print_tripods V print_velocity E $EMITTER d $DEBUG!done attach getModule a $DETECTOR typeset Z STRING typeset Z FLOOR set_variable FIRST_FLOOR for STRING in $STRINGS[*]
Definition: JCanberra.sh:55
esac done
Definition: JAddHDE.sh:21
#define DEBUG(A)
Message macros.
Definition: JMessage.hh:62