def _set_event_PSPL(self): """ Sets internal variable to MulensModel.Event instance that uses PSPL model. """ datasets = [] for band in self.bands: mask = (self._filter == band) if not np.any(mask): continue times = self._JD[mask] flux = self._simulated_flux[band] sigma_flux = self._sigma_flux[band] plot = {} data = MM.MulensData([times, flux, sigma_flux], phot_fmt='flux', bandpass=band, plot_properties={'label': 'LSST ' + band}) datasets.append(data) if self._follow_up_Chilean is not None: if self._follow_up_Chilean.n_epochs > 0: datasets.append(self._follow_up_Chilean) if self._follow_up_nonChilean is not None: if self._follow_up_nonChilean.n_epochs > 0: datasets.append(self._follow_up_nonChilean) model = MM.Model( {p: self._parameters[p] for p in ['t_0', 'u_0', 't_E']}) self._event_PSPL = MM.Event(datasets, model)