ThirdHitPredictionFromInvLine Class Reference

#include <ThirdHitPredictionFromInvLine.h>

Classes
class	MappedPoint

Public Types
typedef PixelRecoRange< float >	Range
typedef TkRotation< double >	Rotation

Public Member Functions
void	add (const GlobalPoint &p, double erroriRPhi=1.)
double	chi2 () const
GlobalPoint	crossing (double radius) const
double	curvature () const
double	errorCurvature () const
void	print () const
void	remove (const GlobalPoint &p, double erroriRPhi=1.)
int	size () const
	ThirdHitPredictionFromInvLine (const GlobalPoint &P1, const GlobalPoint &P2, double errorRPhiP1=1., double errorRPhiP2=1.)

Private Types
typedef MappedPoint< double >	PointUV

Private Member Functions
void	add (const ThirdHitPredictionFromInvLine::PointUV &point, double weight)
void	check ()

Private Attributes
bool	hasParameters
int	nPoints
double	theChi2
double	theCurvatureError
double	theCurvatureValue
Rotation	theRotation
long double	theSum
long double	theSumU
long double	theSumUU
long double	theSumUV
long double	theSumV
long double	theSumVV

Detailed Description

Definition at line 19 of file ThirdHitPredictionFromInvLine.h.

Member Typedef Documentation

typedef MappedPoint<double> ThirdHitPredictionFromInvLine::PointUV

private

Definition at line 71 of file ThirdHitPredictionFromInvLine.h.

typedef PixelRecoRange<float> ThirdHitPredictionFromInvLine::Range

Definition at line 24 of file ThirdHitPredictionFromInvLine.h.

typedef TkRotation<double> ThirdHitPredictionFromInvLine::Rotation

Definition at line 23 of file ThirdHitPredictionFromInvLine.h.

Constructor & Destructor Documentation

ThirdHitPredictionFromInvLine::ThirdHitPredictionFromInvLine	(	const GlobalPoint &	P1,
		const GlobalPoint &	P2,
		double	errorRPhiP1 = 1.,
		double	errorRPhiP2 = 1.
	)

Definition at line 20 of file ThirdHitPredictionFromInvLine.cc.

References add(), theRotation, Vector3DBase< T, FrameTag >::unit(), PV3DBase< T, PVType, FrameType >::x(), and PV3DBase< T, PVType, FrameType >::y().

    : nPoints(0), theSum(0.),  theSumU(0.), theSumUU(0.), theSumV(0.), theSumUV(0.), theSumVV(0.),
      hasParameters(false),  theCurvatureValue(0.), theCurvatureError(0.), theChi2(0.)
{
  GlobalVector aX = GlobalVector( P1.x(), P1.y(), 0.).unit();
  GlobalVector aY( -aX.y(), aX.x(), 0.);
  GlobalVector aZ( 0., 0., 1.);
  theRotation = Rotation(aX,aY,aZ);

  add(P1, errorRPhiP1);
  add(P2, errorRPhiP2);
}

Member Function Documentation

void ThirdHitPredictionFromInvLine::add	(	const GlobalPoint &	p,
		double	erroriRPhi = 1.
	)

Definition at line 48 of file ThirdHitPredictionFromInvLine.cc.

References PV3DBase< T, PVType, FrameType >::perp(), sqr(), theRotation, PV3DBase< T, PVType, FrameType >::x(), and PV3DBase< T, PVType, FrameType >::y().

Referenced by counter.Counter::register(), size(), and ThirdHitPredictionFromInvLine().

{
  double weigth = sqr(sqr(p.perp())/errorRPhi);
  add(PointUV(Point2D(p.x(),p.y()), &theRotation), weigth);
}

void ThirdHitPredictionFromInvLine::add ( const ThirdHitPredictionFromInvLine::PointUV &  point, double  weight  )

private

Definition at line 54 of file ThirdHitPredictionFromInvLine.cc.

References check(), hasParameters, nPoints, sqr(), theSum, theSumU, theSumUU, theSumUV, theSumV, theSumVV, ThirdHitPredictionFromInvLine::MappedPoint< T >::u(), and ThirdHitPredictionFromInvLine::MappedPoint< T >::v().

Referenced by counter.Counter::register().

{
  hasParameters = false;
  nPoints++;
  theSum += weigth;
  theSumU += point.u()*weigth;
  theSumUU += sqr(point.u())*weigth;
  theSumV += point.v()*weigth;
  theSumUV +=  point.u()*point.v()*weigth; 
  theSumVV += sqr(point.v())*weigth;
  check();
}

void ThirdHitPredictionFromInvLine::check ( void  )

private

Definition at line 88 of file ThirdHitPredictionFromInvLine.cc.

References patCaloMETCorrections_cff::A, TtFullHadDaughter::B, hasParameters, rho, sqr(), mathSSE::sqrt(), theChi2, theCurvatureError, theCurvatureValue, theSum, theSumU, theSumUU, theSumUV, theSumV, and theSumVV.

Referenced by add(), and remove().

{
  if (hasParameters) return;

  long double D = theSumUU*theSum - theSumU*theSumU;
  long double A = (theSumUV*theSum-theSumU*theSumV)/D;
  long double B = (theSumUU*theSumV-theSumUV*theSumU)/D;
  double rho =  2.*fabs(B)/sqrt(1+sqr(A));
  double sigmaA2 = theSum/D;
  double sigmaB2 = theSumUU/D;

  hasParameters = true;
  theCurvatureError = sqrt( sqr(rho/B)*sigmaB2+ sqr(rho/(1+sqr(A)))*sigmaA2); 
  theCurvatureValue = 2.*fabs(B)/sqrt(1+sqr(A)); 
  theChi2 = theSumVV - 2*A*theSumUV - 2*B*theSumV + 2*A*B*theSumU + B*B*theSum + A*A*theSumUU;
}

double ThirdHitPredictionFromInvLine::chi2 ( void  ) const

inline

Definition at line 40 of file ThirdHitPredictionFromInvLine.h.

References hasParameters, print(), and theChi2.

{  assert(hasParameters); return theChi2; }

GlobalPoint ThirdHitPredictionFromInvLine::crossing

(

double

radius

)

const

Definition at line 35 of file ThirdHitPredictionFromInvLine.cc.

References patCaloMETCorrections_cff::A, TtFullHadDaughter::B, delta, sqr(), mathSSE::sqrt(), theRotation, theSum, theSumU, theSumUU, theSumUV, theSumV, tmp, MetAnalyzer::u1, Basic2DVector< T >::x(), and Basic2DVector< T >::y().

Referenced by size().

{
  double A = -(theSum*theSumUV-theSumU*theSumV)/(sqr(theSumU)-theSum*theSumUU);
  double B = (theSumU*theSumUV - theSumUU*theSumV)/(sqr(theSumU)-theSum*theSumUU);
  double delta =  sqr(2.*A*B) - 4*(1+sqr(A))*(sqr(B)-sqr(1/radius));
  double sqrtdelta = (delta > 0.) ? sqrt(delta) : 0.; 
  double u1 = (-2.*A*B + sqrtdelta)/2./(1+sqr(A)); 
  double v1 = A*u1+B;
  Point2D tmp = PointUV(u1,v1, &theRotation).unmap();
  return GlobalPoint(tmp.x(), tmp.y(), 0.);
}

double ThirdHitPredictionFromInvLine::curvature ( ) const

inline

Definition at line 36 of file ThirdHitPredictionFromInvLine.h.

References hasParameters, and theCurvatureValue.

{ assert(hasParameters); return theCurvatureValue; }

double ThirdHitPredictionFromInvLine::errorCurvature ( ) const

inline

Definition at line 38 of file ThirdHitPredictionFromInvLine.h.

References hasParameters, and theCurvatureError.

{ assert(hasParameters); return theCurvatureError; }

void ThirdHitPredictionFromInvLine::print

(

void

)

const

Definition at line 82 of file ThirdHitPredictionFromInvLine.cc.

References gather_cfg::cout, nPoints, theSumU, theSumUU, theSumUV, and theSumV.

Referenced by chi2().

{
  std::cout <<" nPoints: " << nPoints <<" theSumU: "<< theSumU <<" theSumUU: "<<theSumUU
                                      <<" theSumV: "<< theSumV <<" theSumUV: "<<theSumUV<<std::endl;
}

void ThirdHitPredictionFromInvLine::remove	(	const GlobalPoint &	p,
		double	erroriRPhi = 1.
	)

Definition at line 67 of file ThirdHitPredictionFromInvLine.cc.

References check(), hasParameters, nPoints, PV3DBase< T, PVType, FrameType >::perp(), point, sqr(), theRotation, theSum, theSumU, theSumUU, theSumUV, theSumV, theSumVV, ThirdHitPredictionFromInvLine::MappedPoint< T >::u(), ThirdHitPredictionFromInvLine::MappedPoint< T >::v(), PV3DBase< T, PVType, FrameType >::x(), and PV3DBase< T, PVType, FrameType >::y().

{
  hasParameters = false;
  PointUV point(Point2D(p.x(),p.y()), &theRotation);
  double weigth = sqr(sqr(p.perp())/errorRPhi);
  nPoints--;
  theSum -= weigth;
  theSumU -= point.u()*weigth;
  theSumUU -= sqr(point.u())*weigth;
  theSumV -= point.v()*weigth;
  theSumUV -=  point.u()*point.v()*weigth;
  theSumVV -= sqr(point.v())*weigth;
  check();
}

int ThirdHitPredictionFromInvLine::size ( void  ) const

inline

Definition at line 28 of file ThirdHitPredictionFromInvLine.h.

References add(), crossing(), nPoints, AlCaHLTBitMon_ParallelJobs::p, and TCMET_cfi::radius.

Referenced by ntupleDataFormat._Collection::__iter__(), and ntupleDataFormat._Collection::__len__().

{ return nPoints;}

Member Data Documentation

bool ThirdHitPredictionFromInvLine::hasParameters

private

Definition at line 79 of file ThirdHitPredictionFromInvLine.h.

Referenced by add(), check(), chi2(), curvature(), errorCurvature(), and remove().

int ThirdHitPredictionFromInvLine::nPoints

private

Definition at line 77 of file ThirdHitPredictionFromInvLine.h.

Referenced by add(), print(), remove(), and size().

double ThirdHitPredictionFromInvLine::theChi2

private

Definition at line 80 of file ThirdHitPredictionFromInvLine.h.

Referenced by check(), and chi2().

double ThirdHitPredictionFromInvLine::theCurvatureError

private

Definition at line 80 of file ThirdHitPredictionFromInvLine.h.

Referenced by check(), and errorCurvature().

double ThirdHitPredictionFromInvLine::theCurvatureValue

private

Definition at line 80 of file ThirdHitPredictionFromInvLine.h.

Referenced by check(), and curvature().

Rotation ThirdHitPredictionFromInvLine::theRotation

private

Definition at line 76 of file ThirdHitPredictionFromInvLine.h.

Referenced by add(), crossing(), remove(), and ThirdHitPredictionFromInvLine().

long double ThirdHitPredictionFromInvLine::theSum

private

Definition at line 78 of file ThirdHitPredictionFromInvLine.h.

Referenced by add(), check(), crossing(), and remove().

long double ThirdHitPredictionFromInvLine::theSumU

private

Definition at line 78 of file ThirdHitPredictionFromInvLine.h.

Referenced by add(), check(), crossing(), print(), and remove().

long double ThirdHitPredictionFromInvLine::theSumUU

private

Definition at line 78 of file ThirdHitPredictionFromInvLine.h.

Referenced by add(), check(), crossing(), print(), and remove().

long double ThirdHitPredictionFromInvLine::theSumUV

private

Definition at line 78 of file ThirdHitPredictionFromInvLine.h.

Referenced by add(), check(), crossing(), print(), and remove().

long double ThirdHitPredictionFromInvLine::theSumV

private

Definition at line 78 of file ThirdHitPredictionFromInvLine.h.

Referenced by add(), check(), crossing(), print(), and remove().

long double ThirdHitPredictionFromInvLine::theSumVV

private

Definition at line 78 of file ThirdHitPredictionFromInvLine.h.

Referenced by add(), check(), and remove().