Jpp  18.5.2
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 FORMULA "[0] + [1] * log(1.0 + 0.5*(x-[2])*(x-[2])/([4]*[4]) + 0.5*(y-[3])*(y-[3])/([5]*[5]))"
14 set_variable FORMULA "[0] + [1] * ((x-[2])*(x-[2])/([4]*[4]) + (y-[3])*(y-[3])/([5]*[5]))"
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 
26 set_variable STRING $argv[1]
27 
28 set_variable DETECTOR_TXT $WORKDIR/detectorxy_${STRING}.txt
29 set_variable DETECTOR_ROOT $WORKDIR/detectorxy_${STRING}.root
30 set_variable G1_TXT $WORKDIR/g1.txt
31 set_variable G1_ROOT $WORKDIR/g1.root
32 
33 if [[ ! -f $DETECTOR_TXT ]] then
34  fatal "No detector data file $DETECTOR_TXT."
35 fi
36 
37 JGraph2D \
38  -f $DETECTOR_TXT \
39  -o $DETECTOR_ROOT
40 
41 sort -gr -k 3,3 $DETECTOR_TXT | tail -1 | read X Y Z RMS
42 
43 echo $X $Y > $G1_TXT
44 
45 JGraph \
46  -f $G1_TXT \
47  -o $G1_ROOT
48 
49 awk '
50  BEGIN {
51  XMIN = +1e3
52  XMAX = -1e3
53  YMIN = +1e3
54  YMAX = -1e3
55  }
56  END { print XMIN, XMAX, YMIN, YMAX }
57  $1 !~ /^#/ {
58  if ( $1 < XMIN ) { XMIN = $1 }
59  if ( $1 > XMAX ) { XMAX = $1 }
60  if ( $2 < YMIN ) { YMIN = $2 }
61  if ( $2 > YMAX ) { YMAX = $2 }
62  }
63  ' $DETECTOR_TXT | read XMIN XMAX YMIN YMAX
64 
65 if (( $X > $XMIN && $X < $XMAX && $Y > $YMIN && $Y < $YMAX )); then
66 
67  let "XMIN = $X - 0.5"
68  let "XMAX = $X + 0.5"
69  let "YMIN = $Y - 0.5"
70  let "YMAX = $Y + 0.5"
71 
72  JFit2D \
73  -f ${DETECTOR_ROOT}:\.\* \
74  -o $WORKDIR/fit.root \
75  -F "$FORMULA" \
76  -x "$XMIN $XMAX" \
77  -y "$YMIN $YMAX" \
78  -@ "p0 = $Z" \
79  -@ "p1 = 2.0" \
80  -@ "p2 = $X" \
81  -@ "p3 = $Y" \
82  -@ "p4 = 0.5" \
83  -@ "p5 = 0.5" \
84  -R "p2 = $XMIN $XMAX" \
85  -R "p3 = $YMIN $YMAX" \
86  -d $DEBUG
87 
88  set_array P2 `JPrintFit -f $WORKDIR/fit.root:user -@p2`
89  set_array P3 `JPrintFit -f $WORKDIR/fit.root:user -@p3`
90 
91  set_variable X $P2[1]
92  set_variable Y $P3[1]
93 
94  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"`
95  printf "optimal %4d %7.3f %7.3f %7.3f\n" $STRING $X $Y `JPrintFit -f $WORKDIR/fit.root:user -x "$X $Y"`
96 
97  mv $WORKDIR/fit.root $DETECTOR_ROOT
98 
99 else
100 
101  echo "Optimum value at limit; -> no fit."
102  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`
103  printf "optimal %4d %7.3f %7.3f %7.3f\n" $STRING $X $Y $Z
104 
105 fi
106 
107 if (( $# == 2 )); then
108 
109  set_variable DETECTOR $argv[2]
110 
111  JEditDetector -a $DETECTOR -S "$STRING add $X $Y 0.0" -o $DETECTOR
112 fi
113 
114 typeset -Z 4 STRING
115 
116 JPlot2D \
117  -f ${DETECTOR_ROOT}:\.\* \
118  -f ${G1_ROOT}:\.\* \
119  -\> "#Deltax [m]" \
120  -\< "#Deltay [m]" \
121  -\^ "#chi^{2}/NDF" \
122  -N "X 505" \
123  -N "Y 505" \
124  -O CONT3 \
125  -T "[$STRING]" \
126  -o $WORKDIR/detectorxy_${STRING}.$FORMAT $BATCH
127 
128 rm -f $DETECTOR_ROOT $G1_TXT $G1_ROOT
129 
130 
$argv
then usage $script< detector identifier >< startdate\"YYYY-MM-DDHH:MM:SS\"><finaldate\"YYYY-MM-DDHH:MM:SS\"><QA/QCfile> fi case set_variable QAQC_TXT $argv[4]
Definition: JDataMonitor.sh:24
k
then fatal No hydrophone data file $HYDROPHONE_TXT fi sort gr k
Definition: hydrophone-phi:fit.sh:35
JCOMPAREHISTOGRAMS::read
std::istream & read(std::istream &in, JTestSummary &summary, const char delimiter= ' ')
Read test summary.
Definition: JTestSummary.hh:222
STRING
clean eval JPrintDetector a $DETECTOR O IDENTIFIER eval JPrintDetector a $DETECTOR O SUMMARY set_variable STRING
Definition: examples/JTrigger/JTriggerEfficiencyRunByRun.sh:64
user
then wget no check certificate user
Definition: archive-put-wiki-detectors.sh:35
f
o $QUALITY_ROOT d $DEBUG!CHECK_EXIT_CODE JPlot1D f
Definition: JDataQuality.sh:76
exit
exit
Definition: JPizza.sh:36
Y
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
Definition: detector-XY:fit1d.sh:33
then
then
Definition: datalogs.sh:31
YMIN
fatal Invalid number of arguments esac let YMIN
Definition: testToT.sh:39
o
then fatal Wrong number of arguments fi JConvertDetectorFormat a o
Definition: software/JDetector/compare-detector.sh:31
rm
then rm
Definition: sftpput.zsh:30
$STRING
then fatal Invalid string $STRING
Definition: detector-XY:run.sh:34
awk
do set_variable STRING_TXT awk
Definition: plot-canberra.sh:195
a
then JCalibrateToT a
Definition: JTuneHV.sh:113
T
do set_variable OUTPUT_DIRECTORY $WORKDIR T
Definition: JCalibrateHeight.sh:61
set_variable
then echo Variable JPP_DIR undefined exit fi source $JPP_DIR setenv sh $JPP_DIR &dev null set_variable
Definition: archive-get-calibration.sh:14
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
FORMULA
&set_variable FORMULA
Definition: JEvtReweightMupageParameterScan.sh:60
p2
p2
Definition: module-Z:fit.sh:74
R
then JCookie sh JDataQuality D $DETECTOR_ID R
Definition: JDataQuality.sh:41
N
then usage $script< input file >[option[primary[working directory]]] nWhere option can be N
Definition: JMuonPostfit.sh:40
$WORKDIR
fi set_variable FILENAME $WORKDIR
Definition: datalogs.sh:34
d
then JMuonMCEvt f $INPUT_FILE o $INTERMEDIATE_FILE d
Definition: JMuonPath.sh:47
$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
X
no fit printf nominal n $STRING awk v X
Definition: detector-XY:fit.sh:102
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
detectorxy_
then fatal Wrong number of arguments fi set_variable STRING $argv[1] set_variable DETECTOR_TXT $WORKDIR detectorxy_
Definition: detector-XY:fit.sh:28
XMIN
then let XMIN
Definition: detector-XY:fit.sh:67
DETECTOR
set_variable DETECTOR
Definition: JARCACalibrateMuon.sh:31
JTOOLS::v
data_type v[N+1][M+1]
Definition: JPolint.hh:866
echo
then echo
Definition: JQAQC.sh:96
p3
p3
Definition: module-Z:fit.sh:74
$G1_TXT
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:35
null
esac $JPP_BIN JLogger sh $LOGGER until pgrep JGetMessage</dev/null > dev null
Definition: examples/JLogger/JLogger.sh:33
DEBUG
#define DEBUG(A)
Message macros.
Definition: JMessage.hh:62
g1
Double_t g1(const Double_t x)
Function.
Definition: JQuantiles.cc:25
Generated by   doxygen 1.8.5