Opened 17 years ago

Closed 16 years ago

#125 closed defect (fixed)

Further to last message

Reported by: m.burton@… Owned by: Malte Marquarding
Priority: normal Milestone: ASAP 2.2
Component: General Version: 2.0
Severity: normal Keywords:
Cc: m.burton@…

Description (last modified by Malte Marquarding)

Malte

Further to my last message, I noticed that there was a directory problem (asap didn't seem to be able to cope with a space in a directory name), so I removed the space (note this wasn't a problem in the old version of asap on my laptop). When I fix this I then get the following error message. Note that reading in the files is taking an absolute age. Again much longer then before??

Michael

ASAP>execfile('071001_asap.txt')
Auto averaging integrations
Importing 2007-10-01_082248_T040.rpf...
Importing 2007-10-01_091752_T040.rpf...
Importing 2007-10-01_094543_T040.rpf...
Importing 2007-10-01_114242_T040.rpf...
Importing 2007-10-01_114813_T040.rpf...
Importing 2007-10-01_123418_T040.rpf...
Importing 2007-10-01_125836_T040.rpf...
Couldn't find any match.
Couldn't find any match.
------------------------------------------------------------
Traceback (most recent call last):
  File "<ipython console>", line 1, in ?
  File "071001_asap.txt", line 33, in ?
    q = quotient(src,ref,preserve=False)
  File "/Library/Frameworks/Python.framework/Versions/2.4/lib/python2.4/site-packages/asap/asapmath.py", line 110, in quotient
    s = scantable(stm._quotient(source, reference, preserve))
ArgumentError: Python argument types in
    stmath._quotient(stmath, NoneType, NoneType, bool)
did not match C++ signature:
    _quotient(asap::STMathWrapper {lvalue}, asap::ScantableWrapper, asap::ScantableWrapper, bool)

*

This is my script:

# Reduction of Tid 34m data for 2007-10-01, Day of Year 274
# Manually adjust relevant source to be reduced (via s=sXY)
# And also change source name (3 times: once in plot header, twice in output files)
#
files = ['2007-10-01_082248_T040.rpf', '2007-10-01_091752_T040.rpf', '2007-10-01_094543_T040.rpf', '2007-10-01_114242_T040.rpf', '
2007-10-01_114813_T040.rpf', '2007-10-01_123418_T040.rpf', '2007-10-01_125836_T040.rpf']

d = scantable(files)

s1=d.get_scan(range(0,16))
s2=d.get_scan(range(16,32))
s3=d.get_scan(range(32,48))
s4=d.get_scan(range(48,64))
s5=d.get_scan(range(66,84))
s6=d.get_scan(range(84,100))
s7=d.get_scan(range(100,116))
s8=d.get_scan(range(116,132))
s9=d.get_scan(range(132,141))

# change this line manually
s=s9

#
sel = selector()
# select the second IF for this day ONLY
sel.set_ifs([1])
s.set_selection(sel)
# Scaling factor to get Tsys approx in K
scaled = s.scale(3.33,tsys=True,insitu=False)
# Separate the on and off-source scans
ref = scaled.get_scan('*[ew]')
src = scaled.get_scan('*[^ew]')
# take the quotient
q = quotient(src,ref,preserve=False)
q.set_freqframe('LSRK')
q.set_unit('km/s')
# calculate antenna Az, El
q.recalc_azel()
# gain-elevation correction for 34m at Ka-band
q.gain_el([0.1757578928574,0.000981543992646,-1.040739732816e-05])
# approx opacity on day 191
#q.opacity(0.05)
q.freq_align()
# HC5N rest freq
q.set_restfreqs(freqs=31951.777e6)
q.set_fluxunit('Jy')
a = q.average_time()
plotter.set_mode('i','t')
# change plot header manually plus two output file names
#plotter.set_title('HC5N G0.84+0.1 DSS-34  2007-10-01',fontsize=18)
#plotter.set_title('HC5N G0.50+0.1 DSS-34  2007-10-01',fontsize=18)
#plotter.set_title('HC5N G5.89-0.3 DSS-34  2007-10-01',fontsize=18)
#plotter.set_title('HC5N G14.99-0. DSS-34  2007-10-01',fontsize=18)
#plotter.set_title('HC5N G16.59-0. DSS-34  2007-10-01',fontsize=18)
#plotter.set_title('HC5N G22.36+0. DSS-34  2007-10-01',fontsize=18)
#plotter.set_title('HC5N G2.54+0.2 DSS-34  2007-10-01',fontsize=18)
#plotter.set_title('HC5N G28.31-0. DSS-34  2007-10-01',fontsize=18)
plotter.set_title('HC5N G28.15-0. DSS-34  2007-10-01',fontsize=18)



plotter.set_range(xstart=-100,xend=150,ystart=-2,yend=4)

# remove residual baseline
#msk = a.create_mask([-100,-10],[30,230],[310,420])
#a.poly_baseline(msk,order=2)
a.auto_poly_baseline(edge=[20],order=5, threshold=8, insitu=True, plot=False)
# smooth
smoo=a.smooth(kernel='hanning',insitu=False)
plotter.plot(smoo)

a.save('9','ASCII','overwrite')
plotter.save('9.png')

Change History (3)

comment:1 by Malte Marquarding, 17 years ago

Description: modified (diff)
Owner: changed from Malte Marquarding to Malte Marquarding
Reporter: changed from m.burton to m.burton
Status: newassigned

comment:2 by Malte Marquarding, 17 years ago

Cc: m.burton@… added
Reporter: changed from m.burton to m.burton@…

This is broken in the release, but will be fixed in the next release.

# Separate the on and off-source scans
ref = scaled.get_scan('*[ew]')
src = scaled.get_scan('*[^ew]')

Please use this instead

ssel = selector()
ssel.set_query("SRCTYPE==0")
scaled.set_selection(ssel)
src = scaled.copy()
ssel.set_query("SRCTYPE==1")
scaled.set_selection(ssel)
ref = scaled.copy()
scaled.set_selection()

I think this (easier) version might also work.

q = scaled.auto_quotient(preserve=False)

comment:3 by Malte Marquarding, 16 years ago

Resolution: fixed
Status: assignedclosed
Note: See TracTickets for help on using tickets.