19 #include "TObjString.h"
30 virtual Int_t
Eval(Int_t npar, Double_t *grad, Double_t &fval, Double_t *pars, Int_t iflag);
38 TMinuit(NumPars), fxFitter(fitter) {
44 Int_t
TMinuitEx::Eval(Int_t npar, Double_t *grad, Double_t &fval, Double_t *pars, Int_t iflag) {
56 TGo4FitterAction(Name,
"Fitter minimization using TMinuit object"), fxCommands(), fxResults() {
71 return ((TObjString*)
fxCommands[n])->GetString().Data();
77 fMinuit.SetPrintLevel(-1);
79 for(Int_t n=0;n<Fitter->
NumPars();n++) {
82 Double_t epsilon = 0, RangeMin = 0, RangeMax = 0;
85 if (!Fitter->
GetParRange(FullName,RangeMin,RangeMax)) { RangeMin = 0; RangeMax = 0; }
87 fMinuit.mnparm(n, FullName, Fitter->
GetParValue(FullName),
88 epsilon, RangeMin, RangeMax, ierflg);
90 if (Fitter->
GetParFixed(FullName)) fMinuit.FixParameter(n);
93 if (
fxCommands.GetLast()<0) fMinuit.Command(
"MIGRAD 500 1");
96 TString cmd ( ((TObjString*)
fxCommands[n])->GetString() );
98 if (cmd.Index(
"result",6,0,TString::kIgnoreCase) == 0) {
100 while ((cmd.Length()>0) && (cmd[0]==
' ')) cmd.Remove(0,1);
101 if (cmd.Length()==0) cmd =
"1000";
102 if (cmd.Length()<4) { std::cerr <<
"invalid result command syntax" << std::endl;
break; }
103 Bool_t getpar = (cmd[0]==
'1');
104 Bool_t geterr = (cmd[1]==
'1');
105 Bool_t getmatr = (cmd[2]==
'1');
106 Bool_t getcontr = (cmd[3]==
'1');
108 while ((cmd.Length()>0) && (cmd[0]==
' ')) cmd.Remove(0,1);
109 if (cmd.Length()==0) cmd=
"Result";
120 }
else fMinuit.Command(cmd);
123 for(Int_t n=0;n<Fitter->
NumPars();n++) {
124 Double_t value,error;
125 fMinuit.GetParameter(n, value, error);
156 std::cout <<
"List of commands:" << std::endl;
158 std::cout <<
" " << ((TObjString*)
fxCommands[n])->String().Data() << std::endl;
160 std::cout <<
"List of stored results:" << std::endl;
Bool_t SetParValue(const char *ParName, Double_t iValue)
virtual Bool_t GetParEpsilon(const char *ParName, Double_t &Epsilon)
void GetContourPlot(TMinuit *fMinuit)
TGo4FitterAbstract * fxFitter
void RemoveResult(TGo4FitMinuitResult *res)
TMinuitEx(Int_t NumPars, TGo4FitterAbstract *fitter)
void AddCommand(const char *iCommand)
virtual Bool_t GetParFixed(const char *ParName)
void CallMNERRS(TMinuit *fMinuit, Int_t nPars)
const char * GetCommand(Int_t n)
void CallMNEMAT(TMinuit *fMinuit, Int_t nPars, Bool_t DoTransform=kTRUE)
Bool_t SetParError(const char *ParName, Double_t iError)
void Print(Option_t *option) const
const char * GetParFullName(Int_t n)
TGo4FitMinuitResult * GetResult(Int_t indx)
Double_t CalculateFitFunction(Double_t *pars=0)
virtual void DoAction(TGo4FitterAbstract *Fitter)
virtual Bool_t GetParRange(const char *ParName, Double_t &RangeMin, Double_t &RangeMax)
virtual Int_t Eval(Int_t npar, Double_t *grad, Double_t &fval, Double_t *pars, Int_t iflag)
virtual void Print(Option_t *option) const
TGo4FitMinuitResult * FindResult(const char *ResName)
void CallMNSTAT(TMinuit *fMinuit)
void AddResult(TGo4FitMinuitResult *res)
void CallMNPOUT(TMinuit *fMinuit, Int_t nPars)
Double_t GetParValue(const char *ParName)