The tutorial uses the data file 2008-03-12_0932-M999.rpf
You can run these commands in batch mode using the batch file tutorial-1.py. This file must be executed in a terminal (and not in ASAP) with the following command:
localhost> python -i tutorial-1.py
Part I - Reduction
To start ASAP type the follwing at the command line prompt
Define the file name
Load the file into memory
Set the plotting mode
Set the doppler convention
Set the rest frame
Set the observed rest frequency in Hz
Define the channel unit
Form the quotient spectra
Average all scans in time, aligning in velocity
Remove the baseline (set to 2 order). Specify the signal-free channels
Average the two polarisations together
Smooth the data with boxcar, full width = 3
Scale the data according to scaling fudge factor
Eg. With beam efficiency of 0.49 at 86 GHz
Make final plot for saving
Save plot as postscript file
Updated 16 May 2008
Any queries should be directed to Kate Brooks (Kate.Brooks @csiro.au)
localhost> asap
file='2008-03-12_0932-M999.rpf'
s = scantable(file, average=True) print s
plotter.set_mode(stacking='i', panelling='t')
s.set_doppler('RADIO')
s.set_freqframe('LSRK')
s.set_restfreqs([86243.37e6])
s.set_unit('km/s')
q=s.auto_quotient()
av = q.average_time(align=True)
msk = av.create_mask([-200,-50],[50,180]) av.poly_baseline(msk,2)
iav = av.average_pol()
siav = iav.smooth(kernel = 'boxcar', width = 3, insitu = False) plotter.plot(siav)
Eg. With beam efficiency of 0.49 at 86 GHz
iav.scale(2)
plotter.set_range(-20,30) plotter.plot(siav) plotter.set_legend(mode=-1) plotter.set_title(['Orion-SiO'], fontsize=18) plotter.text(10,95,"SiO (2-1 v=1) at 86243.440 MHz", fontsize=12) plotter.text(-19,95,"2008/03/12", fontsize=12) plotter.text(-19,90,"Zoom Mode", fontsize=12)
plotter.save('Orion-SiO.ps')
Updated 16 May 2008
Any queries should be directed to Kate Brooks (Kate.Brooks @csiro.au)