Jpp  master_rocky
the software that should make you happy
plot-height.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 
11 set_variable DEBUG 2
12 set_variable WORKDIR ./
13 set_variable APP gandalf
14 set_variable GRAPH_TXT graph-height.txt
15 set_variable GRAPH_ROOT graph-height.root
16 set_variable T T=+20.4
17 set_variable M M=\*
18 
19 if do_usage $*; then
20  usage "$script [working directory]"
21 fi
22 
23 case $# in
24  1) set_variable WORKDIR $1;;
25  0) ;;
26  *) fatal "Wrong number of arguments."
27 esac
28 
29 if ( ! 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
44 fi
45 
46 if (( 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 
69 fi