Opened 13 years ago
Closed 13 years ago
#259 closed defect (fixed)
parallactify flag isn't preserved when scantable is copied.
Reported by: | Malte Marquarding | Owned by: | Malte Marquarding |
---|---|---|---|
Priority: | normal | Milestone: | Unified development |
Component: | General | Version: | 2.0 |
Severity: | normal | Keywords: | |
Cc: | James.Green@… |
Description
Hi Malte, So I've just tried a few tests in asap and I am fairly confident that when you use get_scan the resultant data file forgets the parallactify option - I tried two simple methods, first setting parallactify on the initial file: ASAP>data_1665 = scantable('2005-10-27_0154-P484.rpf') Cross polarization present Found ATPKSHOH data. Auto averaging integrations Importing 2005-10-27_0154-P484.rpf... ASAP>data_1665.parallactify(True) ASAP>sourcedata=data_1665.get_scan(55) ASAP>sourcedata.set_unit('km/s') ASAP>dmsk2=sourcedata.create_mask([-30,-20],[20,30]) The current mask window unit is km/s ASAP>sourcedata.poly_baseline(dmsk2,1) ASAP>plotter.plot(sourcedata) A new scantable is set to the plotter. The masks and data selections are reset. ASAP>stokesp = sourcedata.convert_pol('stokes') ASAP>plotter.plot(stokesp) A new scantable is set to the plotter. The masks and data selections are reset. ASAP>sourcedata.parallactify(True) ASAP>stokesp = sourcedata.convert_pol('stokes') ASAP>plotter.plot(stokesp) A new scantable is set to the plotter. The masks and data selections are reset. *This gives the wrong stokes q and u. The second has parallactify after the get_scan and gives the correct q and u: ASAP>data_1665 = scantable('2005-10-27_0154-P484.rpf') Cross polarization present Found ATPKSHOH data. Auto averaging integrations Importing 2005-10-27_0154-P484.rpf... ASAP>sourcedata=data_1665.get_scan(55) ASAP>sourcedata.parallactify(True) ASAP>sourcedata.set_unit('km/s') ASAP>dmsk2=sourcedata.create_mask([-30,-20],[20,30]) The current mask window unit is km/s ASAP>sourcedata.poly_baseline(dmsk2,1) ASAP>plotter.plot(sourcedata) A new scantable is set to the plotter. The masks and data selections are reset. ASAP>stokesp = sourcedata.convert_pol('stokes') ASAP>plotter.plot(stokesp) A new scantable is set to the plotter. The masks and data selections are reset. Cheers, Jimi
Note:
See TracTickets
for help on using tickets.
fixed in trunk. for existing asap-4.0.0 use
parallactify(True)
explicitly on every scantable or copy.