4 if [ -z $JPP_DIR ];
then
5 echo "Variable JPP_DIR undefined."
9 source $JPP_DIR/setenv.sh $JPP_DIR >& /dev/
null
13 set_array MULS -2.0e-2 -1.7e-2 -1.6e-2 -1.5e-2 -1.4e-2 -1.3e-2 -1.2e-2 -1.1e-2 -1.0e-2 -0.9e-2 -0.9e-2 -0.8e-2 -0.5e-2 +0.0e-2 +0.5e-2
19 usage "$script <detector specific acoustics-fit script> <option>"\
20 "\nAuxiliary script to make scan of pre-stretching of detector strings (see JEditDetector)."\
21 "\nPossible options: run, plot or fit."
24 if (( $# < 1 || $# > 2));
then
25 fatal
"Wrong number of arguments."
34 fatal
"Invalid option $OPTION."
39 for MUL
in $MULS[*];
do
46 # set varibales in script
51 nohup $SCRIPT >& all.log &
55 if [[ $OPTION ==
"plot" || $OPTION ==
"fit" ]];
then
58 set_variable ACOUSTICS_ZMUL_ROOT /tmp/acoustics-zmul.root
60 rm -
f $ACOUSTICS_ZMUL_TXT >& /dev/
null
62 for MUL
in $MULS[*];
do
67 tail -1 $FILE | grep
"Elapsed time" >& /dev/
null
69 if (( $? == 0 ));
then
71 set_array RESULT `grep -
v "warning" $FILE | tail -2 | head -1`
73 printf
"%4.1e %6.3f\n" $MUL $RESULT[-1] >> $ACOUSTICS_ZMUL_TXT
78 -
f $ACOUSTICS_ZMUL_TXT \
79 -
o $ACOUSTICS_ZMUL_ROOT
81 XMIN=`sort -k1g $ACOUSTICS_ZMUL_TXT | head -1 | awk
'{ print $1 }'`
82 XMAX=`sort -k1g $ACOUSTICS_ZMUL_TXT | tail -1 | awk
'{ print $1 }'`
84 let XMIN=
"$XMIN - 1.0e-3"
85 let XMAX=
"$XMAX + 1.0e-3"
87 if [[ $OPTION ==
"fit" ]];
then
89 sort -gr -k2gr $ACOUSTICS_ZMUL_TXT | tail -1 |
read X Y
92 -
f ${ACOUSTICS_ZMUL_ROOT}:\.\* \
101 set_array P2 `JPrintFit -
f /tmp/fit.root:\.\* -@
p2`
103 printf
"Optimal stretching factor %8.2e\n" $P2[1]
105 mv /tmp/fit.root $ACOUSTICS_ZMUL_ROOT
110 -
f $ACOUSTICS_ZMUL_ROOT:\.\* \
115 -
o acoustics-zmul.$FORMAT $BATCH
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 usage $script< detector >< run >< outputfile > fi case set_variable OPTION
then cat $TRIPOD_INITIAL<< EOF1 256877.5 4743716.7-2438.42 256815.5 4743395.0-2435.53 257096.2 4743636.0-2439.5EOFfiJEditDetector-a $DETECTOR_INITIAL-@"z = -2441.5"-s"-1 mul $MUL"-o $DETECTORcp-p $TRIPOD_INITIAL $TRIPODJAcoustics.sh $DETECTOR_IDJAcousticsEventBuilder.sh $DETECTOR $RUNS[*]cd $WORKDIRif[!$HOMEDIR-ef $WORKDIR];then cp-p $HOMEDIR/$TRIPOD $WORKDIR cp-p $HOMEDIR/$DETECTOR $WORKDIR for RUN in $RUNS[*];do cp-p $HOMEDIR/KM3NeT_ ${(l:8::0::0:) DETECTOR_ID}_ ${(l:8::0::0:) RUN}_event.root $WORKDIR done JAcoustics.sh $DETECTOR_IDfifor STAGE in`echo $STAGES`;do $JPP_DIR/examples/JAcoustics/acoustics-fit.sh $DETECTOR $TRIPOD $STAGE $WORKDIR/KM3NeT_ ${(l:8::0::0:) DETECTOR_ID}_ *event.root > $HOMEDIR stage
std::istream & read(std::istream &in, JTestSummary &summary, const char delimiter= ' ')
Read test summary.
then usage $script< detector file >< inputfile > nUtility script to create PDF and CDF of transition time distribution fi case set_variable WORKDIR
do echo Generating $dir eval DIR
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
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 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 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
*fatal Wrong option $OPTION
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