for row in splat: closest = np.argmin(np.abs(barton['Freq']-row['Freq'])) closest_inds.append(closest) barton_ = barton[np.array(closest_inds)] result = table.hstack([splat, barton_]) result.add_column(table.Column(data=result['Freq_1'] - result['Freq_2'], name='Splat-Barton')) result.add_column(table.Column(data=(result['Freq_1'] - result['Freq_2'])/result['Freq_2']*3e5, name='Splat-Barton_kms')) return result KCls = mt(Splatalogue.query_lines(8*u.GHz, 40000*u.GHz, chemical_name=' KCl', line_lists=['SLAIM'])) KCl = load_barton('39K-35Cl') KCl_diff = match_splat_barton(KCls, KCl[KCl['vu']<4]) KCl_offset_model = get_offset_model(KCl_diff) NaCls = mt(Splatalogue.query_lines(8*u.GHz, 40000*u.GHz, chemical_name=' NaCl', line_lists=['SLAIM'])) NaCl = load_barton('23Na-35Cl') NaCl_diff = match_splat_barton(NaCls, NaCl[NaCl['vu']<4]) NaCl_offset_model = get_offset_model(NaCl_diff) cabezasNaCl = Vizier(row_limit=1e9).get_catalogs('J/ApJ/825/150')[1] cabezasNaCl['nuCalc'][cabezasNaCl['x_nuCalc'] == 'cm-1'] = u.Quantity(cabezasNaCl['nuCalc'][cabezasNaCl['x_nuCalc'] == 'cm-1'], u.cm**-1).to(u.MHz, u.spectral()) cabezasNaCl['x_nuCalc'][cabezasNaCl['x_nuCalc'] == 'cm-1'] = 'MHz' cabezasNaCl.add_column(Column(name='Freq', data=u.Quantity(cabezasNaCl['nuCalc'], u.MHz).to(u.GHz))) cabezas_match = match_splat_barton(cabezasNaCl, NaCl) splat_cabezas_match = match_splat_barton(NaCls, cabezasNaCl)
NaClbarton = load_barton('23Na-35Cl') NaCl_diff = match_splat_barton(NaCls, NaClbarton[NaClbarton['vu'] < 4]) NaCl_offset_model = get_offset_model(NaCl_diff) NaClbarton['Freq'] = NaClbarton['Freq'] + NaCl_offset_model( NaClbarton['vu'], NaClbarton['Ju']) # use Cabezas instead of Barton cabezasNaCl = Vizier(row_limit=1e9).get_catalogs('J/ApJ/825/150')[1] cabezasNaCl['nuCalc'][cabezasNaCl['x_nuCalc'] == 'cm-1'] = u.Quantity( cabezasNaCl['nuCalc'][cabezasNaCl['x_nuCalc'] == 'cm-1'], u.cm**-1).to(u.MHz, u.spectral()) cabezasNaCl['x_nuCalc'][cabezasNaCl['x_nuCalc'] == 'cm-1'] = 'MHz' cabezasNaCl.add_column( Column(name='Freq', data=u.Quantity(cabezasNaCl['nuCalc'], u.MHz).to(u.GHz))) cabezasNaCl.rename_column('J1', 'Ju') cabezasNaCl.rename_column('J0', 'Jl') cabezasNaCl.add_column(Column(name='gu', data=cabezasNaCl['Ju'] * 32 + 16)) cabezasNaCl.add_column(Column(name='gl', data=cabezasNaCl['Jl'] * 32 + 16)) cabezasNaCl.rename_column('V1', 'vu') cabezasNaCl.rename_column('V0', 'vl') cabezasNaCl.rename_column('Eup', 'E_U') cabezasNaCl.add_column( Column(name='E_L', data=cabezasNaCl['E_U'] - cabezasNaCl['Freq'].quantity.to( u.eV, u.spectral()).to(u.K, u.temperature_energy()))) cabezasNaCl.add_column( Column(name='QNs', data=[