source:
tags/junk/test/mopra.py
@
2041
Last change on this file since 2041 was 1273, checked in by , 18 years ago | |
---|---|
|
|
File size: 999 bytes |
Line | |
---|---|
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/mopra_dualif.rpf") |
7 | q = s.auto_quotient() |
8 | del s |
9 | restfreqs = [110.201,86.243] # 13CO-1/0, SiO the two IF |
10 | q.set_restfreqs(restfreqs,"GHz") # set the restfrequencies, as not in data |
11 | q.set_unit("km/s") # set the unit to be used from now on |
12 | q.set_freqframe("LSRK") # set frequency frame |
13 | q.auto_poly_baseline() # determine and subtract a poly baseline automatically |
14 | q.convert_flux() # K -> Jy |
15 | msk = q.create_mask([-70,20], [20,70]) # line free region - two windows |
16 | rms = q.stats("rms",msk) |
17 | med = q.stats("median",msk) |
18 | |
19 | print "Fitter test..." |
20 | # select IF1 |
21 | selection = selector() |
22 | selection.set_ifs(1) |
23 | q.set_selection(selection) |
24 | |
25 | rcParams['plotter.gui'] = 0 |
26 | f = fitter() |
27 | f.set_scan(q) |
28 | f.set_function(gauss=2) # fit two gaussians |
29 | f.fit() |
30 | f.plot(filename='output/moprafit.png') |
31 | fp = f.get_parameters() |
32 | print "Mopra Test successful" |
Note: See TracBrowser
for help on using the repository browser.