GlobalDerivatives Class Reference

Class for easier manipulation with global derivatives (and their labels) More...

#include <GlobalDerivatives.h>

Public Member Functions
	GlobalDerivatives (int dim=2)
	Constructor for empty derivative matrix and label vector.
	GlobalDerivatives (const std::pair< std::vector< int >, TMatrixD > &globals)
	Constructor from pair of the vector and the matrix.
	GlobalDerivatives (const std::vector< int > &labels, const TMatrixD &derivs)
	constructor from label vector and derivative matrix
	operator std::pair< std::vector< int >, TMatrixD > ()
	Convenient operator to allow to pass the object directly from RecoHit to globalDerivatives() WARNING: causes zero labels to be removed.
	operator std::vector< int > ()
	Convenient operator to pass only labels (non-zero)
	operator TMatrixD ()
	Covenient operator to pass only derivatives (for non-zero labels)
const std::vector< int > &	getLabels () const
	Get stored lables (includes zeros)
const TMatrixD &	getDerivatives () const
	Return the derivative matrix (includes columns with zero labels)
void	add (const std::pair< std::vector< int >, TMatrixD > &globals)
	Add another set of global labels and derivatives.
void	add (int paramLabel, std::vector< double > dResiduals_dParam)
	Add one parameter - label and the corresponding residual derivative.
void	add (int paramLabel, double drudp)
	Add derivative of local U residual w.r.t.

Static Public Member Functions
static std::pair< std::vector< int >, TMatrixD >	passGlobals (std::pair< std::vector< int >, TMatrixD > globals)
	Static convenient function to remove columns with zero labels (make error in Pede btw.) TODO: refactore interface: return value, params <->

Private Attributes
std::pair< std::vector< int >, TMatrixD >	m_globals {{}, TMatrixD()}
	The global labels and derivatives matrix.

Detailed Description

Class for easier manipulation with global derivatives (and their labels)

Allows to simply add, merge and pass global derivatives matrix and global label vector in RecoHit. For lables = 0 removes given columns from derivative matrix (e.g. save disc space for parameters you do not want to calibrate)

Definition at line 27 of file GlobalDerivatives.h.

Constructor & Destructor Documentation

GlobalDerivatives() [1/3]

GlobalDerivatives ( int  dim = 2 )

inlineexplicit

Constructor for empty derivative matrix and label vector.

Parameters

dim	number of matrix rows (number of local residuals in virtual plane) Usually 2 (also for 1 dim measurement, u/v selected via precision matrix)

Definition at line 32 of file GlobalDerivatives.h.

{m_globals.second.ResizeTo(dim, 0);}

GlobalDerivatives() [2/3]

GlobalDerivatives ( const std::pair< std::vector< int >, TMatrixD > &  globals )

explicit

Constructor from pair of the vector and the matrix.

Definition at line 17 of file GlobalDerivatives.cc.

    {
      m_globals.second.ResizeTo(globals.second);
      m_globals.first = globals.first;
      m_globals.second = globals.second;
    }

GlobalDerivatives() [3/3]

GlobalDerivatives	(	const std::vector< int > &	labels,
		const TMatrixD &	derivs
	)

constructor from label vector and derivative matrix

Definition at line 23 of file GlobalDerivatives.cc.

    {
      m_globals.second.ResizeTo(derivs);
      m_globals.first = labels;
      m_globals.second = derivs;
    }

Member Function Documentation

add() [1/3]

void add

(

const std::pair< std::vector< int >, TMatrixD > &

globals

)

Add another set of global labels and derivatives.

Definition at line 29 of file GlobalDerivatives.cc.

    {
      if (globals.first.empty())
        return;
 
      auto& main = m_globals;
 
      // Create composed matrix of derivatives
      //TODO: check main and globals matrix has the same number of rows
      TMatrixD allDerivatives(main.second.GetNrows(), main.second.GetNcols() + globals.second.GetNcols());
      allDerivatives.Zero();
      allDerivatives.SetSub(0, 0, main.second);
      allDerivatives.SetSub(0, main.second.GetNcols(), globals.second);
 
      // Merge labels
      main.first.insert(main.first.end(), globals.first.begin(), globals.first.end());
      // Update matrix
      main.second.ResizeTo(allDerivatives);
      main.second = allDerivatives;
 
    }

add() [2/3]

void add	(	int	paramLabel,
		double	drudp
	)

Add derivative of local U residual w.r.t.

global parameter Global derivative versus V residual is set to zero. For (u, v) measurement use add(p, {dru/dp, drv/dp}.

Definition at line 61 of file GlobalDerivatives.cc.

    {
      std::vector<double> dResiduals_dParam(m_globals.second.GetNrows(), 0.);
      dResiduals_dParam.at(0) = drudp;
      add(paramLabel, dResiduals_dParam);
    }

add() [3/3]

void add	(	int	paramLabel,
		std::vector< double >	dResiduals_dParam
	)

Add one parameter - label and the corresponding residual derivative.

Parameters

paramLabel	label of the global parameter to calibrate
dResiduals_dParam	vector od derivatives of local residual (U, v) versus global parameter

Definition at line 50 of file GlobalDerivatives.cc.

    {
      int nRows = m_globals.second.GetNrows();
      int nCols = m_globals.second.GetNcols();
 
      m_globals.first.push_back(paramLabel);
      m_globals.second.ResizeTo(nRows, nCols + 1);
      for (int iRow = 0; iRow < nRows; ++iRow) {
        m_globals.second(iRow, nCols) = dResiduals_dParam.at(iRow);
      }
    }

getDerivatives()

const TMatrixD & getDerivatives ( ) const

inline

Return the derivative matrix (includes columns with zero labels)

Definition at line 47 of file GlobalDerivatives.h.

{return m_globals.second;}

getLabels()

const std::vector< int > & getLabels ( ) const

inline

Get stored lables (includes zeros)

Definition at line 45 of file GlobalDerivatives.h.

{return m_globals.first;}

operator std::pair< std::vector< int >, TMatrixD >()

operator std::pair< std::vector< int >, TMatrixD > ( )

inline

Convenient operator to allow to pass the object directly from RecoHit to globalDerivatives() WARNING: causes zero labels to be removed.

Definition at line 39 of file GlobalDerivatives.h.

{return passGlobals(m_globals);}

operator std::vector< int >()

operator std::vector< int > ( )

inline

Convenient operator to pass only labels (non-zero)

Definition at line 41 of file GlobalDerivatives.h.

{return passGlobals(m_globals).first;}

operator TMatrixD()

operator TMatrixD ( )

inline

Covenient operator to pass only derivatives (for non-zero labels)

Definition at line 43 of file GlobalDerivatives.h.

{return passGlobals(m_globals).second;}

passGlobals()

std::pair< std::vector< int >, TMatrixD > passGlobals ( std::pair< std::vector< int >, TMatrixD >  globals )

static

Static convenient function to remove columns with zero labels (make error in Pede btw.) TODO: refactore interface: return value, params <->

Definition at line 67 of file GlobalDerivatives.cc.

    {
      TMatrixD newMatrix(globals.second.GetNrows(), 0);
      std::vector<int> newLabels;
 
      for (unsigned int iOldCol = 0; iOldCol < globals.first.size(); ++iOldCol) {
        auto label = globals.first.at(iOldCol);
        if (label == 0)
          continue;
 
        newLabels.push_back(label);
        newMatrix.ResizeTo(globals.second.GetNrows(), newMatrix.GetNcols() + 1);
        for (int iRow = 0; iRow < globals.second.GetNrows(); ++iRow) {
          newMatrix(iRow, newMatrix.GetNcols() - 1) = globals.second(iRow, iOldCol);
        }
      }
      return {newLabels, newMatrix};
    }

Member Data Documentation

m_globals

std::pair<std::vector<int>, TMatrixD> m_globals {{}, TMatrixD()}

private

The global labels and derivatives matrix.

Definition at line 63 of file GlobalDerivatives.h.

---

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

• alignment/include/GlobalDerivatives.h
• alignment/src/GlobalDerivatives.cc