#include <ParabolaFit.h>

Classes
struct	Column

struct	Point

struct	Result

Public Member Functions
void	addPoint (double x, double y)

void	addPoint (double x, double y, double weight)

double	chi2 () const

int	dof () const

void	fixParC (double val)

double	parA () const

	ParabolaFit ()

double	parB () const

double	parC () const

const Result &	result (bool doErrors) const

void	skipErrorCalculationByDefault ()

double	varAA () const

double	varAB () const

double	varAC () const

double	varBB () const

double	varBC () const

double	varCC () const

Private Member Functions
double	det (const Column &c1, const Column &c2) const

double	det (const Column &c1, const Column &c2, const Column &c3) const

double	fun (double x) const

Private Attributes
bool	doErr

bool	hasErrors

bool	hasFixedParC

bool	hasValues

bool	hasWeights

std::vector< Point >	points

Result	theResult

Detailed Description

The parabola fit y = parA + parB * x + parC * x*x see: R.L. Gluckstern, NIM 24 (1963) 381

Definition at line 11 of file ParabolaFit.h.

Constructor & Destructor Documentation

◆ ParabolaFit()

ParabolaFit::ParabolaFit ( )

inline

Definition at line 17 of file ParabolaFit.h.

17 : doErr(true), hasFixedParC(false), hasValues(false), hasErrors(false), hasWeights(true) {}

Member Function Documentation

◆ addPoint() [1/2]

void ParabolaFit::addPoint	(	double	x,
		double	y
	)

Definition at line 9 of file ParabolaFit.cc.

                                              {
   hasWeights = false;
   addPoint(x, y, 1.);
 }

Referenced by ConformalMappingFit::ConformalMappingFit().

◆ addPoint() [2/2]

void ParabolaFit::addPoint	(	double	x,
		double	y,
		double	weight
	)

Definition at line 14 of file ParabolaFit.cc.

                                                        {
   hasValues = false;
   hasErrors = false;
   Point p = {x, y, w};
   points.push_back(p);
 }

References AlCaHLTBitMon_ParallelJobs::p, HLT_FULL_cff::points, and w.

◆ chi2()

double ParabolaFit::chi2 ( void ) const

Definition at line 111 of file ParabolaFit.cc.

                                {
   double mychi2 = 0.;
   for (vector<Point>::const_iterator ip = points.begin(); ip != points.end(); ip++) {
     mychi2 += ip->w * sqr(ip->y - fun(ip->x));
   }
   return mychi2;
 }

References HLT_FULL_cff::points, and sqr().

Referenced by ConformalMappingFit::chi2().

◆ det() [1/2]

double ParabolaFit::det	(	const Column &	c1,
		const Column &	c2
	)		const

private

Definition at line 134 of file ParabolaFit.cc.

134 { return c1.r1 * c2.r2 - c1.r2 * c2.r1; }

References alignmentValidation::c1, ParabolaFit::Column::r1, and ParabolaFit::Column::r2.

◆ det() [2/2]

double ParabolaFit::det	(	const Column &	c1,
		const Column &	c2,
		const Column &	c3
	)		const

private

Definition at line 129 of file ParabolaFit.cc.

                                                                                   {
   return c1.r1 * c2.r2 * c3.r3 + c2.r1 * c3.r2 * c1.r3 + c3.r1 * c1.r2 * c2.r3 - c1.r3 * c2.r2 * c3.r1 -
          c2.r3 * c3.r2 * c1.r1 - c3.r3 * c1.r2 * c2.r1;
 }

References alignmentValidation::c1, ParabolaFit::Column::r1, ParabolaFit::Column::r2, and ParabolaFit::Column::r3.

◆ dof()

int ParabolaFit::dof ( ) const

Definition at line 121 of file ParabolaFit.cc.

                            {
   int nPar = 3;
   if (hasFixedParC)
     nPar--;
   int nDof = points.size() - nPar;
   return (nDof > 0) ? nDof : 0;
 }

References fftjetvertexadder_cfi::nDof, HcalResponse_cfi::nPar, and HLT_FULL_cff::points.

◆ fixParC()

void ParabolaFit::fixParC ( double val )

inline

Definition at line 23 of file ParabolaFit.h.

                            {
     hasFixedParC = true;
     theResult.parC = val;
   }

References hasFixedParC, ParabolaFit::Result::parC, theResult, and heppy_batch::val.

Referenced by ConformalMappingFit::fixImpactParmaeter().

◆ fun()

double ParabolaFit::fun ( double x ) const

private

Definition at line 119 of file ParabolaFit.cc.

119 { return theResult.parA + theResult.parB * x + theResult.parC * x * x; }

◆ parA()

double ParabolaFit::parA ( ) const

inline

Definition at line 30 of file ParabolaFit.h.

                       {
     if (!hasValues)
       result(doErr);
     return theResult.parA;
   }

References doErr, hasValues, ParabolaFit::Result::parA, result(), and theResult.

Referenced by ConformalMappingFit::charge(), and ConformalMappingFit::curvature().

◆ parB()

double ParabolaFit::parB ( ) const

inline

Definition at line 35 of file ParabolaFit.h.

                       {
     if (!hasValues)
       result(doErr);
     return theResult.parB;
   }

References doErr, hasValues, ParabolaFit::Result::parB, result(), and theResult.

Referenced by ConformalMappingFit::directionPhi().

◆ parC()

double ParabolaFit::parC ( ) const

inline

Definition at line 40 of file ParabolaFit.h.

                       {
     if (!hasValues)
       result(doErr);
     return theResult.parC;
   }

References doErr, hasValues, ParabolaFit::Result::parC, result(), and theResult.

Referenced by ConformalMappingFit::impactParameter().

◆ result()

const ParabolaFit::Result & ParabolaFit::result ( bool doErrors ) const

Definition at line 21 of file ParabolaFit.cc.

                                                                 {
   if (hasErrors)
     return theResult;
   if (hasValues && !doErrors)
     return theResult;
  
   double F0, F1, F2, F3, F4, F0y, F1y, F2y;
   F0 = F1 = F2 = F3 = F4 = F0y = F1y = F2y = 0.;
  
   typedef vector<Point>::const_iterator IT;
   for (IT ip = points.begin(); ip != points.end(); ip++) {
     double pow;
     double x = ip->x;
     double y = ip->y;
     double w = ip->w;
  
     F0 += w;
     F0y += w * y;
     F1 += w * x;
     F1y += w * x * y;
     pow = x * x;
     F2 += w * pow;
     F2y += w * pow * y;  // x^2
     pow *= x;
     F3 += w * pow;  // x^3
     pow *= x;
     F4 += w * pow;  // x^4
   }
  
   Column cA = {F0, F1, F2};
   Column cB = {F1, F2, F3};
   Column cC = {F2, F3, F4};
   Column cY = {F0y, F1y, F2y};
  
   double det0 = det(cA, cB, cC);
  
   if (!hasFixedParC) {
     theResult.parA = det(cY, cB, cC) / det0;
     theResult.parB = det(cA, cY, cC) / det0;
     theResult.parC = det(cA, cB, cY) / det0;
   } else {
     Column cCY = {F0y - theResult.parC * F2, F1y - theResult.parC * F3, F2y - theResult.parC * F4};
     double det0C = det(cA, cB);
     theResult.parA = det(cCY, cB) / det0C;
     theResult.parB = det(cA, cCY) / det0C;
   }
  
   //  std::cout <<" result: A="<<theResult.parA<<" B="<<theResult.parB<<" C="<<theResult.parC<<endl;
   double vAA, vBB, vCC, vAB, vAC, vBC;
   vAA = vBB = vCC = vAB = vAC = vBC = 0.;
  
   hasValues = true;
   if (!doErrors)
     return theResult;
  
   if (!hasWeights && dof() > 0) {
     //     cout <<" CHI2: " << chi2() <<" DOF: " << dof() << endl;
     double scale_w = 1. / chi2() / dof();
     for (IT ip = points.begin(); ip != points.end(); ip++)
       ip->w *= scale_w;
     //     cout <<" CHI2: " << chi2() <<" DOF: " << dof() << endl;
   }
  
   for (IT ip = points.begin(); ip != points.end(); ip++) {
     double w = ip->w;
     Column cX = {1., ip->x, sqr(ip->x)};
  
     double dXBC = det(cX, cB, cC);
     double dAXC = det(cA, cX, cC);
     double dABX = det(cA, cB, cX);
  
     vAA += w * sqr(dXBC);
     vBB += w * sqr(dAXC);
     vCC += w * sqr(dABX);
     vAB += w * dXBC * dAXC;
     vAC += w * dXBC * dABX;
     vBC += w * dAXC * dABX;
   }
  
   theResult.varAA = vAA / sqr(det0);
   theResult.varBB = vBB / sqr(det0);
   theResult.varCC = vCC / sqr(det0);
   theResult.varAB = vAB / sqr(det0);
   theResult.varAC = vAC / sqr(det0);
   theResult.varBC = vBC / sqr(det0);
  
   hasErrors = true;
   return theResult;
 }

References hltPixelTracks_cff::chi2, HLT_FULL_cff::points, funct::pow(), sqr(), and w.

Referenced by parA(), parB(), parC(), varAA(), varAB(), varAC(), varBB(), varBC(), and varCC().

◆ skipErrorCalculationByDefault()

void ParabolaFit::skipErrorCalculationByDefault ( )

inline

Definition at line 22 of file ParabolaFit.h.

22 { doErr = false; }

References doErr.

Referenced by ConformalMappingFit::skipErrorCalculation().

◆ varAA()

double ParabolaFit::varAA ( ) const

inline

Definition at line 45 of file ParabolaFit.h.

                        {
     if (!hasErrors)
       result(true);
     return theResult.varAA;
   }

References hasErrors, result(), theResult, and ParabolaFit::Result::varAA.

Referenced by ConformalMappingFit::curvature().

◆ varAB()

double ParabolaFit::varAB ( ) const

inline

Definition at line 60 of file ParabolaFit.h.

                        {
     if (!hasErrors)
       result(true);
     return theResult.varAB;
   }

References hasErrors, result(), theResult, and ParabolaFit::Result::varAB.

◆ varAC()

double ParabolaFit::varAC ( ) const

inline

Definition at line 65 of file ParabolaFit.h.

                        {
     if (!hasErrors)
       result(true);
     return theResult.varAC;
   }

References hasErrors, result(), theResult, and ParabolaFit::Result::varAC.

◆ varBB()

double ParabolaFit::varBB ( ) const

inline

Definition at line 50 of file ParabolaFit.h.

                        {
     if (!hasErrors)
       result(true);
     return theResult.varBB;
   }

References hasErrors, result(), theResult, and ParabolaFit::Result::varBB.

Referenced by ConformalMappingFit::directionPhi().

◆ varBC()

double ParabolaFit::varBC ( ) const

inline

Definition at line 70 of file ParabolaFit.h.

                        {
     if (!hasErrors)
       result(true);
     return theResult.varBC;
   }

References hasErrors, result(), theResult, and ParabolaFit::Result::varBC.

◆ varCC()

double ParabolaFit::varCC ( ) const

inline

Definition at line 55 of file ParabolaFit.h.

                        {
     if (!hasErrors)
       result(true);
     return theResult.varCC;
   }

References hasErrors, result(), theResult, and ParabolaFit::Result::varCC.

Referenced by ConformalMappingFit::impactParameter().

Member Data Documentation

◆ doErr

bool ParabolaFit::doErr

private

Definition at line 97 of file ParabolaFit.h.

Referenced by parA(), parB(), parC(), and skipErrorCalculationByDefault().

◆ hasErrors

bool ParabolaFit::hasErrors

mutableprivate

Definition at line 98 of file ParabolaFit.h.

Referenced by varAA(), varAB(), varAC(), varBB(), varBC(), and varCC().

◆ hasFixedParC

bool ParabolaFit::hasFixedParC

private

Definition at line 97 of file ParabolaFit.h.

Referenced by fixParC().

◆ hasValues

bool ParabolaFit::hasValues

mutableprivate

Definition at line 98 of file ParabolaFit.h.

Referenced by parA(), parB(), and parC().

◆ hasWeights

bool ParabolaFit::hasWeights

private

Definition at line 99 of file ParabolaFit.h.

◆ points

std::vector<Point> ParabolaFit::points

private

Definition at line 96 of file ParabolaFit.h.

◆ theResult

Result ParabolaFit::theResult

mutableprivate

Definition at line 100 of file ParabolaFit.h.

Referenced by fixParC(), parA(), parB(), parC(), varAA(), varAB(), varAC(), varBB(), varBC(), and varCC().

Classes

Public Member Functions

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ ParabolaFit()

Member Function Documentation

◆ addPoint() [1/2]

◆ addPoint() [2/2]

◆ chi2()

◆ det() [1/2]

◆ det() [2/2]

◆ dof()

◆ fixParC()

◆ fun()

◆ parA()

◆ parB()

◆ parC()

◆ result()

◆ skipErrorCalculationByDefault()

◆ varAA()

◆ varAB()

◆ varAC()

◆ varBB()

◆ varBC()

◆ varCC()

Member Data Documentation

◆ doErr

◆ hasErrors

◆ hasFixedParC

◆ hasValues

◆ hasWeights

◆ points

◆ theResult