MINFIT DT mSR Fitting Programme: Quick Reference
Dr Ivan D Reid, mSR Facility PSI
Chapter 1 Where to get this information
Publicly available copies of the documentation for MINFIT and MINUIT
are available in the directory DISK_244_DAT0:[REID.MINFIT] on the PSI
cluster, and in PRV41:[IVAN.MINFIT] at TRIUMF (the TRIUMF copies aren't
updated as often as those at PSI). All files are available as LATEX source
files (.tex), PostScript files (.ps), and as DEC LN03 files
(.ln3). The files are
- MINUIT: A LATEX version of the CERN MINUIT documentation.
- MINFIT: A description of the MINFIT programme, including commands which
have been added to the basic MINUIT commands.
- MINFIT_HELP: This file, originally based on the VMS HELP file used by
MINFIT. It is a quick reference describing all MINUIT and MINFIT commands,
complete with a table of contents useful as a ``cheat sheet''.
- FCN: A brief mathematical treatment of c2 fitting as it pertains
to mSR data. Mainly useful for those contemplating writing additional
fitting subroutines.
Picture Omitted Figure 2.0: How to prepare for MINFIT
Picture Omitted Figure 2.0: How to load run data into MINFIT
MINFIT is currently being re-arranged to allow it to handle fitting
tasks other than just time-differential positive-muon data. The initial
prompt for MINFIT is now
MINFIT option:
Available options are ASYMM (not implemented), DELTAT (formerly just MINFIT),
FFT (in progress), HELP, MUMINUS (not implemented), and QUIT.
ASYMM
The ASYMM option is not yet implemented.
DELTAT
This option provides the functionality of the previous version of
MINFIT. The initial prompt is DELTAT> until a function is defined,
at which time the function name is used as the prompt.
DELTAT COMMANDS
DELTAT accepts commands in the following format:
<command > <arg1 > [arg2] etc.
<command > ¯ is one of the commands listed below,
<argi > are required arguments, if any.
[argi] are optional arguments, if any.
The arguments (if any) are separated from each other and from
the command by one or more blanks or a comma. Commands may be given in
upper or lower case, and may be abbreviated, usually to three
characters. The shortest recognized abbreviations are indicated by the
capitalized part of the commands listed below. Parameter numbers may be
replaced by their full names (without trailing blanks, case-insensitive);
see the full documentation for restrictions.
!
[DELTAT command] ! [comment string] MINFIT command
Anything in a command to the right of an exclamation mark (!),
including the exclamation mark, will be ignored. In particular, a line
starting with an exclamation mark will be totally ignored. This allows
comments to be entered in command-files for batch runs. The COMMENT
command is used for inserting comments in the output.
?
? MINFIT command
This command is an abbreviation for HELP.
@
@[filename] MINFIT command
The @-command is used to read in further commands from an
indirect command file given in the command. Default extension is `.DAT'.
The command file may include further indirect command file calls, but
the files are not ``stacked'' - command eventually returns directly to
SYS$INPUT.
Indirect files are terminated when the end is reached, or when a
call to a further indirect file cannot be satisfied (e.g. the command
``@''will cause termination).
CALlfcn
CALlfcn <iflag > MINUIT command
Instructs Minuit to call the current ``fcn'' subroutine with the
value of IFLAG= <iflag >. This command is not normally required by a
MINFIT user; in any event, a call with <iflag >=3 is to be avoided. A
call with <iflag > greater than five leads MINUIT to ``forget'' the result
of the fit to date.
CLEar
CLEar MINUIT command
Resets all parameter names and values to undefined. Not
normally required by a MINFIT user.
COMment
COMment <comment string > MINFIT command
This allows comments, similar to the !-command, the difference
being that COMMENT echoes the commentline to the terminal or log-file.
Thus ! comments the input, COMMENT comments the output. If a COMMENT
command contains an exclamation mark, the command will be truncated
before being echoed!
CONtour
CONtour <par1 > <par2 > [devs] [ngrid]
MINUIT command
Traces contour lines of the user function with respect to the
two external parameters <par1 > and <par2 >. Other variable parameters of
the function, if any, will have their values fixed at the current values
during the contour tracing. The optional parameter [devs] (default 2.0)
gives the number of standard deviations in each parameter which should
lie entirely within the plotting area. Optional parameter [ngrid]
(default value 25) determines the plot resolution.
This procedure gives a static picture of FCN as function of the
two parameters specified by the user. All the other variable parameters
(if any) are temporarily fixed at their current values. Contour zero is
defined as the current best function value Fmin, and then the i-th
contour is defined as where FCN has the value Fmin+ i2 *UP. The
method, although not very efficient or precise, is very robust, and
capable of revealing unexpected multiple valleys.
FIX
FIX <parno > MINUIT command
Causes parameter <parno > to be removed from the list of
variable parameters, and its value will remain constant (at the current value)
during subsequent minimizations, etc., until another command changes its value
or its status (e.g. RELEASE or RESTORE). A FIXed parameter is
identified in a .PAR file by a negative value for the error estimate (as
opposed to a constant parameter, which has 0.0 for the error estimate).
HELp
HELp [SET] [SHOw]
Causes MINFIT to list the available commands. The list of SET
and SHOw commands must be requested separately. This VMS help file may
also be accessed by using the DCL command
$ HELP/LIBRARY=MINFIT_HELP:MINFIT_HELP
HESse
HESse [maxcalls] MINUIT command
Instructs Minuit to calculate, by finite differences, the
Hessian or error matrix. That is, it calculates the full matrix of
second derivatives of the function with respect to the currently
variable parameters, and inverts it, printing out the resulting error
matrix. The optional argument [maxcalls] specifies the (approximate)
maximum number of function calls after which the calculation will be
stopped.
IMProve
IMProve [maxcalls] MINUIT command
If a previous minimization has converged, and the current values
of the parameters therefore correspond to a local minimum of the
function, this command requests a search for additional distinct local
minima. The optional argument [maxcalls] specifies the (approximate)
maximum number of function calls after which the calculation will be
stopped.
LOOsen
LOOsen [factor] MINFIT command
Expands error estimate on all parameters by [factor]. Default
value of [factor] is 5.0. If the new error would overlap a limit on
a bounded parameter, the error is reduced to avoid the overlap.
MIGrad
MIGrad [maxcalls] [tolerance] MINUIT command
Causes minimization of the function by the method of Migrad, the
most efficient and complete single method, recommended for general
functions (see also MINImize). The minimization produces as a
by-product the error matrix of the parameters, which is usually reliable
unless warning messages are produced. The optional argument [maxcalls]
specifies the (approximate) maximum number of function calls after which
the calculation will be stopped even if it has not yet converged. The
optional argument [tolerance] specifies required tolerance on the
function value at the minimum. The default tolerance is 0.1, and the
minimization will stop when the estimated vertical distance to the
minimum (EDM) is less than 0.001*[tolerance]*UP (see SET ERR).
MINImize
MINImize [maxcalls] [tolerance] MINUIT command
Causes minimization of the function by the method of Migrad, as
does the MIGrad command, but switches to the SIMplex method if Migrad
fails to converge. Arguments are as for MIGrad. Note that command
requires four characters to be unambiguous with MINOS.
MINOs
MINOs [maxcalls] [parno] [parno] ... MINUIT command
Causes a Minos error analysis to be performed on the parameters
whose numbers [parno] are specified. If none are specified, Minos errors
are calculated for all variable parameters. Minos errors may be
expensive to calculate, but are very reliable since they take account of
non-linearities in the problem as well as parameter correlations, and
are in general asymmetric. The optional argument [maxcalls] specifies
the (approximate) maximum number of function calls per parameter
requested, after which the calculation will be stopped for that
parameter. In MINFIT, a default value of 1000 is taken for [maxcalls] if
the first position is left blank. When parameter numbers are explicitly
given, a maximum of six may be entered for each MINOS call.
MNContour
MNContour <par1 > <par2 > [ngrid]
MINUIT command
This command was omitted in the original documentation. It seems
similar to the CONTOUR command, but from the source it appears that the
only optional parameter is [ngrid], the number of points on the contour
(default 20).
The contour calculated by MNContour is dynamic, and represents
the contour Fmin +UP with respect to all the other NPAR-2 parameters
(if any). In statistical terms, this means that MNContour takes account
of the correlations between the two parameters being plotted, and all
the other variable parameters. (If this feature is not wanted, then the
other parameters must be FIXed before calling MNContour.) MNContour
provides the actual coordinates of the points around the contour,
suitable for plotting with a graphics routine or by hand. As a
by-product, MNContour provides the MINOS errors of the two parameters in
question, since these are just the extreme points of the contour (Use
SHOw MINos to see them).
RELease
RELease <parno > MINUIT command
If <parno > is the number of a previously variable parameter
which has been fixed by a command FIX <parno >, then that parameter will
return to variable status. Otherwise a warning message is printed and
the command is ignored. Note that this command operates only on
parameters which were at one time variable and have been FIXed. It
cannot make constant parameters variable; that must be done by
redefining the parameter with a PARAMETER command.
REStore
REStore [code] MINUIT command
If no [code] is specified, this command restores all previously
FIXed parameters to variable status. If [code]=1, then only the last
parameter FIXed is restored to variable status.
SAVe
SAVe MINUIT command
Causes the current parameter values to be saved in a file which
may be edited into a .CON or .PAR file. If the covariance matrix exists,
it is also output in such a format. The unit number is by default 7, and
the user is prompted for the file name (default FOR007.DAT). This is a
pure MINUIT command; in general PARFILE is to be preferred when creating
.PAR files.
SCAn
SCAn [parno] [numpts] [from] [to] MINUIT command
Scans the value of the user function by varying parameter number
[parno], leaving all other parameters fixed at the current value. If
[parno] is not specified, all variable parameters are scanned in
sequence. The number of points [numpts] in the scan is 40 by default,
and cannot exceed 100. The range of the scan is by default 2 standard
deviations on each side of the current best value, but can be specified
as from [from] to [to]. After each scan, if a new minimum is found, the
best parameter values are retained as start values for future scans or
minimizations. The curve resulting from each scan is plotted on the
output unit in order to show the approximate behaviour of the function.
This command is not intended for minimization, but is sometimes useful
finding a reasonable starting point.
SEEk
SEEk [maxcalls] [devs] MINUIT command
Causes a Monte Carlo minimization of the function, by choosing
random values of the variable parameters, chosen uniformly over a
hypercube centered at the current best value. The region size is by
default 3 standard deviations on each side, but can be changed by
specifying the value of [devs].
SET
SET
Allows many options to be set, affecting MINUIT's strategy,
performance, print-out level, etc.
SET BATch
SET BATch MINUIT command
Informs Minuit that it is running in batch mode. Not necessary
in MINFIT.
SET EPSmachine
SET EPSmachine <accuracy > MINUIT command
Informs Minuit that the relative floating point arithmetic
precision is <accuracy >. Minuit determines the nominal precision itself,
but the SET EPS command can be used to override Minuit's own
determination, when the user knows that the FCN function value is not
calculated to the nominal machine accuracy. Typical values of <accuracy >
are between 1.0E-5 and 1.0E-14. This should not be necessary in MINFIT.
SET ERRordef
SET ERRordef <up > MINUIT command
Sets the value of UP (default value=1.0), defining parameter
errors. Minuit defines parameter errors as the change in parameter value
required to change the function value by UP. Normally, for chi-squared
fits UP=1.0, and for negative log likelihood, UP=0.5. MINFIT users
should normally leave this at 1.0.
SET GRAdient
SET GRAdient [force] MINUIT command
Informs Minuit that the user function is prepared to calculate
its own first derivatives and return their values in the array GRAD when
IFLAG=2. If [force] is not specified, Minuit will calculate the FCN
derivatives by finite differences at the current point and compare with
the user's calculation at that point, accepting the user's values only
if they agree. If [force]=1, Minuit does not do its own derivative
calculation, and uses the derivatives calculated in FCN. This is more a
debugging command, and is not normally needed as MINFIT forces MINUIT to
accept gradients whenever they are calculable (all functions, at
present).
SET INTeractive
SET INTeractive MINUIT command
Informs Minuit that it is running interactively. Not necessary
in MINFIT.
SET LIMits
SET LIMits [parno] [lolim] [uplim] MINUIT command
Allows the user to change the limits on one or all parameters.
If no arguments are specified, all limits are removed from all
parameters. If [parno] alone is specified, limits are removed from
parameter [parno]. If all arguments are specified, then parameter
[parno] will be bounded between [lolim] and [uplim]. Limits can be
specified in either order, Minuit will take the smaller as [lolim] and
the larger as [uplim]. However, if [lolim] is equal to [uplim], an error
condition results. See also SET PARAMETER and CHANGE.
SET LINesperpage
SET LINesperpage MINUIT command
Sets the number of lines that Minuit thinks will fit on one page
of output. The default value is 24 for interactive mode and 56 for
batch.
SET NOGradient
SET NOGradient MINUIT command
The inverse of SET GRAdient, instructs Minuit not to use the
first derivatives calculated by the user in FCN.
SET NOWarnings
SET NOWarnings MINUIT command
Supresses Minuit warning messages. SET WARnings is the default.
SET OUTputfile
SET OUTputfile <unitno > MINUIT command
Instructs Minuit to write further output to unit <unitno > (i.e.
the VMS file FOR <unitno >.DAT). This may be useful in some cases where
the progress of a fit is to be recorded, but in general the LIST command
is more help (being a temporary SET OUT command). To return to terminal
output use SET OUT 6 (SET OUT 5 also works).
SET PAGethrow
SET PAGethrow <integer > MINUIT command
Sets the carriage control character for ``new page'' to <integer >.
Thus the value 1 produces a new page, and 0 produces a blank line, on
some output devices (see TOPofpage command).
SET PARameter
SET PARameter <parno > <value >
MINUIT command
Sets the value of parameter <parno > to <value >. The parameter
in question may be variable, fixed, or constant, but must be defined.
Note that it is not possible to set a bounded parameter outside its
bounds (see also CHANGE).
SET PRIntout
SET PRIntout <level > MINUIT command
Sets the print level, determining how much output Minuit will
produce. The allowed values and their meanings are displayed after a
SHOwPRInt command, and are currently <level >=:
- no output except from SHOW commands (MINFIT default)
- minimum output (no starting values or intermediate results)
- normal output
- additional output giving intermediate results
- maximum output, showing progress of minimizations
Note: See also SET WARnings command.
SET PLOtlevel
SET PLOtlevel <level > MINFIT command
Sets the plot level, determining the form of the graph produced by
the PLOT and HISTPLOT commands. The allowable levels are
- plots specified/default histogram (default)
- plots all fitted histograms
- plots single histogram with parameter list
- plots all fitted histograms with parameter list
SET RANdomgenerator
SET RANdomgenerator <seed > MINUIT command
Sets the seed of the random number generator used in SEEK. This
can be any integer between 10,000 and 900,000,000; for example one
which was output from a SHOw RANdom command of a previous run.
SET STRategy
SET STRategy <level > MINUIT command
Sets the strategy to be used in calculating first and second
derivatives and in certain minimization methods. In general, low values
of <level > mean fewer function calls and high values mean more reliable
minimization. Currently allowed values are 0, 1 (default), and 2.
SET WARnings
SET WARnings MINUIT command
Instructs Minuit to output warning messages when suspicious
conditions arise which may indicate unreliable results. This is the
default.
SET WIDthpage
SET WIDthpage MINUIT command
Informs Minuit of the output page width. Default values are 80
for interactive jobs and 120 for batch.
SHOw
SHOw
All SET commands have a corresponding SHOw command. In
addition, the SHOw commands listed below have no corresponding SET
command for obvious reasons. The full list of SHOw commands is printed
in response to the command HELP SHOw. The LIST command is identical to
the SHOW command, except that output is directed to a user-chosen file.
SHOw FUNCTION
SHOw FUNction MINFIT command
Prints the function currently being used. If this is the COMPOUND
function, a listing of the component signals is also given.
SHOw ALL
SHOw ALL MINFIT command
Prints a fuller listing than SHOW RUN, including the chi-squared
values and the current parameter values and relations if any. The output
from SHOW FUNCTION is also given.
SHOw CORrelations
SHOw CORrelations MINUIT command
Calculates and prints the parameter correlations from the error
matrix.
SHOw COVariance
SHOw COVariance MINUIT command
Prints the (external) covariance (error) matrix.
SHOw EIGenvalues
SHOw EIGenvalues MINUIT command
Calculates and prints the eigenvalues of the covariance matrix.
SHOw FCNvalue
SHOw FCNvalue MINUIT command
Prints the current value of FCN, i.e. the chi-squared sum, plus
information on the fit status and number of calls. See also SHOW
CHISQUARE.
SHOw MINoserrors
SHOw MINoserrors MINUIT command
Prints the current minos errors.
SHOw CHIsquare
SHOw CHIsquare MINFIT command
Shows the current chi-squared value, as well as the overall
reduced chi-squared value and that for each histogram.
SHOw ESTimate
SHOw ESTimate MINFIT command
Estimates and displays the normalisation and
background for each histogram. Binning must be defined. The estimate is
based on a least-squares fit of muon decay to the data, so it may be
invalid when there are large, slowly varying features in the histograms.
(See also NORMALIS).
SHOw RELation
SHOw RELation MINFIT command
Shows the parameters currently defined through relations,
plus the defining relation.
SHOw RUN
SHOw RUN MINFIT command
Shows the current run number, title, and subtitles, as well as
details of the histograms being fitted.
SHOw TIMe
SHOw TIMe MINFIT command
Shows the elapsed CPU and real time since MINFIT was started.
SIMplex
SIMplex [maxcalls] [tolerance] MINUIT command
Performs a function minimization using the simplex method of
Nelder and Mead. Minimization terminates either when the function has
been called (approximately) [maxcalls] times, or when the estimated
vertical distance to minimum (EDM) is less than [tolerance]. The
default value of [tolerance] is 0.1*UP (see SET ERR).
STAndard
STAndard MINUIT command
Causes Minuit to execute the Fortran instruction CALL STAND
where STAND is a subroutine supplied by the user. Does nothing in
MINFIT.
TOPofpage
TOPofpage MINUIT command
Causes Minuit to write the character specified in a SET
PAGethrow command (default = ``1'') to column 1 of the output file, which
may or may not position your output medium to the top of a page
depending on the device and system.
BINNING
BINning MINFIT command
Allows the user to enter the initial binning for a fitting session.
Use REBIN for subsequent changes.
CHAnge
CHAnge [parameter] [value] [error] [lower bound] [upper bound]
MINFIT command
Entered with no parameters, CHANGE enters an ``interactive'' mode,
displaying all the parameters and prompting for a parameter to change.
Once a parameter is selected the value, error, and limits can be
entered. Any option left out will be left unchanged. Note that to
retain the value as unchanged, you must enter <,error... >, as entering
<0,error... > will set the value to zero. Note also that it's not
possible to set the error to zero! The interactive mode is exited by
entering a null value for the parameter to be changed. Entered with
just a parameter number, CHANGE does not display the parameter list, but
goes straight to the prompt for new values.
When entered with values as well as a parameter number, CHANGE
will set the parameter and return directly to the command prompt. Note
that you must enter a new value in this mode, although error and limits
can retain their old values as defaults.
FUNction
FUNction [functname] MINFIT command
Allows the user to choose a new function to be fitted to the
existing data file.
HISTPLOT
HIStplot [histogram] [mode] [start (raw) bin] [end (raw) bin] [packing]
MINFIT command
Plots the raw histogram data.
The start and end bins and the packing
refer to the raw data, not the bins being fitted. All parameters default to
the last value. Initial defaults are the lowest numbered histogram being
fitted, mode 1, the whole histogram and binning by one. If the specified end
bin is too large, it is truncated to the last bin.
The modes are similar to PLOT:
- Raw data: The data are plotted.
- Asymmetry: This is not valid for HISTPLOT.
- Residuals: This is not valid for HISTPLOT.
- Logarithmic: The data are plotted with a logarithmic count axis.
LOAD
LOAd [filename] MINFIT command
This command loads a .PAR file. If the name is not given, it
will be prompted for. This is usually used at the start of a fit session.
NEWFile
NEWfile MINFIT command
Allows the user to read in a new data file. The default option
is to retain the same histograms and binning, and the current fit
parameters when fitting the new file. This facilitates analysing a
series of runs with similar properties.
NORmalis
NORmalis [histograms] MINFIT command
Instructs DELTAT to insert its estimate of the normalisation N0 and
background B for each histogram (see SHOW ESTIMATE into the appropriate
parameters. Binning and parameters must be defined. Default behaviour is
to set N0 and B for all histograms.
NEWPARS
NEWPars MINFIT command
Reads in new parameters from a .CON file, or prompts for new parameters
in the case of the COMPOUND function. FUNCTION must be defined.
REBin
REBin MINFIT command
Allows the user to choose new binning for the current histogram(s).
Current fit parameters are retained. Use BINNING for the initial choice.
PARfile
PARfile [filename] MINFIT command
Allows the user to write out the current fit settings into a
.PAR parameter file, to be used as input into a later fit session. The
default file-name is the current function, and the default extension is
.PAR.
PLOT
PLOt [histogram] [mode] [start (fit) bin] [end (fit) bin] [packing]
[binflag] MINFIT command
Plots the data and the fitted function.
The start and end bins and the packing refer to the bins being
fitted; [binflag], if non-zero, plots raw histogram bin numbers instead of
time on the X-axis. All parameters default to the last value. Initial
defaults are the lowest numbered histogram being fitted, mode 1, first
and last fitted bin, binning by one and X=time. If the specified end bin
is too large, it is truncated to the last fitted bin.
There are four plot modes:
- Raw data: The data is plotted, along with the current fit.
- Asymmetry: The current background and the muon decay are extracted
from the data and the current normalisation used to get
the asymmetry.
- Residuals: The relative difference between the data and the current
fit is plotted.
- Logarithmic: This is a raw data plot with logarithmic count axis.
LISt
LISt <something >
This is an extension of the SHOW command, but the data is
written to a file instead of displayed on the screen. The initial use
of LIST will prompt for a file name; the default extension is .FIT.
WRITE
WRIte MINFIT command
This command allows the user to write the current parameters out
to a MuFit-style .TXT file. It functions almost identically to the
corresponding option in MuFit. MINOS errors are given where they exist;
this has necessitated some slight extensions to the .TXT file format (see
the full documentation for details).
QUIT
QUIT MINFIT command
Leaves DELTAT.
FUNCTIONS
There are several pre-defined full functions defined below, as
well as the COMPOUND function which allows the user to build up his own
customised function from discrete signals. Refer to the full
documentation for the precise mathematical form of the functions.
STANDARD
This is a ``standard'' low-field function consisting of a
Lorentzian relaxing muonium signal plus a non-relaxing muon signal.
LORENTZ
A single Lorentzian relaxing muon signal.
GAUSSIAN
A single Gaussian relaxing muon signal.
EXPROOT
A single muon signal with exponential-square-root relaxation.
POWER
A single muon signal with exponential-power relaxation.
ABRAGAM
A single muon signal with Abragam relaxation.
COMPOUND
This is the generic function built up from the signals described
below.
Y(t) = b + N e-t/t |
é ë
|
1 + |
å
i
|
am Gm(t) |
ù û
|
|
|
The signals include the TF functions, plus several longitudinal
and zero-field signals. They can be referenced by the two-letter codes
given to the signals in MuFit. In general, those with codes ending in
`M' are muon signals where the oscillatory frequency is described by the
applied field F, rather than directly in MHz. Refer to the full
documentation for the precise mathematical forms of the signals.
MU
The low-field muonium-plus-muon signal as in the STANDARD
function.
EM
A muon signal with Lorentzian relaxation.
GM
A muon signal with Gaussian relaxation.
RM
A muon signal with exponential-square-root relaxation
PM
A muon signal with exponential-power relaxation.
AM
A muon signal with Abragam relaxation.
E
An oscillation (in MHz) with Lorentzian relaxation.
G
An oscillation (in MHz) with Gaussian relaxation.
R
An oscillation (in MHz) with exponential-square-root relaxation.
P
An oscillation (in MHz) with exponential-power relaxation.
A
An oscillation (in MHz) with Abragam relaxation.
KL
A ``Lorentzian'' relaxing LF/ZF signal.
KG
A Kubo-Toyabe Gaussian relaxing LF/ZF signal,
E0
An oscillating LF/ZF signal with Lorentzian relaxation.
G0
An oscillating LF/ZF signal with Gaussian relaxation.
Relations
It is possible to define relations for various signal amplitudes
ONLY. This allows you to specify, e.g. that all amplitudes in one
histogram add up to some constant. Relations are specified as a string
containing one to four coefficient/parameter pairs in the format
<sign>[Coefficient][*][<P>n]
where the sign is optional for the first pair (if +!). A missing
parameter (Pn) indicates a constant; a missing coefficient is taken as
1.0. The multipliciation sign is optional. See the topics below to see
how to pass the strings to DELTAT.
Manual entry
In manual entry, relations are entered at the value prompt. They
are signified by entering the character 'R' as the first character in
the entry. Only parameters of the type ``sig_hist'' will prompt for and
accept relations; only amplitudes will successfully use them.
Example:
Value [...] or R<relation> for parameter AsyM\_S2H02: r0.2-p3
PAR-files
In a .PAR file, a relation is specified by a character 'R' in
the first column of the line. The relation string starts in column 30
of the line. This allows the PARFILE command to write a .PAR file with
the derived signal name and calculated value in the same columns as for
a normal parameter. Fitting the ``worst-case'' relation into 50 characters
is the main reason a relation is limited to four coefficients.
Only signal amplitudes can be specified as relations.
Example:
...
GM
7 Feld_S2 85.923455 , 0.86745568
8 Sigm_S2 0.45148085 , 3.28593062E-02
R AsyM_S2H01 5.29725650E-02+.170 -1.00P03
R AsyM_S2H02 4.53263546E-02+.170 -1.00P04
...
FFT
The FFT option is in preparation.
HELP
This option leads the user into the VMS HELP file available for MINFIT.
MUMINUS
The MUMINUS option is not yet available.
QUIT
This option quits MINFIT.
MINFIT may be run as a batch process. There are two ways of
passing commands to MINFIT within a batch job:
- Include the commands directly in the batch command file,
after the command $ RUN MINFIT. Don't start MINFIT commands
with a dollar sign, however!
- Put the commands into a separate file and define this to be
MINFIT input with the command
$ define/user sys$input myfile.name
before RUNning MINFIT.
The use of indirect command files (``@filename'') may be
advantageous in batch runs.
Several VMS logicals are used to help locate the data files,
determine the default file-name convention, and to locate the help
library. These are most conveniently defined in your LOGIN.COM file
to avoid having to type them in each time you want to run MINFIT.
MINFIT_HELP
MINFIT_HELP
This logical should be set to the directory containing the help
library MINFIT_HELP.HLB. e.g. on the PSI cluster:
$ DEFINE MINFIT_HELP DISK_244_DAT0:[REID.MINFIT]
MINFIT_DATA
MINFIT_DATA
This logical should be set to the directory containing the data
files to be analysed. The search strategy is first to look in the
current directory, and then if the data file is not found there, it will
be sought in the directory indicated by MINFIT_DATA.
Note: You MUST have the NETMBX privilege to access files on a
different node, even though they are accessible to DCL commands. Ask
the System Manager to grant you this privilege.
Examples:
$ DEFINE MINFIT_DATA USR_SCROOT:[myaccount]
$ DEFINE MINFIT_DATA PSM020::DUB0:[MUSR.EXP.TD_MUSR.DLOG]
MINFIT_NAME
MINFIT_NAME
This logical should contain a pattern which MINFIT uses to
determine the data file's name from the run number. It should be of the
form <string1 >NNNN <string2 >, where the Ns are replaced by MINFIT with
the zero-padded run number. There should be as many Ns as there are
digits in the file name (minimum 4). The file extension should be
explicitly given in <string2 > as no default is assumed.
Examples:
$ DEFINE MINFIT_NAME RNNNN.TX
$ DEFINE MINFIT_NAME deltat_imp_NNNN.dat
File translated from
TEX
by
TTH,
version 2.67. On 21 Aug 2000, 11:18.
|