Hi Dimitri,
Could you reproduce this problem with a stand alone script?
If you can, send me the file(s)
Rene Brun
Dimitri Peressounko wrote:
>
> Hi Rooters,
> I use TMinuit in my code. Normally everething is OK, but sometimes
> I recieve Arithmetic exception:
>
> Program received signal SIGFPE, Arithmetic exception.
> 0x43d817c3 in TMinuit::mnline ()
> from
> /usr/remote/prive/alice/root/root_v2.25-03.Linux.2.2.14.gcc2952/lib/libMinuit.so
> Current language: auto; currently c
> (gdb) bt
> #0 0x43d817c3 in TMinuit::mnline ()
> from
> /usr/remote/prive/alice/root/root_v2.25-03.Linux.2.2.14.gcc2952/lib/libMinuit.so
> #1 0x43d83616 in TMinuit::mnmigr ()
> from
> /usr/remote/prive/alice/root/root_v2.25-03.Linux.2.2.14.gcc2952/lib/libMinuit.so
> #2 0x43d7ad9c in TMinuit::mnexcm ()
> from
> /usr/remote/prive/alice/root/root_v2.25-03.Linux.2.2.14.gcc2952/lib/libMinuit.so
> #3 0x41b8ea45 in AliPHOSTrackSegmentMakerv1::FindFit (this=0x97d8d70,
> emcRP=0xa210028, maxAt=0x9b38b50, maxAtEnergy=0x9b38b88, nPar=6,
> fitparameters=0xa2c43f0) at AliPHOSTrackSegmentMakerv1.cxx:132
>
> Does someone knows, what is going wrong (or how to check this): TMinuit or
> my program? The data, which should be fitted, looks OK...
>
> Sincerely Dmitri.
>
> Below one can find my code:
> --------------------------------------------------------------------
> AliPHOSTrackSegmentMakerv1::AliPHOSTrackSegmentMakerv1() :
> AliPHOSTrackSegmentMaker()
> {
> // ctor
> if(!gMinuit) gMinuit = new TMinuit(100) ;
> }
>
> //____________________________________________________________________________
> AliPHOSTrackSegmentMakerv1::~AliPHOSTrackSegmentMakerv1()
> {
> // dtor
> delete gMinuit ;
> gMinuit = 0 ;
> }
> //____________________________________________________________________________
> Bool_t AliPHOSTrackSegmentMakerv1::FindFit(AliPHOSEmcRecPoint * emcRP,
> int * maxAt, Float_t * maxAtEnergy,
> Int_t nPar, Float_t * fitparameters)
> {
> // Calls TMinuit to fit the energy distribution of a cluster with
> several maxima
>
> gMinuit->SetPrintLevel(-1) ; // No Printout
> gMinuit->SetFCN(UnfoldingChiSquare) ; // To set the address of the
> //minimization function
> gMinuit->SetObjectFit(emcRP) ; // To tranfer pointer to
> //UnfoldingChiSquare
> for(all parameters)
> gMinuit->mnparm(index, "x", x, 0.1, 0, 0, ierflg) ;
>
> gMinuit->mnexcm("SET STR", &p2, 0, ierflg) ; // force TMinuit to
> // reduce function calls
> gMinuit->mnexcm("SET GRA", &p1, 1, ierflg) ; // force TMinuit to use
> // my gradient
> gMinuit->SetMaxIterations(5);
> gMinuit->mnexcm("SET NOW", &p2 , 0, ierflg) ; // No Warnings
>
> gMinuit->mnexcm("MIGRAD", &p0, 0, ierflg) ; // minimize
>
> }
> //______________________________________________________________________________
> void UnfoldingChiSquare(Int_t & nPar, Double_t * Grad, Double_t & fret,
> Double_t * x, Int_t iflag)
> {
> // Calculates th Chi square for the cluster unfolding minimization
> // Number of parameters, Gradient, Chi squared, parameters, what to do
>
> AliPHOSEmcRecPoint * emcRP = (AliPHOSEmcRecPoint *)gMinuit->GetObjectFit() ;
> // EmcRecPoint to fit
> Grad = ...
> fret = ...
>
> }
This archive was generated by hypermail 2b29 : Tue Jan 01 2002 - 17:50:33 MET