def mapSFR(galaxy, mom0cut, regridDir, outDir):
'''
import sfr map from W4+FUV
Date Programmer Description of Changes
----------------------------------------------------------------------
10/29/2020 Yiqing Song Original Code
12/3/2020 A.A. Kepley Added comments and clarified inputs
'''
sfrhdu = fits.open(
os.path.join(regridDir,
galaxy['NAME'] + '_sfr_fuvw4_gauss15_regrid.fits'))[0]
sfr = sfrhdu.data
sfr[np.isnan(mom0cut)] = np.nan
w = WCS(sfrhdu.header)
## TODO CHECK UNITS HERE
sfrhdu.header[
'BUNIT'] = 'MSUN/YR/KPC^2' #might change when getting new files from Sarah
sfrmap = Projection(sfr, header=sfrhdu.header, wcs=w)
sfrmap.quicklook()
# save plot of map
plt.savefig(os.path.join(outDir, galaxy['NAME'] + '_sfr.png'))
plt.clf()
plt.close()
return sfrmap
def mapLTIR(galaxy, mom0cut, regridDir, outDir):
'''
make LTIR map
galaxy: line from degas_base.fits table with galaxy information
mom0cut: S/N cut on mom0
regridDir: input data directory with all regridded data
outDir: output data directory
Date Programmer Description of Changes
----------------------------------------------------------------------
12/10/2020 A.A. Kepley Original code based on mapStellar
'''
hdu = fits.open(
os.path.join(regridDir,
galaxy['NAME'] + '_LTIR_gauss15_regrid.fits'))[0]
data = hdu.data
data[np.isnan(mom0cut)] = np.nan #apply SN mask (SN >3)
LTIRmap = Projection(data, header=hdu.header, wcs=WCS(hdu.header))
LTIRmap.quicklook()
plt.savefig(os.path.join(outDir, galaxy['NAME'] + '_LTIR.png'))
plt.clf()
plt.close()
return LTIRmap
def mapStellar(galaxy,mom0cut,plotdir='./',stellardatadir='./'):
stellarhdu=fits.open(stellardatadir+galaxy+'_w1_stellarmass_regrid.fits')[0]
starmask=fits.getdata(stellardatadir+galaxy+'_w1_gauss15_stars_regrid.fits')
stellar=stellarhdu.data
stellar[starmask==1.0]=np.nan
stellar[np.isnan(mom0cut)]=np.nan
w=WCS(stellarhdu.header)
stellarhdu.header['BUNIT']='Msun/pc^2' #not the correct unit!!
stellarmap=Projection(stellar,header=stellarhdu.header,wcs=w)
stellarmap.quicklook()
plt.savefig(plotdir+galaxy+'_stellarmass.png')
plt.clf()
plt.close()
return stellarmap
def mapSFR(galaxy,mom0cut,plotdir='./',sfrdatadir='./'):
sfrhdu=fits.open(sfrdatadir+galaxy+'_w4fuv_sfr_regrid.fits')[0]
starmask_fuv=fits.getdata(sfrdatadir+galaxy+'_fuv_gauss15_stars_regrid.fits')
starmask_nuv=fits.getdata(sfrdatadir+galaxy+'_nuv_gauss15_stars_regrid.fits')
sfr=sfrhdu.data
sfr[starmask_fuv==1.0]=np.nan #mask out stars from fuv and nuv
sfr[starmask_nuv==1.0]=np.nan
sfr[np.isnan(mom0cut)]=np.nan
w=WCS(sfrhdu.header)
sfrhdu.header['BUNIT']='MSUN/YR/KPC^2'
sfrmap=Projection(sfr,header=sfrhdu.header,wcs=w)
sfrmap.quicklook()
plt.savefig(plotdir+galaxy+'_sfr.png')
plt.clf()
plt.close()
return sfrmap
def mapStellar(galaxy, mom0cut, regridDir, outDir):
'''
make stellarmass map
galaxy: line from degas_base.fits table with galaxy information
mom0cut: S/N cut on mom0
regridDir: input data directory with all regridded data
outDir: output data directory
Date Programmer Description of Changes
----------------------------------------------------------------------
10/29/2020 Yiqing Song Original Code
12/3/2020 A.A. Kepley Added comments and clarified inputs
'''
# open the stellar mass
## TODO : NEED TO UPDATE WITH NEW ANCILLAR DATA FROM SARAH
stellarhdu = fits.open(
os.path.join(regridDir,
galaxy['NAME'] + '_mstar_gauss15_regrid.fits'))[0]
stellar = stellarhdu.data
#stellar[starmask==1.0]=np.nan #apply star mask ## AAK: I think I can skip this.
stellar[np.isnan(mom0cut)] = np.nan #apply SN mask (SN >3)
w = WCS(stellarhdu.header)
# stellarhdu.header['BUNIT']='Msun/pc^2' #not the correct unit!! ## AAK: I think units are okay for the newdata
stellarmap = Projection(stellar, header=stellarhdu.header, wcs=w)
stellarmap.quicklook()
plt.savefig(os.path.join(outDir, galaxy['NAME'] + '_stellarmass.png'))
plt.clf()
plt.close()
return stellarmap
