Jpp test-rotations-old
the software that should make you happy
Loading...
Searching...
No Matches
plot-height.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-height.txt
15set_variable GRAPH_ROOT graph-height.root
16set_variable T T=+20.4
17set_variable M M=\*
18
19if do_usage $*; then
20 usage "$script [working directory]"
21fi
22
23case $# in
24 1) set_variable WORKDIR $1;;
25 0) ;;
26 *) fatal "Wrong number of arguments."
27esac
28
29if ( ! reuse_file $GRAPH_TXT ); then
30
31 for OUTPUT_DIRECTORY in `eval echo $WORKDIR/${T}:${M}`; do
32
33 INPUT_FILES=($OUTPUT_DIRECTORY/*_${APP}\.root)
34
35 echo "Processing:"
36 echo "$INPUT_FILES[*]" | tr ' ' '\n'
37
38 U=`echo ${OUTPUT_DIRECTORY##.*/} | tr ':' '\n' | grep -vF $T | grep -vF $M | sed 's/.*=//'`
39 V=(`JQuality -f "$INPUT_FILES[*]" |& sed -n "s/Total\ quality\ *//p"`)
40
41 printf "%5.3f %12.1f %12.1f\n" $U $V >> $GRAPH_TXT
42
43 done
44fi
45
46if (( 1 )); then
47
48 JGraph \
49 -f $GRAPH_TXT \
50 -o $GRAPH_ROOT
51
52 JFit \
53 -f ${GRAPH_ROOT}:\.\* \
54 -o fit.root \
55 -F "[0] - [1]*(sqrt(1.0 + 0.5*(x-[2])*(x-[2])/([3]*[3])) - 1.0)" \
56 -@"p0 = GetMaximum" \
57 -@"p1 = GetMaximum * 0.01" \
58 -@"p2 = 1.0" \
59 -@"p3 = 1.0e-3" \
60 -d3
61
62 JPlot1D \
63 -f fit.root:\.\* \
64 -x "0.90 1.10" \
65 -> "height scaling" \
66 -\^ "Likelihood" \
67 -T ""
68
69fi