Документ взят из кэша поисковой машины. Адрес
оригинального документа
: http://www.atnf.csiro.au/computing/software/asap/tutorials/tutorial.html
Дата изменения: Unknown Дата индексирования: Wed Dec 26 02:05:31 2007 Кодировка: Поисковые слова: п п п п п п п п п п п п п п п п п п п п п п п п п |
localhost> asap
commands()
help(scantable) help(scantable.summary)
s = scantable("2005-05-08_0350.rpf") s.summary()
plotter.set_mode("IF","time") plotter.plot(s) # plot s
q = s.auto_quotient() plotter.plot(q) # plot q
restfreqs = [110.201,86.243] # 13CO-1/0, SiO the two IF q.set_restfreqs(restfreqs,"GHz") # set the restfrequencies, as not in data q.set_unit("km/s") # set the unit to be used from now on q.set_freqframe("LSRK") # set frequency frame plotter.plot() # replot, should show velocity now
q.auto_poly_baseline() # determine and subtract a poly baseline automatically plotter.plot() # replot
plotter.set_range(-35,35) # zoom in bit plotter.set_legend([r"$^{13}CO(1\leftarrow 0)$",r"$SiO$"]) # make nice latex IF labels plotter.set_title(['Mopra Tutorial 2005']) # set the title
q.convert_flux() # K -> Jy plotter.plot()
plotter.save("tutorial.png") plotter.save("tutorial.eps")
msk = q.create_mask([-70,20], [20,70]) # line free region - two windows rms = q.stats("rms",msk) med = q.stats("median",msk)
f = fitter()
sel = selector() sel.set_ifs(0) q.set_selection(sel) f.set_scan(q) f.set_function(gauss=2) # fit two gaussians f.fit() f.plot(
sel.set_ifs(1) q.set_selection(sel) f.set_function(gauss=7) f.fit() f.plot(residual=True) f.plot(components=[0,1,2,3,4,5,6,-1]) f.plot(components=[0,1,2,3,4,5,6,-1],plotparms=True)