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))