source: tags/release-1.1.11b/InputComplete

Last change on this file was 801, checked in by MatthewWhiting, 13 years ago

Updating documentation & default inputs for recent changes to beam calcs

File size: 11.4 KB
Line 
1# FULL PARAMETER LIST
2# -------------------
3#
4# This file contains the full list of user-definable parameters. They are
5#  grouped here according to their use, and the values quoted are the defaults:
6#  if a given parameter is not included in your parameter file, this is the
7#  value it will take.
8#
9#
10# INPUT RELATED
11# imageFile -- the FITS image.  THERE IS NO DEFAULT FOR THIS.
12# flagSubsection -- whether to get a subsection of that image.
13# subsection -- the subsection to read in, in format [x1:x2,y1:y2,z1:z2],
14#               or * to indicate the full range
15# flagReconExists -- set to true if the reconstructed cube exists as a FITS
16#                    file (ie. it has been saved from a previous Duchamp run)
17# reconFile -- the location of the FITS file containing the reconstructed cube
18# flagSmoothExists -- set to true if the smoothed cube exists as a FITS
19#                    file (ie. it has been saved from a previous Duchamp run)
20# smoothFile -- the location of the FITS file containing the smoothed cube
21# usePrevious -- whether to read the list of detected objects from a previously-
22#                created log file
23# objectList -- the list of objects to plot in individual plots: used with usePrevious=true
24
25imageFile       your-file-here
26flagSubsection  0
27subsection      [*,*,*]
28flagReconExists 0
29reconFile       your-reconstructed-file-here
30flagSmoothExists 1
31smoothFile      your-smoothed-file-here
32usePrevious     false
33objectList      *
34
35# OUTPUT RELATED
36# flagLog -- log the intermediate results?   
37# logfile-- the file to put that in.
38# outFile -- the final output list.
39# flagSeparateHeader -- whether to write the header information
40#                        (i.e. parameters, statistics and number of detections)
41#                        to a separate file from the outFile (so that the
42#                        outFile just contains the table of detections).
43# headerFile -- the file to write the header information to.
44# spectraFile -- the postscript file of spectra
45# flagTextSpectra -- whether to save a text file detailing the spectra
46#                    of each source
47# spectraTextFile -- the file in which to save the text-formatted spectra
48# flagOutputMomentMap -- whether to save a FITS file containing the moment-0 map
49# fileOutputMomentMap -- the filename for the moment-0 FITS file. If not given, it defaults
50#                   to the format shown (for input file image.fits)
51# flagOutputMask -- whether to save a FITS file containing a mask
52#                   array, showing the locations of detected objects
53# fileOutputMask -- the filename for the mask file. If not given, it defaults
54#                   to the format shown (for input file image.fits)
55# flagMaskWithObjectNum -- if true, the detected pixels in the mask array have
56#                          the object ID number, else they all have the value 1.
57# flagOutputSmooth -- whether to save the smoothed arrays as a FITS file
58# fileOutputSmooth -- the filename for the smoothed array file. If not given, it defaults
59#                     to the format shown (for input file image.fits). See documentation for
60#                     explanation of format.
61# flagOutputRecon/Resid -- whether to save the reconstruction & residual arrays
62#                          as FITS files
63# fileOutputRecon/Resid -- the filename for the reconstruction or residual files.
64#                          If not given, it defaults to the format shown (for input
65#                          file image.fits). See documentation for explanation of format.
66# flagVOT -- make a VOTable file of the results? 
67# votFile -- the file to put it in
68# flagKarma -- make a Karma annotation file of the results? 
69# karmaFile -- the file to put it in
70# flagMaps -- save postscript versions of the detection and 0th moment maps?
71# detectionMap, momentMap -- the postscript files produced.
72# flagXOutput -- display the moment map in a pgplot window
73# newFluxUnits -- brightness units to convert the array values to
74#                 (eg. from Jy/beam to mJy/beam)   
75# precFlux -- Desired precision for Flux value in outputs.
76# precVel -- Desired precision for Velocity/Frequency value in outputs.
77# precSNR -- Desired precision for peak SNR value in ouputs.
78
79flagLog         0
80logFile         duchamp-Logfile.txt
81outFile         duchamp-Results.txt
82flagSeparateHeader 0
83headerFile      duchamp-Results.hdr
84spectraFile     duchamp-Spectra.ps
85flagTextSpectra 0
86spectraTextFile duchamp-Spectra.txt
87flagOutputMomentMap  0
88fileOutputMomentMap  image.MOM0.fits
89flagOutputMask  0
90fileOutputMask  image.MASK.fits
91flagMaskWithObjectNum 0
92flagOutputSmooth 0
93fileOutputSmooth image.SMOOTH-ND-X.fits
94flagOutputRecon 0
95fileOutputRecon image.RECON-A-B-C-D.fits
96flagOutputResid 0
97fileOutputResid image.RESID-A-B-C-D.fits
98flagVOT         0
99votFile         duchamp-Results.xml
100flagKarma       0
101karmaFile       duchamp-Results.ann
102flagMaps        1
103detectionMap    duchamp-DetectionMap.ps
104momentMap       duchamp-MomentMap.ps
105flagXOutput     1
106newFluxUnits    no-default
107precFlux        3
108precVel         3
109precSNR         2
110
111# FIXING UP THE CUBE
112# flagTrim -- whether to trim blank pixels from the edges of the cube.
113# flagBaseline -- whether to subtract spectral baselines before searching
114# flagMW -- ignore a range of channels (to remove Milky Way emission)?
115# minMW, maxMW -- the first and last channels to be ignored
116
117flagTrim        0
118flagMW          0
119minMW           75
120maxMW           112
121flagBaseline    0
122
123# GENERAL DETECTION RELATED
124# flagStatSec -- Whether to only use a subsection of the cube to calculate
125#                the statistics.
126# StatSec -- The subsection used for statistics calculations. It has the
127#            same format as the pixel subsection.
128# flagRobustStats -- Shall we use robust statistics to characterise
129#                    the noise in the image?
130# flagNegative -- Are the features being searched for negative (set to true)
131#                 or positive (false -- the default)?
132# snrCut --  How many sigma above mean is a detection when sigma-clipping
133# threshold -- The threshold flux dividing source and non-source. Used instead
134#              of calculating it from the cube's statistics. If not specified,
135#              it will be calculated.
136# flagGrowth -- Should the detections be "grown" to a lower significance value?
137# growthCut -- The lower threshold used when growing detections
138# growthThreshold -- The lower threshold, used in conjunction with "threshold"
139# beamarea -- the area of the beam in pixels (equivalent to the old beamsize
140#             parameter). This value is overridden by the BMAJ, BMIN, BPA header
141#             parameters if present.
142# beamFWHM -- the full-width at half-maximum of the beam, in pixels. Where given,
143#             it overrides beamarea, but is overridden by the BMAJ, BMIN, BPA
144#             headers.
145# searchType -- how the searching is done. Either "spatial", where each 2D map is
146#               searched then detections combined, or "spectral", where each 1D
147#               spectrum is searched then detections combined.
148
149flagStatSec     0
150StatSec         [*,*,*]
151flagRobustStats 1
152flagNegative    0
153snrCut          3.
154threshold       0.
155flagGrowth      0
156growthCut       2.
157growthThreshold 0.
158beamArea        0.
159beamFWHM        0.
160searchType      spatial
161
162# RECONSTRUCTION RELATED
163# flagATrous -- Whether or not to do the reconstruction before searching
164# reconDim -- The number of dimensions in which to perform the reconstruction.
165# scaleMin -- The minimum scale (starts at 1) to be included in the
166#              reconstruction
167# scaleMax -- The maximum scale to be included in the reconstruction. If it
168#              is <=0 then the maximum scale is calculated from the size of
169#              the array being reconstructed.
170# snrRecon -- The threshold used in filtering the wavelet coefficient arrays.
171# filterCode -- The code number for the choice of filter to be used in the
172#               reconstruction:  1 = B3-spline filter, 2 = Triangle function,
173#                                3 = Haar wavelet. Other numbers default to 2.
174
175flagATrous      false
176reconDim        1
177scaleMin        1
178scaleMax        0
179snrRecon        4.
180filterCode      1
181
182# SMOOTHING
183# flagSmooth -- whether to smooth the cube.
184# smoothType -- either "spectral" or "spatial"
185# hanningWidth -- the width parameter of the Hanning (spectral smoothing)
186#                 function
187# kernMaj -- the FWHM of the major axis of the 2D spatial smoothing gaussian
188# kernMin -- the FWHM of the minor axis of the 2D spatial smoothing gaussian
189# kernPA  -- the position angle of the major axis of the 2D spatial
190#            smoothing gaussian
191
192flagSmooth      false
193smoothType      spectral
194hanningWidth    5
195kernMaj         3.
196kernMin         3.
197kernPA          0.
198
199# FALSE DISCOVERY RATE METHOD
200# flagFDR -- Whether or not to use the false discovery rate method instead of
201#   simple sigma clipping.
202# alphaFDR -- The "alpha" parameter for the FDR method -- what desired
203#             percentage of discoveries will be false. Expressed as a decimal.
204# FDRnumCorChan -- The number of neighbouring channels that are
205#                  correlated. This is used in the FDR algorithm
206
207flagFDR         0
208alphaFDR        0.01
209FDRnumCorChan   2
210
211# MERGING PARAMETERS
212# flagAdjacent -- Whether to use the "adjacent criterion" to judge if objects
213#                  are to be merged.
214# threshSpatial -- If flagAdjacent=false, this is the maximum spatial
215#                   separation between objects for them to be merged.
216# threshVelocity -- The maximum channel separation between objects for them
217#                   to be merged.
218# minChannels -- The minimum number of consecutive channels an object must have
219#                 for it to be accepted.
220# minPix -- Minimum number of spatial pixels a detected object must have for
221#            it to be counted
222# minVoxels -- Minimum number of voxels a detected object must have for it to be
223#              counted
224# flagRejectBeforeMerge -- If true, reject sources for not meeting the above
225#                          criteria before the merging stage
226# flagTwoStageMerging -- If true, do an initial merge of a new source to the list
227#                        (that is, merge it with the first source in the list that
228#                        is close, but stop there)
229
230flagAdjacent           1
231threshSpatial          3
232threshVelocity         7
233minChannels            3
234minPix                 2
235minVoxels              4
236flagRejectBeforeMerge  0
237flagTwoStageMerging    1
238
239# OTHER PARAMETERS
240# verbose -- whether to provide progress indicators on the terminal during
241#             execution
242# drawBorders -- whether to draw borders around the detections on the moment
243#                 maps in the output.
244# drawBlankEdges -- whether to draw an outline around any BLANK regions in the
245#                    moment maps and image cutouts.
246# spectralMethod -- how to plot the spectra in the output -- the spectrum of
247#                    the peak pixel ("peak" -- the default), or integrated over
248#                    all spatial pixels present ("sum")
249# spectralUnits -- what units you want the spectral axis to be plotted in and
250#                   values quoted in (including the units of integrated flux)
251# pixelCentre -- which option to use for quoting the centre of the detection.
252#                 Options are: centroid, average, peak.
253# sortingParam -- which parameter to sort the final detection list on.
254#                  Options are: x-value, y-value, z-value, ra, dec, vel, w50,
255#                                iflux,  pflux, snr.
256
257verbose         1
258drawBorders     1
259drawBlankEdges  1
260spectralMethod  peak
261spectralUnits   km/s
262pixelCentre     centroid
263sortingParam    vel
Note: See TracBrowser for help on using the repository browser.