Inheritance diagram for BaseFitObject:

Public Member Functions
	BaseFitObject ()
	Default constructor.

	BaseFitObject (const BaseFitObject &rhs)
	Copy constructor.

BaseFitObject &	operator= (const BaseFitObject &rhs)
	Assignment.

virtual	~BaseFitObject ()
	Virtual destructor.

virtual BaseFitObject *	copy () const =0
	Return a new copy of itself.

virtual BaseFitObject &	assign (const BaseFitObject &source)
	Assign from anther object, if of same type.

virtual bool	setParam (int ilocal, double par_, bool measured_, bool fixed_=false)
	Set value and measured flag of parameter i; return: significant change.

virtual bool	setParam (int ilocal, double par_)
	Set value of parameter ilocal; return: significant change.

virtual bool	updateParams (double p[], int idim)
	Read values from global vector, readjust vector; return: significant change.

virtual bool	setMParam (int ilocal, double mpar_)
	Set measured value of parameter ilocal; return: success.

virtual bool	setError (int ilocal, double err_)
	Set error of parameter ilocal; return: success.

virtual bool	setCov (int ilocal, int jlocal, double cov_)
	Set covariance of parameters ilocal and jlocal; return: success.

virtual bool	setGlobalParNum (int ilocal, int iglobal)
	Set number of parameter ilocal in global list return true signals OK.

virtual bool	fixParam (int ilocal, bool fix=true)
	Fix a parameter (fix=true), or release it (fix=false)

virtual bool	releaseParam (int ilocal)
	Release a parameter.

virtual bool	isParamFixed (int ilocal) const
	Returns whether parameter is fixed.

virtual double	getParam (int ilocal) const
	Get current value of parameter ilocal.

virtual const char *	getParamName (int ilocal) const =0
	Get name of parameter ilocal.

virtual const char *	getName () const
	Get object's name.

void	setName (const char *name_)
	Set object's name.

virtual double	getMParam (int ilocal) const
	Get measured value of parameter ilocal.

virtual double	getError (int ilocal) const
	Get error of parameter ilocal.

virtual double	getCov (int ilocal, int jlocal) const
	Get covariance between parameters ilocal and jlocal.

virtual double	getRho (int ilocal, int jlocal) const
	Get correlation coefficient between parameters ilocal and jlocal.

virtual bool	isParamMeasured (int ilocal) const
	Get measured flag for parameter ilocal.

virtual int	getGlobalParNum (int ilocal) const
	Get global parameter number of parameter ilocal.

virtual int	getNPar () const =0
	Get total number of parameters of this FitObject.

virtual int	getNMeasured () const
	Get number of measured parameters of this FitObject.

virtual int	getNUnmeasured () const
	Get number of unmeasured parameters of this FitObject.

virtual int	getNFree () const
	Get number of free parameters of this FitObject.

virtual int	getNFixed () const
	Get number of fixed parameters of this FitObject.

virtual double	getChi2 () const
	Get chi squared from measured and fitted parameters.

virtual double	getDChi2DParam (int ilocal) const
	Get derivative of chi squared w.r.t. parameter ilocal.

virtual double	getD2Chi2DParam2 (int ilocal, int jlocal) const
	Get second derivative of chi squared w.r.t. parameters ilocal1 and ilocal2.

virtual std::ostream &	printParams (std::ostream &os) const
	print the parameters and errors

virtual std::ostream &	printRhoValues (std::ostream &os) const
	print the correlation coefficients

virtual std::ostream &	print1stDerivatives (std::ostream &os) const
	print the 1st derivatives wrt metaSet 0 (E, px, py, pz)

virtual std::ostream &	print2ndDerivatives (std::ostream &os) const
	print the 2nd derivatives wrt metaSet 0 (E, px, py, pz)

virtual std::ostream &	print (std::ostream &os) const =0
	print object to ostream

void	invalidateCache () const
	invalidate any cached quantities

virtual void	updateCache () const =0

virtual void	addToGlobCov (double *glcov, int idim) const
	Add covariance matrix elements to global covariance matrix of size idim x idim.

virtual int	addToGlobalChi2DerVector (double *y, int idim) const
	Add derivatives of chi squared to global derivative vector.

virtual void	addToGlobalChi2DerMatrix (double *M, int idim) const
	Add 2nd derivatives of chi squared to global derivative matrix.

virtual void	addToGlobalChi2DerVector (double *y, int idim, double lambda, double der[], int metaSet) const
	Add derivatives of momentum vector to global derivative vector.

virtual void	addTo1stDerivatives (double M[], int idim, double der[], int kglobal, int metaSet) const

virtual void	addTo2ndDerivatives (double der2[], int idim, double factor[], int metaSet) const

virtual void	addTo2ndDerivatives (double M[], int idim, double lambda, double der[], int metaSet) const

virtual double	getFirstDerivative_Meta_Local (int iMeta, int ilocal, int metaSet) const =0

virtual double	getSecondDerivative_Meta_Local (int iMeta, int ilocal, int jlocal, int metaSet) const =0

virtual void	initCov ()

virtual double	getError2 (double der[], int metaset) const

virtual void	getDerivatives (double der[], int idim) const =0

Protected Member Functions
virtual bool	calculateCovInv () const
	Calculate the inverse of the covariance matrix.

Protected Attributes
char *	name

double	par [BaseDefs::MAXPAR]
	fit parameters

double	mpar [BaseDefs::MAXPAR]
	measured parameters

bool	measured [BaseDefs::MAXPAR]
	measured flag

bool	fixed [BaseDefs::MAXPAR]
	fixed flag

int	globalParNum [BaseDefs::MAXPAR]
	global parameter number for each parameter

double	cov [BaseDefs::MAXPAR][BaseDefs::MAXPAR]
	local covariance matrix

double	covinv [BaseDefs::MAXPAR][BaseDefs::MAXPAR]
	inverse pf local covariance matrix

bool	covinvvalid
	flag for valid inverse covariance matrix

bool	cachevalid
	flag for valid cache

Static Protected Attributes
static const double	eps2 = 0.0001

Related Functions
(Note that these are not member functions.)
std::ostream &	operator<< (std::ostream &os, const BaseFitObject &bfo)
	Prints out a BaseFitObject, using its print method.

Detailed Description

Definition at line 91 of file BaseFitObject.h.

Constructor & Destructor Documentation

◆ BaseFitObject() [1/2]

BaseFitObject ( )

Default constructor.

Definition at line 38 of file BaseFitObject.cc.

                                : name(nullptr), par{}, mpar{}, measured{}, fixed{},  globalParNum{}, cov{}, covinv{},  covinvvalid(
        false),
      cachevalid(false)
    {
      setName("???");
      invalidateCache();
 
      for (int ilocal = 0; ilocal < BaseDefs::MAXPAR; ++ilocal) {
        globalParNum[ilocal] = -1;
        fixed[ilocal] = false;
        for (int jlocal = 0; jlocal < BaseDefs::MAXPAR; ++jlocal)
          cov[ilocal][jlocal] = 0;
      }
 
    }

◆ BaseFitObject() [2/2]

BaseFitObject ( const BaseFitObject & rhs )

Copy constructor.

Parameters

rhs	right hand side

Definition at line 54 of file BaseFitObject.cc.

      : name(nullptr), par{}, mpar{}, measured{}, fixed{},  globalParNum{}, cov{}, covinv{}, covinvvalid(false), cachevalid(false)
    {
      BaseFitObject::assign(rhs);
    }

◆ ~BaseFitObject()

~BaseFitObject ( )

virtual

Virtual destructor.

Definition at line 91 of file BaseFitObject.cc.

    {
      delete[] name;
    }

Member Function Documentation

◆ addTo1stDerivatives()

void addTo1stDerivatives	(	double	M[],
		int	idim,
		double	der[],
		int	kglobal,
		int	metaSet
	)		const

virtual

Definition at line 515 of file BaseFitObject.cc.

    {
      if (!cachevalid) updateCache();
      for (int ilocal = 0; ilocal < getNPar(); ilocal++) {
        int iglobal = globalParNum[ilocal];
        if (iglobal >= 0) {
          for (int j = 0; j < BaseDefs::nMetaVars[metaSet]; j++) {
            double x = der[j] * getFirstDerivative_Meta_Local(j, ilocal, metaSet);
            M[idim * kglobal + iglobal] += x;
            M[idim * iglobal + kglobal] += x;
          }
        }
      }
      return;
    }

◆ addTo2ndDerivatives() [1/2]

void addTo2ndDerivatives	(	double	der2[],
		int	idim,
		double	factor[],
		int	metaSet
	)		const

virtual

Definition at line 534 of file BaseFitObject.cc.

    {
      if (!cachevalid) updateCache();
      for (int ilocal = 0; ilocal < getNPar(); ilocal++) {
        int iglobal = getGlobalParNum(ilocal);
        if (iglobal < 0) continue;
        for (int jlocal = ilocal; jlocal < getNPar(); jlocal++) {
          int jglobal = getGlobalParNum(jlocal);
          if (jglobal < 0) continue;
          double sum(0);
          for (int imeta = 0; imeta < BaseDefs::nMetaVars[metaSet]; imeta++) {
            sum += factor[imeta] * getSecondDerivative_Meta_Local(imeta, ilocal, jlocal, metaSet);
          }
          der2[idim * iglobal + jglobal] += sum;
          if (iglobal != jglobal) der2[idim * jglobal + iglobal] += sum;
        }
      }
      return;
    }

◆ addTo2ndDerivatives() [2/2]

void addTo2ndDerivatives	(	double	M[],
		int	idim,
		double	lambda,
		double	der[],
		int	metaSet
	)		const

virtual

Definition at line 559 of file BaseFitObject.cc.

    {
      double factor[BaseDefs::MAXINTERVARS];
      for (int i = 0; i < BaseDefs::nMetaVars[metaSet]; i++) factor[i] = lambda * der[i];
      addTo2ndDerivatives(M, idim, factor, metaSet);
      return;
    }

◆ addToGlobalChi2DerMatrix()

void addToGlobalChi2DerMatrix	(	double *	M,
		int	idim
	)		const

virtual

Add 2nd derivatives of chi squared to global derivative matrix.

Parameters

M	Global derivative matrix
idim	First dimension of global derivative matrix

Definition at line 462 of file BaseFitObject.cc.

    {
      if (!covinvvalid) calculateCovInv();
      for (int ilocal = 0; ilocal < getNPar(); ++ilocal) {
        if (!isParamFixed(ilocal) && isParamMeasured(ilocal)) {
          int iglobal = getGlobalParNum(ilocal);
          assert(iglobal >= 0 && iglobal < idim);
          int ioffs = idim * iglobal;
          for (int jlocal = 0; jlocal < getNPar(); ++jlocal) {
            if (!isParamFixed(jlocal) && isParamMeasured(jlocal)) {
              int jglobal = getGlobalParNum(jlocal);
              assert(jglobal >= 0 && jglobal < idim);
              M[ioffs + jglobal] += getD2Chi2DParam2(ilocal, jlocal);
            }
          }
        }
      }
    }

◆ addToGlobalChi2DerVector() [1/2]

int addToGlobalChi2DerVector	(	double *	y,
		int	idim
	)		const

virtual

Add derivatives of chi squared to global derivative vector.

Parameters

y	Vector of chi2 derivatives
idim	Vector size

Definition at line 482 of file BaseFitObject.cc.

    {
      // This adds the dChi2/dpar piece
      assert(getNPar() <= BaseDefs::MAXPAR);
      if (not covinvvalid and not calculateCovInv()) return 0;
      for (int ilocal = 0; ilocal < getNPar(); ++ilocal) {
        if (!isParamFixed(ilocal) && isParamMeasured(ilocal)) {
          int iglobal = getGlobalParNum(ilocal);
          assert(iglobal >= 0 && iglobal < idim);
          y[iglobal] += getDChi2DParam(ilocal);
        }
      }
      return 0;
    }

◆ addToGlobalChi2DerVector() [2/2]

void addToGlobalChi2DerVector	(	double *	y,
		int	idim,
		double	lambda,
		double	der[],
		int	metaSet
	)		const

virtual

Add derivatives of momentum vector to global derivative vector.

Parameters

y	Vector of chi2 derivatives
idim	Vector size
lambda	The lambda value
der	derivatives of constraint wrt intermediate variables (e.g. 4-vector with dg/dE, dg/dpx, dg/dpy, dg/dpz)
metaSet	which set of intermediate variables

Definition at line 498 of file BaseFitObject.cc.

    {
      // This adds the lambda * dConst/dpar piece
      if (!cachevalid) updateCache();
      for (int ilocal = 0; ilocal < getNPar(); ilocal++) {
        int iglobal = globalParNum[ilocal];
        if (iglobal >= 0) {
          assert(iglobal < idim);
          for (int j = 0; j < BaseDefs::nMetaVars[metaSet]; j++) {
            y[iglobal] += lambda * der[j] * getFirstDerivative_Meta_Local(j, ilocal, metaSet);
          }
        }
      }
    }

◆ addToGlobCov()

void addToGlobCov	(	double *	glcov,
		int	idim
	)		const

virtual

Add covariance matrix elements to global covariance matrix of size idim x idim.

Parameters

glcov	Global covariance matrix
idim	First dimension of global derivative matrix

Definition at line 217 of file BaseFitObject.cc.

    {
      for (int ilocal = 0; ilocal < getNPar(); ++ilocal) {
        if (!isParamFixed(ilocal) && isParamMeasured(ilocal)) {
          int iglobal = getGlobalParNum(ilocal);
          assert(iglobal >= 0 && iglobal < idim);
          int ioffs = idim * iglobal;
          for (int jlocal = 0; jlocal < getNPar(); ++jlocal) {
            if (!isParamFixed(jlocal) && isParamMeasured(jlocal)) {
              int jglobal = getGlobalParNum(jlocal);
              assert(jglobal >= 0 && jglobal < idim);
              globCov[ioffs + jglobal] += getCov(ilocal, jlocal);
            }
          }
        }
      }
    }

◆ assign()

BaseFitObject & assign ( const BaseFitObject & source )

virtual

Assign from anther object, if of same type.

Parameters

source The source object

Reimplemented in ISRPhotonFitObject, JetFitObject, NeutrinoFitObject, ParticleFitObject, and PxPyPzMFitObject.

Definition at line 68 of file BaseFitObject.cc.

    {
      if (&source != this) {
        name = nullptr;
        setName(source.name);
        for (int i = 0; i < BaseDefs::MAXPAR; ++i) {
          par[i]          = source.par[i];
          mpar[i]         = source.mpar[i];
          measured[i]     = source.measured[i];
          fixed[i]        = source.fixed[i];
          globalParNum[i] = source.globalParNum[i];
          for (int j = 0; j < BaseDefs::MAXPAR; ++j) {
            cov[i][j] = source.cov[i][j];
            covinv[i][j] = source.covinv[i][j];
          }
        }
        covinvvalid = false;
        cachevalid = false;
      }
      return *this;
    }

◆ calculateCovInv()

bool calculateCovInv ( ) const

protectedvirtual

Calculate the inverse of the covariance matrix.

Definition at line 235 of file BaseFitObject.cc.

    {
 
      // DANIEL added
 
      int n = getNPar();
 
      gsl_matrix* covm = gsl_matrix_alloc(n, n);
      gsl_matrix_set_identity(covm);
 
      for (int i = 0; i < n; ++i) {
        if (isParamMeasured(i)) {
          for (int j = 0; j < n; ++j) {
            if (isParamMeasured(j)) {
              gsl_matrix_set(covm, i, j, cov[i][j]);
            }
          }
        }
      }
      gsl_error_handler_t* e = gsl_set_error_handler_off();
      int result = gsl_linalg_cholesky_decomp(covm);
      if (result == 0) result = gsl_linalg_cholesky_invert(covm);
      gsl_set_error_handler(e);
 
      for (int i = 0; i < n; ++i) {
        for (int j = 0; j < i; ++j) {
          covinv[i][j] = covinv[j][i] = gsl_matrix_get(covm, i, j);
        }
        covinv[i][i] = gsl_matrix_get(covm, i, i);
      }
 
      gsl_matrix_free(covm);
      covinvvalid = (result == 0);
 
      if (!covinvvalid) {
        B2WARNING("ERROR, COULD NOT INVERT COV MATR!");
 
        B2INFO("COV ");
        for (int i = 0; i < n; ++i) {
          if (isParamMeasured(i)) {
            for (int j = 0; j < n; ++j) {
              if (isParamMeasured(j)) {
                B2INFO(cov[i][j] << " ");
              }
            }
          }
        }
 
        B2INFO("CORREL ");
        for (int i = 0; i < n; ++i) {
          if (isParamMeasured(i)) {
            for (int j = 0; j < n; ++j) {
              if (isParamMeasured(j)) {
                B2INFO(cov[i][j] / sqrt(cov[i][i]*cov[j][j]) << " ");
              }
            }
          }
        }
 
      }
 
      return covinvvalid;
    }

◆ copy()

virtual BaseFitObject * copy ( ) const

pure virtual

Return a new copy of itself.

Implemented in ISRPhotonFitObject, JetFitObject, NeutrinoFitObject, and PxPyPzMFitObject.

◆ fixParam()

bool fixParam	(	int	ilocal,
		bool	fix = `true`
	)

virtual

Fix a parameter (fix=true), or release it (fix=false)

Parameters

ilocal	Local parameter number
fix	fix if true, release if false

Definition at line 352 of file BaseFitObject.cc.

    {
      assert(ilocal >= 0 && ilocal < getNPar());
      return fixed [ilocal] = fix;
    }

◆ getChi2()

double getChi2 ( ) const

virtual

Get chi squared from measured and fitted parameters.

Reimplemented in ParticleFitObject.

Definition at line 414 of file BaseFitObject.cc.

    {
      if (not covinvvalid and not calculateCovInv()) return -1;
      double chi2 = 0;
      static double resid[BaseDefs::MAXPAR];
      static bool chi2contr[BaseDefs::MAXPAR];
      for (int i = 0; i < getNPar(); ++i) {
        resid[i] = par[i] - mpar[i];
 
        B2INFO(" BaseFitObject::getChi2() " << i << " " << resid[i] << " " << covinv[i][i]);
 
        if ((chi2contr[i] = (isParamMeasured(i) && !isParamFixed(i)))) {
          chi2 += resid[i] * covinv[i][i] * resid[i];
          for (int j = 0; j < getNPar(); ++j) {
            if (j < i && chi2contr[j])  chi2 += 2 * resid[i] * covinv[i][j] * resid[j];
          }
        }
      }
      return chi2;
    }

◆ getCov()

double getCov	(	int	ilocal,
		int	jlocal
	)		const

virtual

Get covariance between parameters ilocal and jlocal.

Parameters

ilocal	Local parameter number i
jlocal	Local parameter number j

Reimplemented in JetFitObject.

Definition at line 389 of file BaseFitObject.cc.

    {
      assert(ilocal >= 0 && ilocal < getNPar());
      assert(jlocal >= 0 && jlocal < getNPar());
      return cov[ilocal][jlocal];
    }

◆ getD2Chi2DParam2()

double getD2Chi2DParam2	(	int	ilocal,
		int	jlocal
	)		const

virtual

Get second derivative of chi squared w.r.t. parameters ilocal1 and ilocal2.

Parameters

ilocal	Local parameter number i
jlocal	Local parameter number j

Definition at line 449 of file BaseFitObject.cc.

    {
      assert(ilocal >= 0 && ilocal < getNPar());
      assert(jlocal >= 0 && jlocal < getNPar());
      if (isParamFixed(ilocal) || !isParamMeasured(ilocal) ||
          isParamFixed(jlocal) || !isParamMeasured(jlocal))
        return 0;
      if (not covinvvalid and not calculateCovInv()) return 0;
      return 2 * covinv[ilocal][jlocal]; // JL: ok in absence of soft constraints
    }

◆ getDChi2DParam()

double getDChi2DParam ( int ilocal ) const

virtual

Get derivative of chi squared w.r.t. parameter ilocal.

Parameters

ilocal Local parameter number

Definition at line 435 of file BaseFitObject.cc.

    {
      assert(ilocal >= 0 && ilocal < getNPar());
      if (isParamFixed(ilocal) || !isParamMeasured(ilocal)) return 0;
 
      if (not covinvvalid and not calculateCovInv()) return 0;
 
      double result = 0;
      for (int jlocal = 0; jlocal < getNPar(); jlocal++)
        if (!isParamFixed(jlocal) && isParamMeasured(jlocal))
          result += covinv[ilocal][jlocal] * (par[jlocal] - mpar[jlocal]);
      return 2 * result;
    }

◆ getError()

double getError ( int ilocal ) const

virtual

Get error of parameter ilocal.

Parameters

ilocal Local parameter number

Reimplemented in JetFitObject.

Definition at line 383 of file BaseFitObject.cc.

    {
      assert(ilocal >= 0 && ilocal < getNPar());
      return std::sqrt(cov[ilocal][ilocal]);
    }

◆ getError2()

double getError2	(	double	der[],
		int	metaset
	)		const

virtual

Definition at line 597 of file BaseFitObject.cc.

    {
      if (!cachevalid) updateCache();
      double totError(0);
      for (int i = 0; i < BaseDefs::nMetaVars[metaSet]; i++) {
        for (int j = 0; j < BaseDefs::nMetaVars[metaSet]; j++) {
          double cov_i_j = 0; // this will become the covariance of intermediate variables i and j
          for (int k = 0; k < getNPar(); k++) {
            for (int l = 0; l < getNPar(); l++) {
              cov_i_j += getFirstDerivative_Meta_Local(i, k, metaSet) * cov[k][l] * getFirstDerivative_Meta_Local(j, l, metaSet);
            }
          }
          totError += der[i] * der[j] * cov_i_j;
        }
      }
      return totError;
    }

◆ getFirstDerivative_Meta_Local()

virtual double getFirstDerivative_Meta_Local	(	int	iMeta,
		int	ilocal,
		int	metaSet
	)		const

pure virtual

Implemented in JetFitObject.

◆ getGlobalParNum()

int getGlobalParNum ( int ilocal ) const

virtual

Get global parameter number of parameter ilocal.

Parameters

ilocal Local parameter number

Definition at line 365 of file BaseFitObject.cc.

    {
      if (ilocal < 0 || ilocal >= getNPar()) return -1;
      return globalParNum[ilocal];
    }

◆ getMParam()

double getMParam ( int ilocal ) const

virtual

Get measured value of parameter ilocal.

Parameters

ilocal Local parameter number

Definition at line 377 of file BaseFitObject.cc.

    {
      assert(ilocal >= 0 && ilocal < getNPar());
      return mpar[ilocal];
    }

◆ getName()

const char * getName ( ) const

virtual

Get object's name.

Definition at line 107 of file BaseFitObject.cc.

    {
      return name ? name : "???";
    }

◆ getNFixed()

int getNFixed ( ) const

virtual

Get number of fixed parameters of this FitObject.

Definition at line 130 of file BaseFitObject.cc.

    {
      int nfixed = 0;
      for (int i = 0; i < getNPar(); ++i) if (isParamFixed(i)) ++nfixed;
      return nfixed;
    }

◆ getNFree()

int getNFree ( ) const

virtual

Get number of free parameters of this FitObject.

Definition at line 124 of file BaseFitObject.cc.

    {
      int nfree = 0;
      for (int i = 0; i < getNPar(); ++i) if (!isParamFixed(i)) ++nfree;
      return nfree;
    }

◆ getNMeasured()

int getNMeasured ( ) const

virtual

Get number of measured parameters of this FitObject.

Definition at line 112 of file BaseFitObject.cc.

    {
      int nmeasured = 0;
      for (int i = 0; i < getNPar(); ++i) if (isParamMeasured(i) && !isParamFixed(i)) ++nmeasured;
      return nmeasured;
    }

◆ getNPar()

virtual int getNPar ( ) const

pure virtual

Get total number of parameters of this FitObject.

Implemented in ISRPhotonFitObject, JetFitObject, NeutrinoFitObject, and PxPyPzMFitObject.

◆ getNUnmeasured()

int getNUnmeasured ( ) const

virtual

Get number of unmeasured parameters of this FitObject.

Definition at line 118 of file BaseFitObject.cc.

    {
      int nunmeasrd = 0;
      for (int i = 0; i < getNPar(); ++i) if (!isParamMeasured(i) && !isParamFixed(i)) ++nunmeasrd;
      return nunmeasrd;
    }

◆ getParam()

double getParam ( int ilocal ) const

virtual

Get current value of parameter ilocal.

Parameters

ilocal Local parameter number

Definition at line 371 of file BaseFitObject.cc.

    {
      assert(ilocal >= 0 && ilocal < getNPar());
      return par[ilocal];
    }

◆ getParamName()

virtual const char * getParamName ( int ilocal ) const

pure virtual

Get name of parameter ilocal.

Parameters

ilocal Local parameter number

Implemented in ISRPhotonFitObject, JetFitObject, NeutrinoFitObject, and PxPyPzMFitObject.

◆ getRho()

double getRho	(	int	ilocal,
		int	jlocal
	)		const

virtual

Get correlation coefficient between parameters ilocal and jlocal.

Parameters

ilocal	Local parameter number i
jlocal	Local parameter number j

Definition at line 396 of file BaseFitObject.cc.

    {
      assert(ilocal >= 0 && ilocal < getNPar());
      assert(jlocal >= 0 && jlocal < getNPar());
      return cov[ilocal][jlocal] / std::sqrt(cov[ilocal][ilocal] * cov[jlocal][jlocal]);
    }

◆ initCov()

void initCov ( )

virtual

Definition at line 587 of file BaseFitObject.cc.

    {
      for (int i = 0; i < getNPar(); ++i) {
        for (int j = 0; j < getNPar(); ++j) {
          cov[i][j] = static_cast<double>(i == j);
        }
      }
    }

◆ invalidateCache()

void invalidateCache ( ) const

inline

invalidate any cached quantities

Definition at line 241 of file BaseFitObject.h.

241{cachevalid = false;};

◆ isParamFixed()

bool isParamFixed ( int ilocal ) const

virtual

Returns whether parameter is fixed.

Parameters

ilocal Local parameter number

Definition at line 408 of file BaseFitObject.cc.

    {
      assert(ilocal >= 0 && ilocal < getNPar());
      return fixed[ilocal];
    }

◆ isParamMeasured()

bool isParamMeasured ( int ilocal ) const

virtual

Get measured flag for parameter ilocal.

Parameters

ilocal Local parameter number

Definition at line 402 of file BaseFitObject.cc.

    {
      assert(ilocal >= 0 && ilocal < getNPar());
      return measured[ilocal];
    }

◆ operator=()

BaseFitObject & operator= ( const BaseFitObject & rhs )

Assignment.

Parameters

rhs	right hand side

Definition at line 60 of file BaseFitObject.cc.

    {
      if (this != &rhs) {
        assign(rhs);  // calls virtual function assign of derived class
      }
      return *this;
    }

◆ print()

virtual std::ostream & print ( std::ostream & os ) const

pure virtual

print object to ostream

Parameters

os	The output stream

Implemented in ParticleFitObject.

◆ print1stDerivatives()

std::ostream & print1stDerivatives ( std::ostream & os ) const

virtual

print the 1st derivatives wrt metaSet 0 (E, px, py, pz)

Parameters

os	The output stream

Definition at line 164 of file BaseFitObject.cc.

    {
      int metaSet = 0;
      os << "#";
      for (int i = 0; i < BaseDefs::nMetaVars[metaSet]; ++i) {
        for (int j = 0; j < getNPar(); ++j) {
          if (i > 0 && j == 0) os << ",";
          os << " " << getFirstDerivative_Meta_Local(i, j, metaSet);
        }
      }
      os << "#" << std::endl;
      return os;
    }

◆ print2ndDerivatives()

std::ostream & print2ndDerivatives ( std::ostream & os ) const

virtual

print the 2nd derivatives wrt metaSet 0 (E, px, py, pz)

Parameters

os	The output stream

Definition at line 178 of file BaseFitObject.cc.

    {
      int metaSet = 0;
      os << "##";
      for (int i = 0; i < BaseDefs::nMetaVars[metaSet]; ++i) {
        if (i > 0) os << std::endl << "  ";
        for (int j = 0; j < getNPar(); ++j) {
          for (int k = 0; k < getNPar(); ++k) {
//              if (i>0 && k==0) os << ",";
            os << " " << getSecondDerivative_Meta_Local(i, j, k, metaSet);
            if (k == getNPar() - 1) os << ",";
          }
        }
      }
      os << "##" << std::endl;
      return os;
    }

◆ printParams()

std::ostream & printParams ( std::ostream & os ) const

virtual

print the parameters and errors

Parameters

os	The output stream

Definition at line 137 of file BaseFitObject.cc.

    {
      os << "(";
      for (int i = 0; i < getNPar(); ++i) {
        if (i > 0) os << ", ";
        os << " " << getParam(i);
        if (isParamFixed(i))  os << " fix";
//    else if (getError(i)>0) os << " \261 " << getError(i);
        else if (getError(i) > 0) os << " +- " << getError(i);   // " \261 " appeared as <B1> in ASCII ...    Graham
      }
      os << ")";
      return os;
    }

◆ printRhoValues()

std::ostream & printRhoValues ( std::ostream & os ) const

virtual

print the correlation coefficients

Parameters

os	The output stream

Definition at line 151 of file BaseFitObject.cc.

    {
      os << "{";
      for (int i = 0; i < getNPar(); ++i) {
        for (int j = 0; j < getNPar(); ++j) {
          if (i > 0 && j == 0) os << ",";
          os << " " << getRho(i, j);
        }
      }
      os << "}" << std::endl;
      return os;
    }

◆ releaseParam()

virtual bool releaseParam ( int ilocal )

inlinevirtual

Release a parameter.

Parameters

ilocal Local parameter number

Definition at line 157 of file BaseFitObject.h.

159 { return fixParam(ilocal, false); }

Belle2::OrcaKinFit::BaseFitObject::fixParam

virtual bool fixParam(int ilocal, bool fix=true)

Fix a parameter (fix=true), or release it (fix=false)

Definition: BaseFitObject.cc:352

◆ setCov()

bool setCov	(	int	ilocal,
		int	jlocal,
		double	cov_
	)

virtual

Set covariance of parameters ilocal and jlocal; return: success.

Parameters

ilocal	Local parameter number
jlocal	Local parameter number
cov_	New error value

Definition at line 340 of file BaseFitObject.cc.

    {
      if (!isfinite(cov_)) return false;
      assert(ilocal >= 0 && ilocal < getNPar());
      assert(jlocal >= 0 && jlocal < getNPar());
      invalidateCache();
      covinvvalid = false;
      cov[ilocal][jlocal] = cov[jlocal][ilocal] = cov_;
      return true;
    }

◆ setError()

bool setError	(	int	ilocal,
		double	err_
	)

virtual

Set error of parameter ilocal; return: success.

Parameters

ilocal	Local parameter number
err_	New error value

Definition at line 330 of file BaseFitObject.cc.

    {
      if (!isfinite(err_)) return false;
      assert(ilocal >= 0 && ilocal < getNPar());
      invalidateCache();
      covinvvalid = false;
      cov[ilocal][ilocal] = err_ * err_;
      return true;
    }

◆ setGlobalParNum()

bool setGlobalParNum	(	int	ilocal,
		int	iglobal
	)

virtual

Set number of parameter ilocal in global list return true signals OK.

Parameters

ilocal	Local parameter number
iglobal	New global parameter number

Definition at line 358 of file BaseFitObject.cc.

    {
      if (ilocal < 0 || ilocal >= getNPar()) return false;
      globalParNum[ilocal] = iglobal;
      return true;
    }

◆ setMParam()

bool setMParam	(	int	ilocal,
		double	mpar_
	)

virtual

Set measured value of parameter ilocal; return: success.

Parameters

ilocal	Local parameter number
mpar_	New measured parameter value

Definition at line 320 of file BaseFitObject.cc.

    {
      if (!isfinite(mpar_)) return false;
      assert(ilocal >= 0 && ilocal < getNPar());
      if (mpar[ilocal] == mpar_) return true;
      invalidateCache();
      mpar[ilocal] = mpar_;
      return true;
    }

◆ setName()

void setName ( const char * name_ )

Set object's name.

Definition at line 98 of file BaseFitObject.cc.

    {
      if (name_ == nullptr) return;
      size_t l = strlen(name_);
      if (name) delete[] name;
      name = new char[l + 1];
      strcpy(name, name_);
    }

◆ setParam() [1/2]

bool setParam	(	int	ilocal,
		double	par_
	)

virtual

Set value of parameter ilocal; return: significant change.

Parameters

ilocal	Local parameter number
par_	New parameter value

Definition at line 309 of file BaseFitObject.cc.

    {
      if (!isfinite(par_)) return true;
      assert(ilocal >= 0 && ilocal < getNPar());
      if (par[ilocal] == par_) return false;
      invalidateCache();
      bool result = pow((par_ - par[ilocal]), 2) > eps2 * cov[ilocal][ilocal];
      par[ilocal] = par_;
      return result;
    }

◆ setParam() [2/2]

bool setParam	(	int	ilocal,
		double	par_,
		bool	measured_,
		bool	fixed_ = `false`
	)

virtual

Set value and measured flag of parameter i; return: significant change.

Parameters

ilocal	Local parameter number
par_	New parameter value
measured_	New "measured" flag
fixed_	New "fixed" flag

Definition at line 299 of file BaseFitObject.cc.

    {
      assert(ilocal >= 0 && ilocal < getNPar());
      if (measured[ilocal] != measured_ || fixed[ilocal] != fixed_) invalidateCache();
      measured[ilocal] = measured_;
      fixed[ilocal] = fixed_;
      return setParam(ilocal, par_);
    }

◆ updateParams()

bool updateParams	(	double	p[],
		int	idim
	)

virtual

Read values from global vector, readjust vector; return: significant change.

Parameters

p	The parameter vector
idim	Length of the vector

Reimplemented in ISRPhotonFitObject, JetFitObject, NeutrinoFitObject, and PxPyPzMFitObject.

Definition at line 196 of file BaseFitObject.cc.

    {
      bool result = false;
      invalidateCache();
      for (int ilocal = 0; ilocal < getNPar(); ++ilocal) {
        if (!isParamFixed(ilocal)) {   // daniel added this
          int iglobal = getGlobalParNum(ilocal);
          assert(iglobal >= 0 && iglobal < idim);
          //      result = result || setParam (ilocal, p[iglobal]); // daniel thinks this is a BUG !!! if first param is successfully updated, no further ones are
          // result = setParam (ilocal, p[iglobal]) || result; // daniel thinks this would be OK
          // but to makes it clearer, does the following:
          bool thisresult = setParam(ilocal, p[iglobal]);
          result = result || thisresult;
          // if illegal value: read back legal value
          if (!thisresult)   // daniel added
            p[iglobal] = getParam(ilocal);
        }
      }
      return result;
    }

Friends And Related Function Documentation

◆ operator<<()

std::ostream & operator<<	(	std::ostream &	os,
		const BaseFitObject &	bfo
	)

Parameters

os	The output stream
bfo	The object to print

Definition at line 341 of file BaseFitObject.h.

    {
      return bfo.print(os);
    }

Member Data Documentation

◆ cachevalid

bool cachevalid

mutableprotected

flag for valid cache

Definition at line 333 of file BaseFitObject.h.

◆ cov

double cov[BaseDefs::MAXPAR][BaseDefs::MAXPAR]

protected

local covariance matrix

Definition at line 327 of file BaseFitObject.h.

◆ covinv

double covinv[BaseDefs::MAXPAR][BaseDefs::MAXPAR]

mutableprotected

inverse pf local covariance matrix

Definition at line 329 of file BaseFitObject.h.

◆ covinvvalid

bool covinvvalid

mutableprotected

flag for valid inverse covariance matrix

Definition at line 331 of file BaseFitObject.h.

◆ eps2

const double eps2 = 0.0001

staticprotected

Definition at line 307 of file BaseFitObject.h.

◆ fixed

bool fixed[BaseDefs::MAXPAR]

protected

fixed flag

Definition at line 323 of file BaseFitObject.h.

◆ globalParNum

int globalParNum[BaseDefs::MAXPAR]

protected

global parameter number for each parameter

Definition at line 325 of file BaseFitObject.h.

◆ measured

bool measured[BaseDefs::MAXPAR]

protected

measured flag

Definition at line 321 of file BaseFitObject.h.

◆ mpar

double mpar[BaseDefs::MAXPAR]

protected

measured parameters

Definition at line 319 of file BaseFitObject.h.

◆ name

char* name

protected

Definition at line 306 of file BaseFitObject.h.

◆ par

double par[BaseDefs::MAXPAR]

protected

fit parameters

Definition at line 317 of file BaseFitObject.h.

The documentation for this class was generated from the following files:

analysis/OrcaKinFit/include/BaseFitObject.h
analysis/OrcaKinFit/src/BaseFitObject.cc

Public Member Functions

Protected Member Functions

Protected Attributes

Static Protected Attributes

Related Functions

Detailed Description

Constructor & Destructor Documentation

◆ BaseFitObject() [1/2]

◆ BaseFitObject() [2/2]

◆ ~BaseFitObject()

Member Function Documentation

◆ addTo1stDerivatives()

◆ addTo2ndDerivatives() [1/2]

◆ addTo2ndDerivatives() [2/2]

◆ addToGlobalChi2DerMatrix()

◆ addToGlobalChi2DerVector() [1/2]

◆ addToGlobalChi2DerVector() [2/2]

◆ addToGlobCov()

◆ assign()

◆ calculateCovInv()

◆ copy()

◆ fixParam()

◆ getChi2()

◆ getCov()

◆ getD2Chi2DParam2()

◆ getDChi2DParam()

◆ getError()

◆ getError2()

◆ getFirstDerivative_Meta_Local()

◆ getGlobalParNum()

◆ getMParam()

◆ getName()

◆ getNFixed()

◆ getNFree()

◆ getNMeasured()

◆ getNPar()

◆ getNUnmeasured()

◆ getParam()

◆ getParamName()

◆ getRho()

◆ initCov()

◆ invalidateCache()

◆ isParamFixed()

◆ isParamMeasured()

◆ operator=()

◆ print()

◆ print1stDerivatives()

◆ print2ndDerivatives()

◆ printParams()

◆ printRhoValues()

◆ releaseParam()

◆ setCov()

◆ setError()

◆ setGlobalParNum()

◆ setMParam()

◆ setName()

◆ setParam() [1/2]

◆ setParam() [2/2]

◆ updateParams()

Friends And Related Function Documentation

◆ operator<<()

Member Data Documentation

◆ cachevalid

◆ cov

◆ covinv

◆ covinvvalid

◆ eps2

◆ fixed

◆ globalParNum

◆ measured

◆ mpar

◆ name

◆ par