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