Jpp 19.3.0-rc.1
the software that should make you happy
Loading...
Searching...
No Matches
plot-position.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 1
12set_variable WORKDIR ./
13set_variable APP gandalf
14set_variable GRAPH_TXT graph-position.txt
15set_variable GRAPH_ROOT graph-position.root
16
17if do_usage $*; then
18 usage "$script [working directory]"
19fi
20
21case $# in
22 1) set_variable WORKDIR $1;;
23 0) ;;
24 *) fatal "Wrong number of arguments."
25esac
26
27if ( ! reuse_file $GRAPH_TXT ); then
28
29 set_variable T T=0
30 set_variable X X=\*
31 set_variable Y Y=\*
32 set_variable Z Z=0
33
34 for OUTPUT_DIRECTORY in `eval echo $WORKDIR/${T}:${X}:${Y}:${Z}`; do
35
36 INPUT_FILES=($OUTPUT_DIRECTORY/*_${APP}\.root)
37
38 echo "Processing:"
39 echo "$INPUT_FILES[*]" | tr ' ' '\n'
40
41 U=(`echo ${OUTPUT_DIRECTORY##.*/} | tr ':' '\n' | grep -vF $T | grep -vF $X | grep -vF $Y | grep -vF $Z | sed 's/.*=//'`)
42 V=(`JQuality -f "$INPUT_FILES[*]" |& sed -n "s/Total\ quality\ *//p"`)
43
44 printf "%5.1f %5.1f %12.1f %12.1f\n" $U $V >> $GRAPH_TXT
45
46 done
47fi
48
49if (( 1 )); then
50
51 JGraph2D \
52 -f $GRAPH_TXT \
53 -o $GRAPH_ROOT
54
55 JFit2D \
56 -f ${GRAPH_ROOT}:\.\* \
57 -o fit.root \
58 -F "[0] - [1]*(sqrt(1.0 + 0.5*(x-[2])*(x-[2])/([4]*[4]) + 0.5*(y-[3])*(y-[3])/([4]*[4])) - 1.0)" \
59 -@"p0 = GetMaximum" \
60 -@"p1 = GetMaximum * 0.01" \
61 -@"p2 = 0.0e0" \
62 -@"p3 = 0.0e0" \
63 -@"p4 = 30.0e0" \
64 -d 3
65
66 JPlot2D \
67 -f ${GRAPH_ROOT}:\.\* \
68 -f fit.root:user \
69 -> "#Deltax [m]" \
70 -< "#Deltay [m]" \
71 -\^ "Likelihood" \
72 -O SURF \
73 -T ""
74
75fi