Skip to content

Commit

Permalink
Fix lower/upper case again (temporarily)
Browse files Browse the repository at this point in the history
  • Loading branch information
erykoff committed Jun 17, 2024
1 parent 96216ac commit 9a98313
Showing 1 changed file with 15 additions and 4 deletions.
19 changes: 15 additions & 4 deletions redmapper/dereddener.py
Original file line number Diff line number Diff line change
Expand Up @@ -32,19 +32,30 @@ def deredden_array(self, galaxy_array):
# Lazy-load the map only when needed.
self._reddening_map = healsparse.HealSparseMap.read(self.mapfile)

if 'RA' in galaxy_array.dtype.names:
vals = self._reddening_map.get_values_pos(galaxy_array['RA'], galaxy_array['DEC'], lonlat=True)
if 'ra' in galaxy_array.dtype.names:
lower = True
else:
lower = False

if lower:
vals = self._reddening_map.get_values_pos(galaxy_array['ra'], galaxy_array['dec'], lonlat=True)
else:
vals = self._reddening_map.get_values_pos(galaxy_array['RA'], galaxy_array['DEC'], lonlat=True)

bad = (vals < -100.0)
if np.any(bad):
vals[bad] = np.median(vals[~bad])

galaxy_array['refmag'] -= self.norm*self.constants[self.ref_ind]*vals
if lower:
galaxy_array['refmag'] -= self.norm*self.constants[self.ref_ind]*vals
else:
galaxy_array['REFMAG'] -= self.norm*self.constants[self.ref_ind]*vals

for i in range(self.nmag):
galaxy_array['mag'][:, i] -= self.norm*self.constants[i]*vals
if lower:
galaxy_array['mag'][:, i] -= self.norm*self.constants[i]*vals
else:
galaxy_array['MAG'][:, i] -= self.norm*self.constants[i]*vals

def deredden_catalog(self, galaxies):
"""Deredden a galaxy catalog in-place.
Expand Down

0 comments on commit 9a98313

Please sign in to comment.