Jpp test-rotations-old
the software that should make you happy
Loading...
Searching...
No Matches
sound-V:fit.sh
Go to the documentation of this file.
1#!/usr/bin/env zsh
2script=${0##*/}
3
4if [ -z $JPP_DIR ]; then
5 echo "Variable JPP_DIR undefined."
6 exit
7fi
8
9source $JPP_DIR/setenv.sh $JPP_DIR >& /dev/null
10
11set_variable: DEBUG ACOUSTICS_DEBUG 2
12set_variable: WORKDIR ACOUSTICS_WORKDIR ./
13set_variable FORMULA "[0] + [1]*(sqrt(1.0 + 0.5*(x-[2])*(x-[2])/([3]*[3])) - 1.0)"
14set_variable: FORMAT GRAPHICS_FORMAT gif
15set_variable+ BATCH GRAPHICS_BATCH -B
16
17if do_usage $*; then
18 usage "$script"
19fi
20
21if (( $# != 0 )); then
22 fatal "Wrong number of arguments."
23fi
24
25set_variable VELOCITY_TXT $WORKDIR/soundv.txt
26set_variable VELOCITY_ROOT $WORKDIR/soundv.root
27
28if [[ ! -f $VELOCITY_TXT ]]; then
29 fatal "No sound velocity file $VELOCITY_TXT."
30fi
31
32JGraph \
33 -f $VELOCITY_TXT \
34 -o $VELOCITY_ROOT
35
36sort -gr -k 2,2 $VELOCITY_TXT | tail -1 | read X Y RMS
37
38JFit \
39 -f ${VELOCITY_ROOT}:\.\* \
40 -o $WORKDIR/fit.root \
41 -F "$FORMULA" \
42 -@"p0 = $Y" \
43 -@"p1 = 1.0" \
44 -@"p2 = $X" \
45 -@"p3 = 1.0" \
46 -d $DEBUG
47
48set_array P2 `JPrintFit -f $WORKDIR/fit.root:\.\* -@p2`
49
50printf "result %10.5f +/- %10.5f\n" $P2[1] $P2[2]
51
52mv $WORKDIR/fit.root $VELOCITY_ROOT
53
54JPlot1D \
55 -f $VELOCITY_ROOT:\.\* \
56 -> "V [m/s]" \
57 -\^ "#chi^{2}/NDF" \
58 -T "" \
59 -o $WORKDIR/velocity.$FORMAT $BATCH
60
61rm -f $VELOCITY_ROOT