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