next up previous contents index
Next: Reading a spectrum from Up: Cookbook Previous: Cookbook   Contents   Index

A CLASSic Session

   1    device image white
   2    set angle seconds                    
   3    set coordinates equatorial           
   4    file in raw.30m                      
   5    set line 13co(1-0)                   
   6    set source ic348                     
   7    set telescope iram-30m-b30           
   8    set observed 15-aug-1984             
   9    find /offset 0  25                   
  10    set weight time                      
  11    average                                  
  12    set unit v f                         
  13    set mode x -1 14                     
  14    set mode y -0.5 7.5                  
  15    set plot histogram                   
  16    plot                                 
  17    hardcopy /print
  18    hardcopy spectrum.ps  /dev ps fast        
  19    hardcopy spectrum.eps /dev eps color
  20    set window 3 6 8 10                  
  21    base 4 /plot                         
  22    plot                                 
  23    lines 0                              
  24    gauss                                
  25    fit                                  
  26    residual                             
  27    plot                                 
  28    sic delete reduced.30m               
  29    file out reduced.30m new             
  30    swap                                 
  31    write                                
  32    save ic348                           
  33    exit
1
Define the output device.
2-3
Select the coordinate system and angle unit.
4
Open the input file.
5-9
Build the index according to various criteria.
10-11
Averaged all spectra with weights $w_i=\Delta t\Delta
\nu/T_{\rm sys}^2$.
12-16
Plot the averaged spectrum in a given velocity interval, with velocity for the lower axis and rest frequency for the upper axis.
17-20
Make hardcopies (directly to the printer, in a ps or eps file).
21-23
Subtract a polynomial baseline. Plot the fitted baseline. Plot the baseline subtracted spectrum. The baseline does not take into account channels corresponding to $x$ in the two specified windows, in current units (here velocity).
24-26
Perform single Gaussian fitting and plot the result.
27-28
Compute the residuals and plot them.
29-30
Open a new output file.
31
Recover the baseline subtracted spectrum.
32
Write it.
33
Write the input commands in file ic348.class.



Gildas manager 2024-03-28