def test_dat_list():
"""JXP format :: Likely to be Deprecated
"""
if os.getenv('DLA') is None:
assert True
return
# Load
dlas = DLASurvey.neeleman13_tree()
# tests
assert dlas.nsys == 100
def test_dat_list():
"""JXP format :: Likely to be Deprecated
"""
if os.getenv('DLA') is None:
assert True
return
# Load
dlas = DLASurvey.neeleman13_tree()
# tests
assert dlas.nsys == 100
def grab_meta():
""" Generates the meta data needed for the IGMSpec build
Returns
-------
meta : Table
spec_files : list
List of spec_file names
"""
# Load DLA
from pyigm.surveys.dlasurvey import DLASurvey
hdla100 = DLASurvey.neeleman13_tree()
# Cut down to unique QSOs
spec_files = []
names = []
ra = []
dec = []
coords = hdla100.coord
cnt = 0
for coord in coords:
# Load
names.append('J{:s}{:s}'.format(
coord.ra.to_string(unit=u.hour, sep='', pad=True, precision=2),
coord.dec.to_string(sep='', pad=True, precision=1)))
# RA/DEC
ra.append(coord.ra.value)
dec.append(coord.dec.value)
# SPEC_FILE
fname = hdla100._abs_sys[cnt]._datdict['hi res file'].split('/')[-1]
spec_files.append(fname)
cnt += 1
uni, uni_idx = np.unique(names, return_index=True)
nqso = len(uni_idx)
#
meta = Table()
meta['RA_GROUP'] = np.array(ra)[uni_idx]
meta['DEC_GROUP'] = np.array(dec)[uni_idx]
meta['zem_GROUP'] = hdla100.zem[uni_idx]
meta['sig_zem'] = [0.] * nqso
meta['flag_zem'] = [str('UNKN')] * nqso
meta['STYPE'] = [str('QSO')] * nqso
meta['SPEC_FILE'] = np.array(spec_files)[uni_idx]
# Check
assert chk_meta(meta, chk_cat_only=True)
return meta
def grab_meta():
""" Generates the meta data needed for the IGMSpec build
Returns
-------
meta : Table
spec_files : list
List of spec_file names
"""
# Load DLA
from pyigm.surveys.dlasurvey import DLASurvey
hdla100 = DLASurvey.neeleman13_tree()
# Cut down to unique QSOs
spec_files = []
names = []
ra = []
dec = []
coords = hdla100.coord
cnt = 0
for coord in coords:
# Load
names.append('J{:s}{:s}'.format(coord.ra.to_string(unit=u.hour, sep='', pad=True, precision=2),
coord.dec.to_string(sep='', pad=True, precision=1)))
# RA/DEC
ra.append(coord.ra.value)
dec.append(coord.dec.value)
# SPEC_FILE
fname = hdla100._abs_sys[cnt]._datdict['hi res file'].split('/')[-1]
spec_files.append(fname)
cnt += 1
uni, uni_idx = np.unique(names, return_index=True)
nqso = len(uni_idx)
#
meta = Table()
meta['RA_GROUP'] = np.array(ra)[uni_idx]
meta['DEC_GROUP'] = np.array(dec)[uni_idx]
meta['zem_GROUP'] = hdla100.zem[uni_idx]
meta['sig_zem'] = [0.]*nqso
meta['flag_zem'] = [str('UNKN')]*nqso
meta['STYPE'] = [str('QSO')]*nqso
meta['SPEC_FILE'] = np.array(spec_files)[uni_idx]
# Check
assert chk_meta(meta, chk_cat_only=True)
return meta
