Opened 14 years ago
Closed 14 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.