source: branches/polybatch/test/mopra.py@ 2378

Last change on this file since 2378 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
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
19rcParams['plotter.gui'] = 0
20
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
35print "Mopra Test successful"
Note: See TracBrowser for help on using the repository browser.