How do I use the MaNGA Catalog DAPALL File?

Astropy primer for working with table data

First open the file and select those entries with successful analysis from the DAP using the HYB10-GAU-MILESHC method.

import numpy as np
from astropy.io import fits

dapall = fits.open('dapall-v2_4_3-2.2.1.fits')
indx = (dapall['DAPALL'].data['DAPTYPE'] == 'HYB10-GAU-MILESHC') & \
       (dapall['DAPALL'].data['DAPDONE'] == 1)
tbdata = dapall[1].data[indx]

We can use this file to make simple plot of stellar vs. gas velocity dispersion within the effective radius, for example:

import matplotlib.pyplot as pyplot
pyplot.scatter(np.ma.log10(tbdata['ha_gsigma_1re']),
               np.ma.log10(tbdata['stellar_sigma_1re']),
               alpha=0.5, marker='.', s=30, lw=0)
pyplot.xlim(1,3.2)
pyplot.ylim(0.8,3)
pyplot.xlabel(r'ionized-gas $\sigma$ within 1 $R_e$ [km/s]')
pyplot.ylabel(r'stellar $\sigma$ within 1 $R_e$ [km/s]')
pyplot.show()

The DAPALL file can be read in using for example MRDFITS. This will put the data into an IDL structure.

t=mrdfits('dapall-v2_4_3-2.2.1.fits',1,hdr)
; Print column names
help,t,/str
; Select the successful HYB10-GAU-MILESHC analysis products
indx = where((t.dapdone eq 'T') and (t.daptype eq 'HYB10-GAU-MILESHC'))

We can use this file to make simple plot of stellar vs. gas velocity dispersion within the effective radius, for example:

plot,alog10(t.ha_gsigma_1re[indx]),alog10(t.stellar_sigma_1re[indx]),$
psym=3,xrange=[1,3.2], yrange=[0.8,3.0],$
xtitle='gas velocity dispersion',ytitle='stellar velocity dispersion',$
/ysty,/xsty

Accessing DAPALL file using query

The DAPall file can be accessed using Marvin query tools. First, we select the stellar and gas velocity dispersion within the effective radius as properties of the datamodel.

from marvin.utils.datamodel.query import datamodel

qdm = datamodel['DR15']
dapall = qdm.groups['DAPall Summary']
dapall.ha_gsigma_1re, dapall.stellar_sigma_1re
>>> <QueryParameter full=dapall.ha_gsigma_1re, name=ha_gsigma_1re, short=ha_gsigma_1re, remote=ha_gsigma_1re, display=Ha_Gsigma_1Re>, 
<QueryParameter full=dapall.stellar_sigma_1re, name=stellar_sigma_1re, short=stellar_sigma_1re, remote=stellar_sigma_1re, display=Stellar_Sigma_1Re>

Now we can use those properties to query the MaNGA sample. Let's select MaNGA galaxies with redshift z < 0.1 and the HYB10 bintype, so we only search in one model LOGCUBE to avoid duplicated results:

from marvin.tools.query import doQuery

myparams = dapall.list_params('full', subset=['ha_gsigma_1re','stellar_sigma_1re'])
q, r = doQuery(search_filter='nsa.z < 0.1 and bintype.name==HYB10', return_params=myparams)

And make a scatter plot of the results using the dapall datamodel for the labels:

from marvin.utils.plot.scatter import plot

x = np.ma.log10(r.getListOf('ha_gsigma_1re', return_all=True))
y = np.ma.log10(r.getListOf('stellar_sigma_1re', return_all=True))

from marvin.utils.datamodel.dap import datamodel
dm = datamodel['DR15']
prop1 = dm['emline_gsigma_ha']
prop2 = dm['stellar_sigma']

plot(x,y,xlabel=prop1,ylabel=prop2)

Accessing DAPALL parameters for individual galaxies

Additionally, the DAPall parameters can be accessed from Maps and ModelCube objects via the dapall attribute. For example, for the MaNGA object (plate-ifu) 7443-12703:

from marvin.tools.maps import Maps
my_map = Maps('7443-12703')
my_map.dapall

For more examples see the Marvin documentation.

The DAPALL file exists as a table in CAS. You can query the DAPALL table under the Query tab in CASjobs. For example, to find all galaxies with a star-formation rate within one effective radii > 5 solar masses per year:

SELECT  
mangaid, objra, objdec, daptype, z, sfr_1re 	
FROM dbo.mangaDapAll	
where sfr_1re > 5 and daptype = 'HYB10-GAU-MILESHC'