def testCosmologyCatalog(self):
"""
Does a catalog get written?
"""
dbObj = myTestGals(database=self.dbName)
cat = cosmologicalGalaxyCatalog(dbObj, obs_metadata=self.obs)
cat.write_catalog(self.catName)
def testCosmologyCatalog(self):
"""
Does a catalog get written?
"""
dbObj = myTestGals(database=self.dbName)
cat = cosmologicalGalaxyCatalog(dbObj, obs_metadata=self.obs)
with lsst.utils.tests.getTempFilePath('.txt') as catName:
cat.write_catalog(catName)
def testCosmologyCatalog(self):
"""
Does a catalog get written?
"""
dbObj = myTestGals(database=self.dbName)
cat = cosmologicalGalaxyCatalog(dbObj, obs_metadata=self.obs)
with lsst.utils.tests.getTempFilePath('.txt') as catName:
cat.write_catalog(catName)
def setUp(self):
self.obs_metadata = ObservationMetaData(mjd=52000.7,
boundType = 'circle',
pointingRA=200.0, pointingDec=-30.0,
boundLength=1.0,
m5=[23.9, 25.0, 24.7, 24.0, 23.3, 22.1],
bandpassName=['u', 'g', 'r', 'i', 'z', 'y'])
self.galaxy = myTestGals(database=self.testDB)
self.star = myTestStars(database=self.testDB)
def setUp(self):
defaults = LSSTdefaults()
bandpassName = ['u', 'g', 'r', 'i', 'z', 'y']
self.obs_metadata = ObservationMetaData(mjd=52000.7,
bandpassName=bandpassName,
m5=[defaults.m5(mm) for mm in bandpassName],
boundType='circle',
pointingRA=200.0, pointingDec=-30.0,
boundLength=1.0)
self.galaxy = myTestGals(database=self.testDB)
self.star = myTestStars(database=self.testDB)
def setUp(self):
self.obs_metadata = ObservationMetaData(
mjd=52000.7,
boundType='circle',
pointingRA=200.0,
pointingDec=-30.0,
boundLength=1.0,
m5=[23.9, 25.0, 24.7, 24.0, 23.3, 22.1],
bandpassName=['u', 'g', 'r', 'i', 'z', 'y'])
self.galaxy = myTestGals(database=self.testDB)
self.star = myTestStars(database=self.testDB)
def setUp(self):
defaults = LSSTdefaults()
bandpassName = ['u', 'g', 'r', 'i', 'z', 'y']
self.obs_metadata = ObservationMetaData(
mjd=52000.7,
bandpassName=bandpassName,
m5=[defaults.m5(mm) for mm in bandpassName],
boundType='circle',
pointingRA=200.0,
pointingDec=-30.0,
boundLength=1.0)
self.galaxy = myTestGals(database=self.testDB)
self.star = myTestStars(database=self.testDB)
def testCatalogDistanceModulus(self):
"""
Does cosmologicalDistanceModulus get properly applied
"""
dbObj = myTestGals(database=self.dbName)
cosmoCat = cosmologicalGalaxyCatalog(dbObj, obs_metadata=self.obs)
controlCat = absoluteGalaxyCatalog(dbObj, obs_metadata=self.obs)
cosmoIter = cosmoCat.iter_catalog(chunk_size=self.dbSize)
controlIter = controlCat.iter_catalog(chunk_size=self.dbSize)
cosmology = CosmologyObject()
for (cosmoRow, controlRow) in zip(cosmoIter, controlIter):
modulus = cosmology.distanceModulus(controlRow[25])
self.assertEqual(cosmoRow[0], controlRow[0])
self.assertEqual(cosmoRow[25], controlRow[25])
self.assertEqual(cosmoRow[26], modulus)
for i in range(1, 25):
self.assertAlmostEqual(cosmoRow[i], controlRow[i] + modulus, 6)
def testCatalogDistanceModulus(self):
"""
Does cosmologicalDistanceModulus get properly applied
"""
dbObj = myTestGals(database=self.dbName)
cosmoCat = cosmologicalGalaxyCatalog(dbObj, obs_metadata=self.obs)
controlCat = absoluteGalaxyCatalog(dbObj, obs_metadata=self.obs)
cosmoIter = cosmoCat.iter_catalog(chunk_size=self.dbSize)
controlIter = controlCat.iter_catalog(chunk_size=self.dbSize)
cosmology = CosmologyObject()
for (cosmoRow, controlRow) in zip(cosmoIter, controlIter):
modulus = cosmology.distanceModulus(controlRow[25])
self.assertEqual(cosmoRow[0], controlRow[0])
self.assertEqual(cosmoRow[25], controlRow[25])
self.assertEqual(cosmoRow[26], modulus)
for i in range(1, 25):
self.assertAlmostEqual(cosmoRow[i], controlRow[i] + modulus, 6)
def setUp(self):
self.starDB = myTestStars(database=self.starDbName)
self.galaxyDB = myTestGals(database=self.galaxyDbName)
def setUp(self):
self.starDB = myTestStars(database=self.starDbName)
self.galaxyDB = myTestGals(database=self.galaxyDbName)
