def load_SNsed(self): """ returns a list of SN seds in `lsst.sims.photUtils.Sed` observed within the spatio-temporal range specified by obs_metadata """ c, x1, x0, t0, _z, ra, dec = self.column_by_name('c'),\ self.column_by_name('x1'),\ self.column_by_name('x0'),\ self.column_by_name('t0'),\ self.column_by_name('redshift'),\ self.column_by_name('raJ2000'),\ self.column_by_name('decJ2000') SNobject = SNObject() raDeg = np.degrees(ra) decDeg = np.degrees(dec) sedlist = [] for i in range(self.numobjs): SNobject.set(z=_z[i], c=c[i], x1=x1[i], t0=t0[i], x0=x0[i]) SNobject.setCoords(ra=raDeg[i], dec=decDeg[i]) SNobject.mwEBVfromMaps() sed = SNobject.SNObjectSED(time=self.mjdobs, bandpass=self.lsstBandpassDict, applyExitinction=True) sedlist.append(sed) return sedlist
def load_SNsed(self): """ returns a list of SN seds in `lsst.sims.photUtils.Sed` observed within the spatio-temporal range specified by obs_metadata """ c, x1, x0, t0, _z, ra, dec = self.column_by_name('c'),\ self.column_by_name('x1'),\ self.column_by_name('x0'),\ self.column_by_name('t0'),\ self.column_by_name('redshift'),\ self.column_by_name('raJ2000'),\ self.column_by_name('decJ2000') SNobject = SNObject() raDeg = np.degrees(ra) decDeg = np.degrees(dec) sedlist = [] for i in range(self.numobjs): SNobject.set(z=_z[i], c=c[i], x1=x1[i], t0=t0[i], x0=x0[i]) SNobject.setCoords(ra=raDeg[i], dec=decDeg[i]) SNobject.mwEBVfromMaps() sed = SNobject.SNObjectSED(time=self.mjdobs, bandpass=self.lsstBandpassDict, applyExitinction=True) sedlist.append(sed) return sedlist
def get_snfluxes(self): c, x1, x0, t0, _z, ra, dec = self.column_by_name('c'),\ self.column_by_name('x1'),\ self.column_by_name('x0'),\ self.column_by_name('t0'),\ self.column_by_name('redshift'),\ self.column_by_name('raJ2000'),\ self.column_by_name('decJ2000') raDeg = np.degrees(ra) decDeg = np.degrees(dec) snobject = SNObject() # Initialize return array vals = np.zeros(shape=(self.numobjs, 19)) for i, _ in enumerate(vals): snobject.set(z=_z[i], c=c[i], x1=x1[i], t0=t0[i], x0=x0[i]) snobject.setCoords(ra=raDeg[i], dec=decDeg[i]) snobject.mwEBVfromMaps() # Calculate fluxes vals[i, :6] = snobject.catsimManyBandFluxes( time=self.mjdobs, bandpassDict=self.lsstBandpassDict, observedBandPassInd=None) # Calculate magnitudes vals[i, 6:12] = snobject.catsimManyBandMags( time=self.mjdobs, bandpassDict=self.lsstBandpassDict, observedBandPassInd=None) vals[i, 12:18] = snobject.catsimManyBandADUs( time=self.mjdobs, bandpassDict=self.lsstBandpassDict, photParams=self.photometricparameters) vals[i, 18] = snobject.ebvofMW return (vals[:, 0], vals[:, 1], vals[:, 2], vals[:, 3], vals[:, 4], vals[:, 5], vals[:, 6], vals[:, 7], vals[:, 8], vals[:, 9], vals[:, 10], vals[:, 11], vals[:, 12], vals[:, 13], vals[:, 14], vals[:, 15], vals[:, 16], vals[:, 17], vals[:, 18])
def get_snfluxes(self): c, x1, x0, t0, _z, ra, dec = self.column_by_name('c'),\ self.column_by_name('x1'),\ self.column_by_name('x0'),\ self.column_by_name('t0'),\ self.column_by_name('redshift'),\ self.column_by_name('raJ2000'),\ self.column_by_name('decJ2000') raDeg = np.degrees(ra) decDeg = np.degrees(dec) snobject = SNObject() # Initialize return array vals = np.zeros(shape=(self.numobjs, 19)) for i, _ in enumerate(vals): snobject.set(z=_z[i], c=c[i], x1=x1[i], t0=t0[i], x0=x0[i]) snobject.setCoords(ra=raDeg[i], dec=decDeg[i]) snobject.mwEBVfromMaps() # Calculate fluxes vals[i, :6] = snobject.catsimManyBandFluxes(time=self.mjdobs, bandpassDict=self.lsstBandpassDict, observedBandPassInd=None) # Calculate magnitudes vals[i, 6:12] = snobject.catsimManyBandMags(time=self.mjdobs, bandpassDict=self.lsstBandpassDict, observedBandPassInd=None) vals[i, 12:18] = snobject.catsimManyBandADUs(time=self.mjdobs, bandpassDict=self.lsstBandpassDict, photParams=self.photometricparameters) vals[i, 18] = snobject.ebvofMW return (vals[:, 0], vals[:, 1], vals[:, 2], vals[:, 3], vals[:, 4], vals[:, 5], vals[:, 6], vals[:, 7], vals[:, 8], vals[:, 9], vals[:, 10], vals[:, 11], vals[:, 12], vals[:, 13], vals[:, 14], vals[:, 15], vals[:, 16], vals[:, 17], vals[:, 18])