Data (block) for independent scalar measurement. More...

#include <GblData.h>

Public Member Functions
	GblData (unsigned int aLabel, double aMeas, double aPrec)
	Create data block.

virtual	~GblData ()

void	addDerivatives (unsigned int iRow, const std::vector< unsigned int > &labDer, const SMatrix55 &matDer, unsigned int iOff, const TMatrixD &derLocal, const std::vector< int > &labGlobal, const TMatrixD &derGlobal, unsigned int nLocal, const TMatrixD &derTrans)
	Add derivatives from measurement.

void	addDerivatives (unsigned int iRow, const std::vector< unsigned int > &labDer, const SMatrix27 &matDer, unsigned int nLocal, const TMatrixD &derTrans)
	Add derivatives from kink.

void	addDerivatives (const std::vector< unsigned int > &index, const std::vector< double > &derivatives)
	Add derivatives from external seed.

void	setPrediction (const VVector &aVector)
	Calculate prediction for data from fit (by GblTrajectory::fit).

double	setDownWeighting (unsigned int aMethod)
	Outlier down weighting with M-estimators (by GblTrajectory::fit).

double	getChi2 () const
	Calculate Chi2 contribution.

void	printData () const
	Print data block.

void	getLocalData (double &aValue, double &aWeight, std::vector< unsigned int > &indLocal, std::vector< double > &derLocal)
	Get Data for local fit.

void	getAllData (double &aValue, double &aErr, std::vector< unsigned int > &indLocal, std::vector< double > &derLocal, std::vector< int > &labGlobal, std::vector< double > &derGlobal)
	Get all Data for MP-II binary record.

void	getResidual (double &aResidual, double &aVariance, double &aDownWeight, std::vector< unsigned int > &indLocal, std::vector< double > &derLocal)
	Get data for residual (and errors).

Private Attributes
unsigned int	theLabel
	Label (of measurements point)

double	theValue
	Value (residual)

double	thePrecision
	Precision (1/sigma**2)

double	theDownWeight
	Down-weighting factor (0-1)

double	thePrediction
	Prediction from fit.

std::vector< unsigned int >	theParameters
	List of fit parameters (with non zero derivatives)

std::vector< double >	theDerivatives
	List of derivatives for fit.

std::vector< int >	globalLabels
	Labels for global derivatives.

std::vector< double >	globalDerivatives
	Global derivatives.

Detailed Description

Data (block) for independent scalar measurement.

Data (block) containing value, precision and derivatives for measurements and kinks. Created from attributes of GblPoints, used to construct linear equation system for track fit.

Definition at line 33 of file GblData.h.

Constructor & Destructor Documentation

◆ GblData()

gbl::GblData::GblData	(	unsigned int	aLabel,
		double	aValue,
		double	aPrec
	)

Create data block.

Parameters

[in]	aLabel	Label of corresponding point
[in]	aValue	Value of (scalar) measurement
[in]	aPrec	Precision of (scalar) measurement

Definition at line 19 of file GblData.cc.

                                                                 :
        theLabel(aLabel), theValue(aValue), thePrecision(aPrec), theDownWeight(
                1.), thePrediction(0.), theParameters(), theDerivatives(), globalLabels(), globalDerivatives() {
 
}

◆ ~GblData()

gbl::GblData::~GblData ( )

virtual

Definition at line 25 of file GblData.cc.

25 {

26}

Member Function Documentation

◆ addDerivatives() [1/3]

void gbl::GblData::addDerivatives	(	const std::vector< unsigned int > &	index,
		const std::vector< double > &	derivatives
	)

Add derivatives from external seed.

Add (non-zero) derivatives to data block. Fill list of labels of used fit parameters.

Parameters

[in]	index	Labels for derivatives
[in]	derivatives	Derivatives (vector)

Definition at line 120 of file GblData.cc.

                                              {
    for (unsigned int i = 0; i < derivatives.size(); ++i) // any derivatives
            {
        if (derivatives[i]) {
            theParameters.push_back(index[i]);
            theDerivatives.push_back(derivatives[i]);
        }
    }
}

◆ addDerivatives() [2/3]

void gbl::GblData::addDerivatives	(	unsigned int	iRow,
		const std::vector< unsigned int > &	labDer,
		const SMatrix27 &	matDer,
		unsigned int	extOff,
		const TMatrixD &	extDer
	)

Add derivatives from kink.

Add (non-zero) derivatives to data block. Fill list of labels of used fit parameters.

Parameters

[in]	iRow	Row index (0-1) in 2D kink
[in]	labDer	Labels for derivatives
[in]	matDer	Derivatives (matrix) 'kink vs track fit parameters'
[in]	extOff	Offset for external parameters
[in]	extDer	Derivatives for external Parameters

Definition at line 89 of file GblData.cc.

                                                     {
 
    unsigned int nParMax = 7 + extDer.GetNcols();
    theParameters.reserve(nParMax); // have to be sorted
    theDerivatives.reserve(nParMax);
 
    for (int i = 0; i < extDer.GetNcols(); ++i) // external derivatives
            {
        if (extDer(iRow, i)) {
            theParameters.push_back(extOff + i + 1);
            theDerivatives.push_back(extDer(iRow, i));
        }
    }
 
    for (unsigned int i = 0; i < 7; ++i) // curvature, offset derivatives
            {
        if (labDer[i] and matDer(iRow, i)) {
            theParameters.push_back(labDer[i]);
            theDerivatives.push_back(matDer(iRow, i));
        }
    }
}

◆ addDerivatives() [3/3]

void gbl::GblData::addDerivatives	(	unsigned int	iRow,
		const std::vector< unsigned int > &	labDer,
		const SMatrix55 &	matDer,
		unsigned int	iOff,
		const TMatrixD &	derLocal,
		const std::vector< int > &	labGlobal,
		const TMatrixD &	derGlobal,
		unsigned int	extOff,
		const TMatrixD &	extDer
	)

Add derivatives from measurement.

Add (non-zero) derivatives to data block. Fill list of labels of used fit parameters.

Parameters

[in]	iRow	Row index (0-4) in up to 5D measurement
[in]	labDer	Labels for derivatives
[in]	matDer	Derivatives (matrix) 'measurement vs track fit parameters'
[in]	iOff	Offset for row index for additional parameters
[in]	derLocal	Derivatives (matrix) for additional local parameters
[in]	labGlobal	Labels for additional global (MP-II) parameters
[in]	derGlobal	Derivatives (matrix) for additional global (MP-II) parameters
[in]	extOff	Offset for external parameters
[in]	extDer	Derivatives for external Parameters

Definition at line 41 of file GblData.cc.

                                                     {
 
    unsigned int nParMax = 5 + derLocal.GetNcols() + extDer.GetNcols();
    theParameters.reserve(nParMax); // have to be sorted
    theDerivatives.reserve(nParMax);
 
    for (int i = 0; i < derLocal.GetNcols(); ++i) // local derivatives
            {
        if (derLocal(iRow - iOff, i)) {
            theParameters.push_back(i + 1);
            theDerivatives.push_back(derLocal(iRow - iOff, i));
        }
    }
 
    for (int i = 0; i < extDer.GetNcols(); ++i) // external derivatives
            {
        if (extDer(iRow - iOff, i)) {
            theParameters.push_back(extOff + i + 1);
            theDerivatives.push_back(extDer(iRow - iOff, i));
        }
    }
 
    for (unsigned int i = 0; i < 5; ++i) // curvature, offset derivatives
            {
        if (labDer[i] and matDer(iRow, i)) {
            theParameters.push_back(labDer[i]);
            theDerivatives.push_back(matDer(iRow, i));
        }
    }
 
    globalLabels = labGlobal;
    for (int i = 0; i < derGlobal.GetNcols(); ++i) // global derivatives
        globalDerivatives.push_back(derGlobal(iRow - iOff, i));
}

◆ getAllData()

void gbl::GblData::getAllData	(	double &	aValue,
		double &	aErr,
		std::vector< unsigned int > *&	indLocal,
		std::vector< double > *&	derLocal,
		std::vector< int > *&	labGlobal,
		std::vector< double > *&	derGlobal
	)

Get all Data for MP-II binary record.

Parameters

[out]	aValue	Value
[out]	aErr	Error
[out]	indLocal	List of labels of local parameters
[out]	derLocal	List of derivatives for local parameters
[out]	labGlobal	List of labels of global parameters
[out]	derGlobal	List of derivatives for global parameters

Definition at line 220 of file GblData.cc.

                                                                   {
    aValue = theValue;
    aErr = 1.0 / sqrt(thePrecision);
    indLocal = &theParameters;
    derLocal = &theDerivatives;
    labGlobal = &globalLabels;
    derGlobal = &globalDerivatives;
}

◆ getChi2()

double gbl::GblData::getChi2 ( ) const

Calculate Chi2 contribution.

Returns: (down-weighted) Chi2

Definition at line 173 of file GblData.cc.

                              {
    double aDiff = theValue - thePrediction;
    return aDiff * aDiff * thePrecision * theDownWeight;
}

◆ getLocalData()

void gbl::GblData::getLocalData	(	double &	aValue,
		double &	aWeight,
		std::vector< unsigned int > *&	indLocal,
		std::vector< double > *&	derLocal
	)

Get Data for local fit.

Parameters

[out]	aValue	Value
[out]	aWeight	Weight
[out]	indLocal	List of labels of used (local) fit parameters
[out]	derLocal	List of derivatives for used (local) fit parameters

Definition at line 202 of file GblData.cc.

                                                                          {
 
    aValue = theValue;
    aWeight = thePrecision * theDownWeight;
    indLocal = &theParameters;
    derLocal = &theDerivatives;
}

◆ getResidual()

void gbl::GblData::getResidual	(	double &	aResidual,
		double &	aVariance,
		double &	aDownWeight,
		std::vector< unsigned int > *&	indLocal,
		std::vector< double > *&	derLocal
	)

Get data for residual (and errors).

Parameters

[out]	aResidual	Measurement-Prediction
[out]	aVariance	Variance (of measurement)
[out]	aDownWeight	Down-weighting factor
[out]	indLocal	List of labels of used (local) fit parameters
[out]	derLocal	List of derivatives for used (local) fit parameters

Definition at line 239 of file GblData.cc.

                                      {
    aResidual = theValue - thePrediction;
    aVariance = 1.0 / thePrecision;
    aDownWeight = theDownWeight;
    indLocal = &theParameters;
    derLocal = &theDerivatives;
}

◆ printData()

void gbl::GblData::printData ( ) const

Print data block.

Definition at line 179 of file GblData.cc.

                              {
 
    std::cout << " measurement at label " << theLabel << ": " << theValue
            << ", " << thePrecision << std::endl;
    std::cout << "  param " << theParameters.size() << ":";
    for (unsigned int i = 0; i < theParameters.size(); ++i) {
        std::cout << " " << theParameters[i];
    }
    std::cout << std::endl;
    std::cout << "  deriv " << theDerivatives.size() << ":";
    for (unsigned int i = 0; i < theDerivatives.size(); ++i) {
        std::cout << " " << theDerivatives[i];
    }
    std::cout << std::endl;
}

◆ setDownWeighting()

double gbl::GblData::setDownWeighting ( unsigned int aMethod )

Outlier down weighting with M-estimators (by GblTrajectory::fit).

Parameters

[in] aMethod M-estimator (1: Tukey, 2:Huber, 3:Cauchy)

Definition at line 144 of file GblData.cc.

                                                     {
 
    double aWeight = 1.;
    double scaledResidual = fabs(theValue - thePrediction) * sqrt(thePrecision);
    if (aMethod == 1) // Tukey
            {
        if (scaledResidual < 4.6851) {
            aWeight = (1.0 - 0.045558 * scaledResidual * scaledResidual);
            aWeight *= aWeight;
        } else {
            aWeight = 0.;
        }
    } else if (aMethod == 2) //Huber
            {
        if (scaledResidual >= 1.345) {
            aWeight = 1.345 / scaledResidual;
        }
    } else if (aMethod == 3) //Cauchy
            {
        aWeight = 1.0 / (1.0 + (scaledResidual * scaledResidual / 5.6877));
    }
    theDownWeight = aWeight;
    return aWeight;
}

◆ setPrediction()

void gbl::GblData::setPrediction ( const VVector & aVector )

Calculate prediction for data from fit (by GblTrajectory::fit).

Definition at line 132 of file GblData.cc.

                                                  {
 
    thePrediction = 0.;
    for (unsigned int i = 0; i < theDerivatives.size(); ++i) {
        thePrediction += theDerivatives[i] * aVector(theParameters[i] - 1);
    }
}

Member Data Documentation

◆ globalDerivatives

std::vector<double> gbl::GblData::globalDerivatives

private

Global derivatives.

Definition at line 72 of file GblData.h.

◆ globalLabels

std::vector<int> gbl::GblData::globalLabels

private

Labels for global derivatives.

Definition at line 71 of file GblData.h.

◆ theDerivatives

std::vector<double> gbl::GblData::theDerivatives

private

List of derivatives for fit.

Definition at line 70 of file GblData.h.

◆ theDownWeight

double gbl::GblData::theDownWeight

private

Down-weighting factor (0-1)

Definition at line 67 of file GblData.h.

◆ theLabel

unsigned int gbl::GblData::theLabel

private

Label (of measurements point)

Definition at line 64 of file GblData.h.

◆ theParameters

std::vector<unsigned int> gbl::GblData::theParameters

private

List of fit parameters (with non zero derivatives)

Definition at line 69 of file GblData.h.

◆ thePrecision

double gbl::GblData::thePrecision

private

Precision (1/sigma**2)

Definition at line 66 of file GblData.h.

◆ thePrediction

double gbl::GblData::thePrediction

private

Prediction from fit.

Definition at line 68 of file GblData.h.

◆ theValue

double gbl::GblData::theValue

private

Value (residual)

Definition at line 65 of file GblData.h.

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

genfit/GBL/include/GblData.h
genfit/GBL/src/GblData.cc

Public Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ GblData()

◆ ~GblData()

Member Function Documentation

◆ addDerivatives() [1/3]

◆ addDerivatives() [2/3]

◆ addDerivatives() [3/3]

◆ getAllData()

◆ getChi2()

◆ getLocalData()

◆ getResidual()

◆ printData()

◆ setDownWeighting()

◆ setPrediction()

Member Data Documentation

◆ globalDerivatives

◆ globalLabels

◆ theDerivatives

◆ theDownWeight

◆ theLabel

◆ theParameters

◆ thePrecision

◆ thePrediction

◆ theValue