source: branches/Release12/test/mopra.py@ 1021

Last change on this file since 1021 was 751, checked in by mar637, 19 years ago

renamed test data

  • Property svn:eol-style set to native
  • Property svn:executable set to *
  • Property svn:keywords set to Author Date Id Revision
File size: 938 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()
8q2 = quotient(s.get_scan("*_S"),s.get_scan("*_R"))
9del q2,s
10restfreqs = [110.201,86.243] # 13CO-1/0, SiO the two IF
11q.set_restfreqs(restfreqs,"GHz") # set the restfrequencies, as not in data
12q.set_unit("km/s") # set the unit to be used from now on
13q.set_freqframe("LSRK") # set frequency frame
14q.auto_poly_baseline() # determine and subtract a poly baseline automatically
15q.convert_flux() # K -> Jy
16msk = q.create_mask([-70,20], [20,70]) # line free region - two windows
17rms = q.stats("rms",msk)
18med = q.stats("median",msk)
19#print q.summary()
20print "Fitter test..."
21f = fitter()
22f.set_scan(q)
23q.set_cursor(IF=0)
24f.set_function(gauss=2) # fit two gaussians
25f.fit()
26fp = f.get_parameters()
27print "Mopra Test successful"
Note: See TracBrowser for help on using the repository browser.