8 # ------------------------------------------------------------------------------------------
10 # Utility script to plot JDomino results.
12 # ------------------------------------------------------------------------------------------
14 if [ -z $JPP_DIR ];
then
15 echo "Variable JPP_DIR undefined."
19 source $JPP_DIR/setenv.sh $JPP_DIR
28 usage "$script (input file)+"
31 if (( $# == 0 ));
then
32 fatal
"No input files."
40 -
f "${INPUT_FILES//\.root/.root:job}" \
45 -
T "Average number of photo-electrons per event by a track_in" \
53 -
f "${INPUT_FILES//\.root/.root:trks$}" \
55 -\^
"number of tracks" \
58 -
T "Average occurrence of a track_in per event" \
66 -
f "${INPUT_FILES//\.root/.root:hits_per_E_in}" \
73 -
T "Average number of photo-elecrons per event inside the instrumented volume" \
80 -
f "${INPUT_FILES//\.root/.root:pmt}" \
84 -\^
"number of hits [a.u.]" \
86 -
T "number of photo-electrons per PMT" \
99 for (( XL = $XMIN; $XL < $XMAX; XL += $DX ));
do
103 LABEL=`printf
'10^{%3.1f} GeV < E_{#nu} < 10^{%3.1f} GeV' $XL $XR`
106 -
f "${INPUT_FILES//\.root/.root:${TH2}}" \
110 -z
"$ZMIN $ZMAX" -Z \
121 $WORKDIR/${TH2}_*.gif \
$WORKDIR/${TH2}.gif >& /dev/
null
125 if [[
"$BACKGROUND" !=
"-B" ]];
then
139 for (( XL = $XMIN; $XL < $XMAX; XL += $DX ));
do
143 LABEL=`printf
'10^{%3.1f} GeV < E_{#nu} < 10^{%3.1f} GeV' $XL $XR`
146 -
f "${INPUT_FILES//\.root/.root:${TH2}}" \
149 -z
"$ZMIN $ZMAX" -Z \
150 -\>
"cos(#theta_{0})" \
160 $WORKDIR/${TH2}_*.gif \
$WORKDIR/${TH2}.gif >& /dev/
null
164 if [[
"$BACKGROUND" !=
"-B" ]];
then
180 for (( XL = $XMIN; $XL < $XMAX; XL += $DX ));
do
184 LABEL=`printf
'%3.0f m < D < %3.0f m' $XL $XR`
186 CD=`printf
'%.0f' ${XL}`
189 -
f "${INPUT_FILES//\.root/.root:${TH2}}" \
192 -z
"$ZMIN $ZMAX" -Z \
193 -\>
"cos(#theta_{0})" \
203 $WORKDIR/${TH2}_*.gif \
$WORKDIR/${TH2}.gif >& /dev/
null
207 if [[
"$BACKGROUND" !=
"-B" ]];
then
223 for (( XL = $XMIN; $XL < $XMAX; XL += $DX ));
do
227 LABEL=`printf
'%3.0f m < D < %3.0f m' $XL $XR`
229 CD=`printf
'%.0f' ${XL}`
232 -
f "${INPUT_FILES//\.root/.root:${TH2}}" \
235 -z
"$ZMIN $ZMAX" -Z \
236 -\>
"cos(#theta_{1})" \
246 $WORKDIR/${TH2}_*.gif \
$WORKDIR/${TH2}.gif >& /dev/
null
250 if [[
"$BACKGROUND" !=
"-B" ]];
then
266 for (( XL = $XMIN; $XL < $XMAX; XL += $DX ));
do
270 LABEL=`printf
'%3.0f m < D < %3.0f m' $XL $XR`
272 CD=`printf
'%.0f' ${XL}`
275 -
f "${INPUT_FILES//\.root/.root:${TH2}}" \
279 -z
"$ZMIN $ZMAX" -Z \
290 $WORKDIR/${TH2}_*.gif \
$WORKDIR/${TH2}.gif >& /dev/
null
294 if [[
"$BACKGROUND" !=
"-B" ]];
then
309 for (( XL = $XMIN; $XL < $XMAX; XL += $DX ));
do
313 LABEL=`printf
'10^{%3.1f} GeV < E_{#mu} < 10^{%3.1f} GeV' $XL $XR`
316 -
f "${INPUT_FILES//\.root/.root:${TH2}}" \
320 -z
"$ZMIN $ZMAX" -Z \
331 $WORKDIR/${TH2}_*.gif \
$WORKDIR/${TH2}.gif >& /dev/
null
335 if [[
"$BACKGROUND" !=
"-B" ]];
then
351 for (( XL = $XMIN; $XL < $XMAX; XL += $DX ));
do
355 LABEL=`printf
'%3.0f m < R < %3.0f m' $XL $XR`
357 CD=`printf
'%.0f' ${XL}`
360 -
f "${INPUT_FILES//\.root/.root:${TH2}}" \
363 -z
"$ZMIN $ZMAX" -Z \
364 -\>
"cos(#theta_{1})" \
374 $WORKDIR/${TH2}_*.gif \
$WORKDIR/${TH2}.gif >& /dev/
null
378 if [[
"$BACKGROUND" !=
"-B" ]];
then
394 for (( XL = $XMIN; $XL < $XMAX; XL += $DX ));
do
398 LABEL=`printf
'%3.0f m < R < %3.0f m' $XL $XR`
400 CD=`printf
'%.0f' ${XL}`
403 -
f "${INPUT_FILES//\.root/.root:${TH2}}" \
407 -z
"$ZMIN $ZMAX" -Z \
418 $WORKDIR/${TH2}_*.gif \
$WORKDIR/${TH2}.gif >& /dev/
null
422 if [[
"$BACKGROUND" !=
"-B" ]];
then
428 void set_variable(const std::string &name, const std::string &value)
Set environment variable.
do set_array DAQHEADER JPrintDAQHeader f
then fatal Wrong number of arguments fi set_variable STRING $argv[1] set_variable DETECTORXY_TXT $WORKDIR $DETECTORXY_TXT tail read X Y CHI2 RMS printf optimum n $X $Y $CHI2 $RMS awk v Y
do set_variable OUTPUT_DIRECTORY $WORKDIR T
then echo Variable JPP_DIR undefined exit fi source $JPP_DIR setenv sh $JPP_DIR zmodload zsh mathfunc set_variable DEBUG SIRENE_DEBUG
then fatal Not enough arguments fi set_variable DETECTOR $argv[1] INPUT_FILES
then fatal Wrong number of arguments fi set_variable STRING $argv[1] set_variable DETECTOR_TXT $WORKDIR $DETECTOR_TXT tail read X Y Z RMS let XMIN
*set_variable DETECTOR $argv[1] set_variable INPUT_FILE $argv[2,-2] set_variable OUTPUT_FILE $argv[-1]
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 usage $script< detector file >< inputfile > fi case set_variable WORKDIR
esac $JPP_BIN JLogger sh $LOGGER until pgrep JGetMessage</dev/null > dev null
#define DEBUG(A)
Message macros.