source: branches/Release2.0/test/mopra.py@ 1319

Last change on this file since 1319 was 949, checked in by mar637, 19 years ago

asap2 version of 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: 932 bytes
Line 
1#!/usr/bin/env python
2from asap import *
3rcParams['verbose'] = 0
4print "Test using Mopra dual IF data..."
5print "Reduction test..."
6s = scantable("data/mopra_dualif.rpf")
7q = s.auto_quotient()
8del s
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
16rms = q.stats("rms",msk)
17med = q.stats("median",msk)
18
19print "Fitter test..."
20# select IF1
21selection = selector()
22selection.set_ifs(1)
23q.set_selection(selection)
24
25f = fitter()
26f.set_scan(q)
27f.set_function(gauss=2) # fit two gaussians
28f.fit()
29fp = f.get_parameters()
30print "Mopra Test successful"
Note: See TracBrowser for help on using the repository browser.