6 if [ -z $JPP_DIR ];
then
7 echo "Variable JPP_DIR undefined."
11 source $JPP_DIR/setenv.sh $JPP_DIR >& /dev/
null
23 usage "$script <detector file> <input file>"\
24 "\nInput file corresponds to the output of JKatoomba[.sh]."
27 if (( $# != 2 ));
then
28 fatal
"Wrong number of arguments."
34 eval `JPrintDetector -
a $DETECTOR -O IDENTIFIER`
37 source JAcoustics.sh $DETECTOR_ID
44 get_tripods
$WORKDIR/tripod.txt TRIPODS
46 for EMITTER
in ${(@
k)TRIPODS};
do
51 -
o "$WORKDIR/canberra[${EMITTER}].root" \
63 -
f "${HISTOGRAMS//\.root/.root:.*toa}" \
70 let
XMIN=
"$X - 3.0*$RMS"
71 let XMAX=
"$X + 3.0*$RMS"
76 -@
"p0 = GetMaximum" \
83 let
XMIN=
"$TMIN_MS * 1.0e-3"
84 let XMAX=
"$TMAX_MS * 1.0e-3"
98 kill_child_processes_at_exit
109 echo -
n "Creating graphics for string $STRING.."
113 for EMITTER
in ${(@
k)TRIPODS};
do
120 for (( FLOOR = $FIRST_FLOOR; $FLOOR <= $LAST_FLOOR; FLOOR += 1 ));
do
127 if [[ -
n "$N" ]] && (( $N >= $NUMBER_OF_ENTRIES ));
then
132 let
XMIN=
"$X - 3.0*$RMS"
133 let XMAX=
"$X + 3.0*$RMS"
136 -
f $WORKDIR/canberra\[$EMITTER\].root:${MODULE} \
138 -@
"p0 = GetMaximum" \
148 set_array P1 `JPrintFit -
f $WORKDIR/fit.root:\.\* -@
p1`
149 set_array P2 `JPrintFit -
f $WORKDIR/fit.root:\.\* -@
p2`
151 ((P1[1] = $P1[1] * 1.0e3))
# [ms]
152 ((P2[1] = abs($P2[1]) * 1.0e3))
# [ms]
154 printf
"%2d %9.6f %9.6f\n" $FLOOR $P1[1] $P2[1] >> $STRING_TXT
170 -
f "$HISTOGRAMS[*]" \
173 -y
"$TMIN_MS $TMAX_MS" \
181 echo " -> $WORKDIR/string_${STRING}.$FORMAT"
then usage $script< detector file >< detectorfile > nIf the range of floors is the first detector file is aligned to the second before the comparison fi case set_variable RANGE $argv[3]
then fatal No hydrophone data file $HYDROPHONE_TXT fi sort gr k
do JCanberra a $DETECTOR f $INPUT_FILE o $WORKDIR canberra[${EMITTER}] root T $WORKDIR tripod txt V $WORKDIR sound_velocity txt M $WORKDIR mechanics txt E $EMITTER d $DEBUG!done HISTOGRAMS
clean eval JPrintDetector a $DETECTOR O IDENTIFIER eval JPrintDetector a $DETECTOR O SUMMARY set_variable STRING
then usage $script< detector file >< inputfile > nUtility script to create PDF and CDF of transition time distribution fi case set_variable WORKDIR
do set_variable STRING_TXT $WORKDIR string_
do JCanberra a $DETECTOR f $INPUT_FILE o $WORKDIR canberra[${EMITTER}\] root T $WORKDIR tripod txt V $WORKDIR sound_velocity txt M $WORKDIR mechanics txt E $EMITTER d $DEBUG!done kill_child_processes_at_exit attach getModule a $DETECTOR typeset Z STRING typeset Z FLOOR set_variable FIRST_FLOOR for STRING in $STRINGS[*]
then JShowerPostfit f $INPUT_FILE o $OUTPUT_FILE N
then fatal Invalid string $STRING
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
then JFileTuna f $INPUT_FILE
then JMuonMCEvt f $INPUT_FILE o $INTERMEDIATE_FILE d
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
then usage $script< input_file >< detector_file > fi set_variable OUTPUT_DIR set_variable SELECTOR JDAQTimesliceL1 set_variable DEBUG case set_variable DEBUG
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 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
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/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
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
const JModule & getModule(const JDetector &detector, const JModuleLocation &location)
find module with a given string and floor number