Changeset 767 for branches/Release12


Ignore:
Timestamp:
12/01/05 14:21:49 (18 years ago)
Author:
phi196
Message:

Small typos and Tid example

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/Release12/doc/cookbook.tex

    r763 r767  
    1212\setlength{\parskip}{1ex}
    1313
    14 
    1514\title{ATNF Spectral Analysis Package\\User Guide }
    1615\author{Chris Phillips}
    17 
    1816
    1917\newcommand{\cmd}[1]{{\tt #1}}
     
    170168syntax, specifically they cannot contain ``special'' characters such
    171169as \@ \$ etc and cannot start with a number (but can contain numbers).
     170Variable (and function) names are case sensitive.
    172171
    173172\subsubsection{Unix Interaction}
     
    355354
    356355Examples are the selection of beam, IF and polarisation,  spectral unit
    357 (e.g. km/s) frequency reference frame (e.g. BARY) and velocity Doppler
     356(e.g. km/s), frequency reference frame (e.g. BARY) and velocity Doppler
    358357type (e.g. RADIO).
    359358
     
    594593Quotients can be computed ``automatically''. This requires the data to
    595594have matching source/reference pairs or one reference for multiple
    596 sources.
     595sources. Auto quotient assumes reference scans have a trailing ``\_R''
     596in the source name for data from Parkes and Mopra, and a trailing
     597``e'' or ``w'' for data fro, Tidbinbilla.
    597598
    598599\begin{verbatim}
     
    658659
    659660The default is to use integration time weighting. The alternative is
    660 to use none, variance , Tsys weighting or Tsys \& integration time.
     661to use none, variance, Tsys weighting or Tsys \& integration time.
    661662
    662663\begin{verbatim}
     
    688689
    689690The function \cmd{auto\_poly\_baseline} can be used to automatically
    690 baseline your data with out having to specify channel ranges for
    691 the line free data. It automatically figures out the line-free
    692 emission and fits a polynomial baseline to that data. The user can use
    693 masks to fix the range of channels or velocity range for the fit as
    694 well as mark the band edge as invalid.
     691baseline your data without having to specify channel ranges for the
     692line free data. It automatically figures out the line-free emission
     693and fits a polynomial baseline to that data. The user can use masks to
     694fix the range of channels or velocity range for the fit as well as
     695mark the band edge as invalid.
    695696
    696697Simple example
     
    740741\subsubsection{Brightness Units}
    741742
    742 RPFITS files to not contain any information as to whether the telescope
     743RPFITS files do not contain any information as to whether the telescope
    743744calibration was in units of Kelvin or Janskys.  On reading the data a
    744745default value is set depending on the telescope and frequency of
     
    13811382
    13821383\subsection{Tidbinbilla}
     1384
     1385The following example is processing of some Tidbinbilla observations
     1386of NH$_3$ at 12~mm. Tidbinbilla has (at the time of observations) a
     1387single polarisation, but can process two IFs simultaneously. In the
     1388example, the first half of the observation was observing the (1,1) and
     1389(2,2) transitions simultaneously). The second half observed only the
     1390(4,4) transition due to bandwidth limitations. The data is position
     1391switched, observing first an reference to the west, then the source
     1392twice and finally reference to the east.
     1393
     1394\begin{verbatim}
     1395
     1396# Load the rpfits file and inspect
     1397d = scantable('2003-03-16_082048_t0002.rpf')
     1398print d
     1399
     1400# Make the quotient spectra
     1401q = d.auto_quotient()
     1402print q
     1403
     1404# Plot/select in velocity
     1405q.set_freqframe('LSRK')
     1406q.set_unit('km/s')
     1407
     1408# Seperate data from the (1,1)&(2,2) and (4,4) transitions
     1409g1 = q.get_scan(range(6))     # Rows 0..5
     1410g2 = q.get_scan(range(6,12))  # Rows 6..11
     1411
     1412# Align data in velocity
     1413g1.freq_align(perif=True)
     1414g2.freq_align(perif=True)
     1415
     1416# Average individual scans
     1417a1 = g1.average_time()
     1418a2 = g2.average_time()
     1419
     1420# Rpfits file only contrains a single rest frequency. Set both
     1421a1.set_restfreqs(freqs= [23694.4700e6,23722.6336e6])
     1422
     1423plotter.plot(a1,a2)
     1424plotter.set_mode('i','s')
     1425x = raw_input()
     1426
     1427a1.auto_poly_baseline()
     1428a2.auto_poly_baseline()
     1429
     1430plotter.plot()
     1431
     1432a1.smooth('gauss',5)
     1433a2.smooth('gauss',5)
     1434plotter.plot()
     1435
     1436\end{verbatim}
    13831437
    13841438\newpage
     
    15941648
    15951649\end{document}
     1650
Note: See TracChangeset for help on using the changeset viewer.