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