Total Initial distance square = 8.65172e+07
Total final distance square 19988
****************************************
Minimizer is Minuit2 / Migrad
MinFCN = 19988
NDf = 0
Edm = 1.04086e-06
NCalls = 248
Par_0 = 10.3555 +/- 0.392118
Par_1 = 19.9602 +/- 0.0688151
Par_2 = 0.999233 +/- 0.0444791
Par_3 = 2.00379 +/- 0.00783849
(int) 0
#include <cassert>
void line(
double t,
const double *
p,
double &
x,
double &
y,
double &z) {
z = t;
}
bool first = true;
struct SumDistance2 {
double distance2(
double x,
double y,
double z,
const double *
p) {
return d2;
}
assert(fGraph != nullptr);
double *
x = fGraph->
GetX();
double *
y = fGraph->
GetY();
double * z = fGraph->
GetZ();
int npoints = fGraph->
GetN();
for (int i = 0; i < npoints; ++i) {
double d = distance2(
x[i],
y[i],z[i],par);
}
if (first) {
std::cout <<
"Total Initial distance square = " <<
sum << std::endl;
}
first = false;
}
};
int line3Dfit()
{
int nd = 10000;
double p0[4] = {10,20,1,2};
for (
int N=0;
N<nd;
N++) {
double err = 1;
}
double pStart[4] = {1,1,1,1};
if (!ok) {
Error(
"line3Dfit",
"Line3D Fit failed");
return 1;
}
std::cout <<
"Total final distance square " <<
result.MinFcnValue() << std::endl;
const double * parFit =
result.GetParams();
double t0 = 0;
double dt = 10;
for (
int i = 0; i <
n;++i) {
}
for (
int i = 0; i <
n;++i) {
}
return 0;
}
void Error(const char *location, const char *msgfmt,...)
Use this function in case an error occurred.
TGLVector3 Cross(const TGLVector3 &v1, const TGLVector3 &v2)
winID h TVirtualViewer3D TVirtualGLPainter p
Option_t Option_t TPoint TPoint const char GetTextMagnitude GetFillStyle GetLineColor GetLineWidth GetMarkerStyle GetTextAlign GetTextColor GetTextSize void char Point_t Rectangle_t WindowAttributes_t Float_t Float_t Float_t Int_t Int_t UInt_t UInt_t Rectangle_t result
Option_t Option_t TPoint TPoint const char x1
TRObject operator()(const T1 &t1) const
R__EXTERN TRandom * gRandom
R__EXTERN TStyle * gStyle
const ParameterSettings & ParSettings(unsigned int i) const
get the parameter settings for the i-th parameter (const method)
class containing the result of the fit and all the related information (fitted parameter values,...
Fitter class, entry point for performing all type of fits.
const FitResult & Result() const
get fit result
bool FitFCN(unsigned int npar, Function &fcn, const double *params=nullptr, unsigned int dataSize=0, int fitType=0)
Fit using the a generic FCN function as a C++ callable object implementing double () (const double *)...
const FitConfig & Config() const
access to the fit configuration (const method)
bool SetFCN(unsigned int npar, Function &fcn, const double *params=nullptr, unsigned int dataSize=0, int fitType=0)
Set a generic FCN function as a C++ callable object implementing double () (const double *) Note that...
void SetStepSize(double err)
set the step size
Documentation for class Functor class.
virtual void SetLineColor(Color_t lcolor)
Set the line color.
Graphics object made of three arrays X, Y and Z with the same number of points each.
virtual void SetPoint(Int_t i, Double_t x, Double_t y)
Set x and y values for point number i.
void Draw(Option_t *chopt="") override
Draw this graph with its current attributes.
virtual void Draw(Option_t *option="")
Default Draw method for all objects.
A 3-dimensional polyline.
virtual void SetPoint(Int_t point, Double_t x, Double_t y, Double_t z)
Set point n to x, y, z.
void Draw(Option_t *option="") override
Draw this 3-D polyline with its current attributes.
virtual Double_t Gaus(Double_t mean=0, Double_t sigma=1)
Samples a random number from the standard Normal (Gaussian) Distribution with the given mean and sigm...
virtual Double_t Uniform(Double_t x1=1)
Returns a uniform deviate on the interval (0, x1).
void SetOptStat(Int_t stat=1)
The type of information printed in the histogram statistics box can be selected via the parameter mod...
void SetOptFit(Int_t fit=1)
The type of information about fit parameters printed in the histogram statistics box can be selected ...
SVector< T, D > Unit(const SVector< T, D > &rhs)
Unit.
T Mag2(const SVector< T, D > &rhs)
Vector magnitude square Template to compute .
static uint64_t sum(uint64_t i)