Jpp  18.6.0-rc.1
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 
*fatal Wrong number of arguments esac source JAcousticsToolkit sh mkdir p $WORKDIR post calibration pushd $WORKDIR post calibration post calibration_D0ARCA021 sh $DIR popd set_variable DETECTOR $WORKDIR post calibration detector datx set_variable TRIPOD $WORKDIR post calibration tripod txt set_variable DETECTOR_INITIAL detector_initial datx expand_array RUNS mkdir p $HOMEDIR cd $HOMEDIR JDetectorDB D $DETECTOR_ID r $RUNS[1] V WW o $DETECTOR_INITIAL JEditDetector a $DETECTOR r o detector_tmp datx!eval JPrintDetector a $DETECTOR_INITIAL O SUMMARY
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 $script< input file >[option[primary[working directory]]] nWhere option can be E
Definition: JMuonPostfit.sh:40
clean eval JPrintDetector a $DETECTOR O IDENTIFIER eval JPrintDetector a $DETECTOR O SUMMARY set_variable STRING
then usage eval JPrintDetector a $DETECTOR O IDENTIFIER eval JPrintDetector a $DETECTOR O SUMMARY JAcoustics sh $DETECTOR_ID JCanberra sh $DETECTOR $INPUT_FILES[*] $WORKDIR canberra[%\] root source JAcousticsToolkit sh CHECK_EXIT_CODE typeset A EMITTERS get_tripods $WORKDIR tripod txt EMITTERS get_transmitters $WORKDIR transmitter txt EMITTERS kill_child_processes_at_exit attach getModule a $DETECTOR typeset Z STRING typeset Z FLOOR for STRING in $STRINGS[*]
version
Definition: JEditTuneHV.sh:5
o $QUALITY_ROOT d $DEBUG!CHECK_EXIT_CODE JPlot1D f
Definition: JDataQuality.sh:76
exit
Definition: JPizza.sh:36
then
Definition: datalogs.sh:34
then fatal Wrong number of parameters fi set_variable DETECTOR $argv[1] set_variable INPUT_FILE $argv[2] eval JPrintDetector a $DETECTOR O IDENTIFIER set_variable DETID getDetector D $DETECTOR_ID O string $DIR JPlotPMTParameters a $DETECTOR P $INPUT_FILE o $OUTPUT_FILE d $DEBUG for KEY in TTS_ns QE gain gainSpread
then echo Variable JPP_DIR undefined exit fi source $JPP_DIR setenv sh $JPP_DIR &dev null set_variable DETECTOR_ID set_array RUNS if do_usage *then usage $script< source directory > nThe source directory corresponds to a pre calibrated D0ARCA009 detector fi case set_variable DIR
then fatal Wrong number of arguments fi JConvertDetectorFormat a o
static const double C
Physics constants.
then echo Variable JPP_DIR undefined exit fi source $JPP_DIR setenv sh $JPP_DIR typeset A Y_RANGE Y_RANGE
then fatal Invalid string $STRING
then JCalibrateToT a
Definition: JTuneHV.sh:107
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
do set_variable OUTPUT_DIRECTORY $WORKDIR T
then echo Variable JPP_DIR undefined exit fi source $JPP_DIR setenv sh $JPP_DIR &dev null set_variable
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
then awk F
* usage
&set_variable OUTPUT_FILE
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: JMuonStart.sh:47
then fatal The output file must have the wildcard in the e g root fi 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:48
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
do sput $STRING $FLOOR sget MODULE INPUT_FILES
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
*fatal Wrong number of arguments esac eval JPrintDetector a $DETECTOR O IDENTIFIER eval JPrintDetector a $DETECTOR O SUMMARY set_variable FIRST_FLOOR
then eval ls $DIR $ID $TYPE _
Definition: getArchive.sh:108
set_variable DETECTOR
then echo
Definition: JQAQC.sh:92
source $JPP_DIR setenv csh $JPP_DIR &dev null eval JShellParser o a A
do echo Generating $dir eval D
Definition: JDrawLED.sh:53
then getFile sh $JPP_ARCHIVE $DETECTOR_ID $RUNS[1] $KEY $VERSION $WORKDIR
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
int debug
debug level
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