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