next up previous contents index
Next: How to get Up: Basic concepts of Previous: Basic concepts of

Basic concepts - The transformation for parameters with limits.

For variable parameters with limits, MINUIT uses the following transformation:

Pint= arcsin( 2 {Pext-ab-a}- 1 ) Pext= a + {b - a2}( sinPint+ 1 )

so that the internal value Pint can take on any value, while the external value Pext 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,1010) 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 Pext/ Pint≈0 . Therefore, it is recommended that:


next up previous contents index
Next: How to get Up: Basic concepts of Previous: Basic concepts of

Janne Saarela
Tue May 16 09:09:27 METDST 1995