9 # ------------------------------------------------------------------------------------------
11 # Utility script to test JCalibrateK40/JFitK40.
13 # ------------------------------------------------------------------------------------------
15 if [ -z $JPP_DIR ]; then
16 echo "Variable JPP_DIR undefined."
20 source $JPP_DIR/setenv.sh $JPP_DIR
22 set_variable DIR $JPP_DIR/examples/JCalibrate
24 set_variable WORKDIR /tmp/
26 if ( do_usage $* ); then
27 usage "$script [working directory]"
31 1) set_variable WORKDIR $1;;
34 set_variable DETECTOR $WORKDIR/detector.detx
35 set_variable BACKGROUND_HZ 5e3 500
37 set_variable numberOfSlices 6000
38 set_variable DAQ_TIMESLICE JDAQTimesliceL1
40 set_variable detectorType 1
41 set_variable numberOfLines 1
42 set_variable numberOfFloors 1
46 -@ "detectorType = $detectorType" \
47 -@ "numberOfLines = $numberOfLines" \
48 -@ "numberOfFloors = $numberOfFloors" \
52 set_variable MODULE `getModule -a $DETECTOR -L "1 1"`
54 if ( ! reuse_file $WORKDIR/timeslice.root ); then
56 JRandomTimesliceWriter \
58 -o $WORKDIR/timeslice.root \
64 if ( ! reuse_file $WORKDIR/trigger_processor.root ); then
68 -f $WORKDIR/timeslice.root \
69 -o $WORKDIR/trigger_processor.root \
79 -f $WORKDIR/trigger_processor.root \
82 -o $WORKDIR/monitor.root \
86 -f $WORKDIR/monitor.root \
87 -o $WORKDIR/merge.root \
92 -f $WORKDIR/merge.root \
93 -o $WORKDIR/fit.root \
100 for (( PMT = 0; PMT < 31; PMT += 1 )); do
102 if ( ! reuse_file $WORKDIR/trigger_processor-\[${PMT}\].root ); then
105 -f $WORKDIR/trigger_processor.root \
106 -o $WORKDIR/trigger_processor-\[${PMT}\].root \
114 -f $WORKDIR/trigger_processor-\[${PMT}\].root \
117 -o $WORKDIR/monitor-\[${PMT}\].root \
121 -f $WORKDIR/monitor-\[${PMT}\].root \
122 -o $WORKDIR/merge-\[${PMT}\].root \
127 -f $WORKDIR/merge-\[${PMT}\].root \
128 -o $WORKDIR/fit-\[${PMT}\].root \
135 typeset -a INPUT_FILES
137 HISTOGRAMS=`ls -U $WORKDIR/fit.root $WORKDIR/fit-\[*\].root | eval sed 's/.root/.root:${MODULE}.1QE/g'`
140 -f "$HISTOGRAMS[*]" \
143 -o $WORKDIR/variance.root
146 -f $WORKDIR/variance.root:h0 \
147 -o $WORKDIR/scale.root \
151 -f $WORKDIR/scale.root:h0 \
154 -\^ "number of PMTs" \
156 -T "" -o QE-summary.gif
162 -f $WORKDIR/fit.root:${MODULE}.1QE \
163 -f $WORKDIR/fit-\[0\].root:${MODULE}.1QE \
164 -f $WORKDIR/fit-\[10\].root:${MODULE}.1QE \
165 -f $WORKDIR/fit-\[20\].root:${MODULE}.1QE \
166 -f $WORKDIR/fit-\[30\].root:${MODULE}.1QE \