source: trunk/test/mopra.py@ 1970

Last change on this file since 1970 was 1838, checked in by Malte Marquarding, 14 years ago

update mopra test script

  • Property svn:eol-style set to native
  • Property svn:executable set to *
  • Property svn:keywords set to Author Date Id Revision
File size: 1.1 KB
RevLine 
[730]1#!/usr/bin/env python
2from asap import *
3rcParams['verbose'] = 0
4print "Test using Mopra dual IF data..."
5print "Reduction test..."
[751]6s = scantable("data/mopra_dualif.rpf")
[730]7q = s.auto_quotient()
[937]8del s
[730]9restfreqs = [110.201,86.243] # 13CO-1/0, SiO the two IF
10q.set_restfreqs(restfreqs,"GHz") # set the restfrequencies, as not in data
11q.set_unit("km/s") # set the unit to be used from now on
12q.set_freqframe("LSRK") # set frequency frame
13q.auto_poly_baseline() # determine and subtract a poly baseline automatically
14q.convert_flux() # K -> Jy
15msk = q.create_mask([-70,20], [20,70]) # line free region - two windows
[1838]16rms = q.stats("rms", msk)
17med = q.stats("median", msk)
[949]18
[1838]19rcParams['plotter.gui'] = 0
[949]20
[1838]21def plotfit(scan, ifno):
22 print "Fitter test for IF%d..." % ifno
23 # select IF
24 scan.set_selection(ifs=ifno)
25 f = fitter()
26 f.set_scan(scan)
27 f.set_function(gauss=ifno+1) # fit gaussian(s)
28 f.fit()
29 f.plot(filename='output/moprafit_if%d.png' % ifno)
30 fp = f.get_parameters()
31
32for ifno in q.getifnos():
33 plotfit(q, ifno)
34
[730]35print "Mopra Test successful"
Note: See TracBrowser for help on using the repository browser.