Custom Query (241 matches)

ASAP>s = scantable(file, average=True) Found ATMOPRA data. Auto averaging integrations Importing MOPS.rpf...

ASAP>print s

RuntimeError? Traceback (most recent call last)

/home/jose/<ipython console> in <module>()

/usr/lib/pymodules/python2.7/asap/scantable.pyc in str(self)

438 def str(self): 439 tempFile = tempfile.NamedTemporaryFile?()

--> 440 Scantable._summary(self, tempFile.name)

441 tempFile.seek(0) 442 asaplog.clear()

RuntimeError?: Couldn't convert frequency frame.

• __add__ returns by reference is operand is empty

• rowselection is not copied in copy constructor/assignment operarator

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