NRAO Home  >  Green Bank  |  Wiki Topic:    GB > Data > IDLDocumentation > GBTIDLPosSwitchDualBeamRecipe
   Changes | Index | Contents | Search | Statistics | Jump to Topic

Recipe for Reducing Total Power Dual Beam data with GBTIDL

Command Comments
%sdfits -scans=107:122 /home/archive/test-data/tape-0002/TREG_111504 The first step in reducing GBT data with the gbtidl program is to convert (or "fill") the data from the raw GBT FITS files into an sdfits file. This step accomplishes the task. The output file will be called TREG_111504.raw.acs.fits. Additional information is available on the sdfits command
%gbtidl Starts gbtidl and loads the appropriate scripts needed to reduce GBT data.
filein,'TREG_111504.raw.acs.fits' This command identifies the input data set. If the data are contained all in one file, the filein command is appropriate. It is also possible to use a directory of sdfits files as the input, in which case the dirin command is used to identify the input directory.
list Lists the records available in the input data set.
summary Show a summary of GBT scans
getrec,0 Retrieves the first record from the input data set. The data are copied into the primary data container (PDC).
header Shows some header information
show Displays the data from the PDC on the plotter. In this case, the data are still not calibrated and so the spectrum shows just a bandpass shape. If the plotter is set with Auto Update ON, then this step is redundant.
getnod,107 Calibrates the total power nod data from scans 107-108, stores the result in the PDC, and displays it (as long as Auto Update is ON).
freeze turn Auto Update OFF to speed up the following accumulation loop
for i=107,121,2 do begin getnod,i & accum & end Retrieves scans 107 through 122 and adds each to the accumulator, Although there are only two scans in this loop, the number of scans can be increased trivially by changing the limits in the loop. An alternative to using a simple loop is to use the stack to generate a list of scans to be averaged. See an example of how to use the stack here
unfreeze turn Auto Update back ON
ave averages data in the accumulator, stores the result in the primary data container, and displays the result.
bdrop, 10  
edrop, 10 Drop the first and last 10 channels from the display
sclear Clears the accumulator
setregion Prompt the user to specify the regions to be used for baseline fitting by clicking on the plotter. The baseline regions are stored in !g.nregion. The baseline regions can be specified explicitly by setting !g.nregion directly.
bshape, nfit=5 Shows the result of a fifth order polynomial fit. !g.nfit is set to 5.
show clear the overlay
showregion,/off turns off the baseline box overlay
hanning Hanning smooth the data
bshape, nfit=4 Shows the result of a fourth order polynomial fit. !g.nfit is set to 4. If no order is given, the previous value of !g.nfit will be applied.
baseline Subtracts the fourth order polynomial baseline from the data. The residual is stored in the PDC, and displayed.
stats Provides statistics on the spectrum, Use the mouse to define the range of interest
fitgauss Fits a Gaussian profile to the data. The user will be prompted to enter the range to be fit as well as initial guesses for peak and width by clicking on the plotter.

-- JimBraatz

Topic GBTIDLPosSwitchDualBeamRecipe . { Edit | Attach | Ref-By | Printable | Diffs | r1.6 | > | r1.5 | > | r1.4 | More }
Revision r1.6 - 25 Apr 2006 - 19:47 GMT - BobGarwood
Parents: IDLDocumentation
Content copyright © 1999-2007 by the contributing authors.
All material on this collaboration platform is the property of the contributing authors.