Jpp  master_rocky-43-ge265d140c
the software that should make you happy
multi-plotk40.sh
Go to the documentation of this file.
1 #!/bin/zsh
2 #
3 # \author mdejong
4 #
5 script=${0##*/}
6 
7 # ------------------------------------------------------------------------------------------
8 #
9 # Example script to plot JFitK40 results.
10 #
11 # ------------------------------------------------------------------------------------------
12 
13 if [ -z $JPP_DIR ]; then
14  echo "Variable JPP_DIR undefined."
15  exit
16 fi
17 
18 source $JPP_DIR/setenv.sh $JPP_DIR >& /dev/null
19 
20 set_variable: FORMAT GRAPHICS_FORMAT gif
21 set_variable+ BATCH GRAPHICS_BATCH -B
22 
23 if do_usage $*; then
24  usage "$script <detector file> <input file>"
25 fi
26 
27 case $# in
28  2) set_variable DETECTOR $1
29  set_variable INPUT_FILE $2;;
30  *) fatal "Invalid number of arguments"
31 esac
32 
33 set_variable TYPE F # histogram type
34 
35 eval `JPrintDetector -a $DETECTOR -O SUMMARY`
36 
37 kill_child_processes_at_exit
38 
39 attach getModule -a $DETECTOR
40 
41 typeset -Z 4 STRING
42 typeset -Z 2 FLOOR
43 
44 for STRING in $STRINGS[*]; do
45 
46  set_variable FIRST_FLOOR 1
47 
48  for (( FLOOR = $FIRST_FLOOR; $FLOOR <= $LAST_FLOOR; FLOOR += 1 )); do
49 
50  print -u $FD_O "$STRING $FLOOR"
51 
52  read MODULE <& $FD_I
53 
54  echo "($STRING,$FLOOR) $MODULE"
55 
56  JPlot2D \
57  -f ${INPUT_FILE}:${MODULE}.2${TYPE} \
58  -z "-0.1 1.5" \
59  -> "PMT pair" \
60  -< "#Deltat [ns]" \
61  -\^ "Rate [Hz]" \
62  -T "($STRING,$FLOOR)" \
63  -O COLZ \
64  -o ${TMPDIR:-/tmp}/${TYPE}_${STRING}_${FLOOR}.$FORMAT $BATCH
65 
66  done
67 
68  montage \
69  -tile 6x3 \
70  -geometry +0+0 \
71  ${TMPDIR:-/tmp}/${TYPE}_${STRING}_*.$FORMAT \
72  ${TYPE}_${STRING}.$FORMAT
73 
74 done
75 
76 detach