Ignore:
Timestamp:
03/27/06 19:36:17 (18 years ago)
Author:
mar637
Message:

asap2 version of test script

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/test/parkes-pol.py

    r744 r937  
    33
    44rcParams['verbose'] = 0
     5rcParams['plotter.ganged'] = 0
    56
    67# Don't plot to the screen...
     
    1617data_1665.set_freqframe('LSRK')
    1718
    18 # Make a copy for the 1667 transition
    19 data_1667 = data_1665.copy()
    20 data_1667.set_restfreqs(lines = ['OH1667'])
    21 data_1667.set_unit('km/s')
    22 data_1667.set_freqframe('LSRK')
     19# Look at the first scan
     20from asap._asap import selector
     21selection = selector()
     22selection._setscans([0])
     23data_1665._setselection(selection)
    2324
     25d1_5 = data_1665.copy()
     26d1_7 = data_1665.copy()
    2427
    25 # Look at the first scan
    26 d1_5 = data_1665.get_scan(0)
    27 d1_7 = data_1667.get_scan(0)
     28d1_7.set_restfreqs([1667.0],'MHz')
     29# merge the two scans back together into a new scantable
     30plotscans = merge(d1_5,d1_7)
     31print plotscans.summary()
     32del d1_5,d1_7
    2833
    2934# Baseline
    30 msk = d1_5.create_mask([-30,-25],[-5,0])
    31 d1_5.poly_baseline(msk,1)
    32 msk = d1_7.create_mask([-30,-25],[-5,0])
    33 d1_7.poly_baseline(msk,1)
    34 
     35msk = plotscans.create_mask([-30,-25],[-5,0])
     36plotscans.poly_baseline(msk,1)
    3537# Plot the results
    36 plotter.plot(d1_5,d1_7)
    3738plotter.set_mode('p','s')
    3839plotter.set_layout(2,1)
    3940plotter.set_range(-30,0)
    40 plotter.set_cursor(pol=['I','Q','U','V'])
    41 plotter.set_cursor(pol=['I','Plinear'])
    42 plotter.set_cursor(pol=['RR','LL'])
    43 plotter.save('output/parkes_rrll.png')
     41selection._reset()
     42selection._setpolstrings(['I','Q', 'U', 'V'])
     43#plotter.set_selection(selection)
     44#selection._setpolstrings(['I','Plinear'])
     45#plotter.set_selection(selection)
     46#selection._setpolstrings(['RR','LL'])
     47plotter.plot(plotscans)
     48plotter.set_selection(selection)
     49plotter.save('/nfs/wwwpeople/mmarquar/parkes.png', dpi=80)
     50#plotter.save('output/parkes_rrll.png')
    4451
    4552print "Parkes-Pol Test successful"
Note: See TracChangeset for help on using the changeset viewer.