Jpp  15.0.3
the software that should make you happy
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
software/JCompass/JCompass.sh
Go to the documentation of this file.
1 #!/bin/zsh
2 #
3 #
4 # \author mdejong
5 #
6 version=1.0
7 script=${0##*/}
8 
9 if [ -z $JPP_DIR ]; then
10  echo "Variable JPP_DIR undefined."
11  exit
12 fi
13 
14 source $JPP_DIR/setenv.sh $JPP_DIR
15 
16 set_variable: DEBUG COMPASS_DEBUG 2
17 set_variable: WORKDIR COMPASS_WORKDIR ./
18 set_variable: TMAX_S COMPASS_TMAX_S 600
19 set_variable: SIGMA_DEG COMPASS_SIGMA_DEG 0.5
20 set_variable: STDEV COMPASS_STDEV 5.0
21 set_variable: MESTIMATOR COMPASS_MESTIMATOR 1
22 set_variable: NUMBER_OF_OUTLIERS COMPASS_NUMBER_OF_OUTLIERS 0
23 
24 if do_usage $*; then
25  usage "$script <detector file> (run[-run])+ <output file>"\
26  "\nAuxiliary script for in situ compass calibration."
27 fi
28 
29 if (( $# < 3 )); then
30  fatal "Wrong number of parameters."
31 fi
32 
34 set_array RUNS $argv[2,-2]
35 set_variable OUTPUT_FILE $argv[-1]
36 
37 set_variable AHRS_CALIBRATION_FILE $WORKDIR/ahrs_calibration.txt
38 
39 if [[ ! ${DETECTOR##*\.} == "datx" ]]; then
40  fatal "Detector file should be binary format."
41 fi
42 
43 eval `JPrintDetector -a $DETECTOR -O IDENTIFIER`
44 
45 typeset -a INPUT_FILES
46 
47 for RANGE in $RUNS[*]; do
48 
49  for (( RUN = ${RANGE%%-*}; $RUN <= ${RANGE##*-}; RUN += 1 )); do
50 
51  set_variable INPUT_FILE $WORKDIR/KM3NeT_${(l:8::0::0:)DETECTOR_ID}_${(l:8::0::0:)RUN}_ahrs.root
52 
53  if [[ ! -f $INPUT_FILE ]]; then
54 
55  echo -n "Writing $INPUT_FILE... "
56 
57  JConvertDB \
58  -q "ahrs" \
59  -@ "detid = $DETECTOR_ID" \
60  -@ "minrun = $RUN" \
61  -@ "maxrun = $RUN" \
62  -o $INPUT_FILE \
63  -d $DEBUG --
64 
65  if (( $? == 0 )); then
66  echo "OK"
67  else
68  continue
69  fi
70  fi
71 
73 
74  done
75 done
76 
77 if [[ ! -f $AHRS_CALIBRATION_FILE ]]; then
78 
79  JAHRSCalibration \
80  -D $DETECTOR_ID \
81  -o $AHRS_CALIBRATION_FILE \
82  -d $DEBUG
83 fi
84 
86 
87 JCompass \
88  -a $DETECTOR \
89  -A \
90  -c $AHRS_CALIBRATION_FILE \
91  -f "$INPUT_FILES[*]" \
92  -o $OUTPUT_FILE \
93  -@ "Tmax_s = $TMAX_S" \
94  -@ "sigma_deg = $SIGMA_DEG" \
95  -@ "stdev = $STDEV" \
96  -@ "mestimator = $MESTIMATOR" \
97  -@ "numberOfOutliers = $NUMBER_OF_OUTLIERS" \
98  -d $DEBUG --!
99 
100 timer_stop
101 timer_print
then usage $script< detector file >< detectorfile > nIf the range of floors is the first detector file is aligned to the second before the comparison nIn only modules within the given range of floors are used fi case set_variable RANGE $argv[3]
set_variable INPUT_FILE
o $QUALITY_ROOT d $DEBUG!JPlot1D f
Definition: JDataQuality.sh:66
exit
Definition: JPizza.sh:36
then echo
const int n
Definition: JPolint.hh:660
then timer_start(JRandomTimesliceWriter-a ${DETECTOR}-o ${PIPE}-n ${NUMBER_OF_SLICES}-P ${PMT_FILE}-d ${DEBUG}-S ${SEED}-B"$BACKGROUND_HZ"-N"$RECYCLING"--!&&for((i=0;$i!=10;++i));do;echo 'EOF' > $PIPE;usleep 100000;done &&rm-f $PIPE)&JTimesliceReprocessor-a $
then echo Variable JPP_DIR undefined exit fi source $JPP_DIR setenv sh $JPP_DIR &dev null if do_usage *then usage $script< detector identifier > fi case set_variable DETECTOR_ID
Definition: JTDC.sh:11
then fatal Detector file should be binary format fi eval JPrintDetector a $DETECTOR O IDENTIFIER typeset a INPUT_FILES for RANGE in $RUNS[*]
then echo Variable JPP_DIR undefined exit fi source $JPP_DIR setenv sh $JPP_DIR &dev null set_variable
Definition: JAcoustics.sh:20
then rm i $OUTPUT_FILE fi let RUN
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 ${TMPDIR:-/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
then echo Variable JPP_DIR undefined exit fi source $JPP_DIR setenv sh $JPP_DIR &dev null set_variable DEBUG set_variable WORKDIR
Definition: JLegolas.sh:20
* usage
then JConvertDetectorFormat a $DETECTOR[1] o
then JFileTuna f $INPUT_FILE
Definition: file-Tuna.sh:66
do sput $STRING $FLOOR sget MODULE INPUT_FILES
Definition: JCanberra.sh:73
&set_variable OUTPUT_FILE
then JMuonMCEvt f $INPUT_FILE o $INTERMEDIATE_FILE d
Definition: JMuonPath.sh:47
then JCalibrateToT a
Definition: JTuneHV.sh:116
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
then display $WORKDIR
Definition: plot-Domino.sh:128
$WORKDIR ev_configure_domsimulator txt echo process $DOM_SIMULATOR $i $SOURCE_HOST[$index] csh c(setenv ROOTSYS $ROOTSYS &&source $JPP_DIR/setenv.csh $JPP_DIR &&($DOM_SIMULATOR\-u\$NAME\$\-H\$SERVER\$\-M\$LOGGER\$\-d $DEBUG</dev/null > &/dev/null &))'
set_variable DETECTOR
do set_variable INPUT_FILE $WORKDIR KM3NeT_
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
Definition: JCanberra.sh:42
source $JPP_DIR setenv csh $JPP_DIR &dev null eval JShellParser o a A
do echo Generating $dir eval D
Definition: JDrawLED.sh:53
version
Definition: JCalibratePMT.sh:7
set_array RUNS
Definition: JDataQuality.sh:23
do if[[!-f $ACOUSTICS_WORKDIR/${KEY}.txt]]
Definition: JAcoustics.sh:39