Jpp test-rotations-old
the software that should make you happy
Loading...
Searching...
No Matches
plot-orientation.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 2
12set_variable WORKDIR ./
13set_variable APP gandalf
14set_variable GRAPH_TXT graph-orientation-%.txt
15set_variable GRAPH_ROOT graph-orientation-%.root
16set_variable ID ID=2
17set_variable R R=\*
18
19if do_usage $*; then
20 usage "$script <string identifier> [working directory]"
21fi
22
23case $# in
24 2) set_variable WORKDIR $2;&
25 1) set_variable ID $1;;
26 *) fatal "Wrong number of arguments."
27esac
28
29set_variable GRAPH_TXT ${GRAPH_TXT/\%/${ID}}
30set_variable GRAPH_ROOT ${GRAPH_ROOT/\%/${ID}}
31
32if ( ! reuse_file $GRAPH_TXT ); then
33
34 for OUTPUT_DIRECTORY in `eval echo $WORKDIR/${ID}:${R}`; 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 $ID | grep -vF $R | sed 's/.*=//'`
42 V=(`JQuality -f "$INPUT_FILES[*]" |& sed -n "s/Total\ quality\ *//p"`)
43
44 printf "%5.1f %12.1f %12.1f\n" $U $V >> $GRAPH_TXT
45
46 done
47fi
48
49if (( 1 )); then
50
51 JGraph \
52 -f $GRAPH_TXT \
53 -o $GRAPH_ROOT
54
55 JFit \
56 -f ${GRAPH_ROOT}:\.\* \
57 -o fit.root \
58 -F "[0] + [1]*cos(x-[2])" \
59 -@"p0 = GetMaximum" \
60 -@"p1 = GetMaximum * 0.01" \
61 -@"p2 = 0.0e0" \
62 -d3
63
64 JPlot1D \
65 -f fit.root:\.\* \
66 -x "-3.14 +3.14" \
67 -> "#phi [rad]" \
68 -\^ "Likelihood" \
69 -T "String ${ID}"
70
71fi