#import atnf_sd module
from atnf_sd import *
# create a reader
r = sdreader()
# open an RPFITS file
r.open('/u/mmarquar/zorro/singledish/data/2001-09-01_0332_P363.rpf') 
# create a vector with numbers [0..109]
integrations = range(110)
r.read(integrations)
# get the data out of  the reader
scans = r.getdata()
# close the reader
r = None

# Test sdwriter.
print 'Begin sdwriter tests...'
# Create an MS2 writer.
w = sdwriter('MS2')
# Change to SDFITS output format (the default).
w.setformat()
# Write out the spectra.
w.write(scans, 'test_SDWriter.sdfits')

# print a short summary of the data
scans.summary()
# get the scan with the number '1'
scan = scans.getscan(1)
# get the scan with the number '1'
ref = scans.getscan(2)
# close the data table
scans = None
# open the math server
m = sdmath()
# average the "on" scan
scanav = m.average(scan)
# get rid of the original scan
scan = None
# print a summary of the scan
scanav.summary()
# average the "off" scan
refav = m.average(ref)
# get rid of the original scan
ref = None
# print a summary of the scan
refav.summary()
# form the quotione spectrum
quot = m.quotient(scanav,refav)
# set the cursor to polarisation 0
quot.setpol(0)
# get the spectrum for polarisation 0
v0 = quot.getspectrum()
#print  the first ten channel
print v0[0:10]
# set the cursor to polarisation 0
quot.setpol(1)
# get the spectrum for polarisation 0
v1 = quot.getspectrum()
#print  the first ten channel
print v1[0:10]
# write it to disk for further use
quot.makepersistent('/tmp/myfirstquotient.tbl')