Class to hold one picture of the BPix survey and the local fit. More...

#include <SurveyPxbImageLocalFit.h>

Inheritance diagram for SurveyPxbImageLocalFit:

Public Types
typedef unsigned int	count_t

typedef std::vector< coord_t >	fidpoint_t

typedef std::vector< value_t >	localpars_t

typedef float	pede_deriv_t

typedef int	pede_label_t

Public Types inherited from SurveyPxbImage
typedef Point3DBase< value_t, LocalTag >	coord_t

typedef unsigned int	count_t

typedef unsigned int	id_t

typedef std::pair< id_t, id_t >	idPair_t

enum	location { ll, ul, lr, ur }
	`enum` to help access specific locations on a picture More...

typedef double	value_t

Public Member Functions
void	doFit (const fidpoint_t &fidpointvec)
	Invoke the fit. More...

void	doFit (const fidpoint_t &fidpointvec, const pede_label_t &label1, const pede_label_t &label2)

void	doFit (value_t x1, value_t y1, value_t g1, value_t x2, value_t y2, value_t g2)

value_t	getChi2 ()
	returns the chi^2 of the fit More...

const pede_label_t *	getGlobalDerivsLabelPtr (count_t i)

const pede_deriv_t *	getGlobalDerivsPtr (count_t i)

pede_label_t	getGlobalDerivsSize ()

const pede_deriv_t *	getLocalDerivsPtr (count_t i)

pede_label_t	getLocalDerivsSize ()

localpars_t	getLocalParameters ()
	returns local parameters after fit More...

pede_deriv_t	getResiduum (count_t i)

pede_deriv_t	getSigma (count_t i)

bool	isFitValid ()
	returns validity flag More...

void	setGlobalDerivsToZero (count_t i)

void	setLocalDerivsToZero (count_t i)

	SurveyPxbImageLocalFit ()

	SurveyPxbImageLocalFit (std::istringstream &iss)
	Constructor from istringstream. More...

Public Member Functions inherited from SurveyPxbImage
void	fill (std::istringstream &iss)

const coord_t	getCoord (count_t m)

id_t	getIdFirst ()
	Get `Id` of first module. More...

const idPair_t	getIdPair ()
	Get `Id` pair. More...

id_t	getIdSecond ()
	Get `Id` of second module. More...

value_t	getSigmaX ()
	Get Gaussian error in u direction. More...

value_t	getSigmaY ()
	Get Gaussian error in u direction. More...

bool	isValid ()
	returns validity flag More...

	SurveyPxbImage ()

	SurveyPxbImage (std::istringstream &iss)

Static Public Attributes
static const count_t	nFidpoints = 4

static const count_t	nGlD = 3

static const count_t	nLcD = 4

static const count_t	nLcPars = 4

static const count_t	nMsrmts = 8

Private Member Functions
value_t	dist (const coord_t &p1, const coord_t &p2)
	Distance. More...

void	initFidPoints ()
	Initialise the fiducial points. More...

Private Attributes
localpars_t	a_
	Local parameters. More...

value_t	chi2_
	chi2 of the local fit More...

bool	derivsValidFlag_

std::vector< coord_t >	fidpoints_
	True position of the fiducial points on a sensor wrt. local frame (u,v) More...

bool	fitValidFlag_
	Validity Flag. More...

ROOT::Math::SMatrix< pede_deriv_t, nMsrmts, nGlD >	globalDerivsMatrix_
	Matrix with global derivs. More...

std::vector< pede_label_t >	labelVec1_
	Vector with labels to global derivs. More...

std::vector< pede_label_t >	labelVec2_

ROOT::Math::SMatrix< pede_deriv_t, nMsrmts, nLcD >	localDerivsMatrix_
	Matrix with local derivs. More...

ROOT::Math::SVector< value_t, nMsrmts >	r
	Vector of residuals. More...

Additional Inherited Members
Protected Attributes inherited from SurveyPxbImage
idPair_t	idPair_

bool	isRotated_
	Flag if the image was rotated or not. More...

bool	isValidFlag_
	Validity Flag. More...

std::vector< coord_t >	measurementVec_
	Vector to hold four measurements. More...

value_t	sigma_x_
	Gaussian errors. More...

value_t	sigma_y_

Detailed Description

Class to hold one picture of the BPix survey and the local fit.

Definition at line 14 of file SurveyPxbImageLocalFit.h.

Member Typedef Documentation

◆ count_t

typedef unsigned int SurveyPxbImageLocalFit::count_t

Definition at line 18 of file SurveyPxbImageLocalFit.h.

◆ fidpoint_t

typedef std::vector<coord_t> SurveyPxbImageLocalFit::fidpoint_t

Definition at line 17 of file SurveyPxbImageLocalFit.h.

◆ localpars_t

typedef std::vector<value_t> SurveyPxbImageLocalFit::localpars_t

Definition at line 16 of file SurveyPxbImageLocalFit.h.

◆ pede_deriv_t

typedef float SurveyPxbImageLocalFit::pede_deriv_t

Definition at line 26 of file SurveyPxbImageLocalFit.h.

◆ pede_label_t

typedef int SurveyPxbImageLocalFit::pede_label_t

Definition at line 25 of file SurveyPxbImageLocalFit.h.

Constructor & Destructor Documentation

◆ SurveyPxbImageLocalFit() [1/2]

SurveyPxbImageLocalFit::SurveyPxbImageLocalFit ( )

inline

Definition at line 29 of file SurveyPxbImageLocalFit.h.

       : SurveyPxbImage(), a_(nLcPars, 0), fidpoints_(nFidpoints), fitValidFlag_(false), derivsValidFlag_(false) {
     initFidPoints();
   };

References initFidPoints().

◆ SurveyPxbImageLocalFit() [2/2]

SurveyPxbImageLocalFit::SurveyPxbImageLocalFit ( std::istringstream & iss )

inline

Constructor from istringstream.

Definition at line 35 of file SurveyPxbImageLocalFit.h.

       : SurveyPxbImage(iss), a_(nLcPars, 0), fidpoints_(nFidpoints), fitValidFlag_(false), derivsValidFlag_(false) {
     initFidPoints();
   };

References initFidPoints().

Member Function Documentation

◆ dist()

value_t SurveyPxbImageLocalFit::dist	(	const coord_t &	p1,
		const coord_t &	p2
	)

inlineprivate

Distance.

Definition at line 104 of file SurveyPxbImageLocalFit.h.

                                                      {
     value_t dx = p1.x() - p2.x();
     value_t dy = p1.y() - p2.y();
     return sqrt(dx * dx + dy * dy);
   }

References PVValHelper::dx, PVValHelper::dy, p1, p2, and mathSSE::sqrt().

◆ doFit() [1/3]

void SurveyPxbImageLocalFit::doFit ( const fidpoint_t & fidpointvec )

Invoke the fit.

Definition at line 29 of file SurveyPxbImageLocalFit.cc.

                                                                 {
   fitValidFlag_ = false;
  
   // Calculate gamma of right module w.r.t left modules' fram
   const value_t dxr = fidpointvec[3].x() - fidpointvec[2].x();
   const value_t dyr = fidpointvec[3].y() - fidpointvec[2].y();
   const value_t gammar = atan(dyr / dxr);
   const value_t dxl = fidpointvec[1].x() - fidpointvec[0].x();
   const value_t dyl = fidpointvec[1].y() - fidpointvec[0].y();
   const value_t gammal = atan(dyl / dxl);
   const value_t gamma = gammar - gammal;
   //const value_t gamma = 0.; // Testhalber
   const value_t sing = sin(gamma);
   const value_t cosg = cos(gamma);
 #ifdef DEBUG
   std::cout << "gamma: " << gamma << " gamma left: " << gammal << " gamma right: " << gammar << std::endl;
 #endif
  
 #ifdef DEBUG
   std::cout << "&fidpointvec: " << std::endl;
   for (count_t i = 0; i != fidpointvec.size(); i++)
     std::cout << i << ": " << fidpointvec[i] << std::endl;
   std::cout << "&fidpoints_: " << std::endl;
   for (count_t i = 0; i != fidpoints_.size(); i++)
     std::cout << i << ": " << fidpoints_[i] << std::endl;
 #endif
  
   // Matrix of the local derivatives
   ROOT::Math::SMatrix<value_t, nMsrmts, nLcD> A;  // 8x4
   A(0, 0) = 1.;
   A(0, 1) = 0;
   A(0, 2) = +fidpointvec[0].x();
   A(0, 3) = +fidpointvec[0].y();
   A(1, 0) = 0.;
   A(1, 1) = 1;
   A(1, 2) = +fidpointvec[0].y();
   A(1, 3) = -fidpointvec[0].x();
   A(2, 0) = 1.;
   A(2, 1) = 0;
   A(2, 2) = +fidpointvec[1].x();
   A(2, 3) = +fidpointvec[1].y();
   A(3, 0) = 0.;
   A(3, 1) = 1;
   A(3, 2) = +fidpointvec[1].y();
   A(3, 3) = -fidpointvec[1].x();
   A(4, 0) = 1.;
   A(4, 1) = 0;
   A(4, 2) = +fidpointvec[2].x();
   A(4, 3) = +fidpointvec[2].y();
   A(5, 0) = 0.;
   A(5, 1) = 1;
   A(5, 2) = +fidpointvec[2].y();
   A(5, 3) = -fidpointvec[2].x();
   A(6, 0) = 1.;
   A(6, 1) = 0;
   A(6, 2) = +fidpointvec[3].x();
   A(6, 3) = +fidpointvec[3].y();
   A(7, 0) = 0.;
   A(7, 1) = 1;
   A(7, 2) = +fidpointvec[3].y();
   A(7, 3) = -fidpointvec[3].x();
 #ifdef DEBUG
   std::cout << "A: \n" << A << std::endl;
 #endif
  
   // Covariance matrix
   ROOT::Math::SMatrix<value_t, nMsrmts, nMsrmts> W;  // 8x8
   //ROOT::Math::MatRepSym<value_t,8> W;
   const value_t sigma_x2inv = 1. / (sigma_x_ * sigma_x_);
   const value_t sigma_y2inv = 1. / (sigma_y_ * sigma_y_);
   W(0, 0) = sigma_x2inv;
   W(1, 1) = sigma_y2inv;
   W(2, 2) = sigma_x2inv;
   W(3, 3) = sigma_y2inv;
   W(4, 4) = sigma_x2inv;
   W(5, 5) = sigma_y2inv;
   W(6, 6) = sigma_x2inv;
   W(7, 7) = sigma_y2inv;
 #ifdef DEBUG
   std::cout << "W: \n" << W << std::endl;
 #endif
  
   // Prepare for the fit
   ROOT::Math::SMatrix<value_t, nLcD, nLcD> ATWA;  // 4x4
   ATWA = ROOT::Math::Transpose(A) * W * A;
   //ATWA = ROOT::Math::SimilarityT(A,W); // W muss symmterisch sein -> aendern.
   //std::cout << "ATWA: \n" << ATWA << std::endl;
   ROOT::Math::SMatrix<value_t, nLcD, nLcD> ATWAi;  // 4x4
   int ifail = 0;
   ATWAi = ATWA.Inverse(ifail);
   if (ifail != 0) {  // TODO: ifail Pruefung auf message logger ausgeben statt cout
     std::cout << "Problem singular - fit impossible." << std::endl;
     fitValidFlag_ = false;
     return;
   }
 #ifdef DEBUG
   std::cout << "ATWA-1: \n" << ATWAi << ifail << std::endl;
 #endif
  
   // Measurements
   ROOT::Math::SVector<value_t, nMsrmts> y;  // 8
   y(0) = measurementVec_[0].x();
   y(1) = measurementVec_[0].y();
   y(2) = measurementVec_[1].x();
   y(3) = measurementVec_[1].y();
   y(4) = measurementVec_[2].x();
   y(5) = measurementVec_[2].y();
   y(6) = measurementVec_[3].x();
   y(7) = measurementVec_[3].y();
 #ifdef DEBUG
   std::cout << "y: " << y << std::endl;
 #endif
  
   // do the fit
   ROOT::Math::SVector<value_t, nLcD> a;  // 4
   a = ATWAi * ROOT::Math::Transpose(A) * W * y;
   chi2_ = ROOT::Math::Dot(y, W * y) - ROOT::Math::Dot(a, ROOT::Math::Transpose(A) * W * y);
 #ifdef DEBUG
   std::cout << "a: " << a << " S= " << sqrt(a[2] * a[2] + a[3] * a[3]) << " phi= " << atan(a[3] / a[2])
             << " chi2= " << chi2_ << std::endl;
   std::cout << "A*a: " << A * a << std::endl;
 #endif
   a_.assign(a.begin(), a.end());
  
   // Calculate vector of residuals
   r = y - A * a;
 #ifdef DEBUG
   std::cout << "r: " << r << std::endl;
 #endif
  
   // Fill matrix for global fit with local derivatives
   localDerivsMatrix_(0, 0) = 1.;
   localDerivsMatrix_(0, 1) = 0;
   localDerivsMatrix_(1, 0) = 0.;
   localDerivsMatrix_(1, 1) = 1;
   localDerivsMatrix_(2, 0) = 1.;
   localDerivsMatrix_(2, 1) = 0;
   localDerivsMatrix_(3, 0) = 0.;
   localDerivsMatrix_(3, 1) = 1;
   localDerivsMatrix_(4, 0) = 1.;
   localDerivsMatrix_(4, 1) = 0;
   localDerivsMatrix_(5, 0) = 0.;
   localDerivsMatrix_(5, 1) = 1;
   localDerivsMatrix_(6, 0) = 1.;
   localDerivsMatrix_(6, 1) = 0;
   localDerivsMatrix_(7, 0) = 0.;
   localDerivsMatrix_(7, 1) = 1;
   localDerivsMatrix_(0, 2) = +fidpointvec[0].x() + cosg * fidpoints_[0].x() - sing * fidpoints_[0].y();
   localDerivsMatrix_(0, 3) = +fidpointvec[0].y() + cosg * fidpoints_[0].y() + sing * fidpoints_[0].x();
   localDerivsMatrix_(1, 2) = +fidpointvec[0].y() + cosg * fidpoints_[0].y() + sing * fidpoints_[0].x();
   localDerivsMatrix_(1, 3) = -fidpointvec[0].x() - cosg * fidpoints_[0].x() + sing * fidpoints_[0].y();
   localDerivsMatrix_(2, 2) = +fidpointvec[1].x() + cosg * fidpoints_[1].x() - sing * fidpoints_[1].y();
   localDerivsMatrix_(2, 3) = +fidpointvec[1].y() + cosg * fidpoints_[1].y() + sing * fidpoints_[1].x();
   localDerivsMatrix_(3, 2) = +fidpointvec[1].y() + cosg * fidpoints_[1].y() + sing * fidpoints_[1].x();
   localDerivsMatrix_(3, 3) = -fidpointvec[1].x() - cosg * fidpoints_[1].x() + sing * fidpoints_[1].y();
   localDerivsMatrix_(4, 2) = +fidpointvec[2].x() + cosg * fidpoints_[2].x() - sing * fidpoints_[2].y();
   localDerivsMatrix_(4, 3) = +fidpointvec[2].y() + cosg * fidpoints_[2].y() + sing * fidpoints_[2].x();
   localDerivsMatrix_(5, 2) = +fidpointvec[2].y() + cosg * fidpoints_[2].y() + sing * fidpoints_[2].x();
   localDerivsMatrix_(5, 3) = -fidpointvec[2].x() - cosg * fidpoints_[2].x() + sing * fidpoints_[2].y();
   localDerivsMatrix_(6, 2) = +fidpointvec[3].x() + cosg * fidpoints_[3].x() - sing * fidpoints_[3].y();
   localDerivsMatrix_(6, 3) = +fidpointvec[3].y() + cosg * fidpoints_[3].y() + sing * fidpoints_[3].x();
   localDerivsMatrix_(7, 2) = +fidpointvec[3].y() + cosg * fidpoints_[3].y() + sing * fidpoints_[3].x();
   localDerivsMatrix_(7, 3) = -fidpointvec[3].x() - cosg * fidpoints_[3].x() + sing * fidpoints_[3].y();
  
   // Fill vector with global derivatives and labels (8x3)
   globalDerivsMatrix_(0, 0) = +a(2);
   globalDerivsMatrix_(0, 1) = +a(3);
   globalDerivsMatrix_(0, 2) = +cosg * (a(3) * fidpoints_[0].x() - a(2) * fidpoints_[0].y()) -
                               sing * (a(2) * fidpoints_[0].x() + a(3) * fidpoints_[0].y());
   globalDerivsMatrix_(1, 0) = -a(3);
   globalDerivsMatrix_(1, 1) = +a(2);
   globalDerivsMatrix_(1, 2) = +cosg * (a(2) * fidpoints_[0].x() + a(3) * fidpoints_[0].y()) -
                               sing * (a(2) * fidpoints_[0].y() - a(3) * fidpoints_[0].x());
   globalDerivsMatrix_(2, 0) = +a(2);
   globalDerivsMatrix_(2, 1) = +a(3);
   globalDerivsMatrix_(2, 2) = +cosg * (a(3) * fidpoints_[1].x() - a(2) * fidpoints_[1].y()) -
                               sing * (a(2) * fidpoints_[1].x() + a(3) * fidpoints_[1].y());
   globalDerivsMatrix_(3, 0) = -a(3);
   globalDerivsMatrix_(3, 1) = +a(2);
   globalDerivsMatrix_(3, 2) = +cosg * (a(2) * fidpoints_[1].x() + a(3) * fidpoints_[1].y()) -
                               sing * (a(2) * fidpoints_[1].y() - a(3) * fidpoints_[1].x());
  
   globalDerivsMatrix_(4, 0) = +a(2);
   globalDerivsMatrix_(4, 1) = +a(3);
   globalDerivsMatrix_(4, 2) = +cosg * (a(3) * fidpoints_[2].x() - a(2) * fidpoints_[2].y()) -
                               sing * (a(2) * fidpoints_[2].x() + a(3) * fidpoints_[2].y());
   globalDerivsMatrix_(5, 0) = -a(3);
   globalDerivsMatrix_(5, 1) = +a(2);
   globalDerivsMatrix_(5, 2) = +cosg * (a(2) * fidpoints_[2].x() + a(3) * fidpoints_[2].y()) -
                               sing * (a(2) * fidpoints_[2].y() - a(3) * fidpoints_[2].x());
   globalDerivsMatrix_(6, 0) = +a(2);
   globalDerivsMatrix_(6, 1) = +a(3);
   globalDerivsMatrix_(6, 2) = +cosg * (a(3) * fidpoints_[3].x() - a(2) * fidpoints_[3].y()) -
                               sing * (a(2) * fidpoints_[3].x() + a(3) * fidpoints_[3].y());
   globalDerivsMatrix_(7, 0) = -a(3);
   globalDerivsMatrix_(7, 1) = +a(2);
   globalDerivsMatrix_(7, 2) = +cosg * (a(2) * fidpoints_[3].x() + a(3) * fidpoints_[3].y()) -
                               sing * (a(2) * fidpoints_[3].y() - a(3) * fidpoints_[3].x());
  
   fitValidFlag_ = true;
 }

References MaterialEffects_cfi::A, a, a_, chi2_, funct::cos(), gather_cfg::cout, fidpoints_, fitValidFlag_, CustomPhysics_cfi::gamma, globalDerivsMatrix_, mps_fire::i, localDerivsMatrix_, SurveyPxbImage::measurementVec_, r, SurveyPxbImage::sigma_x_, SurveyPxbImage::sigma_y_, funct::sin(), mathSSE::sqrt(), and y.

Referenced by doFit().

◆ doFit() [2/3]

void SurveyPxbImageLocalFit::doFit	(	const fidpoint_t &	fidpointvec,
		const pede_label_t &	label1,
		const pede_label_t &	label2
	)

Definition at line 15 of file SurveyPxbImageLocalFit.cc.

                                                                {
   labelVec1_.clear();
   labelVec1_.push_back(label1 + 0);
   labelVec1_.push_back(label1 + 1);
   labelVec1_.push_back(label1 + 5);
   labelVec2_.clear();
   labelVec2_.push_back(label2 + 0);
   labelVec2_.push_back(label2 + 1);
   labelVec2_.push_back(label2 + 5);
   doFit(fidpointvec);
 }

References doFit(), bTagCommon_cff::label1, bTagCommon_cff::label2, labelVec1_, and labelVec2_.

◆ doFit() [3/3]

void SurveyPxbImageLocalFit::doFit	(	value_t	x1,
		value_t	y1,
		value_t	g1,
		value_t	x2,
		value_t	y2,
		value_t	g2
	)

Definition at line 231 of file SurveyPxbImageLocalFit.cc.

                                                                                                          {
   // Creating vectors with the global parameters of the two modules
   ROOT::Math::SVector<value_t, 4> mod1, mod2;
   mod1(0) = u1;
   mod1(1) = v1;
   mod1(2) = cos(g1);
   mod1(3) = sin(g1);
   mod2(0) = u2;
   mod2(1) = v2;
   mod2(2) = cos(g2);
   mod2(3) = sin(g2);
   //std::cout << "mod1: " << mod1 << std::endl;
   //std::cout << "mod2: " << mod2 << std::endl;
  
   // Create a matrix for the transformed position of the fidpoints
   ROOT::Math::SMatrix<value_t, 4, 4> M1, M2;
   M1(0, 0) = 1.;
   M1(0, 1) = 0.;
   M1(0, 2) = +fidpoints_[0].x();
   M1(0, 3) = -fidpoints_[0].y();
   M1(1, 0) = 0.;
   M1(1, 1) = 1.;
   M1(1, 2) = +fidpoints_[0].y();
   M1(1, 3) = +fidpoints_[0].x();
   M1(2, 0) = 1.;
   M1(2, 1) = 0.;
   M1(2, 2) = +fidpoints_[1].x();
   M1(2, 3) = -fidpoints_[1].y();
   M1(3, 0) = 0.;
   M1(3, 1) = 1.;
   M1(3, 2) = +fidpoints_[1].y();
   M1(3, 3) = +fidpoints_[1].x();
   M2(0, 0) = 1.;
   M2(0, 1) = 0.;
   M2(0, 2) = +fidpoints_[2].x();
   M2(0, 3) = -fidpoints_[2].y();
   M2(1, 0) = 0.;
   M2(1, 1) = 1.;
   M2(1, 2) = +fidpoints_[2].y();
   M2(1, 3) = +fidpoints_[2].x();
   M2(2, 0) = 1.;
   M2(2, 1) = 0.;
   M2(2, 2) = +fidpoints_[3].x();
   M2(2, 3) = -fidpoints_[3].y();
   M2(3, 0) = 0.;
   M2(3, 1) = 1.;
   M2(3, 2) = +fidpoints_[3].y();
   M2(3, 3) = +fidpoints_[3].x();
  
   //std::cout << "M1:\n" << M1 << std::endl;
   //std::cout << "M2:\n" << M2 << std::endl;
  
   ROOT::Math::SVector<value_t, 4> mod_tr1, mod_tr2;
   mod_tr1 = M1 * mod2;
   mod_tr2 = M2 * mod1;
   //std::cout << "mod_tr1: " << mod_tr1 << std::endl;
   //std::cout << "mod_tr2: " << mod_tr2 << std::endl;
  
   fidpoint_t fidpointvec;
   fidpointvec.push_back(coord_t(mod_tr1(0), mod_tr1(1)));
   fidpointvec.push_back(coord_t(mod_tr1(2), mod_tr1(3)));
   fidpointvec.push_back(coord_t(mod_tr2(0), mod_tr2(1)));
   fidpointvec.push_back(coord_t(mod_tr2(2), mod_tr2(3)));
  
   doFit(fidpointvec);
 }

References funct::cos(), doFit(), fidpoints_, diffTwoXMLs::g1, diffTwoXMLs::g2, funct::sin(), testProducerWithPsetDescEmpty_cfi::u1, and MetAnalyzer::u2.

◆ getChi2()

SurveyPxbImageLocalFit::value_t SurveyPxbImageLocalFit::getChi2 ( )

returns the chi^2 of the fit

Definition at line 306 of file SurveyPxbImageLocalFit.cc.

                                                               {
   if (!fitValidFlag_)
     throw std::logic_error(
         "SurveyPxbImageLocalFit::getChi2(): Fit is not valid. Call doFit(...) before calling this function.");
   return chi2_;
 }

References chi2_, and fitValidFlag_.

◆ getGlobalDerivsLabelPtr()

const pede_label_t* SurveyPxbImageLocalFit::getGlobalDerivsLabelPtr ( count_t i )

inline

Definition at line 58 of file SurveyPxbImageLocalFit.h.

58 { return i < 4 ? &labelVec1_[0] : &labelVec2_[0]; };

References mps_fire::i, labelVec1_, and labelVec2_.

◆ getGlobalDerivsPtr()

const pede_deriv_t* SurveyPxbImageLocalFit::getGlobalDerivsPtr ( count_t i )

inline

Definition at line 57 of file SurveyPxbImageLocalFit.h.

57 { return globalDerivsMatrix_.Array() + i * nGlD; };

References globalDerivsMatrix_, mps_fire::i, and nGlD.

◆ getGlobalDerivsSize()

pede_label_t SurveyPxbImageLocalFit::getGlobalDerivsSize ( )

inline

Definition at line 55 of file SurveyPxbImageLocalFit.h.

55 { return nGlD; };

References nGlD.

◆ getLocalDerivsPtr()

const pede_deriv_t* SurveyPxbImageLocalFit::getLocalDerivsPtr ( count_t i )

inline

Definition at line 56 of file SurveyPxbImageLocalFit.h.

56 { return localDerivsMatrix_.Array() + i * nLcD; };

References mps_fire::i, localDerivsMatrix_, and nLcD.

◆ getLocalDerivsSize()

pede_label_t SurveyPxbImageLocalFit::getLocalDerivsSize ( )

inline

Definition at line 54 of file SurveyPxbImageLocalFit.h.

54 { return nLcD; };

References nLcD.

◆ getLocalParameters()

SurveyPxbImageLocalFit::localpars_t SurveyPxbImageLocalFit::getLocalParameters ( )

returns local parameters after fit

Definition at line 298 of file SurveyPxbImageLocalFit.cc.

                                                                              {
   if (!fitValidFlag_)
     throw std::logic_error(
         "SurveyPxbImageLocalFit::getLocalParameters(): Fit is not valid. Call doFit(...) before calling this "
         "function.");
   return a_;
 }

References a_, and fitValidFlag_.

◆ getResiduum()

pede_deriv_t SurveyPxbImageLocalFit::getResiduum ( count_t i )

inline

Definition at line 59 of file SurveyPxbImageLocalFit.h.

59 { return (pede_deriv_t)r(i); };

References mps_fire::i, and r.

◆ getSigma()

pede_deriv_t SurveyPxbImageLocalFit::getSigma ( count_t i )

inline

Definition at line 60 of file SurveyPxbImageLocalFit.h.

60 { return i % 2 ? sigma_x_ : sigma_y_; };

References mps_fire::i, SurveyPxbImage::sigma_x_, and SurveyPxbImage::sigma_y_.

◆ initFidPoints()

void SurveyPxbImageLocalFit::initFidPoints ( )

inlineprivate

Initialise the fiducial points.

Definition at line 96 of file SurveyPxbImageLocalFit.h.

                        {
     fidpoints_[0] = coord_t(-0.91, -3.30);
     fidpoints_[1] = coord_t(+0.91, -3.30);
     fidpoints_[2] = coord_t(-0.91, +3.30);
     fidpoints_[3] = coord_t(+0.91, +3.30);
   }

References fidpoints_.

Referenced by SurveyPxbImageLocalFit().

◆ isFitValid()

bool SurveyPxbImageLocalFit::isFitValid ( )

inline

returns validity flag

Definition at line 46 of file SurveyPxbImageLocalFit.h.

46 { return fitValidFlag_; };

References fitValidFlag_.

◆ setGlobalDerivsToZero()

void SurveyPxbImageLocalFit::setGlobalDerivsToZero ( count_t i )

Definition at line 320 of file SurveyPxbImageLocalFit.cc.

                                                             {
   if (!(j < nGlD))
     throw std::range_error("SurveyPxbImageLocalFit::setLocalDerivsToZero(j): j out of range.");
   for (count_t i = 0; i != nMsrmts; i++)
     globalDerivsMatrix_(i, j) = 0;
 }

References globalDerivsMatrix_, mps_fire::i, dqmiolumiharvest::j, nGlD, and nMsrmts.

◆ setLocalDerivsToZero()

void SurveyPxbImageLocalFit::setLocalDerivsToZero ( count_t i )

Definition at line 313 of file SurveyPxbImageLocalFit.cc.

                                                            {
   if (!(j < nLcD))
     throw std::range_error("SurveyPxbImageLocalFit::setLocalDerivsToZero(j): j out of range.");
   for (count_t i = 0; i != nMsrmts; i++)
     localDerivsMatrix_(i, j) = 0;
 }

References mps_fire::i, dqmiolumiharvest::j, localDerivsMatrix_, nLcD, and nMsrmts.

Member Data Documentation

◆ a_

localpars_t SurveyPxbImageLocalFit::a_

private

Local parameters.

Definition at line 67 of file SurveyPxbImageLocalFit.h.

Referenced by doFit(), and getLocalParameters().

◆ chi2_

value_t SurveyPxbImageLocalFit::chi2_

private

chi2 of the local fit

Definition at line 87 of file SurveyPxbImageLocalFit.h.

Referenced by doFit(), and getChi2().

◆ derivsValidFlag_

bool SurveyPxbImageLocalFit::derivsValidFlag_

private

Definition at line 93 of file SurveyPxbImageLocalFit.h.

◆ fidpoints_

std::vector<coord_t> SurveyPxbImageLocalFit::fidpoints_

private

True position of the fiducial points on a sensor wrt. local frame (u,v)

Definition at line 84 of file SurveyPxbImageLocalFit.h.

Referenced by doFit(), and initFidPoints().

◆ fitValidFlag_

bool SurveyPxbImageLocalFit::fitValidFlag_

private

Validity Flag.

Definition at line 90 of file SurveyPxbImageLocalFit.h.

Referenced by doFit(), getChi2(), getLocalParameters(), and isFitValid().

◆ globalDerivsMatrix_

ROOT::Math::SMatrix<pede_deriv_t, nMsrmts, nGlD> SurveyPxbImageLocalFit::globalDerivsMatrix_

private

Matrix with global derivs.

Definition at line 74 of file SurveyPxbImageLocalFit.h.

Referenced by doFit(), getGlobalDerivsPtr(), and setGlobalDerivsToZero().

◆ labelVec1_

std::vector<pede_label_t> SurveyPxbImageLocalFit::labelVec1_

private

Vector with labels to global derivs.

Definition at line 81 of file SurveyPxbImageLocalFit.h.

Referenced by doFit(), and getGlobalDerivsLabelPtr().

◆ labelVec2_

std::vector<pede_label_t> SurveyPxbImageLocalFit::labelVec2_

private

Definition at line 81 of file SurveyPxbImageLocalFit.h.

Referenced by doFit(), and getGlobalDerivsLabelPtr().

◆ localDerivsMatrix_

ROOT::Math::SMatrix<pede_deriv_t, nMsrmts, nLcD> SurveyPxbImageLocalFit::localDerivsMatrix_

private

Matrix with local derivs.

Definition at line 78 of file SurveyPxbImageLocalFit.h.

Referenced by doFit(), getLocalDerivsPtr(), and setLocalDerivsToZero().

◆ nFidpoints

const count_t SurveyPxbImageLocalFit::nFidpoints = 4

static

Definition at line 23 of file SurveyPxbImageLocalFit.h.

◆ nGlD

const count_t SurveyPxbImageLocalFit::nGlD = 3

static

Definition at line 19 of file SurveyPxbImageLocalFit.h.

Referenced by getGlobalDerivsPtr(), getGlobalDerivsSize(), and setGlobalDerivsToZero().

◆ nLcD

const count_t SurveyPxbImageLocalFit::nLcD = 4

static

Definition at line 20 of file SurveyPxbImageLocalFit.h.

Referenced by getLocalDerivsPtr(), getLocalDerivsSize(), and setLocalDerivsToZero().

◆ nLcPars

const count_t SurveyPxbImageLocalFit::nLcPars = 4

static

Definition at line 22 of file SurveyPxbImageLocalFit.h.

◆ nMsrmts

const count_t SurveyPxbImageLocalFit::nMsrmts = 8

static

Definition at line 21 of file SurveyPxbImageLocalFit.h.

Referenced by MillePedeAlignmentAlgorithm::addPxbSurvey(), setGlobalDerivsToZero(), and setLocalDerivsToZero().

◆ r

ROOT::Math::SVector<value_t, nMsrmts> SurveyPxbImageLocalFit::r

private

Vector of residuals.

Definition at line 70 of file SurveyPxbImageLocalFit.h.

Referenced by doFit(), and getResiduum().

Public Types

Public Member Functions

Static Public Attributes

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Member Typedef Documentation

◆ count_t

◆ fidpoint_t

◆ localpars_t

◆ pede_deriv_t

◆ pede_label_t

Constructor & Destructor Documentation

◆ SurveyPxbImageLocalFit() [1/2]

◆ SurveyPxbImageLocalFit() [2/2]

Member Function Documentation

◆ dist()

◆ doFit() [1/3]

◆ doFit() [2/3]

◆ doFit() [3/3]

◆ getChi2()

◆ getGlobalDerivsLabelPtr()

◆ getGlobalDerivsPtr()

◆ getGlobalDerivsSize()

◆ getLocalDerivsPtr()

◆ getLocalDerivsSize()

◆ getLocalParameters()

◆ getResiduum()

◆ getSigma()

◆ initFidPoints()

◆ isFitValid()

◆ setGlobalDerivsToZero()

◆ setLocalDerivsToZero()

Member Data Documentation

◆ a_

◆ chi2_

◆ derivsValidFlag_

◆ fidpoints_

◆ fitValidFlag_

◆ globalDerivsMatrix_

◆ labelVec1_

◆ labelVec2_

◆ localDerivsMatrix_

◆ nFidpoints

◆ nGlD

◆ nLcD

◆ nLcPars

◆ nMsrmts

◆ r