LMU/musrfit
LMU/musrfit
5
0
Fork 0
Code Issues 6 Pull Requests Actions Packages Projects Releases 4 Activity

added Minuit2 Hesse

Browse Source
This commit is contained in:
nemu 2008-09-01 11:32:08 +00:00
parent 4f66c78e4a
commit b96ac14a18
2 changed files with 19 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

1
src/ToDo.txt
View File

@ -92,6 +92,7 @@ short term:
**DONE** 08-09-01 **DONE** 08-09-01
* PFitter.cpp: implement HESSE * PFitter.cpp: implement HESSE
**DONE** 08-09-01
* implement FFT with msr-interface * implement FFT with msr-interface

18
src/classes/PFitter.cpp
View File

@ -36,6 +36,7 @@ using namespace std;
#include <math.h> #include <math.h>
#include "Minuit2/FunctionMinimum.h" #include "Minuit2/FunctionMinimum.h"
#include "Minuit2/MnHesse.h"
#include "Minuit2/MnMinimize.h" #include "Minuit2/MnMinimize.h"
#include "Minuit2/MnMigrad.h" #include "Minuit2/MnMigrad.h"
#include "Minuit2/MnMinos.h" #include "Minuit2/MnMinos.h"
@ -311,6 +312,23 @@ bool PFitter::ExecuteHesse()
if (fFcnMin) if (fFcnMin)
fMnUserParams = fFcnMin->UserParameters(); fMnUserParams = fFcnMin->UserParameters();
// create the hesse object
ROOT::Minuit2::MnHesse hesse;
// call hesse
ROOT::Minuit2::MnUserParameterState mnState = hesse((*fFitterFcn), fMnUserParams);
if (!mnState.IsValid()) {
cout << endl << "**WARNING** PFitter::ExecuteHesse(): Hesse encountered some problems!";
cout << endl;
return false;
}
// fill parabolic errors
for (unsigned int i=0; i<fParams.size(); i++) {
fRunInfo->SetMsrParamStep(i, mnState.Error(i));
}
return true; return true;
} }