test_ctatools.sh
2.25 KB
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
#! /bin/bash
#
# ctatools test script
# ====================================================================
#
# Print Header
#
echo "*****************************"
echo "* Test ctatools executables *"
echo "*****************************"
#
# Remove any existing result files
# ================================
rm -rf *.fits *.log *.xml
#
# Test ctobssim
# ============
echo -n "Test ctobssim: "
ctobssim infile="data/crab.xml" \
outfile="events.fits" \
caldb="irf" \
irf="kb_E_50h_v3" \
ra=83.63 \
dec=22.01 \
rad=10.0 \
tmin=0.0 \
tmax=1800.0 \
emin=0.1 \
emax=100.0
echo -n "."
if [ -s "events.fits" ]
then
echo -n "."
else
echo " events.fits file is not valid"
exit 1
fi
echo " ok"
#
# Test ctbin
# ==========
echo -n "Test ctbin: "
ctbin evfile="events.fits" \
outfile="cntmap.fits" \
emin=0.1 \
emax=100.0 \
enumbins=20 \
nxpix=200 \
nypix=200 \
binsz=0.02 \
coordsys="CEL" \
xref=83.63 \
yref=22.01 \
proj="CAR"
echo -n "."
if [ -s "cntmap.fits" ]
then
echo -n "."
else
echo " cntmap.fits file is not valid"
exit 1
fi
echo " ok"
#
# Test ctselect
# =============
echo -n "Test ctselect: "
ctselect infile="events.fits" \
outfile="selected_events.fits" \
ra=83.63 \
dec=22.01 \
rad=10.0 \
tmin=0.0 \
tmax=1800.0 \
emin=0.1 \
emax=100.0
echo -n "."
if [ -s "selected_events.fits" ]
then
echo -n "."
else
echo " selected_events.fits file is not valid"
exit 1
fi
echo " ok"
#
# Test ctlike
# ===========
echo -n "Test ctlike: "
ctlike cntmap="cntmap.fits" \
srcmdl="data/crab.xml" \
outmdl="results_binned.xml" \
method="BINNED" \
caldb="irf" \
irf="kb_E_50h_v3"
echo -n "."
if [ -s "results_binned.xml" ]
then
echo -n "."
else
echo " results_binned.xml file is not valid"
exit 1
fi
ctlike evfile="selected_events.fits" \
srcmdl="data/crab.xml" \
outmdl="results_unbinned.xml" \
method="UNBINNED" \
caldb="irf" \
irf="kb_E_50h_v3"
echo -n "."
if [ -s "results_unbinned.xml" ]
then
echo -n "."
else
echo " results_unbinned.xml file is not valid"
exit 1
fi
echo " ok"