test_csviscube.py
4.04 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
130
131
132
133
134
135
136
137
138
139
140
#! /usr/bin/env python
# ==========================================================================
# This scripts performs unit tests for the csviscube script
#
# Copyright (C) 2016-2017 Juergen Knoedlseder
#
# This program is free software: you can redistribute it and/or modify
# it under the terms of the GNU General Public License as published by
# the Free Software Foundation, either version 3 of the License, or
# (at your option) any later version.
#
# This program is distributed in the hope that it will be useful,
# but WITHOUT ANY WARRANTY; without even the implied warranty of
# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
# GNU General Public License for more details.
#
# You should have received a copy of the GNU General Public License
# along with this program. If not, see <http://www.gnu.org/licenses/>.
#
# ==========================================================================
import gammalib
import cscripts
from testing import test
# =============================== #
# Test class for csviscube script #
# =============================== #
class Test(test):
"""
Test class for csviscube script
"""
# Constructor
def __init__(self):
"""
Constructor
"""
# Call base class constructor
test.__init__(self)
# Return
return
# Set test functions
def set(self):
"""
Set all test functions
"""
# Set test name
self.name('csviscube')
# Append tests
self.append(self._test_cmd, 'Test csviscube on command line')
self.append(self._test_python, 'Test csviscube from Python')
# Return
return
# Test csviscube on command line
def _test_cmd(self):
"""
Test csviscube on the command line
"""
# Set script name
csviscube = self._script('csviscube')
# Setup csviscube command
cmd = csviscube+' tmin=0.0 tmax=10000.0 geolon=79.4041'+\
' geolat=-24.6272 binsz=10.0 outfile=csviscube_cmd1.fits'+\
' logfile="csviscube_cmd1.log" chatter=1'
# Check if execution of wrong command fails
self.test_assert(self._execute('command_that_does_not_exist') != 0,
'Self test of test script')
# Check if execution was successful
self.test_assert(self._execute(cmd) == 0,
'Check successful execution from command line')
# Check visibility cube
self._check_viscube('csviscube_cmd1.fits')
# Check csviscube --help
self._check_help(csviscube)
# Return
return
# Test csviscube from Python
def _test_python(self):
"""
Test csviscube from Python
"""
# Set-up csviscube
viscube = cscripts.csviscube()
viscube['tmin'] = 0.0
viscube['tmax'] = 10000.0
viscube['geolon'] = 79.4041
viscube['geolat'] = -24.6272
viscube['binsz'] = 10.0
viscube['outfile'] = 'csviscube_py1.fits'
viscube['logfile'] = 'csviscube_py1.log'
viscube['chatter'] = 2
viscube['publish'] = True
# Execute script
viscube.logFileOpen() # Make sure we get a log file
viscube.execute()
# Check visibility cube
self._check_viscube('csviscube_py1.fits')
# Return
return
# Check visibility cube
def _check_viscube(self, filename, nx=36, ny=18, nz=60):
"""
Check visibility cube
Parameters
----------
filename : str
Visibility cube file name
nx : int, optional
Number of X pixels
ny : int, optional
Number of Y pixels
"""
# Open result file
viscube = gammalib.GSkyMap(filename)
# Check dimensions
self.test_value(viscube.nmaps(), nz, 'Check number of maps')
self.test_value(viscube.nx(), nx, 'Check for number of X pixels')
self.test_value(viscube.ny(), ny, 'Check for number of Y pixels')
# Return
return