logger.info("Took %.2f seconds to read the base (%d files)" % (time.clock() - t1, base.sspfile.size))
wl_norm_window = (l_ssp < 5680.0) & (l_ssp > 5590.0)
################################################################################
##########
########## Morphology model setup
##########
################################################################################
logger.info("Creating original B-D model.")
norm_model = get_model(args.trueModel, with_default=True)
norm_params = np.array(norm_model.getParams(), dtype=norm_model.dtype)
logger.info("Original model at normalization window:\n%s\n" % str(norm_model))
logger.info('Creating PSF (FWHM = %.2f ", beta = %.2f)' % (args.psfFWHM, args.psfBeta))
PSF = moffat_psf(args.psfFWHM, args.psfBeta, size=max(args.Nx, args.Ny))
# Pad image to avoid artifacts at the borders when convolving.
Ny_psf = PSF.shape[0]
Nx_psf = PSF.shape[1]
imshape_pad = (args.Ny + 2 * Ny_psf, args.Nx + 2 * Nx_psf)
ifsshape_pad = (len(l_ssp),) + imshape_pad
# Fix the model center for padded images.
norm_y0 = norm_model.y0.value + Ny_psf
norm_x0 = norm_model.x0.value + Nx_psf
norm_model.y0.setValue(norm_y0)
norm_model.x0.setValue(norm_x0)
t1 = time.clock()
logger.info("Creating bulge spectra (tau proportional to distance).")
bulge_image_pad, disk_image_pad = create_model_images(norm_model, imshape_pad, PSF=None)
def setSynthPSF(self, FWHM=0.0, beta=None, size=15):
if beta is not None:
self.setPSF(moffat_psf(FWHM, beta, size=size))
else:
self.setPSF(gaussian_psf(FWHM, size=size))
