Commit 5c3a38fa authored Oct 09, 2026 by Wei Chengliang

Set galaxy size limits and Sersic truncation

parent ea4c8191

catalog/newCatalog.py

+3 −1

Original line number	Diff line number	Diff line
		@@ -390,7 +390,9 @@ class Catalog(CatalogBase):
		param["disksize"] = np.rad2deg(gals["diskradius"][igals] / angdist) * 3600.
		param["hlr_bulge"] = param["bulgesize"]
		param["hlr_disk"] = param["disksize"]
		if param["hlr_bulge"] < 0.001 and param["hlr_disk"] < 0.001:
		if param["hlr_bulge"] < 0.01 and param["hlr_disk"] < 0.01:
		continue
		if param["hlr_bulge"] > 5.0 or param["hlr_disk"] > 5.0:
		continue

		# SED coefficients

observation_sim/mock_objects/Galaxy.py

+4 −0

Original line number	Diff line number	Diff line
		@@ -171,6 +171,8 @@ class Galaxy(MockObject):
		half_light_radius=self.hlr_disk,
		flux=1.0,
		gsparams=gsp,
		trunc=self.hlr_disk*10,
		flux_untruncated=True,
		)
		disk_shape = galsim.Shear(g1=self.e1_disk, g2=self.e2_disk)
		disk = disk.shear(disk_shape)
		@@ -179,6 +181,8 @@ class Galaxy(MockObject):
		half_light_radius=self.hlr_bulge,
		flux=1.0,
		gsparams=gsp,
		trunc=self.hlr_bulge*50,
		flux_untruncated=True,
		)
		bulge_shape = galsim.Shear(g1=self.e1_bulge, g2=self.e2_bulge)
		bulge = bulge.shear(bulge_shape)