4 if [ -z $JPP_DIR ];
then
5 echo "Variable JPP_DIR undefined."
9 source $JPP_DIR/setenv.sh $JPP_DIR >& /dev/
null
14 set_variable FORMULA
"[0] + [1] * log(1.0 + 0.5*(x-[2])*(x-[2])/([3]*[3]))"
17 usage "$script [<tripod file>]"
21 fatal
"Wrong number of arguments."
24 source JAcoustics.sh --
30 if [[ ! -
f $TRIPOD3Z_TXT ]]
then
31 fatal
"No detector data file $TRIPOD3Z_TXT."
34 grep -
v '^#' $TRIPOD3Z_TXT | sort -gr -
k 4,4 | tail -1 |
read X Y Z
CHI2 RMS
36 printf
"optimum %7.3f %7.3f %7.3f %7.3f %7.3f\n" $X $Y $Z $CHI2 $RMS
38 awk
'$1 ~ /^#/ { for (i = 2; i <= NF; ++i ) { printf " %s", $i }}' $TRIPOD3Z_TXT |
read A B
C
40 awk -
v Y=$Y -
v Z=$Z
'$1 !~ /^#/ { if ( $2 == Y && $3 == Z ) { printf "%7.3f %7.3f %7.3f\n", $1, $4, $5 } }' $TRIPOD3Z_TXT >!
$WORKDIR/
tripod3z_${
A}.txt
41 awk -
v X=$X -
v Z=$Z
'$1 !~ /^#/ { if ( $1 == X && $3 == Z ) { printf "%7.3f %7.3f %7.3f\n", $2, $4, $5 } }' $TRIPOD3Z_TXT >!
$WORKDIR/
tripod3z_${B}.txt
42 awk -
v X=$X -
v Y=$Y
'$1 !~ /^#/ { if ( $1 == X && $2 == Y ) { printf "%7.3f %7.3f %7.3f\n", $3, $4, $5 } }' $TRIPOD3Z_TXT >!
$WORKDIR/
tripod3z_${
C}.txt
44 for TRIPOD
in $A $B $C;
do
68 END { print XMIN, XMAX }
70 if ( $1 < XMIN ) { XMIN = $1 }
71 if ( $1 > XMAX ) { XMAX = $1 }
75 if (( $X > $XMIN && $X < $XMAX ));
then
78 -
f ${TRIPOD_ROOT}:\.\* \
87 set_array P2 `JPrintFit -
f $WORKDIR/fit.root:\.\* -@
p2`
91 printf
"nominal %2d %7.3f %7.3f\n" $TRIPOD 0.0 `JPrintFit -
f $WORKDIR/fit.root:\.\* -x
"0.0"`
92 printf
"optimal %2d %7.3f %7.3f\n" $TRIPOD $X `JPrintFit -
f $WORKDIR/fit.root:\.\* -x
"$X"`
98 echo "Optimum value at limit; -> no fit."
99 printf
"nominal %2d %7.3f %7.3f\n" $TRIPOD 0.0 `awk -
v X=0.0
'{ if ($1 == X) { print $2 }}' $TRIPOD_TXT`
100 printf
"optimal %2d %7.3f %7.3f\n" $TRIPOD $X $Y
104 if (( $# == 1 ));
then
105 JEditTripod -
f $argv[1] -
T "$TRIPOD add 0.0 0.0 $X"
109 -
f ${TRIPOD_ROOT}:\.\* \
116 rm -
f $TRIPOD_TXT $TRIPOD_ROOT
$G1_TXT $G1_ROOT
then fatal No hydrophone data file $HYDROPHONE_TXT fi sort gr k
bool read(Vec &v, std::istream &is)
Read a Vec(tor) from a stream.
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 TRIPOD_TXT $WORKDIR tripod3z_
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
do set_variable OUTPUT_DIRECTORY $WORKDIR T
then break fi done getCenter read X Y Z let X
*set_variable DETECTOR $argv[1] set_variable INPUT_FILE $argv[2,-2] set_variable OUTPUT_FILE $argv[-1]
then JMuonMCEvt f $INPUT_FILE o $INTERMEDIATE_FILE d
then usage $script< input_file >< detector_file > fi set_variable OUTPUT_DIR set_variable SELECTOR JDAQTimesliceL1 set_variable DEBUG case set_variable DEBUG
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 TRIPOD_FILE $argv[2] set_variable STAGE $argv[3] set_array INPUT_FILES $argv[4,-1] eval JPrintDetector a $DETECTOR O IDENTIFIER eval JPrintDetector a $DETECTOR O SUMMARY source JAcoustics sh $DETECTOR_ID typeset A TRIPODS get_tripods $TRIPOD_FILE TRIPODS typeset a CHI2 CHI2
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
do typeset Z STRING typeset Z FLOOR set_variable MODULE_TXT $WORKDIR $MODULE_TXT tail read X ZMAX RMS sort gr $MODULE_TXT tail read X ZMIN RMS echo $X $ZMIN $G1_TXT
then fatal Invalid tripod $TRIPOD
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
Double_t g1(const Double_t x)
Function.