Changeset 767 for branches/Release12
- Timestamp:
- 12/01/05 14:21:49 (19 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/Release12/doc/cookbook.tex
r763 r767 12 12 \setlength{\parskip}{1ex} 13 13 14 15 14 \title{ATNF Spectral Analysis Package\\User Guide } 16 15 \author{Chris Phillips} 17 18 16 19 17 \newcommand{\cmd}[1]{{\tt #1}} … … 170 168 syntax, specifically they cannot contain ``special'' characters such 171 169 as \@ \$ etc and cannot start with a number (but can contain numbers). 170 Variable (and function) names are case sensitive. 172 171 173 172 \subsubsection{Unix Interaction} … … 355 354 356 355 Examples are the selection of beam, IF and polarisation, spectral unit 357 (e.g. km/s) frequency reference frame (e.g. BARY) and velocity Doppler356 (e.g. km/s), frequency reference frame (e.g. BARY) and velocity Doppler 358 357 type (e.g. RADIO). 359 358 … … 594 593 Quotients can be computed ``automatically''. This requires the data to 595 594 have matching source/reference pairs or one reference for multiple 596 sources. 595 sources. Auto quotient assumes reference scans have a trailing ``\_R'' 596 in the source name for data from Parkes and Mopra, and a trailing 597 ``e'' or ``w'' for data fro, Tidbinbilla. 597 598 598 599 \begin{verbatim} … … 658 659 659 660 The default is to use integration time weighting. The alternative is 660 to use none, variance 661 to use none, variance, Tsys weighting or Tsys \& integration time. 661 662 662 663 \begin{verbatim} … … 688 689 689 690 The function \cmd{auto\_poly\_baseline} can be used to automatically 690 baseline your data with out having to specify channel ranges for691 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 fitas694 well asmark the band edge as invalid.691 baseline your data without having to specify channel ranges for the 692 line free data. It automatically figures out the line-free emission 693 and fits a polynomial baseline to that data. The user can use masks to 694 fix the range of channels or velocity range for the fit as well as 695 mark the band edge as invalid. 695 696 696 697 Simple example … … 740 741 \subsubsection{Brightness Units} 741 742 742 RPFITS files to not contain any information as to whether the telescope743 RPFITS files do not contain any information as to whether the telescope 743 744 calibration was in units of Kelvin or Janskys. On reading the data a 744 745 default value is set depending on the telescope and frequency of … … 1381 1382 1382 1383 \subsection{Tidbinbilla} 1384 1385 The following example is processing of some Tidbinbilla observations 1386 of NH$_3$ at 12~mm. Tidbinbilla has (at the time of observations) a 1387 single polarisation, but can process two IFs simultaneously. In the 1388 example, 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 1391 switched, observing first an reference to the west, then the source 1392 twice and finally reference to the east. 1393 1394 \begin{verbatim} 1395 1396 # Load the rpfits file and inspect 1397 d = scantable('2003-03-16_082048_t0002.rpf') 1398 print d 1399 1400 # Make the quotient spectra 1401 q = d.auto_quotient() 1402 print q 1403 1404 # Plot/select in velocity 1405 q.set_freqframe('LSRK') 1406 q.set_unit('km/s') 1407 1408 # Seperate data from the (1,1)&(2,2) and (4,4) transitions 1409 g1 = q.get_scan(range(6)) # Rows 0..5 1410 g2 = q.get_scan(range(6,12)) # Rows 6..11 1411 1412 # Align data in velocity 1413 g1.freq_align(perif=True) 1414 g2.freq_align(perif=True) 1415 1416 # Average individual scans 1417 a1 = g1.average_time() 1418 a2 = g2.average_time() 1419 1420 # Rpfits file only contrains a single rest frequency. Set both 1421 a1.set_restfreqs(freqs= [23694.4700e6,23722.6336e6]) 1422 1423 plotter.plot(a1,a2) 1424 plotter.set_mode('i','s') 1425 x = raw_input() 1426 1427 a1.auto_poly_baseline() 1428 a2.auto_poly_baseline() 1429 1430 plotter.plot() 1431 1432 a1.smooth('gauss',5) 1433 a2.smooth('gauss',5) 1434 plotter.plot() 1435 1436 \end{verbatim} 1383 1437 1384 1438 \newpage … … 1594 1648 1595 1649 \end{document} 1650
Note:
See TracChangeset
for help on using the changeset viewer.