Jpp
plot-PMTParameters2D.sh
Go to the documentation of this file.
1 #!/bin/zsh
2 #
3 #
4 # \author mdejong
5 #
6 version=1.0
7 script=${0##*/}
8 
9 if [ -z $JPP_DIR ]; then
10  echo "Variable JPP_DIR undefined."
11  exit
12 fi
13 
14 source $JPP_DIR/setenv.sh $JPP_DIR
15 
16 
17 set_variable DEBUG 3
18 set_variable WORKDIR `pwd`
19 set_variable OUTPUT_FILE $WORKDIR/pmt_parameters.root
20 set_variable DIR $JPP_DIR/examples/JDetector/
21 set_variable NUMBER_OF_PMTS 31
22 
23 if ( do_usage $* ); then
24  usage "$script <detector file> (PMT parameters file)+"
25 fi
26 
27 if (( $# < 2 )); then
28  fatal "Wrong number of parameters."
29 fi
30 
31 set_variable DETECTOR $argv[1]
32 
33 INPUT_FILES=($argv[2,-1])
34 
35 $DIR/JPlotPMTParameters2D \
36  -a $DETECTOR \
37  -P "$INPUT_FILES" \
38  -o $OUTPUT_FILE \
39  -d $DEBUG --!
40 
41 # parameters and plot ranges
42 
43 typeset -A Y_RANGE
44 
45 Y_RANGE+=(EFFICIENCY " 0.0 1.5")
46 Y_RANGE+=(QE " 0.0 1.5")
47 Y_RANGE+=(gain " 0.0 1.5")
48 Y_RANGE+=(gainSpread " 0.0 1.0")
49 Y_RANGE+=(riseTime_ns " 0.0 15.0")
50 Y_RANGE+=(TTS_ns "-5.0 +5.0")
51 Y_RANGE+=(threshold " 0.0 1.0")
52 Y_RANGE+=(slope " 0.0 15.0")
53 Y_RANGE+=(saturation " 0.0 250.0")
54 Y_RANGE+=(slewing " 0.0 1.5")
55 
56 
57 eval `JPrintDetector -a $DETECTOR -O SUMMARY`
58 
59 attach getModule -a $DETECTOR
60 
61 typeset -Z 3 STRING
62 typeset -Z 2 FLOOR
63 
64 let XMIN=-0.5
65 let XMAX=${NUMBER_OF_PMTS}-0.5
66 
67 for STRING in $STRINGS[*]; do
68 
69  for (( FLOOR = $FIRST_FLOOR; $FLOOR <= $LAST_FLOOR; FLOOR += 1 )); do
70 
71  print -u $FD_O "$STRING $FLOOR"
72 
73  read MODULE <& $FD_I
74 
75  for KEY in ${(k)Y_RANGE}; do
76 
77  debug "($STRING,$FLOOR) $KEY"
78 
79  JPlot2D \
80  -f ${OUTPUT_FILE}:${MODULE}\.${KEY}\$ \
81  -z $Y_RANGE[${KEY}] \
82  -O COLZ \
83  -> "index" \
84  -< "PMT" \
85  -\^ "$KEY" \
86  -T "(${STRING},${FLOOR})" \
87  -o /tmp/${KEY}_${STRING}_${FLOOR}.gif -B
88 
89  JProject2D \
90  -f ${OUTPUT_FILE}:${MODULE}\.${KEY}\$ \
91  -P X -y "$XMIN $XMAX" \
92  -o /tmp/project.root
93 
94  JScale1D \
95  -f /tmp/project.root:\.\* \
96  -o /tmp/scale.root \
97  -F "1.0/${NUMBER_OF_PMTS}"
98 
99  JPlot1D \
100  -f /tmp/scale.root:\.\* \
101  -y $Y_RANGE[${KEY}] \
102  -> "index" \
103  -\^ "<$KEY>" \
104  -T "(${STRING},${FLOOR})" \
105  -o /tmp/P_${KEY}_${STRING}_${FLOOR}.gif -B
106 
107 
108  done
109  done
110 
111  for KEY in ${(k)Y_RANGE}; do
112 
113  montage \
114  -tile 6x3 \
115  -geometry +0+0 \
116  /tmp/${KEY}_${STRING}_*.gif \
117  $WORKDIR/${KEY}_${STRING}.gif
118  done
119 
120  for KEY in ${(k)Y_RANGE}; do
121 
122  montage \
123  -tile 6x3 \
124  -geometry +0+0 \
125  /tmp/P_${KEY}_${STRING}_*.gif \
126  $WORKDIR/P_${KEY}_${STRING}.gif
127  done
128 done
129 
130 detach