Jpp  16.0.0-rc.2
the software that should make you happy
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
JMakePDF.sh
Go to the documentation of this file.
1 #!/bin/zsh
2 #
3 # \author mdejong
4 #
5 # utility script to produce PDF and CDF tables.
6 script=${0##*/}
7 
8 
9 # ------------------------------------------------------------------------------------------
10 #
11 # Utility script to produce PDF and CDF tables.
12 #
13 # ------------------------------------------------------------------------------------------
14 
15 if [ -z $JPP_DIR ]; then
16  echo "Variable JPP_DIR undefined."
17  exit
18 fi
19 
20 source $JPP_DIR/setenv.sh $JPP_DIR >& /dev/null
21 
22 set_variable DEBUG 2
23 set_variable WORKDIR $JPP_DATA
24 set_variable MAKE_PDF 0
25 set_variable MERGE_PDF 0
26 set_variable MAKE_CDF 0
27 set_variable CLEAN 0
28 set_variable BACKGROUND
29 set_variable WGET
30 set_variable RMAX_M 999
31 set_variable absorptionLengthFactor 1.0
32 set_variable scatteringLengthFactor 1.0
33 
34 # check whether script is already running
35 
36 BUFFER=`get_pid $script`
37 BUFFER=${BUFFER/${$}/}
38 
39 if [[ -n "${BUFFER/ */}" ]]; then
40  notice "${script} already running <${BUFFER}>"
41  exit 0
42 fi
43 
44 #
45 # Method to print usage of script.
46 #
47 function usage()
48 {
49  print "$script \n\
50  -d <debug level> \n\
51  -W <working directory PDF and CDF tables> \n\
52  -P create PDF tables \n\
53  -M merge PDF tables \n\
54  -C convert PDF tables to CDF tables \n\
55  -c clean PDF tables \n\
56  -B parallel running of jobs in background \n\
57  -w download files from git server \n\
58  -R <maximal distance between muon and PMT> \
59  "
60  exit
61 }
62 
63 #
64 # Parse command line options.
65 #
66 while getopts "d:W:PMCcBw:hA:S:R:" option; do
67  case "$option" in
68  d) set_variable DEBUG $OPTARG;;
69  W) set_variable WORKDIR $OPTARG;;
70  P) set_variable MAKE_PDF 1;;
71  M) set_variable MERGE_PDF 1;;
72  C) set_variable MAKE_CDF 1;;
73  c) set_variable CLEAN 1;;
74  B) set_variable BACKGROUND \&;;
75  w) set_variable WGET $OPTARG;;
76  h) usage; exit;;
77  A) set_variable absorptionLengthFactor $OPTARG;;
78  S) set_variable scatteringLengthFactor $OPTARG;;
79  R) set_variable RMAX_M $OPTARG;;
80  esac
81 done
82 
83 check_variable ROOTSYS JPP_DIR WORKDIR
84 
85 mkdir -p $WORKDIR
86 
87 status "Start $script"
88 
89 
90 if [[ -n $WGET ]]; then
91  cd $WORKDIR
92  rm -f *.dat
93  wget -A dat,txt,root,evt.gz -r -l 1 -nd -q $WGET
94 fi
95 
96 
97 if (( $MAKE_PDF )); then
98 
99  check_process JMakePDF JMakePDG
100 
101  if (( 1 )); then
102 
103  if [[ ! -f $WORKDIR/J1p.dat ]]; then
104  eval nohup JMakePDF -F1 -A $absorptionLengthFactor -S $scatteringLengthFactor -o $WORKDIR/J1p.dat -d $DEBUG $BACKGROUND
105  fi
106  if [[ ! -f $WORKDIR/J2p.dat ]]; then
107  eval nohup JMakePDF -F2 -A $absorptionLengthFactor -S $scatteringLengthFactor -o $WORKDIR/J2p.dat -d $DEBUG $BACKGROUND
108  fi
109 
110  for R in 0.1 0.3 0.5 1.0 2.0 3.0 4.0 5.0 6.0 7.0 8.0 9.0 10.0 \
111  11.0 12.0 13.0 14.0 15.0 16.0 17.0 18.0 19.0 20.0 \
112  22.0 24.0 26.0 28.0 30.0 32.0 34.0 36.0 38.0 40.0 42.0 44.0 46.0 48.0 50.0 \
113  55.0 60.0 65.0 70.0 75.0 80.0 85.0 90.0 95.0 100.0 110.0 120.0 130.0 140.0 150.0 \
114  170.0 190.0 210.0 230.0 250.0 270.0 290.0 310.0; do
115 
116  if (( $R <= $RMAX_M )); then
117 
118  if [[ ! -f $WORKDIR/J3p.dat ]]; then
119  eval nohup JMakePDF -F3 -A $absorptionLengthFactor -S $scatteringLengthFactor -R $R -o $WORKDIR/J3p-$R.dat -d $DEBUG $BACKGROUND
120  fi
121  if [[ ! -f $WORKDIR/J4p.dat ]]; then
122  eval nohup JMakePDF -F4 -A $absorptionLengthFactor -S $scatteringLengthFactor -R $R -o $WORKDIR/J4p-$R.dat -d $DEBUG $BACKGROUND
123  fi
124  fi
125  done
126 
127  if [[ ! -f $WORKDIR/J5p.dat ]]; then
128  eval nohup JMakePDF -F5 -A $absorptionLengthFactor -S $scatteringLengthFactor -o $WORKDIR/J5p.dat -d $DEBUG $BACKGROUND
129  fi
130 
131  for R in 0.1 0.3 0.5 1.0 2.0 3.0 4.0 5.0 6.0 7.0 8.0 9.0 10.0 \
132  11.0 12.0 13.0 14.0 15.0 16.0 17.0 18.0 19.0 20.0 \
133  22.0 24.0 26.0 28.0 30.0 32.0 34.0 36.0 38.0 40.0 42.0 44.0 46.0 48.0 50.0 \
134  55.0 60.0 65.0 70.0 75.0 80.0 85.0 90.0 95.0 100.0 110.0 120.0 130.0 140.0 150.0 \
135  170.0 190.0 210.0 230.0 250.0; do
136 
137  if (( $R <= $RMAX_M )); then
138 
139  if [[ ! -f $WORKDIR/J6p.dat ]]; then
140  eval nohup JMakePDF -F6 -A $absorptionLengthFactor -S $scatteringLengthFactor -R $R -o $WORKDIR/J6p-$R.dat -d $DEBUG $BACKGROUND
141  fi
142  fi
143  done
144 
145  let DMAX_M="${RMAX_M}*1.5"
146 
147  for D in 0.1 0.5 1.0 5.0 10.0 20.0 30.0 40.0 50.0 60.0 70.0 80.0 90.0 100.0 \
148  120.0 150.0 170.0 190.0 210.0 230.0 250.0 270.0 290.0 310.0; do
149 
150  if (( $D <= $DMAX_M )); then
151  if [[ ! -f $WORKDIR/J12p.dat ]]; then
152  eval nohup JMakePDG -F12 -A $absorptionLengthFactor -S $scatteringLengthFactor -D $D -o $WORKDIR/J12p-$D.dat -d $DEBUG $BACKGROUND
153  fi
154  fi
155  done
156 
157  if [[ ! -f $WORKDIR/J13p.dat ]]; then
158  eval nohup JMakePDG -F13 -A $absorptionLengthFactor -S $scatteringLengthFactor -o $WORKDIR/J13p.dat -d $DEBUG $BACKGROUND
159  fi
160 
161  let DMAX_M="${RMAX_M}*1.5"
162 
163  for D in 0.1 0.5 1.0 5.0 10.0 20.0 30.0 40.0 50.0 60.0 70.0 80.0 90.0 100.0 \
164  120.0 150.0 170.0 190.0 210.0 230.0 250.0 270.0 290.0 310.0 \
165  340.0 370.0 400.0 430.0 460.0 490.0 520.0 550.0 600.0 650.0 700.0 750.0 800.0; do
166 
167  if (( $D <= $DMAX_M )); then
168  if [[ ! -f $WORKDIR/J14p.dat ]]; then
169  eval nohup JMakePDG -F14 -A $absorptionLengthFactor -S $scatteringLengthFactor -D $D -o $WORKDIR/J14p-$D.dat -d $DEBUG $BACKGROUND
170  fi
171  fi
172  done
173 
174  if [[ ! -f $WORKDIR/J23p.dat ]]; then
175  eval nohup JMakePD0 -F23 -A $absorptionLengthFactor -S $scatteringLengthFactor -o $WORKDIR/J23p.dat -d $DEBUG $BACKGROUND
176  fi
177  if [[ ! -f $WORKDIR/J24p.dat ]]; then
178  eval nohup JMakePD0 -F24 -A $absorptionLengthFactor -S $scatteringLengthFactor -o $WORKDIR/J24p.dat -d $DEBUG $BACKGROUND
179  fi
180  fi
181 
182 fi
183 
184 
185 if (( $MERGE_PDF )); then
186 
187  if [[ ! -f $WORKDIR/J3p.dat ]]; then
188 
189  FILE_LIST=`ls $WORKDIR/J3p-*.dat`
190 
191  JMergePDF -f $FILE_LIST -o $WORKDIR/J3p.dat -d $DEBUG
192  fi
193 
194  if [[ ! -f $WORKDIR/J4p.dat ]]; then
195 
196  FILE_LIST=`ls $WORKDIR/J4p-*.dat`
197 
198  JMergePDF -f $FILE_LIST -o $WORKDIR/J4p.dat -d $DEBUG
199  fi
200 
201  if [[ ! -f $WORKDIR/J6p.dat ]]; then
202 
203  FILE_LIST=`ls $WORKDIR/J6p-*.dat`
204 
205  JMergePDF -f $FILE_LIST -o $WORKDIR/J6p.dat -d $DEBUG
206  fi
207 
208  if [[ ! -f $WORKDIR/J12p.dat ]]; then
209 
210  FILE_LIST=`ls $WORKDIR/J12p-*.dat`
211 
212  JMergePDG -f $FILE_LIST -o $WORKDIR/J12p.dat -d $DEBUG
213  fi
214 
215  if [[ ! -f $WORKDIR/J14p.dat ]]; then
216 
217  FILE_LIST=`ls $WORKDIR/J14p-*.dat`
218 
219  JMergePDG -f $FILE_LIST -o $WORKDIR/J14p.dat -d $DEBUG
220  fi
221 fi
222 
223 
224 if (( $MAKE_CDF )); then
225 
226  check_process JMakeCDF JMakeCDG
227 
228  if [[ ! -f $WORKDIR/I1p.dat ]]; then
229  JMakeCDF -f $WORKDIR/J1p.dat -e 1e-5 -o $WORKDIR/I1p.dat -d $DEBUG
230  fi
231  if [[ ! -f $WORKDIR/I2p.dat ]]; then
232  JMakeCDF -f $WORKDIR/J2p.dat -e 1e-4 -o $WORKDIR/I2p.dat -d $DEBUG
233  fi
234  if [[ ! -f $WORKDIR/I5p.dat ]]; then
235  JMakeCDF -f $WORKDIR/J5p.dat -e 1e-5 -o $WORKDIR/I5p.dat -d $DEBUG
236  fi
237  if [[ ! -f $WORKDIR/I6p.dat ]]; then
238  JMakeCDF -f $WORKDIR/J6p.dat -e 1e-5 -o $WORKDIR/I6p.dat -d $DEBUG
239  fi
240  if [[ ! -f $WORKDIR/I12p.dat ]]; then
241  JMakeCDG -f $WORKDIR/J12p.dat -e 1e-5 -o $WORKDIR/I12p.dat -d $DEBUG
242  fi
243  if [[ ! -f $WORKDIR/I13p.dat ]]; then
244  JMakeCDG -f $WORKDIR/J13p.dat -e 1e-5 -o $WORKDIR/I13p.dat -d $DEBUG
245  fi
246  if [[ ! -f $WORKDIR/I14p.dat ]]; then
247  JMakeCDG -f $WORKDIR/J14p.dat -e 1e-5 -o $WORKDIR/I14p.dat -d $DEBUG
248  fi
249 fi
250 
251 if (( $CLEAN )); then
252 
253  if ( rm -f $WORKDIR/J*p-*.dat ); then; fi
254 
255 fi
256 
257 status "End $script"
JTOOLS::w
data_type w[N+1][M+1]
Definition: JPolint.hh:757
M
do $JPP JMEstimator M
Definition: JMEstimator.sh:37
txt
then usage $script< directory A >< directoryB > nIn each there should be a file detector datx and tripod txt
Definition: JAlignDetector.sh:25
f
o $QUALITY_ROOT d $DEBUG!JPlot1D f
Definition: JDataQuality.sh:66
absorptionLengthFactor
double absorptionLengthFactor
Scaling of absorption and scattering length.
Definition: JDrawPD0.cc:24
exit
exit
Definition: JPizza.sh:36
JTOOLS::r
data_type r[M+1]
Definition: JPolint.hh:758
echo
then echo
Definition: JDOMDAQDriver.sh:56
JPHYSICS::C
static const double C
Physics constants.
Definition: JPhysics/JConstants.hh:25
JTOOLS::n
const int n
Definition: JPolint.hh:676
set_variable
then echo Variable JPP_DIR undefined exit fi source $JPP_DIR setenv sh $JPP_DIR &dev null set_variable
Definition: JAcoustics.sh:20
done
set_variable NUMBER_OF_ITERATIONS set_variable EPSILON cat acoustics_fit_parameters txt<< EOF $CONFIGURATION[*]Nmin=3;sigma_s=100.0e-6;stdev=10.0;mestimator=0;fixStrings=0;EOF for STRING in $STRINGS[*];do#fit stretching and(z) position of given string set_variable DETECTOR_TMP ${TMPDIR:-/tmp}/detector_A.datx JEditDetector-a $DETECTOR-o $DETECTOR_TMP-r $STRING JEditDetector-a $DETECTOR-o $DETECTOR-k $STRING for MUL in 0.005 0.001;do DX_M=0.2 for((N=0;$N< $NUMBER_OF_ITERATIONS;++N));do CHI2[3]=$CHI2[1] fitPositionOfString $STRING Z $DX_M fitStretchingOfString $STRING $MUL if(($CHI2[3]-$CHI2[1]< $EPSILON));then break fi done if(($N >=$NUMBER_OF_ITERATIONS));then printf"warning: reached maximum number of iterations %d - converenge %7.3f\n"$N $(($CHI2[3]-$CHI2[1])) fi done JMergeDetector-a $DETECTOR-a $DETECTOR_TMP-o $DETECTOR rm-f $DETECTOR_TMP JConvertDetectorFormat-a $DETECTOR-o $DETECTOR-r-d 0 > &dev null done
Definition: acoustics-fit.sh:563
WORKDIR
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
usage
* usage
Definition: JLigier-local.sh:34
o
then JConvertDetectorFormat a $DETECTOR[1] o
Definition: JCompareDetector.sh:45
print
print
Definition: JConvertDusj.sh:44
R
then usage $script[distance] fi case set_variable R
Definition: JDrawLED.sh:43
S
do set_variable SIGMA_NS set_variable OUTLIERS set_variable OUTPUT_FILE matrix[${ALPHA_DEG}\deg\] root $JPP JMatrixNZ a $DETECTOR f $INPUT_FILE o $OUTPUT_FILE S
Definition: JMatrixNZ.sh:58
root
do JPlot2D f $WORKDIR detector root
Definition: JCompareDetector.sh:95
d
then JMuonMCEvt f $INPUT_FILE o $INTERMEDIATE_FILE d
Definition: JMuonPath.sh:47
$script
then echo Variable JPP_DIR undefined exit fi source $JPP_DIR setenv sh $JPP_DIR set_variable DEBUG set_variable WORKDIR if do_usage *then usage $script(input file)+" fi if (( $ fatal "Invalid number of arguments." fi JCookie.sh set_array INPUT_FILES $argv[1
DEBUG
then usage $script< input_file >< detector_file > fi set_variable OUTPUT_DIR set_variable SELECTOR JDAQTimesliceL1 set_variable DEBUG case set_variable DEBUG
Definition: JCalibratePMT.sh:21
scatteringLengthFactor
double scatteringLengthFactor
Definition: JDrawPD0.cc:25
$WORKDIR
then display $WORKDIR
Definition: plot-Domino.sh:128
BUFFER
then echo Variable JPP_DIR undefined exit fi source $JPP_DIR setenv sh $JPP_DIR &dev null set_variable DEBUG set_variable WORKDIR $JPP_DATA set_variable MAKE_PDF set_variable MERGE_PDF set_variable MAKE_CDF set_variable CLEAN set_variable BACKGROUND set_variable WGET set_variable RMAX_M set_variable absorptionLengthFactor set_variable scatteringLengthFactor BUFFER
Definition: JMakePDF.sh:36
c
$WORKDIR ev_configure_domsimulator txt echo process $DOM_SIMULATOR $i $SOURCE_HOST[$index] csh c(setenv ROOTSYS $ROOTSYS &&source $JPP_DIR/setenv.csh $JPP_DIR &&($DOM_SIMULATOR\-u\$NAME\$\-H\$SERVER\$\-M\$LOGGER\$\-d $DEBUG</dev/null > &/dev/null &))'
in
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 source JAcoustics sh $DETECTOR_ID CHECK_EXIT_CODE typeset A TRIPODS get_tripods $WORKDIR tripod txt TRIPODS for EMITTER in
Definition: JCanberra.sh:42
A
source $JPP_DIR setenv csh $JPP_DIR &dev null eval JShellParser o a A
Definition: JShellParser.csh:15
script
script
Definition: JAcoustics.sh:2
D
do echo Generating $dir eval D
Definition: JDrawLED.sh:53
null
esac $JPP_BIN JLogger sh $LOGGER until pgrep JGetMessage</dev/null > dev null
Definition: examples/JLogger/JLogger.sh:33
P
then $DIR JPlotNPE PDG P
Definition: JPlotNPE-PDG.sh:62
Generated by   doxygen 1.8.5