Re: [ROOT] Interactive minuit access from ROOT? -FixParameter functionmissing?

From: Arthur E. Snyder (snyder@SLAC.stanford.edu)
Date: Mon Dec 03 2001 - 04:58:15 MET


Also, I noticed that the implementation (given in the documentation) of
TF1::FixParameter seems to set the upper and lower limits equal using
SetParLimits. This would normally correspond to 'constant' in minuit not
'fixed'. If I recall correctly this procedure would not allow me to step a
fixed parameter and refit (as in mapping out a likelihood function) as the
value after a step would be outside the allowed limits. I guess this could
be dodged by using SetParLimits as well as well as SetParameter when making
such a step.

Art S.


----- Original Message -----
From: "Arthur E. Snyder" <snyder@SLAC.Stanford.EDU>
To: "Rene Brun" <brun@pcbrun.cern.ch>
Cc: <roottalk@pcroot.cern.ch>; "Arthur E. Snyder"
<snyder@MAILBOX.SLAC.Stanford.EDU>
Sent: Sunday, December 02, 2001 7:44 PM
Subject: Re: [ROOT] Interactive minuit access from ROOT? -FixParameter
functionmissing?


> Hello Rene,
>
> I'm trying to use the function TF1::FixParameter(Int_t,Double_t). It does
> not seem to be there. I see it in the write up, but when I try to execute
it
> complains that it does not exist. Also <tab> completion does not find it,
> but does find other functions, such as GetParameter.
>
> ?
>
> Thanks, Art S.
>
> ----- Original Message -----
> From: "Rene Brun" <brun@pcbrun.cern.ch>
> To: "Arthur E. Snyder" <snyder@SLAC.Stanford.EDU>
> Cc: <roottalk@pcroot.cern.ch>
> Sent: Sunday, December 02, 2001 2:49 PM
> Subject: Re: [ROOT] Interactive minuit access from ROOT?
>
>
> > Art,
> >
> > Even in PAW, people included the Minuit commands in a KUIP macro.
> > You should do the same. However, currently, it is not easy to get
> > the control after each iteration. You get the control only at the end
> > of the fit, unless you implement your own fitting model (like
> > H1FitChisquare or likelihood). One could imagine a new option
> > to force a call to a function or macro after each iteration.
> > If there are more suggestions, they are welcome.
> >
> > Rene Brun
> >
> > On Sun, 2 Dec 2001, Arthur E. Snyder wrote:
> >
> > > Hi Rene,
> > >
> > > I was able to grab gMinuit and use mnexcm to execute minuit commands,
> refit,
> > > etc. It is pretty awkward interactively, but works. However, the
results
> do
> > > not propagate back to the function attached to my histogram, so when I
> > > redraw it I see the result of the original default fit. Also, if I
refit
> it
> > > seems to start again from scratch and forget any 'fix' or 'set par'
I've
> > > done to gMinuit object.
> > >
> > > Perhaps, I would be better off using the get/set functions provided in
> TF1?
> > > It looks pretty easy to write my own little macro to manipulate these
> > > functions with minimal typing ...
> > >
> > > Art
> > >
> > >
> > > ----- Original Message -----
> > > From: "Rene Brun" <brun@pcbrun.cern.ch>
> > > To: "Arthur E. Snyder" <snyder@SLAC.Stanford.EDU>
> > > Cc: <roottalk@pcroot.cern.ch>
> > > Sent: Sunday, December 02, 2001 1:39 PM
> > > Subject: Re: [ROOT] Interactive minuit access from ROOT?
> > >
> > >
> > > > Use the global gMinuit
> > > >
> > > > Rene Brun
> > > >
> > > > On Sun, 2 Dec 2001, Arthur E. Snyder wrote:
> > > >
> > > > > Thanks, Rene. Can you tell me how I access TMinuit object when I'm
> > > fitting a
> > > > > histogram? I don't find an assessor either in TF1 or TH1. I've
also
> > > looked
> > > > > up the inheritance tree but don't find anything that looks like
the
> > > right
> > > > > thing.
> > > > >
> > > > > Art
> > > > >
> > > > >
> > > > >
> > > > > ----- Original Message -----
> > > > > From: "Rene Brun" <brun@pcbrun.cern.ch>
> > > > > To: "Arthur E. Snyder" <snyder@SLAC.Stanford.EDU>
> > > > > Cc: <roottalk@pcroot.cern.ch>
> > > > > Sent: Sunday, December 02, 2001 9:49 AM
> > > > > Subject: Re: [ROOT] Interactive minuit access from ROOT?
> > > > >
> > > > >
> > > > > > Hi Art,
> > > > > >
> > > > > > The direct equivalent of PAW option "M" does not exist.
> > > > > > However, you can invoke any function from TMinuit directly.
> > > > > > This should provide even more flexibility at the expense
> > > > > > of slightly more typing.
> > > > > >
> > > > > > Rene Brun
> > > > > >
> > > > > > On Sat, 1 Dec 2001, Arthur E. Snyder wrote:
> > > > > >
> > > > > > > Try again to post this with more sensible subject line!
> > > > > > >
> > > > > > > ----- Original Message -----
> > > > > > > From: "Arthur E. Snyder" <snyder@SLAC.Stanford.EDU>
> > > > > > > To: <roottalk@pcroot.cern.ch>
> > > > > > > Sent: Saturday, December 01, 2001 12:29 AM
> > > > > > > Subject: Re: [ROOT] tiny Bug in root
> > > > > > >
> > > > > > >
> > > > > > > > Is it possible to invoke interactive minuit from within
ROOT?
> > > I.e.,
> > > > > the
> > > > > > > > equivalent of 'h/fit 100 myfunc.f m' in paw? The ROOT fit
> provides
> > > > > most of
> > > > > > > > the options of paw, but not this one. The 'm' option in ROOT
> seems
> > > to
> > > > > be
> > > > > > > > something different.
> > > > > > > >
> > > > > > > > I've looked in TMinuit, TF, and TH1 classes, but have not
been
> > > able to
> > > > > > > find
> > > > > > > > anyway to get to the interactive interface of minuit or some
> > > > > alternative.
> > > > > > > >
> > > > > > > > Art Snyder at SLAC
> > > > > > > >
> > > > > > > >
> > > > > > > >
> > > > > > >
> > > > > > >
> > > > >
> > >



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