Jpp
JPostfit2F.sh
Go to the documentation of this file.
1 #!/bin/zsh
2 #
3 # \author mdejong
4 #
5 version=1.0
6 script=${0##*/}
7 
8 # ------------------------------------------------------------------------------------------
9 #
10 # Utility script to optimise sorting of fits.
11 #
12 # ------------------------------------------------------------------------------------------
13 
14 
15 if [ -z $JPP_DIR ]; then
16  echo "Variable JPP_DIR undefined."
17  exit
18 fi
19 
20 
21 source $JPP_DIR/setenv.sh $JPP_DIR
22 
23 set_variable DEBUG 2
24 set_variable INPUT_FILE_A $JPP_DATA/gandalf.root
25 set_variable INPUT_FILE_B $JPP_DATA/gandalf-mc.root
26 set_variable OUTPUT_FILE postfit.root
27 set_variable ANGLE_DEG 2.0
28 set_variable JPP $JPP_DIR/examples/JFit/
29 
30 
31 if ( do_usage $* ); then
32  usage "$script [input file]"
33 fi
34 
35 case $# in
36  2) INPUT_FILE_B=$2;&
37  1) INPUT_FILE_A=$1;;
38 esac
39 
40 
41 if (( 1 )); then
42 
43  if ( ! reuse_file $OUTPUT_FILE); then
44 
45  $JPP/JPostfit2F \
46  -a $INPUT_FILE_A \
47  -b $INPUT_FILE_B \
48  -o $OUTPUT_FILE \
49  -A $ANGLE_DEG \
50  -d ${DEBUG} --!
51  fi
52 fi
53 
54 
55 if (( 1 )); then
56 
57  JPlot1D \
58  -f postfit.root:h.A \
59  -f postfit.root:h.B \
60  -> "angle [deg]" \
61  -\^ "number of events [a.u]" \
62  -XX \
63  -L TR -T "" \
64  -o angle.gif
65 
66 fi
67 
68 
69 if (( 1 )); then
70 
71  JPlot2D \
72  -f postfit.root:h2 \
73  -> "quality" \
74  -< "angle [deg]" \
75  -O COLZ \
76  -T "" \
77  -o angle2d.gif
78 
79 fi
80 
81 
82 for PARAMETER in \
83  JQUALITY \
84  JGANDALF_BETA0_RAD \
85  JGANDALF_BETA1_RAD \
86  JGANDALF_NUMBER_OF_HITS \
87  JSTART_NPE_MIP \
88  JSTART_NPE_MIP_TOTAL \
89  JSTART_LENGTH_METRES \
90  JENERGY_ENERGY \
91  JENERGY_CHI2; \
92  do
93 
94  if (( 1 )); then
95 
96  JPlot1D \
97  -f postfit.root:'\^.*C.'${PARAMETER}'.$' \
98  -\^ "number of events [a.u]" \
99  -L TR -T ""
100 
101  fi
102 
103 done
Generated by   doxygen 1.8.16