Jpp  16.0.2
the software that should make you happy
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
JGizmo/JHistogram2D.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 test JHistogram2D, JFit2D and JPlot2D.
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 
24 set_variable DEBUG 2
25 set_variable WORKDIR ${TMPDIR:-/tmp}/
26 
27 if do_usage $*; then
28  usage "$script [working directory]"
29 fi
30 
31 case $# in
32  1) set_variable WORKDIR $1;;
33 esac
34 
35 
36 if (( 1 )); then
37 
38  set_variable FORMULA "[0] * exp(-0.5 * (x-[1])*(x-[1]) / ([2]*[2])) * exp(-0.5 * (y-[1])*(y-[1]) / ([2]*[2]))"
39 
40  JHistogram2D \
41  -o $WORKDIR/histogram.root \
42  -F "$FORMULA" \
43  -@ "p0 = 1.0; p1 = 0.0; p2 = 2.0" \
44  -x "100 -10 +10" \
45  -y "100 -10 +10" \
46  -n 100000
47 
48 
49  JFit2D \
50  -f $WORKDIR/histogram.root:h0 \
51  -o $WORKDIR/fit.root \
52  -F "$FORMULA" \
53  -@ "p0 = GetMaximum" \
54  -@ "p1 = GetMean" \
55  -@ "p2 = GetRMS" \
56  -d $DEBUG
57 
58  echo -n "chi2/NDF: "
59  JPrintFit \
60  -f $WORKDIR/fit.root:h0 \
61  -F "GetChisquare/GetNDF"
62 
63  JPlot2D \
64  -f $WORKDIR/fit.root:h0 \
65  -O COLZ
66 
67 fi
68 
69 
70 if (( 1 )); then
71 
72  set_variable FORMULA "[0] * pow(10,[1]*x)"
73 
74  $JPP_DIR/examples/JGizmo/JHistogram2D \
75  -o $WORKDIR/histogram.root \
76  -F "$FORMULA" \
77  -@ "p0 = 1.0; p1 = -2" \
78  -x "100 0 7" \
79  -y "100 -1 +1" \
80  -n 10000000
81 
82 
83  JFit2D \
84  -f $WORKDIR/histogram.root:h0 \
85  -o $WORKDIR/fit.root \
86  -F "$FORMULA" \
87  -@ "p0 = GetMaximum" \
88  -@ "p1 = -2" \
89  -d3
90 
91  echo -n "chi2/NDF: "
92  JPrintFit \
93  -f $WORKDIR/fit.root:h0 \
94  -F "GetChisquare/GetNDF"
95 
96  JPlot2D \
97  -f $WORKDIR/fit.root:h0 \
98  -z "1e-1 1e5" \
99  -XXZ \
100  -O COLZ
101 
102 fi
103 
104 
105 if (( 1 )); then
106 
107  set_variable FORMULA "[0] * pow(10,[1]*y)"
108 
109  $JPP_DIR/examples/JGizmo/JHistogram2D \
110  -o $WORKDIR/histogram.root \
111  -F "$FORMULA" \
112  -@ "p0 = 1.0; p1 = -2" \
113  -x "100 -1 +1" \
114  -y "100 0 7" \
115  -n 10000000
116 
117 
118  JFit2D \
119  -f $WORKDIR/histogram.root:h0 \
120  -o $WORKDIR/fit.root \
121  -F "$FORMULA" \
122  -@ "p0 = GetMaximum" \
123  -@ "p1 = -2" \
124  -d3
125 
126  echo -n "chi2/NDF: "
127  JPrintFit \
128  -f $WORKDIR/fit.root:h0 \
129  -F "GetChisquare/GetNDF"
130 
131  JPlot2D \
132  -f $WORKDIR/fit.root:h0 \
133  -z "1e-1 1e5" \
134  -YYZ \
135  -O COLZ
136 
137 fi
f
o $QUALITY_ROOT d $DEBUG!JPlot1D f
Definition: JDataQuality.sh:66
exit
exit
Definition: JPizza.sh:36
echo
then echo
Definition: JDOMDAQDriver.sh:56
JTOOLS::n
const int n
Definition: JPolint.hh:676
set_variable
then echo Variable JPP_DIR undefined exit fi source $JPP_DIR setenv sh $JPP_DIR &dev null set_variable
Definition: JAcoustics.sh:20
TMPDIR
then echo Variable JPP_DIR undefined exit fi source $JPP_DIR setenv sh $JPP_DIR if do_usage *then usage $script[(input file)+] fi set_variable DEBUG set_variable WORKDIR TMPDIR
Definition: JEnergyCorrection.sh:21
WORKDIR
then echo Variable JPP_DIR undefined exit fi source $JPP_DIR setenv sh $JPP_DIR &dev null set_variable DEBUG set_variable WORKDIR
Definition: JLegolas.sh:20
F
then awk F
Definition: getMIIseacurrent.sh:86
usage
* usage
Definition: JLigier-local.sh:34
o
then JConvertDetectorFormat a $DETECTOR[1] o
Definition: JCompareDetector.sh:45
FORMULA
&set_variable FORMULA
Definition: JEvtReweightMupageParameterScan.sh:54
d
then JMuonMCEvt f $INPUT_FILE o $INTERMEDIATE_FILE d
Definition: JMuonPath.sh:47
d3
then echo Launching message logger in new window xterm geometry e JGetMessage H $HOST T $TAG d3
Definition: software/JLogger/JLogger.sh:39
DEBUG
then usage $script< input_file >< detector_file > fi set_variable OUTPUT_DIR set_variable SELECTOR JDAQTimesliceL1 set_variable DEBUG case set_variable DEBUG
Definition: JCalibratePMT.sh:21
$WORKDIR
then display $WORKDIR
Definition: plot-Domino.sh:128
in
then fatal Wrong number of arguments fi set_variable DETECTOR $argv[1] set_variable INPUT_FILE $argv[2] eval JPrintDetector a $DETECTOR O IDENTIFIER eval JPrintDetector a $DETECTOR O SUMMARY source JAcoustics sh $DETECTOR_ID CHECK_EXIT_CODE typeset A TRIPODS get_tripods $WORKDIR tripod txt TRIPODS for EMITTER in
Definition: JCanberra.sh:42
version
version
Definition: JCalibratePMT.sh:7
Generated by   doxygen 1.8.5