Implementation in C++ of the Minuit package written by Fred James.
This is a straightforward conversion of the original Fortran version.
The main changes are:
The MINUIT package acts on a multiparameter Fortran function to which one must give the generic name FCN
. In the ROOT implementation, the function FCN
is defined via the MINUIT SetFCN member function when an Histogram.Fit command is invoked. The value of FCN
will in general depend on one or more variable parameters.
To take a simple example, in case of ROOT histograms (classes TH1C,TH1S,TH1F,TH1D) the Fit function defines the Minuit fitting function as being H1FitChisquare or H1FitLikelihood depending on the options selected. H1FitChisquare calculates the chisquare between the user fitting function (gaussian, polynomial, user defined,etc) and the data for given values of the parameters. It is the task of MINUIT to find those values of the parameters which give the lowest value of chisquare.
For variable parameters with limits, MINUIT uses the following transformation:
\[ P_{\mathrm{int}} = \arcsin \left( 2\: \frac{P_{\mathrm{ext}}-a}{b-a} - 1 \right) P_{\mathrm{ext}} = a + \frac{b - a}{2} \left( \sin P_{\mathrm{int}} + 1 \right) \]
so that the internal value \(P_{\mathrm{int}}\) can take on any value, while the external value \(P_{\mathrm{ext}}\) can take on values only between the lower limit \(a\) and the upper limit \(b\). Since the transformation is necessarily non-linear, it would transform a nice linear problem into a nasty non-linear one, which is the reason why limits should be avoided if not necessary. In addition, the transformation does require some computer time, so it slows down the computation a little bit, and more importantly, it introduces additional numerical inaccuracy into the problem in addition to what is introduced in the numerical calculation of the FCN value. The effects of non-linearity and numerical roundoff both become more important as the external value gets closer to one of the limits (expressed as the distance to nearest limit divided by distance between limits). The user must therefore be aware of the fact that, for example, if he puts limits of \((0,10^{10})\) on a parameter, then the values \(0.0\) and \(1.0\) will be indistinguishable to the accuracy of most machines.
The transformation also affects the parameter error matrix, of course, so Minuit does a transformation of the error matrix (and the `‘parabolic’' parameter errors) when there are parameter limits. Users should however realize that the transformation is only a linear approximation, and that it cannot give a meaningful result if one or more parameters is very close to a limit, where \(\partial P_{\mathrm{ext}} / \partial P_{\mathrm{int}} \approx 0\). Therefore, it is recommended that:
MINUIT offers the user a choice of several minimization algorithms. The MIGRAD algorithm is in general the best minimizer for nearly all functions. It is a variable-metric method with inexact line search, a stable metric updating scheme, and checks for positive-definiteness. Its main weakness is that it depends heavily on knowledge of the first derivatives, and fails miserably if they are very inaccurate.
If parameter limits are needed, in spite of the side effects, then the user should be aware of the following techniques to alleviate problems caused by limits:
If MIGRAD converges normally to a point where no parameter is near one of its limits, then the existence of limits has probably not prevented MINUIT from finding the right minimum. On the other hand, if one or more parameters is near its limit at the minimum, this may be because the true minimum is indeed at a limit, or it may be because the minimizer has become `‘blocked’' at a limit. This may normally happen only if the parameter is so close to a limit (internal value at an odd multiple of \(\pm \frac{\pi}{2}\) that MINUIT prints a warning to this effect when it prints the parameter values.
The minimizer can become blocked at a limit, because at a limit the derivative seen by the minimizer \(\partial F / \partial P_{\mathrm{int}}\) is zero no matter what the real derivative \(\partial F / \partial P_{\mathrm{ext}}\) is.
\[ \frac{\partial F}{\partial P_{\mathrm{int}}} = \frac{\partial F}{\partial P_{\mathrm{ext}}} \frac{\partial P_{\mathrm{ext}}}{\partial P_{\mathrm{int}}} = \frac{\partial F}{\partial P_{\mathrm{ext}}} = 0 \]
In the best case, where the minimum is far from any limits, MINUIT will correctly transform the error matrix, and the parameter errors it reports should be accurate and very close to those you would have got without limits. In other cases (which should be more common, since otherwise you wouldn't need limits), the very meaning of parameter errors becomes problematic. Mathematically, since the limit is an absolute constraint on the parameter, a parameter at its limit has no error, at least in one direction. The error matrix, which can assign only symmetric errors, then becomes essentially meaningless.
There are two kinds of problems that can arise: the reliability of MINUIT's error estimates, and their statistical interpretation, assuming they are accurate.
For discussion of basic concepts, such as the meaning of the elements of the error matrix, or setting of exact confidence levels see:
MINUIT always carries around its own current estimates of the parameter errors, which it will print out on request, no matter how accurate they are at any given point in the execution. For example, at initialization, these estimates are just the starting step sizes as specified by the user. After a HESSE step, the errors are usually quite accurate, unless there has been a problem. MINUIT, when it prints out error values, also gives some indication of how reliable it thinks they are. For example, those marked CURRENT GUESS ERROR
are only working values not to be believed, and APPROXIMATE ERROR
means that they have been calculated but there is reason to believe that they may not be accurate.
If no mitigating adjective is given, then at least MINUIT believes the errors are accurate, although there is always a small chance that MINUIT has been fooled. Some visible signs that MINUIT may have been fooled are:
EDM
too big (estimated Distance to Minimum).The best way to be absolutely sure of the errors, is to use independent'' calculations and compare them, or compare the calculated errors with a picture of the function. Theoretically, the covariance matrix for a
physical'' function must be positive-definite at the minimum, although it may not be so for all points far away from the minimum, even for a well-determined physical problem. Therefore, if MIGRAD reports that it has found a non-positive-definite covariance matrix, this may be a sign of one or more of the following:
On its way to the minimum, MIGRAD may have traversed a region which has unphysical behaviour, which is of course not a serious problem as long as it recovers and leaves such a region.
If the matrix is not positive-definite even at the minimum, this may mean that the solution is not well-defined, for example that there are more unknowns than there are data points, or that the parameterisation of the fit contains a linear dependence. If this is the case, then MINUIT (or any other program) cannot solve your problem uniquely, and the error matrix will necessarily be largely meaningless, so the user must remove the under-determinedness by reformulating the parameterisation. MINUIT cannot do this itself.
It is possible that the apparent lack of positive-definiteness is in fact only due to excessive roundoff errors in numerical calculations in the user function or not enough precision. This is unlikely in general, but becomes more likely if the number of free parameters is very large, or if
the parameters are badly scaled (not all of the same order of magnitude), and correlations are also large. In any case, whether the non-positive-definiteness is real or only numerical is largely irrelevant, since in both cases the error matrix will be unreliable and the minimum suspicious.
For questions of parameter dependence, see the discussion above on positive-definiteness.
Possible other mathematical problems are the following:
Be especially careful of exponential and factorial functions which get big very quickly and lose accuracy.
The function may have unphysical local minima, especially at infinity in some variables.
This concerns the way Minuit reports the symmetric (or parabolic) errors on parameters. It does not apply to the errors reported from Minos, which are in general asymmetric.
The symmetric errors reported by Minuit are always calculated from the covariance matrix, assuming that this matrix has been calculated, usually as the result of a Migrad minimization or a direct calculation by Hesse which inverts the second derivative matrix.
When there are no limits on the parameter in question, the error reported by Minuit should therefore be exactly equal to the square root of the corresponding diagonal element of the error matrix reported by Minuit.
However, when there are limits on the parameter, there is a transformation between the internal parameter values seen by Minuit (which are unbounded) and the external parameter values seen by the user in FCN (which remain inside the desired limits). Therefore the internal error matrix kept by Minuit must be transformed to an external error matrix for the user. This is done by multiplying the (I,J)th element by DEXDIN(I)*DEXDIN(J), where DEXDIN is the derivative of the external value with respect to the internal value at the minimum. This is a linearisation of the transformation, and is the only way to produce an error matrix in external coordinates meaningful to the user. But when reporting the individual parabolic errors for limited parameters, Minuit can do a little better, so it does. In this case, Minuit actually transforms the ends of the internal "error bar" to external coordinates and reports the length of this transformed interval. Strictly speaking, it is now asymmetric, but since the origin of the asymmetry is only an artificial transformation it does not make much sense, so the transformed errors are symmetrized.
The result of all the above is that for parameters with limits, the error reported by Minuit is not exactly equal to the square root of the diagonal element of the error matrix. The difference is a measure of how much the limits deform the problem. If possible, it is suggested not to use limits on parameters, and the problem goes away. If for some reason limits are necessary, and you are sensitive to the difference between the two ways of calculating the errors, it is suggested to use Minos errors which take into account the non-linearities much more precisely.
Public Types | |
enum | { kMAXWARN =100 } |
Public Types inherited from TObject | |
enum | { kIsOnHeap = 0x01000000 , kNotDeleted = 0x02000000 , kZombie = 0x04000000 , kInconsistent = 0x08000000 , kBitMask = 0x00ffffff } |
enum | { kSingleKey = BIT(0) , kOverwrite = BIT(1) , kWriteDelete = BIT(2) } |
enum | EDeprecatedStatusBits { kObjInCanvas = BIT(3) } |
enum | EStatusBits { kCanDelete = BIT(0) , kMustCleanup = BIT(3) , kIsReferenced = BIT(4) , kHasUUID = BIT(5) , kCannotPick = BIT(6) , kNoContextMenu = BIT(8) , kInvalidObject = BIT(13) } |
Public Member Functions | |
TMinuit () | |
Minuit normal constructor. | |
TMinuit (Int_t maxpar) | |
Minuit normal constructor. | |
virtual | ~TMinuit () |
Minuit default destructor. | |
virtual void | BuildArrays (Int_t maxpar=15) |
Create internal Minuit arrays for the maxpar parameters. | |
virtual TObject * | Clone (const char *newname="") const |
Make a clone of an object using the Streamer facility. | |
virtual Int_t | Command (const char *command) |
Execute a Minuit command. | |
virtual TObject * | Contour (Int_t npoints=10, Int_t pa1=0, Int_t pa2=1) |
Creates a TGraph object describing the n-sigma contour of a TMinuit fit. | |
virtual Int_t | DefineParameter (Int_t parNo, const char *name, Double_t initVal, Double_t initErr, Double_t lowerLimit, Double_t upperLimit) |
Define a parameter. | |
virtual void | DeleteArrays () |
Delete internal Minuit arrays. | |
virtual Int_t | Eval (Int_t npar, Double_t *grad, Double_t &fval, Double_t *par, Int_t flag) |
Evaluate the minimisation function Input parameters: | |
virtual Int_t | FixParameter (Int_t parNo) |
fix a parameter | |
Int_t | GetMaxIterations () const |
TMethodCall * | GetMethodCall () const |
virtual Int_t | GetNumFixedPars () const |
returns the number of currently fixed parameters | |
virtual Int_t | GetNumFreePars () const |
returns the number of currently free parameters | |
virtual Int_t | GetNumPars () const |
returns the total number of parameters that have been defined as fixed or free. | |
TObject * | GetObjectFit () const |
virtual Int_t | GetParameter (Int_t parNo, Double_t ¤tValue, Double_t ¤tError) const |
return parameter value and error | |
virtual TObject * | GetPlot () const |
Int_t | GetStatus () const |
virtual Int_t | Migrad () |
invokes the MIGRAD minimizer | |
virtual void | mnamin () |
Initialize AMIN. | |
virtual void | mnbins (Double_t a1, Double_t a2, Int_t naa, Double_t &bl, Double_t &bh, Int_t &nb, Double_t &bwid) |
Compute reasonable histogram intervals. | |
virtual void | mncalf (Double_t *pvec, Double_t &ycalf) |
Transform FCN to find further minima. | |
virtual void | mncler () |
Resets the parameter list to UNDEFINED. | |
virtual void | mncntr (Int_t ke1, Int_t ke2, Int_t &ierrf) |
Print function contours in two variables, on line printer. | |
virtual void | mncomd (const char *crdbin, Int_t &icondn) |
Reads a command string and executes. | |
virtual void | mncont (Int_t ke1, Int_t ke2, Int_t nptu, Double_t *xptu, Double_t *yptu, Int_t &ierrf) |
Find points along a contour where FCN is minimum. | |
virtual void | mncrck (TString crdbuf, Int_t maxcwd, TString &comand, Int_t &lnc, Int_t mxp, Double_t *plist, Int_t &llist, Int_t &ierr, Int_t isyswr) |
Cracks the free-format input. | |
virtual void | mncros (Double_t &aopt, Int_t &iercr) |
Find point where MNEVAL=AMIN+UP. | |
virtual void | mncuve () |
Makes sure that the current point is a local minimum. | |
virtual void | mnderi () |
Calculates the first derivatives of FCN (GRD) | |
virtual void | mndxdi (Double_t pint, Int_t ipar, Double_t &dxdi) |
Calculates the transformation factor between ext/internal values. | |
virtual void | mneig (Double_t *a, Int_t ndima, Int_t n, Int_t mits, Double_t *work, Double_t precis, Int_t &ifault) |
Compute matrix eigen values. | |
virtual void | mnemat (Double_t *emat, Int_t ndim) |
Calculates the external error matrix from the internal matrix. | |
virtual void | mnerrs (Int_t number, Double_t &eplus, Double_t &eminus, Double_t &eparab, Double_t &gcc) |
Utility routine to get MINOS errors. | |
virtual void | mneval (Double_t anext, Double_t &fnext, Int_t &ierev) |
Evaluates the function being analysed by MNCROS. | |
virtual void | mnexcm (const char *comand, Double_t *plist, Int_t llist, Int_t &ierflg) |
Interprets a command and takes appropriate action. | |
virtual void | mnexin (Double_t *pint) |
Transforms the external parameter values U to internal values. | |
virtual void | mnfixp (Int_t iint, Int_t &ierr) |
Removes parameter IINT from the internal parameter list. | |
virtual void | mnfree (Int_t k) |
Restores one or more fixed parameter(s) to variable status. | |
virtual void | mngrad () |
Interprets the SET GRAD command. | |
virtual void | mnhelp (const char *command="") |
interface to Minuit help | |
virtual void | mnhelp (TString comd) |
HELP routine for MINUIT interactive commands. | |
virtual void | mnhes1 () |
Calculate first derivatives (GRD) and uncertainties (DGRD) | |
virtual void | mnhess () |
Calculates the full second-derivative matrix of FCN. | |
virtual void | mnimpr () |
Attempts to improve on a good local minimum. | |
virtual void | mninex (Double_t *pint) |
Transforms from internal coordinates (PINT) to external (U) | |
virtual void | mninit (Int_t i1, Int_t i2, Int_t i3) |
Main initialization member function for MINUIT. | |
virtual void | mnlims () |
Interprets the SET LIM command, to reset the parameter limits. | |
virtual void | mnline (Double_t *start, Double_t fstart, Double_t *step, Double_t slope, Double_t toler) |
Perform a line search from position START. | |
virtual void | mnmatu (Int_t kode) |
Prints the covariance matrix v when KODE=1. | |
virtual void | mnmigr () |
Performs a local function minimization. | |
virtual void | mnmnos () |
Performs a MINOS error analysis. | |
virtual void | mnmnot (Int_t ilax, Int_t ilax2, Double_t &val2pl, Double_t &val2mi) |
Performs a MINOS error analysis on one parameter. | |
virtual void | mnparm (Int_t k, TString cnamj, Double_t uk, Double_t wk, Double_t a, Double_t b, Int_t &ierflg) |
Implements one parameter definition. | |
virtual void | mnpars (TString &crdbuf, Int_t &icondn) |
Implements one parameter definition. | |
virtual void | mnpfit (Double_t *parx2p, Double_t *pary2p, Int_t npar2p, Double_t *coef2p, Double_t &sdev2p) |
To fit a parabola to npar2p points. | |
virtual void | mnpint (Double_t &pexti, Int_t i, Double_t &pinti) |
Calculates the internal parameter value PINTI. | |
virtual void | mnplot (Double_t *xpt, Double_t *ypt, char *chpt, Int_t nxypt, Int_t npagwd, Int_t npagln) |
Plots points in array xypt onto one page with labelled axes. | |
virtual void | mnpout (Int_t iuext, TString &chnam, Double_t &val, Double_t &err, Double_t &xlolim, Double_t &xuplim, Int_t &iuint) const |
Provides the user with information concerning the current status. | |
virtual void | mnprin (Int_t inkode, Double_t fval) |
Prints the values of the parameters at the time of the call. | |
virtual void | mnpsdf () |
Calculates the eigenvalues of v to see if positive-def. | |
virtual void | mnrazz (Double_t ynew, Double_t *pnew, Double_t *y, Int_t &jh, Int_t &jl) |
Called only by MNSIMP (and MNIMPR) to add a new point. | |
virtual void | mnrn15 (Double_t &val, Int_t &inseed) |
This is a super-portable random number generator. | |
virtual void | mnrset (Int_t iopt) |
Resets function value and errors to UNDEFINED. | |
virtual void | mnsave () |
Writes current parameter values and step sizes onto file ISYSSA. | |
virtual void | mnscan () |
Scans the values of FCN as a function of one parameter. | |
virtual void | mnseek () |
Performs a rough (but global) minimization by monte carlo search. | |
virtual void | mnset () |
Interprets the commands that start with SET and SHOW. | |
virtual void | mnsimp () |
Minimization using the simplex method of Nelder and Mead. | |
virtual void | mnstat (Double_t &fmin, Double_t &fedm, Double_t &errdef, Int_t &npari, Int_t &nparx, Int_t &istat) |
Returns concerning the current status of the minimization. | |
virtual void | mntiny (volatile Double_t epsp1, Double_t &epsbak) |
To find the machine precision. | |
Bool_t | mnunpt (TString &cfname) |
Returns .TRUE. | |
virtual void | mnvert (Double_t *a, Int_t l, Int_t m, Int_t n, Int_t &ifail) |
Inverts a symmetric matrix. | |
virtual void | mnwarn (const char *copt, const char *corg, const char *cmes) |
Prints Warning messages. | |
virtual void | mnwerr () |
Calculates the WERR, external parameter errors. | |
virtual Int_t | Release (Int_t parNo) |
release a parameter | |
virtual Int_t | SetErrorDef (Double_t up) |
To get the n-sigma contour the error def parameter "up" has to set to n^2. | |
virtual void | SetFCN (void(*fcn)(Int_t &, Double_t *, Double_t &f, Double_t *, Int_t)) |
To set the address of the minimization function. | |
virtual void | SetGraphicsMode (Bool_t mode=kTRUE) |
virtual void | SetMaxIterations (Int_t maxiter=500) |
virtual void | SetObjectFit (TObject *obj) |
virtual Int_t | SetPrintLevel (Int_t printLevel=0) |
set Minuit print level. | |
Public Member Functions inherited from TNamed | |
TNamed () | |
TNamed (const char *name, const char *title) | |
TNamed (const TNamed &named) | |
TNamed copy ctor. | |
TNamed (const TString &name, const TString &title) | |
virtual | ~TNamed () |
TNamed destructor. | |
virtual void | Clear (Option_t *option="") |
Set name and title to empty strings (""). | |
virtual Int_t | Compare (const TObject *obj) const |
Compare two TNamed objects. | |
virtual void | Copy (TObject &named) const |
Copy this to obj. | |
virtual void | FillBuffer (char *&buffer) |
Encode TNamed into output buffer. | |
virtual const char * | GetName () const |
Returns name of object. | |
virtual const char * | GetTitle () const |
Returns title of object. | |
virtual ULong_t | Hash () const |
Return hash value for this object. | |
virtual Bool_t | IsSortable () const |
virtual void | ls (Option_t *option="") const |
List TNamed name and title. | |
TNamed & | operator= (const TNamed &rhs) |
TNamed assignment operator. | |
virtual void | Print (Option_t *option="") const |
Print TNamed name and title. | |
virtual void | SetName (const char *name) |
Set the name of the TNamed. | |
virtual void | SetNameTitle (const char *name, const char *title) |
Set all the TNamed parameters (name and title). | |
virtual void | SetTitle (const char *title="") |
Set the title of the TNamed. | |
virtual Int_t | Sizeof () const |
Return size of the TNamed part of the TObject. | |
Public Member Functions inherited from TObject | |
TObject () | |
TObject constructor. | |
TObject (const TObject &object) | |
TObject copy ctor. | |
virtual | ~TObject () |
TObject destructor. | |
void | AbstractMethod (const char *method) const |
Use this method to implement an "abstract" method that you don't want to leave purely abstract. | |
virtual void | AppendPad (Option_t *option="") |
Append graphics object to current pad. | |
virtual void | Browse (TBrowser *b) |
Browse object. May be overridden for another default action. | |
ULong_t | CheckedHash () |
Check and record whether this class has a consistent Hash/RecursiveRemove setup (*) and then return the regular Hash value for this object. | |
virtual const char * | ClassName () const |
Returns name of class to which the object belongs. | |
virtual void | Delete (Option_t *option="") |
Delete this object. | |
virtual Int_t | DistancetoPrimitive (Int_t px, Int_t py) |
Computes distance from point (px,py) to the object. | |
virtual void | Draw (Option_t *option="") |
Default Draw method for all objects. | |
virtual void | DrawClass () const |
Draw class inheritance tree of the class to which this object belongs. | |
virtual TObject * | DrawClone (Option_t *option="") const |
Draw a clone of this object in the current selected pad for instance with: gROOT->SetSelectedPad(gPad) . | |
virtual void | Dump () const |
Dump contents of object on stdout. | |
virtual void | Error (const char *method, const char *msgfmt,...) const |
Issue error message. | |
virtual void | Execute (const char *method, const char *params, Int_t *error=0) |
Execute method on this object with the given parameter string, e.g. | |
virtual void | Execute (TMethod *method, TObjArray *params, Int_t *error=0) |
Execute method on this object with parameters stored in the TObjArray. | |
virtual void | ExecuteEvent (Int_t event, Int_t px, Int_t py) |
Execute action corresponding to an event at (px,py). | |
virtual void | Fatal (const char *method, const char *msgfmt,...) const |
Issue fatal error message. | |
virtual TObject * | FindObject (const char *name) const |
Must be redefined in derived classes. | |
virtual TObject * | FindObject (const TObject *obj) const |
Must be redefined in derived classes. | |
virtual Option_t * | GetDrawOption () const |
Get option used by the graphics system to draw this object. | |
virtual const char * | GetIconName () const |
Returns mime type name of object. | |
virtual char * | GetObjectInfo (Int_t px, Int_t py) const |
Returns string containing info about the object at position (px,py). | |
virtual Option_t * | GetOption () const |
virtual UInt_t | GetUniqueID () const |
Return the unique object id. | |
virtual Bool_t | HandleTimer (TTimer *timer) |
Execute action in response of a timer timing out. | |
Bool_t | HasInconsistentHash () const |
Return true is the type of this object is known to have an inconsistent setup for Hash and RecursiveRemove (i.e. | |
virtual void | Info (const char *method, const char *msgfmt,...) const |
Issue info message. | |
virtual Bool_t | InheritsFrom (const char *classname) const |
Returns kTRUE if object inherits from class "classname". | |
virtual Bool_t | InheritsFrom (const TClass *cl) const |
Returns kTRUE if object inherits from TClass cl. | |
virtual void | Inspect () const |
Dump contents of this object in a graphics canvas. | |
void | InvertBit (UInt_t f) |
virtual Bool_t | IsEqual (const TObject *obj) const |
Default equal comparison (objects are equal if they have the same address in memory). | |
virtual Bool_t | IsFolder () const |
Returns kTRUE in case object contains browsable objects (like containers or lists of other objects). | |
R__ALWAYS_INLINE Bool_t | IsOnHeap () const |
R__ALWAYS_INLINE Bool_t | IsZombie () const |
void | MayNotUse (const char *method) const |
Use this method to signal that a method (defined in a base class) may not be called in a derived class (in principle against good design since a child class should not provide less functionality than its parent, however, sometimes it is necessary). | |
virtual Bool_t | Notify () |
This method must be overridden to handle object notification. | |
void | Obsolete (const char *method, const char *asOfVers, const char *removedFromVers) const |
Use this method to declare a method obsolete. | |
void | operator delete (void *ptr) |
Operator delete. | |
void | operator delete[] (void *ptr) |
Operator delete []. | |
void * | operator new (size_t sz) |
void * | operator new (size_t sz, void *vp) |
void * | operator new[] (size_t sz) |
void * | operator new[] (size_t sz, void *vp) |
TObject & | operator= (const TObject &rhs) |
TObject assignment operator. | |
virtual void | Paint (Option_t *option="") |
This method must be overridden if a class wants to paint itself. | |
virtual void | Pop () |
Pop on object drawn in a pad to the top of the display list. | |
virtual Int_t | Read (const char *name) |
Read contents of object with specified name from the current directory. | |
virtual void | RecursiveRemove (TObject *obj) |
Recursively remove this object from a list. | |
void | ResetBit (UInt_t f) |
virtual void | SaveAs (const char *filename="", Option_t *option="") const |
Save this object in the file specified by filename. | |
virtual void | SavePrimitive (std::ostream &out, Option_t *option="") |
Save a primitive as a C++ statement(s) on output stream "out". | |
void | SetBit (UInt_t f) |
void | SetBit (UInt_t f, Bool_t set) |
Set or unset the user status bits as specified in f. | |
virtual void | SetDrawOption (Option_t *option="") |
Set drawing option for object. | |
virtual void | SetUniqueID (UInt_t uid) |
Set the unique object id. | |
virtual void | SysError (const char *method, const char *msgfmt,...) const |
Issue system error message. | |
R__ALWAYS_INLINE Bool_t | TestBit (UInt_t f) const |
Int_t | TestBits (UInt_t f) const |
virtual void | UseCurrentStyle () |
Set current style settings in this object This function is called when either TCanvas::UseCurrentStyle or TROOT::ForceStyle have been invoked. | |
virtual void | Warning (const char *method, const char *msgfmt,...) const |
Issue warning message. | |
virtual Int_t | Write (const char *name=0, Int_t option=0, Int_t bufsize=0) |
Write this object to the current directory. | |
virtual Int_t | Write (const char *name=0, Int_t option=0, Int_t bufsize=0) const |
Write this object to the current directory. | |
Private Member Functions | |
TMinuit (const TMinuit &m) | |
Private TMinuit copy ctor. TMinuit can not be copied. | |
TMinuit & | operator= (const TMinuit &m) |
Additional Inherited Members | |
Static Public Member Functions inherited from TObject | |
static Long_t | GetDtorOnly () |
Return destructor only flag. | |
static Bool_t | GetObjectStat () |
Get status of object stat flag. | |
static void | SetDtorOnly (void *obj) |
Set destructor only flag. | |
static void | SetObjectStat (Bool_t stat) |
Turn on/off tracking of objects in the TObjectTable. | |
Protected Types inherited from TObject | |
enum | { kOnlyPrepStep = BIT(3) } |
Protected Member Functions inherited from TObject | |
virtual void | DoError (int level, const char *location, const char *fmt, va_list va) const |
Interface to ErrorHandler (protected). | |
void | MakeZombie () |
Protected Attributes inherited from TNamed | |
TString | fName |
TString | fTitle |
#include <TMinuit.h>
|
private |
Private TMinuit copy ctor. TMinuit can not be copied.
Definition at line 488 of file TMinuit.cxx.
TMinuit::TMinuit | ( | ) |
Minuit normal constructor.
Definition at line 346 of file TMinuit.cxx.
TMinuit::TMinuit | ( | Int_t | maxpar | ) |
Minuit normal constructor.
maxpar is the maximum number of parameters used with this TMinuit object.
Definition at line 463 of file TMinuit.cxx.
|
virtual |
Minuit default destructor.
Definition at line 496 of file TMinuit.cxx.
Create internal Minuit arrays for the maxpar parameters.
Definition at line 511 of file TMinuit.cxx.
|
virtual |
Make a clone of an object using the Streamer facility.
Function pointer is copied to Clone
Reimplemented from TNamed.
Reimplemented in TMVA::MinuitWrapper.
Definition at line 595 of file TMinuit.cxx.
|
virtual |
Execute a Minuit command.
Equivalent to MNEXCM except that the command is given as a character string. See TMinuit::mnhelp for the full list of available commands See also the complete documentation of all the available commands
Returns the status of the execution:
Definition at line 625 of file TMinuit.cxx.
Creates a TGraph object describing the n-sigma contour of a TMinuit fit.
The contour of the parameters pa1 and pa2 is calculated using npoints (>=4) points. The TMinuit status will be
To get the n-sigma contour the ERRDEF parameter in Minuit has to set to n^2. The fcn function has to be set before the routine is called.
The TGraph object is created via the interpreter. The user must cast it to a TGraph*. Note that the TGraph is created with npoints+1 in order to close the contour (setting last point equal to first point).
You can find an example in $ROOTSYS/tutorials/fit/fitcont.C
Definition at line 652 of file TMinuit.cxx.
|
virtual |
Define a parameter.
Definition at line 694 of file TMinuit.cxx.
|
virtual |
Delete internal Minuit arrays.
Definition at line 707 of file TMinuit.cxx.
|
virtual |
Evaluate the minimisation function Input parameters:
The meaning of the parameters par is of course defined by the user, who uses the values of those parameters to calculate their function value. The starting values must be specified by the user. Later values are determined by Minuit as it searches for the minimum or performs whatever analysis is requested by the user.
Note that this virtual function may be redefined in a class derived from TMinuit. The default function calls the function specified in SetFCN
Example of Minimisation function:
Reimplemented in TMVA::MinuitWrapper.
Definition at line 799 of file TMinuit.cxx.
fix a parameter
Definition at line 826 of file TMinuit.cxx.
|
inline |
|
virtual |
returns the number of currently fixed parameters
Definition at line 854 of file TMinuit.cxx.
|
virtual |
returns the number of currently free parameters
Definition at line 862 of file TMinuit.cxx.
|
virtual |
returns the total number of parameters that have been defined as fixed or free.
The constant parameters are not counted.
Definition at line 871 of file TMinuit.cxx.
|
virtual |
return parameter value and error
Definition at line 840 of file TMinuit.cxx.
|
inlinevirtual |
|
virtual |
invokes the MIGRAD minimizer
Definition at line 879 of file TMinuit.cxx.
|
virtual |
Initialize AMIN.
Called from many places. Initializes the value of AMIN by calling the user function. Prints out the function value and parameter values if Print Flag value is high enough.
Definition at line 971 of file TMinuit.cxx.
|
virtual |
Compute reasonable histogram intervals.
Function TO DETERMINE REASONABLE HISTOGRAM INTERVALS GIVEN ABSOLUTE UPPER AND LOWER BOUNDS A1 AND A2 AND DESIRED MAXIMUM NUMBER OF BINS NAA PROGRAM MAKES REASONABLE BINNING FROM BL TO BH OF WIDTH BWID F. JAMES, AUGUST, 1974 , stolen for Minuit, 1988
Definition at line 996 of file TMinuit.cxx.
Transform FCN to find further minima.
Called only from MNIMPR. Transforms the function FCN by dividing out the quadratic part in order to find further minima. Calculates ycalf = (f-fmin)/(x-xmin)*v*(x-xmin)
Definition at line 1069 of file TMinuit.cxx.
|
virtual |
Resets the parameter list to UNDEFINED.
Called from MINUIT and by option from MNEXCM
Definition at line 1102 of file TMinuit.cxx.
Print function contours in two variables, on line printer.
input arguments: parx, pary, devs, ngrid
Definition at line 1131 of file TMinuit.cxx.
Reads a command string and executes.
Called by user. 'Reads' a command string and executes. Equivalent to MNEXCM except that the command is given as a character string. ICONDN = - 0: command executed normally - 1: command is blank, ignored - 2: command line unreadable, ignored - 3: unknown command, ignored - 4: abnormal termination (e.g., MIGRAD not converged) - 5: command is a request to read PARAMETER definitions - 6: 'SET INPUT' command - 7: 'SET TITLE' command - 8: 'SET COVAR' command - 9: reserved - 10: END command - 11: EXIT or STOP command - 12: RETURN command
Definition at line 1309 of file TMinuit.cxx.
|
virtual |
Find points along a contour where FCN is minimum.
Find NPTU points along a contour where the function
FMIN (X(KE1),X(KE2)) = AMIN+UP where FMIN is the minimum of FCN with respect to all the other NPAR-2 variable parameters (if any).
IERRF on return will be equal to the number of points found:
input arguments: parx, pary, devs, ngrid
Definition at line 1394 of file TMinuit.cxx.
|
virtual |
Cracks the free-format input.
Cracks the free-format input, expecting zero or more alphanumeric fields (which it joins into COMAND(1:LNC)) followed by one or more numeric fields separated by blanks and/or one comma. The numeric fields are put into the LLIST (but at most MXP) elements of PLIST. IERR : - = 0 if no errors, - = 1 if error(s).
Definition at line 1676 of file TMinuit.cxx.
Find point where MNEVAL=AMIN+UP.
Find point where MNEVAL=AMIN+UP, along the line through XMIDCR,YMIDCR with direction XDIRCR,YDIRCR, where X and Y are parameters KE1CR and KE2CR. If KE2CR=0 (from MINOS), only KE1CR is varied. From MNCONT, both are varied. Crossing point is at
(U(KE1),U(KE2)) = (XMID,YMID) + AOPT*(XDIR,YDIR)
Definition at line 1797 of file TMinuit.cxx.
|
virtual |
Makes sure that the current point is a local minimum.
Makes sure that the current point is a local minimum and that the error matrix exists, or at least something good enough for MINOS and MNCONT
Definition at line 2129 of file TMinuit.cxx.
|
virtual |
Calculates the first derivatives of FCN (GRD)
Calculates the first derivatives of FCN (GRD), either by finite differences or by transforming the user- supplied derivatives to internal coordinates, according to whether fISW[2] is zero or one.
Definition at line 2177 of file TMinuit.cxx.
Calculates the transformation factor between ext/internal values.
calculates the transformation factor between external and internal parameter values. this factor is one for parameters which are not limited. called from MNEMAT.
Definition at line 2292 of file TMinuit.cxx.
|
virtual |
Compute matrix eigen values.
Definition at line 2304 of file TMinuit.cxx.
Calculates the external error matrix from the internal matrix.
Note that if the matrix is declared like Double_t matrix[5][5] in the calling program, one has to call mnemat with, eg
gMinuit->mnemat(&matrix[0][0],5);
Definition at line 2500 of file TMinuit.cxx.
|
virtual |
Utility routine to get MINOS errors.
Called by user.
NUMBER is the parameter number
values returned by MNERRS:
Definition at line 2577 of file TMinuit.cxx.
Evaluates the function being analysed by MNCROS.
Evaluates the function being analysed by MNCROS, which is generally the minimum of FCN with respect to all remaining variable parameters. The class data members contains the data necessary to know the values of U(KE1CR) and U(KE2CR) to be used, namely U(KE1CR) = XMIDCR + ANEXT*XDIRCR and (if KE2CR .NE. 0) U(KE2CR) = YMIDCR + ANEXT*YDIRCR
Definition at line 2619 of file TMinuit.cxx.
|
virtual |
Interprets a command and takes appropriate action.
either directly by skipping to the corresponding code in MNEXCM, or by setting up a call to a function
recognized MINUIT commands: obsolete commands: IERFLG is now (94.5) defined the same as ICONDN in MNCOMD =
see also the possible list of all Minuit commands.
Definition at line 2663 of file TMinuit.cxx.
Transforms the external parameter values U to internal values.
Transforms the external parameter values U to internal values in the dense array PINT.
Definition at line 3150 of file TMinuit.cxx.
Removes parameter IINT from the internal parameter list.
and arranges the rest of the list to fill the hole.
Definition at line 3168 of file TMinuit.cxx.
Restores one or more fixed parameter(s) to variable status.
Restores one or more fixed parameter(s) to variable status by inserting it into the internal parameter list at the appropriate place.
- K = 0 means restore all parameters - K = 1 means restore the last parameter fixed - K = -I means restore external parameter I (if possible) - IQ = fix-location where internal parameters were stored - IR = external number of parameter being restored - IS = internal number of parameter being restored
Definition at line 3255 of file TMinuit.cxx.
|
virtual |
Interprets the SET GRAD command.
- Called from MNSET - Interprets the SET GRAD command, which informs MINUIT whether - the first derivatives of FCN will be calculated by the user - inside FCN. It can check the user derivative calculation - by comparing it with a finite difference approximation.
Definition at line 3361 of file TMinuit.cxx.
|
virtual |
interface to Minuit help
Definition at line 3421 of file TMinuit.cxx.
HELP routine for MINUIT interactive commands.
COMD =Command_name: print detailed help for one command. Note that at least 3 characters must be given for the command name.
Author: Rene Brun comments extracted from the MINUIT documentation file.
Definition at line 3438 of file TMinuit.cxx.
|
virtual |
Calculate first derivatives (GRD) and uncertainties (DGRD)
and appropriate step sizes GSTEP Called from MNHESS and MNGRAD
Definition at line 4217 of file TMinuit.cxx.
|
virtual |
Calculates the full second-derivative matrix of FCN.
by taking finite differences. When calculating diagonal elements, it may iterate so that step size is nearly that which gives function change= UP/10. The first derivatives of course come as a free side effect, but with a smaller step size in order to obtain a known accuracy.
Definition at line 3992 of file TMinuit.cxx.
|
virtual |
Attempts to improve on a good local minimum.
Attempts to improve on a good local minimum by finding a better one. The quadratic part of FCN is removed by MNCALF and this transformed function is minimised using the simplex method from several random starting points.
ref. – Goldstein and Price, Math.Comp. 25, 569 (1971)
Definition at line 4294 of file TMinuit.cxx.
Transforms from internal coordinates (PINT) to external (U)
The minimising routines which work in internal coordinates call this routine before calling FCN.
Definition at line 4505 of file TMinuit.cxx.
Main initialization member function for MINUIT.
It initializes some constants (including the logical I/O unit nos.),
Definition at line 4525 of file TMinuit.cxx.
|
virtual |
Interprets the SET LIM command, to reset the parameter limits.
Called from MNSET
Definition at line 4615 of file TMinuit.cxx.
|
virtual |
Perform a line search from position START.
along direction STEP, where the length of vector STEP gives the expected position of minimum.
SLAMBG and ALPHA control the maximum individual steps allowed. The first step is always =1. The max length of second step is SLAMBG. The max size of subsequent steps is the maximum previous successful step multiplied by ALPHA + the size of most recent successful step, but cannot be smaller than SLAMBG.
Definition at line 4735 of file TMinuit.cxx.
Prints the covariance matrix v when KODE=1.
always prints the global correlations, and calculates and prints the individual correlation coefficients
Definition at line 4969 of file TMinuit.cxx.
|
virtual |
Performs a local function minimization.
Performs a local function minimization using basically the method of Davidon-Fletcher-Powell as modified by Fletcher
ref. – Fletcher, Comp.J. 13,317 (1970) "switching method"
Definition at line 5049 of file TMinuit.cxx.
|
virtual |
Performs a MINOS error analysis.
Performs a MINOS error analysis on those parameters for which it is requested on the MINOS command by calling MNMNOT for each parameter requested.
Definition at line 5389 of file TMinuit.cxx.
Performs a MINOS error analysis on one parameter.
The parameter ILAX is varied, and the minimum of the function with respect to the other parameters is followed until it crosses the value FMIN+UP.
Definition at line 5464 of file TMinuit.cxx.
|
virtual |
Implements one parameter definition.
Called from MNPARS and user-callable Implements one parameter definition, that is:
Definition at line 5666 of file TMinuit.cxx.
Implements one parameter definition.
Called from MNREAD and user-callable Implements one parameter definition, that is: parses the string CRDBUF and calls MNPARM
output conditions:
Definition at line 5868 of file TMinuit.cxx.
|
virtual |
To fit a parabola to npar2p points.
Definition at line 5957 of file TMinuit.cxx.
Calculates the internal parameter value PINTI.
corresponding to the external value PEXTI for parameter I.
Definition at line 6018 of file TMinuit.cxx.
|
virtual |
Plots points in array xypt onto one page with labelled axes.
CHPT(I) = character to be plotted at this position the input point arrays XPT, YPT, CHPT are destroyed.
If fGraphicsmode is true (default), a TGraph object is produced via the Plug-in handler. To get the plot, you can do:
Definition at line 6069 of file TMinuit.cxx.
|
virtual |
Provides the user with information concerning the current status.
of parameter number IUEXT. Namely, it returns:
Note also: If IUEXT is negative, then it is -internal parameter number, and IUINT is returned as the EXTERNAL number. Except for IUINT, this is exactly the inverse of MNPARM User-called
Definition at line 6248 of file TMinuit.cxx.
Prints the values of the parameters at the time of the call.
also prints other relevant information such as function value, estimated distance to minimum, parameter errors, step sizes.
According to the value of IKODE, the printout is: IKODE=INKODE=
when INKODE=5, MNPRIN chooses IKODE=1,2, or 3, according to fISW[1]
Definition at line 6305 of file TMinuit.cxx.
|
virtual |
Calculates the eigenvalues of v to see if positive-def.
if not, adds constant along diagonal to make positive.
Definition at line 6495 of file TMinuit.cxx.
|
virtual |
Called only by MNSIMP (and MNIMPR) to add a new point.
and remove an old one from the current simplex, and get the estimated distance to minimum.
Definition at line 6569 of file TMinuit.cxx.
This is a super-portable random number generator.
It should not overflow on any 32-bit machine. The cycle is only ~10**9, so use with care! Note especially that VAL must not be undefined on input.
Set Default Starting Seed
Definition at line 6618 of file TMinuit.cxx.
Resets function value and errors to UNDEFINED.
Definition at line 6660 of file TMinuit.cxx.
|
virtual |
Writes current parameter values and step sizes onto file ISYSSA.
in format which can be reread by Minuit for restarting. The covariance matrix is also output if it exists.
Definition at line 6694 of file TMinuit.cxx.
|
virtual |
Scans the values of FCN as a function of one parameter.
and plots the resulting values as a curve using MNPLOT. It may be called to scan one parameter or all parameters. retains the best function and parameter values found.
Definition at line 6707 of file TMinuit.cxx.
|
virtual |
Performs a rough (but global) minimization by monte carlo search.
Each time a new minimum is found, the search area is shifted to be centered at the best value. Random points are chosen uniformly over a hypercube determined by current step sizes. The Metropolis algorithm accepts a worse point with probability exp(-d/UP), where d is the degradation. Improved points are of course always accepted. Actual steps are random multiples of the nominal steps (DIRIN).
Definition at line 6817 of file TMinuit.cxx.
|
virtual |
Interprets the commands that start with SET and SHOW.
Called from MNEXCM file characteristics for SET INPUT 'SET ' or 'SHOW', 'ON ' or 'OFF', 'SUPPRESSED' or 'REPORTED ' explanation of print level numbers -1:3 and strategies 0:2 identification of debug options things that can be set or shown options not intended for normal users
Definition at line 6912 of file TMinuit.cxx.
|
virtual |
Minimization using the simplex method of Nelder and Mead.
Performs a minimization using the simplex method of Nelder and Mead (ref. – Comp. J. 7,308 (1965)).
Definition at line 7430 of file TMinuit.cxx.
|
virtual |
Returns concerning the current status of the minimization.
User-called Namely, it returns:
Definition at line 7637 of file TMinuit.cxx.
To find the machine precision.
Compares its argument with the value 1.0, and returns the value .TRUE. if they are equal. To find EPSMAC safely by foiling the Fortran optimiser
Definition at line 7660 of file TMinuit.cxx.
Returns .TRUE.
if CFNAME contains unprintable characters
Definition at line 7669 of file TMinuit.cxx.
Inverts a symmetric matrix.
inverts a symmetric matrix. matrix is first scaled to have all ones on the diagonal (equivalent to change of units) but no pivoting is done since matrix is positive-definite.
Definition at line 7695 of file TMinuit.cxx.
|
virtual |
Prints Warning messages.
Definition at line 7783 of file TMinuit.cxx.
|
virtual |
Calculates the WERR, external parameter errors.
and the global correlation coefficients, to be called whenever a new covariance matrix is available.
Definition at line 7860 of file TMinuit.cxx.
release a parameter
Definition at line 893 of file TMinuit.cxx.
To get the n-sigma contour the error def parameter "up" has to set to n^2.
Definition at line 907 of file TMinuit.cxx.
To set the address of the minimization function.
Definition at line 919 of file TMinuit.cxx.
set Minuit print level.
printlevel:
Definition at line 951 of file TMinuit.cxx.
TString* TMinuit::fCpnam |
TMethodCall* TMinuit::fMethodCall |