Jpp  master_rocky-43-ge265d140c
the software that should make you happy
JDAQTimesliceSelector.sh
Go to the documentation of this file.
1 #!/bin/zsh
2 #
3 # \author mdejong
4 #
5 version=1.0
6 script=${0##*/}
7 
8 # ------------------------------------------------------------------------------------------
9 #
10 # Example script for JDAQTimeslice.
11 #
12 # ------------------------------------------------------------------------------------------
13 
14 
15 if [ -z $JPP_DIR ]; then
16  echo "Variable JPP_DIR undefined."
17  exit
18 fi
19 
20 
21 source $JPP_DIR/setenv.sh $JPP_DIR
22 
23 
24 set_variable DEBUG 2
25 set_variable WORKDIR ${TMPDIR:-/tmp}/
26 set_variable OUTPUT_FILE histogram\[%\].root
27 
28 if do_usage $*; then
29  usage "$script <input file> [output file]"
30 fi
31 
32 case $# in
33  2) set_variable OUTPUT_FILE $2;&
34  1) set_variable INPUT_FILE $1;;
35  *) fatal "Wrong number of arguments."
36 esac
37 
38 check_input_file $INPUT_FILE
39 
40 for TIMESLICE in JDAQTimeslice JDAQTimesliceL0 JDAQTimesliceL1 JDAQTimesliceL2 JDAQTimesliceSN; do
41 
42  echo Processing $INPUT_FILE with selection $TIMESLICE.
43 
44  $JPP_DIR/examples/JDAQ/JDAQTimesliceSelector \
45  -f ${INPUT_FILE} \
46  -o ${OUTPUT_FILE/\%/${TIMESLICE}} \
47  -C $TIMESLICE \
48  -d $DEBUG
49 done
50 
51 JPlot1D \
52  -f ${OUTPUT_FILE/\%/JDAQTimeslice}:h0 \
53  -f ${OUTPUT_FILE/\%/JDAQTimesliceL0}:h0 \
54  -f ${OUTPUT_FILE/\%/JDAQTimesliceL1}:h0 \
55  -f ${OUTPUT_FILE/\%/JDAQTimesliceL2}:h0 \
56  -f ${OUTPUT_FILE/\%/JDAQTimesliceSN}:h0 \
57  -> "PMT" \
58  -\^ "number of events [a.u.]" \
59  -T "" -L TR
60 
61 JPlot1D \
62  -f ${OUTPUT_FILE/\%/JDAQTimeslice}:h1 \
63  -f ${OUTPUT_FILE/\%/JDAQTimesliceL0}:h1 \
64  -f ${OUTPUT_FILE/\%/JDAQTimesliceL1}:h1 \
65  -f ${OUTPUT_FILE/\%/JDAQTimesliceL2}:h1 \
66  -f ${OUTPUT_FILE/\%/JDAQTimesliceSN}:h1 \
67  -> "time-over-threshold [ns]" \
68  -\^ "number of events [a.u.]" \
69  -T "" -L TR
70 
71 JPlot1D \
72  -f ${OUTPUT_FILE/\%/JDAQTimeslice}:h2 \
73  -f ${OUTPUT_FILE/\%/JDAQTimesliceL0}:h2 \
74  -f ${OUTPUT_FILE/\%/JDAQTimesliceL1}:h2 \
75  -f ${OUTPUT_FILE/\%/JDAQTimesliceL2}:h2 \
76  -f ${OUTPUT_FILE/\%/JDAQTimesliceSN}:h2 \
77  -> "time [ns]" \
78  -\^ "number of events [a.u.]" \
79  -w 800x400 \
80  -T "" -L TR