Re: [ROOT] [root] TMinuit Arithmetic exception.

From: Rene Brun (Rene.Brun@cern.ch)
Date: Tue Jan 09 2001 - 17:51:00 MET


Hi Dimitri,
If you change the problem in the same session, you must reset the fitter.
If you use TMinuit directly, you should call TMinuit::mncler.
If you use the TVirtualFitter interface (eg TFitter) use TFitter::Clear
See an example in TH1::Fit


Rene Brun

Dimitri Peressounko wrote:
> 
>    Dear Rene,
> 
> I meneged to reproduce this error with root
> compiled in debug mode, and obtain the following bt:
> 
> Program received signal SIGFPE, Arithmetic exception.
> 0x434342cb in TMath::Abs () at include/TUnixSystem.h:173
> 173     include/TUnixSystem.h: No such file or directory.
> (gdb) bt
> #0  0x434342cb in TMath::Abs () at include/TUnixSystem.h:173
> #1  0x44099e08 in TMinuit::mnline (this=0x97d55d8, start=0x440d2ce0,
> fstart=0.031689250648864867,
>     step=0x440d1980, slope=-1.757748811030452e-07,
> toler=0.050000000000000003) at minuit/src/TMinuit.cxx:4530
> #2  0x4409c1ab in TMinuit::mnmigr (this=0x97d55d8) at
> minuit/src/TMinuit.cxx:4959
> #3  0x44091127 in TMinuit::mnexcm (this=0x97d55d8, command=0x41c89ae3
> "MIGRAD", plist=0xbffe3bd8, llist=0,
>     ierflg=@0xbffe3bf4) at minuit/src/TMinuit.cxx:2570
> #4  0x41c081b9 in AliPHOSTrackSegmentMakerv1::FindFit (this=0x97d55b0,
> emcRP=0xa247400, maxAt=0x9ce2e00,
>     maxAtEnergy=0x9ce2e18, nPar=6, fitparameters=0xa0dde98) at
> AliPHOSTrackSegmentMakerv1.cxx:130
> .............
> 
> This takes place because variable fNpar (the number of free parameters)
> is not reseted each new call of minuit, and the resting positions of
> arrais are filled with some garbage...
> 
> So, could you add something like TMinuit::Reset() ?
> or check of the number of free parameters before minimization?
> 
>            Sincerely yours, Dmitri.
> 
> On Mon, 8 Jan 2001, Rene Brun wrote:
> 
> > Hi Dimitri,
> > Could you reproduce this problem with a stand alone script?
> > If you can, send me the file(s)
> >
> > Rene Brun
> >



This archive was generated by hypermail 2b29 : Tue Jan 01 2002 - 17:50:33 MET