Jpp test-rotations-new
the software that should make you happy
Loading...
Searching...
No Matches
JPostfit2F.sh
Go to the documentation of this file.
1#!/usr/bin/env zsh
2#
3# \author mdejong
4#
5version=1.0
6script=${0##*/}
7
8# ------------------------------------------------------------------------------------------
9#
10# Utility script to optimise sorting of fits.
11#
12# ------------------------------------------------------------------------------------------
13
14
15if [ -z $JPP_DIR ]; then
16 echo "Variable JPP_DIR undefined."
17 exit
18fi
19
20
21source $JPP_DIR/setenv.sh $JPP_DIR
22
23set_variable DEBUG 2
24set_variable WORKDIR ${TMPDIR:-/tmp}/
25set_variable ANGLE_DEG 0.0
26set_variable JPP $JPP_DIR/examples/JReconstruction/
27set_variable: FORMAT GRAPHICS_FORMAT gif
28set_variable+ BATCH GRAPHICS_BATCH -B
29
30if do_usage $*; then
31 usage "$script <input file> <input file> [<angle [deg]>]"
32fi
33
34case $# in
35 3) set_variable ANGLE_DEG $3;&
36 2) set_variable INPUT_FILE_B $2;
37 set_variable INPUT_FILE_A $1;;
38 *) fatal "Wrong number of arguments."
39esac
40
41$JPP/JPostfit2F \
42 -a $INPUT_FILE_A \
43 -b $INPUT_FILE_B \
44 -o $WORKDIR/postfit.root \
45 -A $ANGLE_DEG \
46 -d ${DEBUG} --!
47
48JPlot1D \
49 -f $WORKDIR/postfit.root:h.A \
50 -f $WORKDIR/postfit.root:h.B \
51 -> "angle [deg]" \
52 -\^ "number of events [a.u.]" \
53 -XX \
54 -L TR -T "" \
55 -o angle1D.$FORMAT $BATCH
56
57JPlot2D \
58 -f $WORKDIR/postfit.root:h2 \
59 -> "quality" \
60 -< "angle [deg]" \
61 -O COLZ \
62 -T "" \
63 -o angle2D.$FORMAT $BATCH
64
65for PARAMETER in \
66 JQUALITY \
67 JGANDALF_BETA0_RAD \
68 JGANDALF_BETA1_RAD \
69 JGANDALF_NUMBER_OF_HITS \
70 JSTART_NPE_MIP \
71 JSTART_NPE_MIP_TOTAL \
72 JSTART_LENGTH_METRES \
73 JENERGY_ENERGY \
74 JENERGY_CHI2; \
75 do
76
77 JPlot1D \
78 -f "$WORKDIR/postfit.root:^\[A\].${PARAMETER}$" \
79 -f "$WORKDIR/postfit.root:^\[B\].${PARAMETER}$" \
80 -\^ "number of events [a.u.]" \
81 -N "X 505" \
82 -L TR \
83 -T "$PARAMETER" \
84 -o AB.${PARAMETER:l}.$FORMAT $BATCH
85
86 JPlot1D \
87 -f "$WORKDIR/postfit.root:^\[C\].${PARAMETER}$" \
88 -f "$WORKDIR/postfit.root:^\[D\].${PARAMETER}$" \
89 -> "#Delta" \
90 -\^ "number of events [a.u.]" \
91 -N "X 505" \
92 -L TR \
93 -T "$PARAMETER" \
94 -o CD.${PARAMETER:l}.$FORMAT $BATCH
95done