4 if [ -z $JPP_DIR ];
then
5 echo "Variable JPP_DIR undefined."
9 source $JPP_DIR/setenv.sh $JPP_DIR >& /dev/
null
21 usage "$script <detector file> <string identifier> [<detector file>]"
24 if (( $# == 0 || $# > 3 ));
then
25 fatal
"Wrong number of arguments."
30 set_array QUANTILES 0.333 0.666
31 set_variable
FORMULA "[0] * exp(-0.5 * (x-[1])*(x-[1]) / ([2]*[2]))"
34 source JAcousticsToolkit.sh
38 get_tripods
$WORKDIR/tripod.txt TRIPODS
44 for TRIPOD
in ${(
k)TRIPODS[*]};
do
47 -
f "hydrophone.root:${MODULE}.${TRIPOD}" \
51 if (( ${#RESULT} == ${#QUANTILES} ));
then
53 P1=$((0.5*($RESULT[-1] + $RESULT[1]))) # average
54 P2=$((0.5*($RESULT[-1] - $RESULT[1]))) # width
57 -
f "hydrophone.root:${MODULE}.${TRIPOD}" \
60 -
@"p0 = GetMaximum * 0.8" \
68 VALUE=$((1.0e+6 * $P1[1])) # [us]
69 TOTAL=$(($TOTAL + 1.0))
71 if (( $TOTAL == 1.0 ));
then
76 STDEV[2]=$(($STDEV[1] + ($VALUE -
$XMEAN[1]) * ($VALUE -
$XMEAN[2])))
81 printf
"string %4d tripod %2d %5.0f %3.0f us\n" $STRING $TRIPOD $(($P1[1] * 1.0e6)) $(($P2[1] * 1.0e6))
85 if (( $TOTAL > 1.0 ));
then
86 printf
"string %4d mean RMS %5.0f %3.0f us\n" $STRING $XMEAN[1] $((sqrt($TOTAL * $STDEV[1] / ($TOTAL * ($TOTAL - 1.0)))))
92 -
f hydrophone.
root:$MODULE \
101 P1=$((0.5*($RESULT[-1] + $RESULT[1]))) # average
102 P2=$((0.5*($RESULT[-1] - $RESULT[1]))) # width
108 -
@"p0 = GetMaximum * 0.8" \
116 printf
"string %4d t0 %6.0f +/- %6.0f us\n" $STRING $(($P1[1] * 1.0e6)) $(($P1[2] * 1.0e6))
117 printf
"string %4d sigma %6.0f +/- %6.0f us\n" $STRING $(($P2[1] * 1.0e6)) $(($P2[2] * 1.0e6))
119 if (( ${#} == 3 ));
then
125 let
"TIME_OFFSET_NS = $P1[1] * 1.0e9"
127 printf
"Add time offset %12.6d to hydrophone.\n" $TIME_OFFSET_NS
132 -M
"$MODULE SUB $TIME_OFFSET_NS" \
140 -\^
"number of events [a.u.]" \
143 -
o hydrophonet0_${
STRING}.$FORMAT $BATCH
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]
then fatal No hydrophone data file $HYDROPHONE_TXT fi sort gr k
Q(UTCMax_s-UTCMin_s)-livetime_s
std::istream & read(std::istream &in, JTestSummary &summary, const char delimiter= ' ')
Read test summary.
clean eval JPrintDetector a $DETECTOR O IDENTIFIER eval JPrintDetector a $DETECTOR O SUMMARY set_variable STRING
o $QUALITY_ROOT d $DEBUG!CHECK_EXIT_CODE JPlot1D f
usr share Modules init zsh export TMPDIR
static const double H
Planck constant [eV s].
then fatal Wrong number of arguments fi JConvertDetectorFormat a o
then fatal Invalid string $STRING
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 echo Variable JPP_DIR undefined exit fi source $JPP_DIR setenv sh $JPP_DIR &dev null set_variable DEBUG set_variable WORKDIR
then usage $script< input file >[option[primary[working directory]]] nWhere option can be N
fi set_variable FILENAME $WORKDIR
then JMuonMCEvt f $INPUT_FILE o $INTERMEDIATE_FILE d
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 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
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 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
then printf string mean RMS us n $STRING $XMEAN((sqrt($TOTAL *$STDEV[1]/($TOTAL *($TOTAL-1.0))))) fi typeset-Z 4 STRING JOpera1D-f hydrophone.root
source $JPP_DIR setenv csh $JPP_DIR &dev null eval JShellParser o a A
then fatal Invalid tripod $TRIPOD
do JPlot2D f $WORKDIR canberra[${EMITTER}] root
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
#define DEBUG(A)
Message macros.