Jpp  15.0.3
the software that should make you happy
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
detector-Zmul+zoom:run.sh
Go to the documentation of this file.
1 #!/bin/zsh
2 script=${0##*/}
3 
4 if [ -z $JPP_DIR ]; then
5  echo "Variable JPP_DIR undefined."
6  exit
7 fi
8 
9 source $JPP_DIR/setenv.sh $JPP_DIR >& /dev/null
10 
11 set_variable: DEBUG ACOUSTICS_DEBUG 2
12 set_variable: WORKDIR ACOUSTICS_WORKDIR ./
13 set_variable: FIT ACOUSTICS_FIT 2
14 
15 if do_usage $*; then
16  usage "$script <string identifier> <detector file> (input file)+"
17 fi
18 
19 if (( $# < 3 )); then
20  fatal "Wrong number of arguments."
21 fi
22 
23 set_variable STRING $argv[1]
24 set_variable DETECTOR $argv[2]
25 set_array INPUT_FILES $argv[3,-1]
26 
27 eval `JPrintDetector -a $DETECTOR -O IDENTIFIER`
28 eval `JPrintDetector -a $DETECTOR -O SUMMARY`
29 
30 set_variable MODULE `getModule -a $DETECTOR -L "$STRING 0"`
31 
32 source JAcoustics.sh $DETECTOR_ID
33 
34 CHECK_EXIT_CODE
35 
36 if [[ -z "${STRINGS[(r)$STRING]}" ]]; then
37  fatal "Invalid string $STRING; "\
38  "possible values: ${STRINGS[*]}"
39 fi
40 
41 set_variable DETECTOR_TXT $WORKDIR/detectorzmul_${STRING}.txt
42 set_variable TMPDIR $WORKDIR/.$$
43 mkdir -p $TMPDIR
44 
45 for (( Z = -2.0; $Z <= +2.0; Z += 0.1 )); do
46  for (( MUL = -0.01; $MUL <= 0.01; MUL += 0.001 )); do
47 
48  rm -f $TMPDIR/detector.detx
49 
50  JEditDetector \
51  -a $DETECTOR \
52  -k "$STRING" \
53  -s "$STRING add 0.0 0.0 $Z" \
54  -s "$STRING mul $MUL" \
55  -o $TMPDIR/detector.detx
56 
57  rm -f $TMPDIR/katoomba.root
58 
59  JKatoomba \
60  -a $TMPDIR/detector.detx \
61  -f "$INPUT_FILES[*]" \
62  -o $TMPDIR/katoomba.root \
63  -T $WORKDIR/tripod.txt \
64  -V $WORKDIR/sound_velocity.txt \
65  -M $WORKDIR/mechanics.txt \
66  -@ $WORKDIR/acoustics_fit_parameters.txt \
67  -! $WORKDIR/disable.txt \
68  -F $FIT \
69  -u \
70  -d $DEBUG --!
71 
72  if (( `JPrintResult -f $TMPDIR/katoomba.root:chi2 -F "GetSumOfWeights"` > 1.0 )); then
73 
74  set_array RESULT `JPrintResult -f $TMPDIR/katoomba.root:chi2 -F "GetMean" -F "GetRMS / sqrt(GetSumOfWeights)"`
75 
76  printf "%7.3f %7.3f %7.3f %7.4f\n" $Z $MUL $RESULT[*] >> $DETECTOR_TXT
77  else
78  error "No result."
79  fi
80  done
81 done
82 
83 rm -rf $TMPDIR
detectorzmul_
then fatal Wrong number of arguments fi set_variable STRING $argv[1] set_variable DETECTOR_TXT $WORKDIR detectorzmul_
Definition: detector-Zmul+zoom:fit.sh:32
$argv
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]
Definition: JCompareDetector.sh:31
k
then fatal No hydrophone data file $HYDROPHONE_TXT fi sort gr k
Definition: hydrophone-phi:fit.sh:35
M
do $JPP JMEstimator M
Definition: JMEstimator.sh:37
STRING
clean eval JPrintDetector a $DETECTOR O IDENTIFIER eval JPrintDetector a $DETECTOR O SUMMARY set_variable STRING
Definition: examples/JTrigger/JTriggerEfficiencyRunByRun.sh:64
f
o $QUALITY_ROOT d $DEBUG!JPlot1D f
Definition: JDataQuality.sh:66
exit
exit
Definition: JPizza.sh:36
V
V(JDAQEvent-JTriggerReprocessor)*1.0/(JDAQEvent+1.0e-10)
echo
then echo
Definition: JDOMDAQDriver.sh:56
set_variable
then echo Variable JPP_DIR undefined exit fi source $JPP_DIR setenv sh $JPP_DIR &dev null set_variable
Definition: JAcoustics.sh:20
T
do set_variable OUTPUT_DIRECTORY $WORKDIR T
Definition: JCalibrateHeight.sh:61
done
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
Definition: acoustics-fit.sh:563
TMPDIR
then echo Variable JPP_DIR undefined exit fi source $JPP_DIR setenv sh $JPP_DIR if do_usage *then usage $script[(input file)+] fi set_variable DEBUG set_variable WORKDIR TMPDIR
Definition: examples/JReconstruction/JEnergyCorrection.sh:21
WORKDIR
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
F
then awk F
Definition: getMIIseacurrent.sh:86
usage
* usage
Definition: JLigier-local.sh:34
o
then JConvertDetectorFormat a $DETECTOR[1] o
Definition: JCompareDetector.sh:45
INPUT_FILES
do sput $STRING $FLOOR sget MODULE INPUT_FILES
Definition: JCanberra.sh:73
d
then JMuonMCEvt f $INPUT_FILE o $INTERMEDIATE_FILE d
Definition: JMuonPath.sh:47
a
then JCalibrateToT a
Definition: JTuneHV.sh:116
$DETECTOR
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
Definition: module-Z:fit.sh:118
DEBUG
then usage $script< input_file >< detector_file > fi set_variable OUTPUT_DIR set_variable SELECTOR JDAQTimesliceL1 set_variable DEBUG case set_variable DEBUG
Definition: JCalibratePMT.sh:21
$WORKDIR
then display $WORKDIR
Definition: plot-Domino.sh:128
Z
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
Definition: JAlignDetector.sh:81
DETECTOR
set_variable DETECTOR
Definition: JARCACalibrateMuon.sh:31
detector
do set_variable DETECTOR_TXT $WORKDIR detector
Definition: detector-XY:fit1d.sh:38
JTOOLS::u
double u[N+1]
Definition: JPolint.hh:739
script
script
Definition: JAcoustics.sh:2
null
esac $JPP_BIN JLogger sh $LOGGER until pgrep JGetMessage</dev/null > dev null
Definition: examples/JLogger/JLogger.sh:33
getModule
const JModule & getModule(const JDetector &detector, const JModuleLocation &location)
find module with a given string and floor number
Definition: JBeaconSimulator.cc:928
Generated by   doxygen 1.8.5