source: tags/release-1.1.1/InputComplete

Last change on this file was 285, checked in by Matthew Whiting, 17 years ago

Changes worth talking about are:

  • New parameters:
    • The trimming is now controlled by flagTrim. The BLANK etc keywords are always read in -- if they aren't there, flagTrim set to false. User access of flagBlankPix and blankPixVal is disabled. For clarity, flagTrimmed changed to hasBeenTrimmed.
    • Default value of kernMin set to -1 -- then if it is found to be negative, it is set to match kernMaj.
    • Slight changes to parameter output.
  • Fixed bug in findAllStats, so that the correct array is used.
  • Number of pixels for FDR setup set to 2*beamsize only if we are 3-dimensional.
  • Added pixel information to VOTable.
  • Simplified some code in drawMomentCutout
  • Added check for recon array allocation at start of ReconSearch?, as well as checks on dimensionality of cube.
  • Updated images and documentation.
File size: 6.8 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.
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
22imageFile       your-file-here
23flagSubsection  0
24subsection      [*,*,*]
25flagReconExists 0
26reconFile       your-reconstructed-file-here
27flagReconExists 1
28reconFile       your-smoothed-file-here
29
30# OUTPUT RELATED
31# flagLog -- log the intermediate results?   
32# logfile-- the file to put that in.
33# outFile -- the final output list.
34# spectraFile -- the postscript file of spectra
35# flagOutputSmooth -- whether to save the smoothed arrays as a FITS file
36# flagOutputRecon/Resid -- whether to save the reconstruction & residual arrays
37#                          as FITS files
38# flagVOT -- make a VOTable file of the results? 
39# votFile -- the file to put it in
40# flagKarma -- make a Karma annotation file of the results? 
41# karmaFile -- the file to put it in
42# flagMaps -- save postscript versions of the detection and 0th moment maps?
43# detectionMap, momentMap -- the postscript files produced.
44# flagXOutput -- display the moment map in a pgplot window
45
46flagLog         0
47logFile         duchamp-Logfile.txt
48outFile         duchamp-Results.txt
49spectraFile     duchamp-Spectra.ps
50flagOutputSmooth 0
51flagOutputRecon 0
52flagOutputResid 0
53flagVOT         0
54votFile         duchamp-Results.xml
55flagKarma       0
56karmaFile       duchamp-Results.ann
57flagMaps        1
58detectionMap    duchamp-DetectionMap.ps
59momentMap       duchamp-MomentMap.ps
60flagXOutput     1
61
62# FIXING UP THE CUBE
63# flagTrim -- whether to trim blank pixels from the edges of the cube.
64# blankPixValue -- what value they should take. This is overridden by
65#                  FITS header information.
66# flagBaseline -- whether to subtract spectral baselines before searching
67# flagMW -- ignore a range of channels (to remove Milky Way emission)?
68# minMW, maxMW -- the first and last channels to be ignored
69
70flagTrim        0
71blankPixValue   -8.00061
72flagMW          0
73minMW           75
74maxMW           112
75flagBaseline    0
76
77# GENERAL DETECTION RELATED
78# flagStatSec -- Whether to only use a subsection of the cube to calculate
79#                the statistics.
80# StatSec -- The subsection used for statistics calculations. It has the
81#            same format as the pixel subsection.
82# flagNegative -- Are the features being searched for negative (set to true)
83#                 or positive (false -- the default)?
84# snrCut --  How many sigma above mean is a detection when sigma-clipping
85# threshold -- The threshold flux dividing source and non-source. Used instead
86#              of calculating it from the cube's statistics. If not specified,
87#              it will be calculated.
88# flagGrowth -- Should the detections be "grown" to a lower significance value?
89# growthCut -- The lower threshold used when growing detections
90
91flagStatSec     0
92StatSec         [*,*,*]
93flagNegative    0
94snrCut          3.
95threshold       0.
96flagGrowth      0
97growthCut       2.
98
99# RECONSTRUCTION RELATED
100# flagATrous -- Whether or not to do the reconstruction before searching
101# reconDim -- The number of dimensions in which to perform the reconstruction.
102# scaleMin -- The minimum scale (starts at 1) to be included in the
103#              reconstruction
104# snrRecon -- The threshold used in filtering the wavelet coefficient arrays.
105# filterCode -- The code number for the choice of filter to be used in the
106#               reconstruction:  1 = B3-spline filter, 2 = Triangle function,
107#                                3 = Haar wavelet. Other numbers default to 2.
108
109flagATrous      false
110reconDim        1
111scaleMin        1
112snrRecon        4.
113filterCode      1
114
115# SMOOTHING
116# flagSmooth -- whether to smooth the cube.
117# smoothType -- either "spectral" or "spatial"
118# hanningWidth -- the width parameter of the Hanning (spectral smoothing)
119#                 function
120# kernMaj -- the FWHM of the major axis of the 2D spatial smoothing gaussian
121# kernMin -- the FWHM of the minor axis of the 2D spatial smoothing gaussian
122# kernPA  -- the position angle of the major axis of the 2D spatial
123#            smoothing gaussian
124
125flagSmooth      false
126smoothType      spectral
127hanningWidth    5
128kernMaj         3.
129kernMin         3.
130kernPA          0.
131
132# FALSE DISCOVERY RATE METHOD
133# flagFDR -- Whether or not to use the false discovery rate method instead of
134#   simple sigma clipping.
135# alphaFDR -- The "alpha" parameter for the FDR method -- what desired
136#             percentage of discoveries will be false. Expressed as a decimal.
137
138flagFDR         0
139alphaFDR        0.01
140
141# MERGING PARAMETERS
142# flagAdjacent -- Whether to use the "adjacent criterion" to judge if objects
143#                  are to be merged.
144# threshSpatial -- If flagAdjacent=false, this is the maximum spatial
145#                   separation between objects for them to be merged.
146# threshVelocity -- The maximum channel separation between objects for them
147#                   to be merged.
148# minChannels -- The minimum number of consecutive channels an object must have
149#                 for it to be accepted.
150# minPix -- Minimum number of spatial pixels a detected object must have for
151#            it to be counted
152
153flagAdjacent    1
154threshSpatial   3
155threshVelocity  7
156minChannels     3
157minPix          2
158
159# OTHER PARAMETERS
160# verbose -- whether to provide progress indicators on the terminal during
161#             execution
162# drawBorders -- whether to draw borders around the detections on the moment
163#                 maps in the output.
164# drawBlankEdges -- whether to draw an outline around any BLANK regions in the
165#                    moment maps and image cutouts.
166# spectralMethod -- how to plot the spectra in the output -- the spectrum of
167#                    the peak pixel ("peak" -- the default), or integrated over
168#                    all spatial pixels present ("sum")
169# spectralUnits -- what units you want the spectral axis to be plotted in and
170#                   values quoted in (including the units of integrated flux)
171# pixelCentre -- which option to use for quoting the centre of the detection.
172#                 Options are: centroid, average, peak.
173
174verbose         1
175drawBorders     1
176drawBlankEdges  1
177spectralMethod  peak
178spectralUnits   km/s
179pixelCentre     centroid
Note: See TracBrowser for help on using the repository browser.