def test_two_dependent_variables(self):
self.addTerm("x = ")
self.addTerm("y = 3 * x + 1", "y = Newtons")
r = Result('equations')
results = r.get()
self.assertEqual(results['x'], (None, None, '', 'x', None))
self.assertEqual(results['y'], (None, 'Newtons', '3 * x + 1', 'y', None))
self.assertEqual(results['foo'], (None, None, "size + time + bar", 'foo', None))
r.set('bar', 2)
r.set('x', 5)
results = r.get()
self.assertEqual(results['foo'], (82, None, "size + time + bar", 'foo', None))
self.assertEqual(results['x'], (5, None, '', 'x', None))
self.assertEqual(results['y'], (16, 'Newtons', '3 * x + 1', 'y', None))
r.__del__()
def test_galactic_model(self):
self.results = Result('equations')
self.results.set('frame', 'topocentric frame')
self.results.set('topocentric_freq', 1420)
self.results.set('declination', 38.43)
self.results.set('right_ascension', 10.25)
self.results.set('galactic', 'model')
value, _, _, _, _ = self.results.get('t_galactic_model')
self.assertTrue(value is not None)
self.assertAlmostEqual(value, 0.5586, 4)
self.results.__del__()
def test_galactic_model(self):
self.results = Result("equations")
self.results.set("frame", "topocentric frame")
self.results.set("topocentric_freq", 1420)
self.results.set("declination", 38.43)
self.results.set("right_ascension", 10.25)
self.results.set("galactic", "model")
value, _, _, _, _ = self.results.get("t_galactic_model")
self.assertTrue(value is not None)
self.assertAlmostEqual(value, 0.5586, 4)
self.results.__del__()
def setUp(self):
super(TestResultCalculations, self).setUp()
self.results = Result('equations')
self.results.set('frame', 'topocentric frame')
self.results.set('topocentric_freq', 1420)
class TestResultCalculations(unittest.TestCase):
def setUp(self):
super(TestResultCalculations, self).setUp()
self.results = Result('equations')
self.results.set('frame', 'topocentric frame')
self.results.set('topocentric_freq', 1420)
def tearDown(self):
super(TestResultCalculations, self).tearDown()
self.results.__del__()
def test_fwhm(self):
value, _, _, _, _ = self.results.get('fwhm')
self.assertAlmostEqual(8.677, value, 3)
def test_source_diameter(self):
self.results.set('source_diameter_slider', 2)
value, _, _, _, _ = self.results.get('source_diameter')
self.assertAlmostEqual(1.735, value, 3)
def test_topocentric_wavelength(self):
value, _, _, _, _ = self.results.get('topocentric_wavelength')
self.assertAlmostEqual(21.112, value, 3)
def test_aperture_eff(self):
self.results.set('declination', 38.43)
self.results.set('min_elevation', 5)
value, _, _, _, _ = self.results.get('aperture_eff')
self.assertAlmostEqual(0.6999, value, 3)
def test_extended_source_eff(self):
self.results.set('declination', 38.43)
self.results.set('min_elevation', 5)
self.results.set('source_diameter_slider', 2)
self.results.set('units', 'flux')
value, _, _, _, _ = self.results.get('extended_source_eff')
self.assertAlmostEqual(0.690, value, 3)
def test_confusion_limit(self):
self.results.set('backend', 'GBT Spectrometer')
self.results.set('declination', 38.43)
self.results.set('min_elevation', 5)
self.results.set('source_diameter_slider', 2)
self.results.set('units', 'flux')
value, _, _, _, _ = self.results.get('confusion_limit')
self.assertAlmostEqual(60.819, value, 3)
self.results.set('units', 'tr')
value, _, _, _, _ = self.results.get('confusion_limit')
self.assertAlmostEqual(173.663, value, 3)
self.results.set('units', 'ta')
value, _, _, _, _ = self.results.get('confusion_limit')
self.assertAlmostEqual(119.792, value, 3)
def test_airmass(self):
def getResults(dec):
self.results.set('declination', dec)
self.results.set('min_elevation', dec2Els(dec)[0])
min_el = float(self.results.get('min_elevation')[0])
max_el = float(self.results.get('max_elevation')[0])
air_mass = float(self.results.get('air_mass')[0])
return min_el, max_el, air_mass
min_el, max_el, air_mass = getResults(-20)
self.assertAlmostEqual(min_el, 5, 2)
self.assertAlmostEqual(max_el, 31.57, 2)
self.assertAlmostEqual(air_mass, 4.03, 2)
min_el, max_el, air_mass = getResults(20)
self.assertAlmostEqual(min_el, 5, 2)
self.assertAlmostEqual(max_el, 71.57, 2)
self.assertAlmostEqual(air_mass, 2.41, 2)
min_el, max_el, air_mass = getResults(38.43)
self.assertAlmostEqual(min_el, 5, 2)
self.assertAlmostEqual(max_el, 90, 2)
self.assertAlmostEqual(air_mass, 2.11, 2)
min_el, max_el, air_mass = getResults(50)
self.assertAlmostEqual(min_el, 5, 2)
self.assertAlmostEqual(max_el, 78.43, 2)
self.assertAlmostEqual(air_mass, 2.28, 2)
min_el, max_el, air_mass = getResults(60)
self.assertAlmostEqual(min_el, 8.43, 2)
self.assertAlmostEqual(max_el, 68.43, 2)
self.assertAlmostEqual(air_mass, 2.12, 2)
min_el, max_el, air_mass = getResults(85)
self.assertAlmostEqual(min_el, 33.43, 2)
self.assertAlmostEqual(max_el, 43.43, 2)
self.assertAlmostEqual(air_mass, 1.62, 2)
def test_est0(self):
self.results.set('backend', 'GBT Spectrometer')
self.results.set('declination', 38.43)
self.results.set('min_elevation', 5)
value = self.results.get('est0')[0]
self.assertTrue(value is not None)
est0 = float(value)
self.assertAlmostEqual(est0, 16.634, 3)
def test_attenuation(self):
self.results.set('backend', 'GBT Spectrometer')
self.results.set('declination', 38.43)
self.results.set('min_elevation', 5)
value = self.results.get('attenuation')[0]
self.assertTrue(value is not None)
attenuation = float(value)
self.assertAlmostEqual(attenuation, 1.015, 3)
def test_estTS(self):
self.results.set('backend', 'GBT Spectrometer')
self.results.set('declination', 38.43)
self.results.set('min_elevation', 5)
self.results.set('right_ascension', 0.0)
self.results.set('galactic', 'no_correction')
self.results.set('estimated_continuum', 3)
value = self.results.get('est_ts')[0]
self.assertTrue(value is not None)
est_ts = float(value)
self.assertAlmostEqual(est_ts, 19.852, 3)
def test_t_sys(self):
self.results.set('backend', 'GBT Spectrometer')
self.results.set('declination', 38.43)
self.results.set('min_elevation', 5)
self.results.set('right_ascension', 0.0)
self.results.set('galactic', 'no_correction')
self.results.set('estimated_continuum', 3)
value = self.results.get('t_sys')[0]
self.assertTrue(value is not None)
t_sys = float(value)
self.assertAlmostEqual(t_sys, 19.558, 3)
def setUp(self):
super(TestResultCalculations, self).setUp()
self.results = Result("equations")
self.results.set("frame", "topocentric frame")
self.results.set("topocentric_freq", 1420)
class TestResultCalculations(unittest.TestCase):
def setUp(self):
super(TestResultCalculations, self).setUp()
self.results = Result("equations")
self.results.set("frame", "topocentric frame")
self.results.set("topocentric_freq", 1420)
def tearDown(self):
super(TestResultCalculations, self).tearDown()
self.results.__del__()
def test_fwhm(self):
value, _, _, _, _ = self.results.get("fwhm")
self.assertAlmostEqual(8.677, value, 3)
def test_source_diameter(self):
self.results.set("units", "tr")
self.results.set("source_diameter_slider", 2)
value, _, _, _, _ = self.results.get("source_diameter")
self.assertAlmostEqual(1.735, value, 3)
def test_topocentric_wavelength(self):
value, _, _, _, _ = self.results.get("topocentric_wavelength")
self.assertAlmostEqual(21.112, value, 3)
def test_aperture_eff(self):
self.results.set("declination", 38.43)
self.results.set("min_elevation", 5)
value, _, _, _, _ = self.results.get("aperture_eff")
self.assertAlmostEqual(0.6999, value, 3)
def test_extended_source_eff(self):
self.results.set("declination", 38.43)
self.results.set("min_elevation", 5)
self.results.set("source_diameter_slider", 2)
self.results.set("units", "flux")
value, _, _, _, _ = self.results.get("extended_source_eff")
self.assertAlmostEqual(0.686957, value, 3)
def test_confusion_limit(self):
self.results.set("backend", "GBT Spectrometer")
self.results.set("declination", 38.43)
self.results.set("min_elevation", 5)
self.results.set("source_diameter_slider", 2)
self.results.set("units", "flux")
value, _, _, _, _ = self.results.get("confusion_limit")
self.assertAlmostEqual(60.819, value, 3)
self.results.set("units", "tr")
value, _, _, _, _ = self.results.get("confusion_limit")
self.assertAlmostEqual(173.663, value, 3)
self.results.set("units", "ta")
value, _, _, _, _ = self.results.get("confusion_limit")
self.assertAlmostEqual(119.792, value, 3)
def test_airmass(self):
def getResults(dec):
self.results.set("declination", dec)
self.results.set("min_elevation", dec2Els(dec)[0])
min_el = float(self.results.get("min_elevation")[0])
max_el = float(self.results.get("max_elevation")[0])
air_mass = float(self.results.get("air_mass")[0])
return min_el, max_el, air_mass
min_el, max_el, air_mass = getResults(-20)
self.assertAlmostEqual(min_el, 5, 2)
self.assertAlmostEqual(max_el, 31.57, 2)
self.assertAlmostEqual(air_mass, 4.03, 2)
min_el, max_el, air_mass = getResults(20)
self.assertAlmostEqual(min_el, 5, 2)
self.assertAlmostEqual(max_el, 71.57, 2)
self.assertAlmostEqual(air_mass, 2.41, 2)
min_el, max_el, air_mass = getResults(38.43)
self.assertAlmostEqual(min_el, 5, 2)
self.assertAlmostEqual(max_el, 90, 2)
self.assertAlmostEqual(air_mass, 2.11, 2)
min_el, max_el, air_mass = getResults(50)
self.assertAlmostEqual(min_el, 5, 2)
self.assertAlmostEqual(max_el, 78.43, 2)
self.assertAlmostEqual(air_mass, 2.28, 2)
min_el, max_el, air_mass = getResults(60)
self.assertAlmostEqual(min_el, 8.43, 2)
self.assertAlmostEqual(max_el, 68.43, 2)
self.assertAlmostEqual(air_mass, 2.12, 2)
min_el, max_el, air_mass = getResults(85)
self.assertAlmostEqual(min_el, 33.43, 2)
self.assertAlmostEqual(max_el, 43.43, 2)
self.assertAlmostEqual(air_mass, 1.62, 2)
def test_est0(self):
self.results.set("backend", "GBT Spectrometer")
self.results.set("declination", 38.43)
self.results.set("min_elevation", 5)
value = self.results.get("est0")[0]
self.assertTrue(value is not None)
est0 = float(value)
self.assertAlmostEqual(est0, 16.634, 3)
def test_attenuation(self):
self.results.set("backend", "GBT Spectrometer")
self.results.set("declination", 38.43)
self.results.set("min_elevation", 5)
value = self.results.get("attenuation")[0]
self.assertTrue(value is not None)
attenuation = float(value)
self.assertAlmostEqual(attenuation, 1.015, 3)
def test_estTS(self):
self.results.set("backend", "GBT Spectrometer")
self.results.set("declination", 38.43)
self.results.set("min_elevation", 5)
self.results.set("right_ascension", 0.0)
self.results.set("galactic", "no_correction")
self.results.set("estimated_continuum", 3)
value = self.results.get("est_ts")[0]
self.assertTrue(value is not None)
est_ts = float(value)
self.assertAlmostEqual(est_ts, 19.852, 3)
def test_t_sys(self):
self.results.set("backend", "GBT Spectrometer")
self.results.set("declination", 38.43)
self.results.set("min_elevation", 5)
self.results.set("right_ascension", 0.0)
self.results.set("galactic", "no_correction")
self.results.set("estimated_continuum", 3)
value = self.results.get("t_sys")[0]
self.assertTrue(value is not None)
t_sys = float(value)
self.assertAlmostEqual(t_sys, 19.558, 3)
def test_t_tot(self):
# general
data0 = {
"units": [u"tr"],
"conversion": [u"Sensitivity to Time"],
"sensitivity": [u"1"],
"semester": [u"A"],
"time": [u"1"],
}
# Hardware Information
data1 = {
"receiver": [u"L (1.15 - 1.73 GHz)"],
"beams": [u"1"],
"switching": [u"In-Band Frequency Switching"],
"bandwidth": [u"12.5"],
"polarization": [u"Dual"],
"backend": [u"GBT Spectrometer"],
"windows": [u"1"],
"switching": [u"In-Band Frequency Switching"],
"polarization": [u"Dual"],
"windows": [u"1"],
"bandwidth": [u"12.5"],
"receiver": [u"L (1.15 - 1.73 GHz)"],
"mode": [u"Spectral Line"], # not 'mode-hidden'
"backend": [u"GBT Spectrometer"],
}
# Source Information
data2 = {
"doppler": [u"Optical"],
"galactic": [u"no_correction"],
"redshift": [u"0"],
"frame": [u"Rest Frame"],
"source_diameter_slider": [u"1"], # to avoid t_tot None
"right_ascension": [u"0"], # needs to not be '' to eval
"rest_freq": [u"1440.0"],
"source_velocity": [u"0"],
"declination": [u"38"],
"estimated_continuum": [u"0"],
"topocentric_freq": [u"1440.0"],
"doppler": [u"Optical"]
# need a value to force execution
,
"effective_bw": [u"0"],
"duty_cycle": [u"0"],
"min_elevation": [u"5"],
}
# Data Reduction Information
data3 = {
"smoothing_resolution": [u"1"],
"diff_signal": [u"true"],
"smoothing": [u"velocity_resolution_rest"],
"r_sig_ref": [u"1"],
"bw": [u"0.004803322970701629"],
"smoothing_factor": [u"1"],
"avg_pol": [u"true"],
"no_avg_ref": [u"1"],
}
# evalute by setting all the terms above
for d in [data0, data1, data2, data3]:
for k, v in d.items():
try:
# send over those strs that are floats
v = float(v[0])
self.results.set(k, v)
except:
self.results.set(k, v[0])
# of course, these values all depend on equations.cfg;
# so expect these to fail when those equations change.
value = self.results.get("t_sys")[0]
self.assertAlmostEqual(16.61232745, value, 5)
value = self.results.get("t_tot")[0]
self.assertAlmostEqual(1667672275.52183, value, 2)
def setUp(self):
super(TestResult, self).setUp()
self.results = Result('equations')
class TestResult(CalculatorTestCase):
def setUp(self):
super(TestResult, self).setUp()
self.results = Result('equations')
def tearDown(self):
self.results.__del__()
super(TestResult, self).tearDown()
def testGet(self):
self.assertEquals(5, self.results.get('size')[0])
self.assertEquals('meters', self.results.get('size')[1])
self.assertEquals(15, self.results.get('sensitivity')[0])
self.assertEquals(75, self.results.get('time')[0])
self.assertEquals('sec', self.results.get('time')[1])
self.assertEquals(95, self.results.get('airmass')[0])
self.results.set('size', None)
self.assertEquals(5, self.results.get('size')[0])
self.assertEquals(None, self.results.get('sensitivity')[0])
self.assertEquals(None, self.results.get('time')[0])
self.assertEquals(None, self.results.get('airmass')[0])
self.results.set('size', 1)
self.assertEquals(1, self.results.get('size')[0])
self.assertEquals(3, self.results.get('sensitivity')[0])
self.assertEquals(3, self.results.get('time')[0])
self.assertEquals(7, self.results.get('airmass')[0])
def test_two_dependent_variables(self):
self.addTerm("x = ")
self.addTerm("y = 3 * x + 1", "y = Newtons")
r = Result('equations')
results = r.get()
self.assertEqual(results['x'], (None, None, '', 'x', None))
self.assertEqual(results['y'], (None, 'Newtons', '3 * x + 1', 'y', None))
self.assertEqual(results['foo'], (None, None, "size + time + bar", 'foo', None))
r.set('bar', 2)
r.set('x', 5)
results = r.get()
self.assertEqual(results['foo'], (82, None, "size + time + bar", 'foo', None))
self.assertEqual(results['x'], (5, None, '', 'x', None))
self.assertEqual(results['y'], (16, 'Newtons', '3 * x + 1', 'y', None))
r.__del__()