Jpp test-rotations-old
the software that should make you happy
Loading...
Searching...
No Matches
JMEstimator.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 compare M-estimators.
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 WORKDIR ${TMPDIR:-/tmp}/
24set_variable DEBUG 2
25
26if do_usage $*; then
27 usage "$script"
28fi
29
30OPTIONS=(normal lorentzian linear tukey normalwithbackground)
31
32for mestimator in $OPTIONS[*]; do
33 $JPP_DIR/examples/JFit/JMEstimator \
34 -M ${mestimator} \
35 -o $WORKDIR/\[${mestimator}\].root
36done
37
38JPlot1D \
39 -f$WORKDIR/\[${^OPTIONS}\].root:rho \
40 -> "sigma" \
41 -\^ "#chi^{2}" \
42 -y "0 15" \
43 -O "][" \
44 -L TR
45
46JPlot1D \
47 -f$WORKDIR/\[${^OPTIONS}\].root:psi \
48 -> "sigma" \
49 -\^ "#partial#chi^{2}/#partialz" \
50 -y "-5 +5" \
51 -O "][" \
52 -L TR
53
54