Jpp - the software that should make you happy
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
tripod-XY: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 
12 set_variable WORKDIR ${ACOUSTICS_WORKDIR:-./}
13 set_variable BACKGROUND ${ACOUSTICS_BACKGROUND:+-B}
14 set_variable FORMULA "[0] + [1] * (sqrt(1.0 + 0.5*(x-[2])*(x-[2])/([5]*[5]) + 0.5*(y-[3])*(y-[3])/([6]*[6]) + [4]*(x-[2])*(y-[3])/([5]*[6])) - 1.0)"
15 
16 if ( do_usage $* ); then
17  usage "$script <tripod identifier> [<tripod file>]"
18 fi
19 
20 if (( $# == 0 || $# > 2 )); then
21  fatal "Wrong number of arguments."
22 fi
23 
24 set_variable TRIPOD $argv[1]
25 
26 source JAcoustics.sh --
27 
28 CHECK_EXIT_CODE
29 
30 set_variable TRIPOD_TXT $WORKDIR/tripodxy_${TRIPOD}.txt
31 set_variable TRIPOD_ROOT $WORKDIR/tripodxy_${TRIPOD}.root
32 set_variable G1_TXT $WORKDIR/g1.txt
33 set_variable G1_ROOT $WORKDIR/g1.root
34 
35 if [[ ! -f $TRIPOD_TXT ]] then
36  fatal "No tripod data file $TRIPOD_TXT."
37 fi
38 
39 JGraph2D \
40  -f $TRIPOD_TXT \
41  -o $TRIPOD_ROOT
42 
43 sort -gr -k 3,3 $TRIPOD_TXT | tail -1 | read X Y Z RMS
44 
45 echo $X $Y > $G1_TXT
46 
47 JGraph \
48  -f $G1_TXT \
49  -o $G1_ROOT
50 
51 awk '
52  BEGIN {
53  XMIN = +1e3
54  XMAX = -1e3
55  YMIN = +1e3
56  YMAX = -1e3
57  }
58  END { print XMIN, XMAX, YMIN, YMAX }
59  $1 !~ /^#/ {
60  if ( $1 < XMIN ) { XMIN = $1 }
61  if ( $1 > XMAX ) { XMAX = $1 }
62  if ( $2 < YMIN ) { YMIN = $2 }
63  if ( $2 > YMAX ) { YMAX = $2 }
64  }
65  ' $TRIPOD_TXT | read XMIN XMAX YMIN YMAX
66 
67 if (( $X > $XMIN && $X < $XMAX && $Y > $YMIN && $Y < $YMAX )); then
68 
69  let XMIN="$X - 0.5"
70  let XMAX="$X + 0.5"
71  let YMIN="$Y - 0.5"
72  let YMAX="$Y + 0.5"
73 
74  JFit2D \
75  -f ${TRIPOD_ROOT}:\.\* \
76  -o $WORKDIR/fit.root \
77  -F "$FORMULA" \
78  -x "$XMIN $XMAX" \
79  -y "$YMIN $YMAX" \
80  -@ "p0 = $Z" \
81  -@ "p1 = 2.0" \
82  -@ "p2 = $X" \
83  -@ "p3 = $Y" \
84  -@ "p4 = 1.0" \
85  -@ "p5 = 5.0" \
86  -@ "p6 = 5.0" \
87  -R "p2 = $XMIN $XMAX" \
88  -R "p3 = $YMIN $YMAX" \
89  -d $DEBUG
90 
91  set_array P2 `JPrintFit -f $WORKDIR/fit.root:user -@p2`
92  set_array P3 `JPrintFit -f $WORKDIR/fit.root:user -@p3`
93 
94  set_variable X $P2[1]
95  set_variable Y $P3[1]
96 
97  printf "nominal %2d %7.3f %7.3f %7.3f\n" $TRIPOD 0.0 0.0 `JPrintFit -f $WORKDIR/fit.root:user -x "0.0 0.0"`
98  printf "optimal %2d %7.3f %7.3f %7.3f\n" $TRIPOD $X $Y `JPrintFit -f $WORKDIR/fit.root:user -x "$X $Y"`
99 
100  mv $WORKDIR/fit.root $TRIPOD_ROOT
101 
102 else
103 
104  echo "Optimum value at limit; -> no fit."
105  printf "nominal %2d %7.3f %7.3f %7.3f\n" $TRIPOD 0.0 0.0 `awk -v X=0.0 -v Y=0.0 '{ if ($1 == X && $2 == Y) { print $3 }}' $TRIPOD_TXT`
106  printf "optimal %2d %7.3f %7.3f %7.3f\n" $TRIPOD $X $Y $Z
107 
108 fi
109 
110 if (( $# == 2 )); then
111  JEditTripod -f $argv[2] -T "$TRIPOD add $X $Y 0.0"
112 fi
113 
114 JPlot2D \
115  -f $TRIPOD_ROOT:\.\* \
116  -f ${G1_ROOT}:\.\* \
117  -\> "#Deltax [m]" \
118  -\< "#Deltay [m]" \
119  -\^ "#chi^{2}/NDF" \
120  -O CONT3 \
121  -T "[$TRIPOD]" \
122  -o $WORKDIR/tripodxy_${TRIPOD}.gif $BACKGROUND
123 
124 rm -f $TRIPOD_ROOT $G1_TXT $G1_ROOT
void set_variable(const std::string &name, const std::string &value)
Set environment variable.
then fatal Wrong number of arguments fi set_variable TRIPOD $argv[1] source JAcoustics sh CHECK_EXIT_CODE set_variable TRIPOD_TXT $WORKDIR tripodxy_
bool read(Vec &v, std::istream &is)
Read a Vec(tor) from a stream.
Definition: io_ascii.hh:141
then fatal No sound hydrophone file $HYDROPHONE_TXT fi JGraph f $HYDROPHONE_TXT o $HYDROPHONE_ROOT sort gr k
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
fatal Invalid number of arguments esac let YMIN
Definition: testToT.sh:39
then JPizza f
Definition: JPizza.sh:46
then echo
do set_variable OUTPUT_DIRECTORY $WORKDIR T
then break fi done getCenter read X Y Z let X
do $DIR JTransitTime o
Definition: JTransitTime.sh:44
then awk F
* usage
then usage $script[distance] fi case set_variable R
Definition: JDrawLED.sh:40
then echo Variable JPP_DIR undefined exit fi source $JPP_DIR setenv sh $JPP_DIR &dev null set_variable ACOUSTICS_DEBUG
Definition: JAcoustics.sh:17
p2
Definition: module-Z:fit.sh:72
*set_variable DETECTOR $argv[1] set_variable INPUT_FILE $argv[2,-2] set_variable OUTPUT_FILE $argv[-1]
Definition: JDomino.sh:39
then JMuonMCEvt f $INPUT_FILE o $INTERMEDIATE_FILE d
Definition: JMuonPath.sh:45
then display $WORKDIR
Definition: plot-Domino.sh:127
then let XMIN
data_type v[N+1][M+1]
Definition: JPolint.hh:740
script
Definition: JAcoustics.sh:2
p3
Definition: module-Z:fit.sh:72
do typeset Z STRING typeset Z FLOOR set_variable MODULE_TXT $WORKDIR $MODULE_TXT tail read X ZMAX RMS sort gr $MODULE_TXT tail read X ZMIN RMS echo $X $ZMIN $G1_TXT
Definition: module-Z:fit.sh:33
then fatal Invalid tripod $TRIPOD
then usage $script< detector file >< inputfile > fi case set_variable WORKDIR
Definition: JLegolas.sh:28
esac $JPP_BIN JLogger sh $LOGGER until pgrep JGetMessage</dev/null > dev null
#define DEBUG(A)
Message macros.
Definition: JMessage.hh:62
Double_t g1(const Double_t x)
Function.
Definition: JQuantiles.cc:25