def test_hist_1d(self): ftmp = NamedTemporaryFile(suffix='.h5', delete=False) try: ftmp.close() h = Histogram(3, [0, 3]) h.data[:] = [-3, 0, 5] h.save(ftmp.name) htmp = Histogram.load(ftmp.name) self.assertTrue(h.isidentical(htmp)) h.axes[0].label = 'x (cm)' h.save(ftmp.name) htmp = Histogram.load(ftmp.name) #print(type(h.axes[0].label), h.axes[0].label) #print(type(htmp.axes[0].label), htmp.axes[0].label) self.assertTrue(h.isidentical(htmp)) h.label = 'counts' h.save(ftmp.name) htmp = Histogram.load(ftmp.name) self.assertTrue(h.isidentical(htmp)) h.title = 'title' h.save(ftmp.name) htmp = Histogram.load(ftmp.name) self.assertTrue(h.isidentical(htmp)) h.uncert = [2, 3, 4] h.save(ftmp.name) htmp = Histogram.load(ftmp.name) self.assertTrue(h.isidentical(htmp)) finally: os.remove(ftmp.name)
def test_hist_4d(self): ftmp = NamedTemporaryFile(suffix='.root', delete=False) try: h = Histogram(3, [0, 1], 3, [0, 1], 3, [0, 1], 3, [0, 1]) with self.assertRaises(ValueError): h.save(ftmp.name) finally: os.remove(ftmp.name)
def test_hist_3d(self): h = Histogram(3, [0, 1], 'xx', 4, [-1, 1], 'yy', 5, [5, 10], 'zz', 'counts', 'data') ftmp = NamedTemporaryFile(suffix='.root', delete=False) try: ftmp.close() with warnings.catch_warnings(record=True) as w: h.save(ftmp.name) self.assertEqual(len(w), 1) self.assertRegex(str(w[-1].message), 'label') finally: os.remove(ftmp.name)
def test_root_2d(self): ftmp = NamedTemporaryFile(suffix='.root', delete=False) try: ftmp.close() h = Histogram(3, [0, 3], 4, [0, 4]) h.data[:] = [[-3, 0, 5, 3], [-2, 0, 4, 2], [-1, 0, 3, 1024]] h.uncert = np.sqrt(np.abs( h.data)) # ROOT will always return uncert h.save(ftmp.name) htmp = Histogram.load(ftmp.name) self.assertTrue(h.isidentical(htmp)) h.axes[0].label = 'x (cm)' h.axes[1].label = 'y (cm)' h.save(ftmp.name) htmp = Histogram.load(ftmp.name) self.assertTrue(h.isidentical(htmp)) h.label = 'counts' h.save(ftmp.name) htmp = Histogram.load(ftmp.name) self.assertTrue(h.isidentical(htmp)) h.title = 'title' h.save(ftmp.name) htmp = Histogram.load(ftmp.name) self.assertTrue(h.isidentical(htmp)) finally: os.remove(ftmp.name)
def test_npz_2d(self): ftmp = NamedTemporaryFile(suffix='.npz', delete=False) try: ftmp.close() h = Histogram(3, [0, 3], 4, [0, 4]) h.data[:] = [[-3, 0, 5, 3], [-2, 0, 4, 2], [-1, 0, 3, 1024]] h.save(ftmp.name) htmp = Histogram.load(ftmp.name) self.assertTrue(h.isidentical(htmp)) h.axes[0].label = 'x (cm)' h.axes[1].label = 'y (cm)' h.save(ftmp.name) htmp = Histogram.load(ftmp.name) self.assertTrue(h.isidentical(htmp)) h.label = 'counts' h.save(ftmp.name) htmp = Histogram.load(ftmp.name) self.assertTrue(h.isidentical(htmp)) h.title = 'title' h.save(ftmp.name) htmp = Histogram.load(ftmp.name) self.assertTrue(h.isidentical(htmp)) finally: os.remove(ftmp.name)
def test_unicode(self): ftmp = NamedTemporaryFile(suffix='.npz', delete=False) try: h = Histogram(3, [0, 3]) h.data[:] = [-3, 0, 5] h.title = 'χ-squared' h.label = 'αβγ' h.axes[0].label = 'θ' h.save(ftmp.name) htmp = Histogram.load(ftmp.name) self.assertTrue(h.isidentical(htmp)) finally: os.remove(ftmp.name)
def test_unicode(self): ftmp = NamedTemporaryFile(suffix='.root', delete=False) try: ftmp.close() h = Histogram(3, [0, 3]) h.data[:] = [-3, 0, 5] h.uncert = np.sqrt(np.abs( h.data)) # ROOT will always return uncert h.title = 'χ-squared' h.label = 'αβγ' h.axes[0].label = 'θ' h.save(ftmp.name) htmp = Histogram.load(ftmp.name) self.assertTrue(h.isidentical(htmp)) finally: os.remove(ftmp.name)