1 GAUSS_1D Multi component one-dimension gaussian fitting
GAUSS_1D
This task makes a multi-component non-linear gaussian fit into col-
umns of a table.
It fits a function Y = f(X) where f is the sum of up to five
gaussian, and X and Y data are taken in specified columns of the input
table. The method is identical to that used by the GAUSS command in
CLASS (simplex followed by conjugate gradient). The format for input of
initial parameters is also similar. All messages are sent to a
specified output file. The fitted profiles can be kept in separate
columns of the input table.
2 IN$
TASK\FILE "Input Table" IN$
The name of the input table.
2 LIST$
TASK\FILE "Output File" LIST$
The name of the output formatted file. All messages are written to this
file, including the parameters of the fitted gaussian(s).
2 COLUMN_IN$
TASK\INTEGER "Column Numbers" COLUMN_IN$[2]
The number of the table column used as X and Y data.
2 NLINE$
TASK\INTEGER "Number of gaussian components" NLINE$
The number of gaussian components used in the fit. NLINE$ should not be
larger than 5. If NLINE$ is 0, the program will attempt to guess
initial values for a single component from the moments of the spectrum.
If NLINE$ > 0, it will need some initial values.
2 LINE_1$
TASK\REAL "Parameters for first component" LINE_1$[6]
The initial parameters for the first gaussian component in the profile
fit. These values are entered in free format as follows :
Code, Value, Code, Value, Code, Value, Code, Value
-Intensity- -Position- -Width-
The code is an integer interpreted as follows:
0 adjustable parameter
1 fixed parameter
2 adjustable parameter (head of group)
3 parameter fixed with respect to parameter coded 2 or 4
4 fixed parameter (head of group)
Codes 2 3 and 4 are used to fit depending lines (e.g. HCN, for which
the displacements are 4.842 and -7.064 km/s, or -1.431 and 2.088 MHz,
and line ratios 1:0.6:0.2).
2 LINE_2$
TASK\REAL "Parameters for second component" LINE_2$[6]
The initial parameters for the second gaussian component in the profile
fit.
2 LINE_3$
TASK\REAL "Parameters for second component" LINE_3$[6]
The initial parameters for the third gaussian component in the profile
fit.
2 LINE_4$
TASK\REAL "Parameters for second component" LINE_4$[6]
The initial parameters for the fourth gaussian component in the profile
fit.
2 LINE_5$
TASK\REAL "Parameters for second component" LINE_5$[6]
The initial parameters for the fifth gaussian component in the profile
fit.
2 COLUMN_OUT$
TASK\INTEGER "Column numbers for fitted profiles" COL-
UMN_OUT$['NLINE$']
These are the numbers of columns where values of the fitted com-
ponents will be kept. Specify one empty (or useless) column of the in-
put table for each component that you want to store, and 0 for each
component that you don't want to keep.
1 ENDOFHELP