[ROOT] Re: [ROOTDEV] How to get TMinuit's fit-results??

From: Rene Brun (Rene.Brun@cern.ch)
Date: Wed Jun 23 2004 - 18:47:08 MEST


Hi Christian,

At the end of your script, do:

    //Try to get the result (alpha1 and alpha2) from minuit 
    Double_t perrors[6], params[6];
    gMinuit->GetParameter(0,params[0], perrors[0]);
    gMinuit->GetParameter(1,params[1], perrors[1]);
    cout<<"params[0] "<<params[0]<<" params[1] "<<params[1]<<endl;


Rene Brun

Christian Hansen wrote:
> 
> Dear Rooters!
> 
> In the attached example I'm doing a fit, and the
> result is printed out correctly:
> 
> NO.   NAME      VALUE            ERROR          SIZE      DERIVATIVE
>  1  alpha1       9.89932e-04   3.87311e-05   3.84553e-09  -1.15652e+00
>  2  alpha2      -5.89321e+05   2.28881e+04   2.28035e+00   1.96090e-09
> 
> But how do I get hold of the results alpha1 and alpha2
> in variables???
> I've made a try in line 64 in the attached file to use
> gMinuit->GetParameters(&params[3]);
> but it does not work.
> 
> How should I do instead?
> --
> / Christian Hansen
> 
>   --------------------------------------------------------------------------------
>                  Name: simpleEx.C
>    simpleEx.C    Type: Plain Text (TEXT/PLAIN)
>              Encoding: BASE64



This archive was generated by hypermail 2b29 : Sun Jan 02 2005 - 05:50:08 MET