Jpp
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
detector-Zmul: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] * log(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]))"
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/detectorzmul_${STRING}.root
28 
29 if [[ ! -f $DETECTOR_TXT ]] then
30  fatal "No detector data file $DETECTOR_TXT."
31 fi
32 
33 JGraph2D \
34  -f $DETECTOR_TXT \
35  -o $DETECTOR_ROOT
36 
37 sort -gr -k 3,3 $DETECTOR_TXT | tail -1 | read X Y Z RMS
38 
39 let XMIN="$X - 1.0"
40 let XMAX="$X + 1.0"
41 let YMIN="$Y - 0.005"
42 let YMAX="$Y + 0.005"
43 
44 JFit2D \
45  -f ${DETECTOR_ROOT}:\.\* \
46  -o $WORKDIR/fit.root \
47  -F "$FORMULA" \
48  -x "$XMIN $XMAX" \
49  -y "$YMIN $YMAX" \
50  -@ "p0 = $Z" \
51  -@ "p1 = 1.0" \
52  -@ "p2 = $X" \
53  -@ "p3 = $Y" \
54  -@ "p4 = 1.0" \
55  -@ "p5 = 1.0" \
56  -@ "p6 = 0.001" \
57  -O LL \
58  -d $DEBUG
59 
60 set_array P2 `JPrintFit -f $WORKDIR/fit.root:user -@p2`
61 set_array P3 `JPrintFit -f $WORKDIR/fit.root:user -@p3`
62 
63 printf "nominal %4d %7.3f %7.4f %7.3f\n" $STRING 0.0 0.0 `JPrintFit -f $WORKDIR/fit.root:user -x " 0.0 0.0 "`
64 printf "optimal %4d %7.3f %7.4f %7.3f\n" $STRING $P2[1] $P3[1] `JPrintFit -f $WORKDIR/fit.root:user -x "$P2[1] $P3[1]"`
65 
66 if (( $# == 2 )); then
67 
69  set_variable MODULE `getModule -a $DETECTOR -L "$STRING 0"`
70 
71  JEditDetector -a $DETECTOR -S "$STRING add 0.0 0.0 $P2[1]" -S "$STRING mul $P3[1]" -M "$MODULE sub 0.0 0.0 $P2[1]" -o $DETECTOR
72 fi
73 
74 mv $WORKDIR/fit.root $DETECTOR_ROOT
75 
76 typeset -Z 4 STRING
77 
78 JPlot2D \
79  -f ${DETECTOR_ROOT}:\.\* \
80  -\> "#Deltaz [m]" \
81  -\< "stretching" \
82  -\^ "#chi^{2}/NDF" \
83  -N "X 505" \
84  -N "Y 505" \
85  -O CONT3 \
86  -T "[$STRING]" \
87  -o $WORKDIR/detectorzmul_${STRING}.gif $BACKGROUND
88 
89 rm -f $DETECTOR_ROOT
90 
91 
set_variable DETECTOR
do $JPP JMEstimator M
Definition: JMEstimator.sh:37
void set_variable(const std::string &name, const std::string &value)
Set environment variable.
clean eval JPrintDetector a $DETECTOR O IDENTIFIER eval JPrintDetector a $DETECTOR O SUMMARY set_variable STRING
then fatal Wrong number of arguments fi set_variable STRING $argv[1] set_variable DETECTOR_TXT $WORKDIR detectorzmul_
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
do set_array DAQHEADER JPrintDAQHeader f
Definition: JTuneHV.sh:74
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 echo
then fatal Invalid string $STRING
do set_variable OUTPUT_DIRECTORY $WORKDIR T
do $DIR JTransitTime o
Definition: JTransitTime.sh:44
then awk F
* usage
then fatal Wrong number of arguments fi set_variable STRING $argv[1] set_variable DETECTOR_TXT $WORKDIR $DETECTOR_TXT tail read X Y Z RMS let XMIN
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:48
*set_variable DETECTOR $argv[1] set_variable INPUT_FILE $argv[2,-2] set_variable OUTPUT_FILE $argv[-1]
Definition: JDomino.sh:39
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:103
then display $WORKDIR
Definition: plot-Domino.sh:127
then usage $script[input file[working directory[option]]] nWhere option can be N
Definition: JMuonPostfit.sh:37
script
Definition: JAcoustics.sh:2
p3
Definition: module-Z:fit.sh:48
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
const JModule & getModule(const JDetector &detector, const JModuleLocation &location)
find module with a given string and floor number
#define DEBUG(A)
Message macros.
Definition: JMessage.hh:62