Commit 4f12b49b authored by Fang Yuedong's avatar Fang Yuedong
Browse files

scale wavelengths by 1+z for qso sed

parent 24fcd003

Catalog/C6_50sqdeg.py

+1 −1
Original line number Original line Diff line number Diff line
@@ -454,7 +454,7 @@ class Catalog(CatalogBase): 
            elif obj.type == 'quasar':

            elif obj.type == 'quasar':

                flux = self.agn_seds[obj.agnsed_file][int(obj.qsoindex)] * 1e-17

                flux = self.agn_seds[obj.agnsed_file][int(obj.qsoindex)] * 1e-17

                flux[flux < 0] = 0.

                flux[flux < 0] = 0.

                wave = self.lamb_gal

                wave = self.lamb_gal * (1.0 + obj.z)

        else:

        else:

            raise ValueError("Object type not known")

            raise ValueError("Object type not known")

        speci = interpolate.interp1d(wave, flux)

        speci = interpolate.interp1d(wave, flux)