source:
trunk/test/mopra.py
@
743
Last change on this file since 743 was 730, checked in by , 19 years ago | |
---|---|
|
|
File size: 941 bytes |
Rev | Line | |
---|---|---|
[730] | 1 | #!/usr/bin/env python |
2 | from asap import * | |
3 | rcParams['verbose'] = 0 | |
4 | print "Test using Mopra dual IF data..." | |
5 | print "Reduction test..." | |
6 | s = scantable("data/2005-05-08_0350.rpf") | |
7 | q = s.auto_quotient() | |
8 | q2 = quotient(s.get_scan("*_S"),s.get_scan("*_R")) | |
9 | del q2,s | |
10 | restfreqs = [110.201,86.243] # 13CO-1/0, SiO the two IF | |
11 | q.set_restfreqs(restfreqs,"GHz") # set the restfrequencies, as not in data | |
12 | q.set_unit("km/s") # set the unit to be used from now on | |
13 | q.set_freqframe("LSRK") # set frequency frame | |
14 | q.auto_poly_baseline() # determine and subtract a poly baseline automatically | |
15 | q.convert_flux() # K -> Jy | |
16 | msk = q.create_mask([-70,20], [20,70]) # line free region - two windows | |
17 | rms = q.stats("rms",msk) | |
18 | med = q.stats("median",msk) | |
19 | #print q.summary() | |
20 | print "Fitter test..." | |
21 | f = fitter() | |
22 | f.set_scan(q) | |
23 | q.set_cursor(IF=0) | |
24 | f.set_function(gauss=2) # fit two gaussians | |
25 | f.fit() | |
26 | fp = f.get_parameters() | |
27 | print "Mopra Test successful" |
Note: See TracBrowser
for help on using the repository browser.