Jpp  17.0.0
the software that should make you happy
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
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 
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: FORMAT GRAPHICS_FORMAT gif
22 set_variable+ BATCH GRAPHICS_BATCH -B
23 
24 if do_usage $*; then
25  usage "$script <detector file> (PMT parameters file)+"
26 fi
27 
28 if (( $# < 2 )); then
29  fatal "Wrong number of parameters."
30 fi
31 
33 set_array INPUT_FILES $argv[2,-1]
34 
35 $DIR/JPlotPMTParameters2D \
36  -a $DETECTOR \
37  -P "$INPUT_FILES" \
38  -o $OUTPUT_FILE \
39  -d $DEBUG -A --!
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 eval `JPrintDetector -a $DETECTOR -O SUMMARY`
57 
59 
60 kill_child_processes_at_exit
61 
62 attach getModule -a $DETECTOR
63 
64 typeset -Z 4 STRING
65 typeset -Z 2 FLOOR
66 
67 typeset -A PMT_RANGE
68 
69 PMT_RANGE[-]="0-30"
70 PMT_RANGE[A]="0-0"
71 PMT_RANGE[B]="1-6"
72 PMT_RANGE[C]="7-12"
73 PMT_RANGE[D]="13-18"
74 PMT_RANGE[E]="19-24"
75 PMT_RANGE[F]="25-30"
76 
77 for STRING in $STRINGS[*]; do
78 
79  for (( FLOOR = $FIRST_FLOOR; $FLOOR <= $LAST_FLOOR; FLOOR += 1 )); do
80 
81  sput $STRING $FLOOR
82  sget MODULE
83 
84  for KEY in ${(k)Y_RANGE}; do
85 
86  debug "($STRING,$FLOOR) $KEY"
87 
88  JPlot2D \
89  -f ${OUTPUT_FILE}:${MODULE}\.${KEY}\$ \
90  -z $Y_RANGE[${KEY}] \
91  -O COLZ \
92  -\> "index" \
93  -\< "PMT" \
94  -\^ "$KEY" \
95  -T "(${STRING},${FLOOR})" \
96  -o ${TMPDIR:-/tmp}/${KEY}_${STRING}_${FLOOR}.$FORMAT $BATCH
97 
98  for RING RANGE in ${(@kv)PMT_RANGE[*]}; do
99 
100  let "YMIN = ${RANGE%%-*} - 0.5"
101  let "YMAX = ${RANGE##*-} + 0.5"
102 
103  JProject2D \
104  -f ${OUTPUT_FILE}:${MODULE}\.${KEY}\$ \
105  -P X -y "$YMIN $YMAX" \
106  -o ${TMPDIR:-/tmp}/project.root
107 
108  JScale1D \
109  -f ${TMPDIR:-/tmp}/project.root:\.\* \
110  -o ${TMPDIR:-/tmp}/scale\[$RING\].root \
111  -F "1.0 / ($YMAX - $YMIN)" \
112  -O "nosw2"
113  done
114 
115  set_array ABC `ls -f ${TMPDIR:-/tmp}/scale\[-\].root ${TMPDIR:-/tmp}/scale\[[A-Z]\].root`
116 
117  JPlot1D \
118  -f${^ABC}:\.\* \
119  -y $Y_RANGE[${KEY}] \
120  -\> "index" \
121  -\^ "<$KEY>" \
122  -L BL \
123  -T "(${STRING},${FLOOR})" \
124  -o ${TMPDIR:-/tmp}/P_${KEY}_${STRING}_${FLOOR}.$FORMAT $BATCH
125  done
126  done
127 
128  for KEY in ${(k)Y_RANGE}; do
129 
130  montage \
131  -tile 6x3 \
132  -geometry +0+0 \
133  ${TMPDIR:-/tmp}/${KEY}_${STRING}_{01..18}.$FORMAT \ $WORKDIR/${KEY}_${STRING}.$FORMAT >& /dev/null
134  done
135 
136  for KEY in ${(k)Y_RANGE}; do
137 
138  montage \
139  -tile 6x3 \
140  -geometry +0+0 \
141  ${TMPDIR:-/tmp}/P_${KEY}_${STRING}_{01..18}.$FORMAT \ $WORKDIR/P_${KEY}_${STRING}.$FORMAT >& /dev/null
142  done
143 done
144 
145 detach
146 
then usage $script< detector identifier >< startdate\"YYYY-MM-DDHH:MM:SS\"><finaldate\"YYYY-MM-DDHH:MM:SS\"><QA/QCfile> fi case set_variable QAQC_TXT $argv[4]
Definition: JDataMonitor.sh:24
then fatal No hydrophone data file $HYDROPHONE_TXT fi sort gr k
then usage E
Definition: JMuonPostfit.sh:35
clean eval JPrintDetector a $DETECTOR O IDENTIFIER eval JPrintDetector a $DETECTOR O SUMMARY set_variable STRING
version
Definition: JEditTuneHV.sh:5
o $QUALITY_ROOT d $DEBUG!JPlot1D f
Definition: JDataQuality.sh:66
exit
Definition: JPizza.sh:36
do echo Generating $dir eval DIR
Definition: JPlotNPE1D.sh:54
then echo
do JCanberra a $DETECTOR f $INPUT_FILE o $WORKDIR canberra[${EMITTER}\] root T $WORKDIR tripod txt V $WORKDIR sound_velocity txt M $WORKDIR mechanics txt H $WORKDIR hydrophone txt E $EMITTER $DISABLE d $DEBUG!done kill_child_processes_at_exit attach getModule a $DETECTOR typeset Z STRING typeset Z FLOOR for STRING in $STRINGS[*]
Definition: JCanberra.sh:68
static const double C
Physics constants.
then fatal Invalid string $STRING
then echo Variable JPP_DIR undefined exit fi source $JPP_DIR setenv sh $JPP_DIR &dev null set_variable
Definition: JAcoustics.sh:21
then echo Variable JPP_DIR undefined exit fi source $JPP_DIR setenv sh $JPP_DIR if do_usage *then usage $script[(input file)+] fi set_variable DEBUG set_variable WORKDIR TMPDIR
*fatal Wrong number of arguments esac eval JPrintDetector a $DETECTOR O IDENTIFIER eval JPrintDetector a $DETECTOR O SUMMARY set_variable FIRST_FLOOR
do set_variable OUTPUT_DIRECTORY $WORKDIR T
then $JPP_DIR examples JDetector JSlewing P slewing
Definition: JSlewing.sh:43
then echo Variable JPP_DIR undefined exit fi source $JPP_DIR setenv sh $JPP_DIR &dev null set_variable DEBUG set_variable WORKDIR
Definition: JLegolas.sh:20
int debug
debug level
Definition: JSirene.cc:66
then awk F
* usage
then JConvertDetectorFormat a $DETECTOR[1] o
&set_variable OUTPUT_FILE
do sput $STRING $FLOOR sget MODULE INPUT_FILES
Definition: JCanberra.sh:77
then fatal Wrong number of parameters fi set_variable DETECTOR $argv[1] set_array INPUT_FILE $argv[2] $DIR JPlotPMTParameters a $DETECTOR P $INPUT_FILE o $OUTPUT_FILE d $DEBUG A typeset A Y_RANGE Y_RANGE
then $JPP_DIR examples JDetector JToT o $OUTPUT_FILE n N $NPE P gain
Definition: JToT.sh:47
then JMuonMCEvt f $INPUT_FILE o $INTERMEDIATE_FILE d
Definition: JMuonPath.sh:47
then JCalibrateToT a
Definition: JTuneHV.sh:116
then set_variable MODULE getModule a $DETECTOR L $STRING $FLOOR JEditDetector a $DETECTOR M $MODULE add $X o $DETECTOR else echo No update of detector $DETECTOR
then display $WORKDIR
Definition: plot-Domino.sh:128
no fit printf nominal n $STRING awk v X
do set_variable MODULE getModule a $WORKDIR detector_a datx L $STRING JEditDetector a $WORKDIR detector_a datx M $MODULE setz o $WORKDIR detector_a datx JEditDetector a $WORKDIR detector_b datx M $MODULE setz o $WORKDIR detector_b datx done echo Output stored at $WORKDIR detector_a datx and $WORKDIR tripod_a txt JDrawDetector2D a $WORKDIR detector_a datx a $WORKDIR detector_b datx L BL o detector $FORMAT $BATCH JDrawDetector2D T $WORKDIR tripod_a txt T $WORKDIR tripod_b txt L BL o tripod $FORMAT $BATCH JCompareDetector a $WORKDIR detector_a datx b $WORKDIR detector_b datx o $WORKDIR abc root &dev null for KEY in X Y Z
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 JAcoustics sh $DETECTOR_ID source JAcousticsToolkit sh CHECK_EXIT_CODE typeset A EMITTERS get_tripods $WORKDIR tripod txt EMITTERS get_transmitters $WORKDIR transmitter txt EMITTERS for EMITTER in
Definition: JCanberra.sh:46
source $JPP_DIR setenv csh $JPP_DIR &dev null eval JShellParser o a A
do echo Generating $dir eval D
Definition: JDrawLED.sh:53
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
then fatal Invalid detector identifier $DETECTOR_ID fi set_variable RUNSETUPID typeset a RANGE RANGE[1]
then $DIR JPlotNPE PDG P
Definition: JPlotNPE-PDG.sh:62
esac done
Definition: JAddHDE.sh:21
#define DEBUG(A)
Message macros.
Definition: JMessage.hh:62
void scale(vector< double > &v, double c)
scale vector content