5 # utility script to produce PDF and CDF tables.
9 # ------------------------------------------------------------------------------------------
11 # Utility script to produce PDF and CDF tables.
13 # ------------------------------------------------------------------------------------------
15 if [ -z $JPP_DIR ];
then
16 echo "Variable JPP_DIR undefined."
20 source $JPP_DIR/setenv.sh $JPP_DIR >& /dev/
null
34 # check whether script is already running
39 if [[ -
n "${BUFFER/ */}" ]];
then
40 notice
"${script} already running <${BUFFER}>"
45 # Method to print usage of script.
49 print -u2
"$script \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> \
64 # Parse command line options.
66 while getopts
"d:W:PMCcBw:hA:S:R:" option;
do
83 check_variable ROOTSYS JPP_DIR
WORKDIR
87 status
"Start $script"
90 if [[ -
n $WGET ]];
then
93 wget -
A dat,txt,
root,evt.gz -
r -l 1 -nd -q $WGET
97 if (( $MAKE_PDF ));
then
99 check_process JMakePDF JMakePDG
104 eval nohup JMakePDF -F1 -
A $absorptionLengthFactor -
S $scatteringLengthFactor -
o $WORKDIR/J1p.dat -
d $DEBUG $BACKGROUND
107 eval nohup JMakePDF -F2 -
A $absorptionLengthFactor -
S $scatteringLengthFactor -
o $WORKDIR/J2p.dat -
d $DEBUG $BACKGROUND
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
116 if (( $R <= $RMAX_M ));
then
119 eval nohup JMakePDF -F3 -
A $absorptionLengthFactor -
S $scatteringLengthFactor -
R $R -
o $WORKDIR/J3p-$R.dat -
d $DEBUG $BACKGROUND
122 eval nohup JMakePDF -F4 -
A $absorptionLengthFactor -
S $scatteringLengthFactor -
R $R -
o $WORKDIR/J4p-$R.dat -
d $DEBUG $BACKGROUND
128 eval nohup JMakePDF -F5 -
A $absorptionLengthFactor -
S $scatteringLengthFactor -
o $WORKDIR/J5p.dat -
d $DEBUG $BACKGROUND
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
137 if (( $R <= $RMAX_M ));
then
140 eval nohup JMakePDF -F6 -
A $absorptionLengthFactor -
S $scatteringLengthFactor -
R $R -
o $WORKDIR/J6p-$R.dat -
d $DEBUG $BACKGROUND
146 eval nohup JMakePDG -F13 -
A $absorptionLengthFactor -
S $scatteringLengthFactor -
o $WORKDIR/J13p.dat -
d $DEBUG $BACKGROUND
149 let DMAX_M=
"${RMAX_M}*1.5"
151 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 \
152 120.0 150.0 170.0 190.0 210.0 230.0 250.0 270.0 290.0 310.0 \
153 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
155 if (( $D <= $DMAX_M ));
then
157 eval nohup JMakePDG -F14 -
A $absorptionLengthFactor -
S $scatteringLengthFactor -
D $D -
o $WORKDIR/J14p-$D.dat -
d $DEBUG $BACKGROUND
163 eval nohup JMakePD0 -F23 -
A $absorptionLengthFactor -
S $scatteringLengthFactor -
o $WORKDIR/J23p.dat -
d $DEBUG $BACKGROUND
166 eval nohup JMakePD0 -F24 -
A $absorptionLengthFactor -
S $scatteringLengthFactor -
o $WORKDIR/J24p.dat -
d $DEBUG $BACKGROUND
173 if (( $MERGE_PDF ));
then
179 JMergePDF -
f $FILE_LIST -
o $WORKDIR/J3p.dat -
d $DEBUG
186 JMergePDF -
f $FILE_LIST -
o $WORKDIR/J4p.dat -
d $DEBUG
193 JMergePDF -
f $FILE_LIST -
o $WORKDIR/J6p.dat -
d $DEBUG
200 JMergePDG -
f $FILE_LIST -
o $WORKDIR/J14p.dat -
d $DEBUG
205 if (( $MAKE_CDF ));
then
207 check_process JMakeCDF JMakeCDG
229 if (( $CLEAN ));
then
do echo Generating $dir eval D
then JPlot1D f $WORKDIR postfit[prefit\] root
double absorptionLengthFactor
Scaling of absorption and scattering length.
static const double C
Physics constants.
then echo Variable JPP_DIR undefined exit fi source $JPP_DIR setenv sh $JPP_DIR &dev null set_variable
then echo Variable JPP_DIR undefined exit fi source $JPP_DIR setenv sh $JPP_DIR set_variable NORTH set_variable EAST set_variable SOUTH set_variable WEST set_variable WORKDIR tmp set_variable R set_variable CT set_variable YMAX set_variable YMIN if do_usage *then usage $script[distance] fi case set_variable R
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
then JMuonMCEvt f $INPUT_FILE o $INTERMEDIATE_FILE d
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
alias put_queue eval echo n
then usage $script< input_file >< detector_file > fi set_variable OUTPUT_DIR set_variable SELECTOR JDAQTimesliceL1 set_variable DEBUG case set_variable DEBUG
double scatteringLengthFactor
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
then echo Variable JPP_DIR undefined exit fi source $JPP_DIR setenv sh $JPP_DIR &dev null set_variable DEBUG set_variable WORKDIR tmp set_variable JDAQ_TIMESLICE JDAQTimeslice if do_usage *then usage $script< detector file >< inputfile > fi case set_variable WORKDIR
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
source $JPP_DIR setenv csh $JPP_DIR eval JShellParser o a A
esac $JPP_BIN JLogger sh $LOGGER until pgrep JGetMessage</dev/null > dev null