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

Change History (1)

comment:1 by Malte Marquarding, 13 years ago

Resolution: fixed
Status: newclosed

fixed in trunk. for existing asap-4.0.0 use parallactify(True) explicitly on every scantable or copy.

Note: See TracTickets for help on using tickets.