Jpp  15.0.3
the software that should make you happy
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
hydrophone-Z:fit.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 FORMULA "[0] + [1]*(sqrt(1.0 + 0.5*(x-[2])*(x-[2])/([3]*[3])) - 1.0)"
14 set_variable FORMULA "[0] + [1]*((x-[2])*(x-[2])/([3]*[3]))"
15 set_variable: FORMAT GRAPHICS_FORMAT gif
16 set_variable+ BATCH GRAPHICS_BATCH -B
17 
18 if do_usage $*; then
19  usage "$script <string identifier> [<detector file>]"
20 fi
21 
22 if (( $# == 0 || $# > 2 )); then
23  fatal "Wrong number of arguments."
24 fi
25 
27 
29 set_variable HYDROPHONE_ROOT $WORKDIR/hydrophonez_${STRING}.root
30 
31 if [[ ! -f $HYDROPHONE_TXT ]]; then
32  fatal "No hydrophone data file $HYDROPHONE_TXT."
33 fi
34 
35 sort -gr -k 2,2 $HYDROPHONE_TXT | tail -1 | read X Y
36 
37 JGraph \
38  -f $HYDROPHONE_TXT \
39  -o $HYDROPHONE_ROOT
40 
41 let XMIN="$X - 2.5"
42 let XMAX="$X + 2.5"
43 
44 JFit \
45  -f ${HYDROPHONE_ROOT}:\.\* \
46  -o $WORKDIR/fit.root \
47  -F "$FORMULA" \
48  -x "$XMIN $XMAX" \
49  -@"p0 = $Y" \
50  -@"p1 = 0.001" \
51  -@"p2 = $X" \
52  -@"p3 = 0.3" \
53  -d $DEBUG
54 
55 set_array P2 `JPrintFit -f $WORKDIR/fit.root:\.\* -@p2`
56 
57 printf "result %7.3f +/- %7.3f [m]\n" $P2[1] $P2[2]
58 
59 if (( $# == 2 )); then
60 
62  set_variable MODULE `getModule -a $DETECTOR -L "$STRING 0"`
63  set_variable Z $P2[1]
64 
65  JEditDetector -a $DETECTOR -o $DETECTOR -M "$MODULE add 0.0 0.0 $Z"
66 fi
67 
68 mv $WORKDIR/fit.root $HYDROPHONE_ROOT
69 
70 typeset -Z 4 STRING
71 
72 JPlot1D \
73  -f ${HYDROPHONE_ROOT}:\.\* \
74  -\> "#Deltaz [m]" \
75  -\^ "RMS [ms]" \
76  -T "[$STRING]" \
77  -o $WORKDIR/hydrophonez_${STRING}.$FORMAT $BATCH
78 
79 rm -f $HYDROPHONE_ROOT
80 rm -f ${TMPDIR:-/tmp}/hydrophone.txt
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]
then fatal No hydrophone data file $HYDROPHONE_TXT fi sort gr k
do $JPP JMEstimator M
Definition: JMEstimator.sh:37
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!JPlot1D f
Definition: JDataQuality.sh:66
exit
Definition: JPizza.sh:36
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 fatal Wrong number of arguments fi set_variable STRING $argv[1] set_variable HYDROPHONE_TXT $WORKDIR hydrophonez_
then echo
then echo Variable JPP_DIR undefined exit fi source $JPP_DIR setenv sh $JPP_DIR &dev null set_variable
Definition: JAcoustics.sh:20
do set_variable OUTPUT_DIRECTORY $WORKDIR T
then break fi done getCenter read X Y Z let X
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
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
then awk F
* usage
then JConvertDetectorFormat a $DETECTOR[1] o
&set_variable FORMULA
p2
Definition: module-Z:fit.sh:74
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
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 let XMIN
set_variable DETECTOR
script
Definition: JAcoustics.sh:2
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