ThirdHitPredictionFromInvParabola Class Reference

#include <ThirdHitPredictionFromInvParabola.h>

Public Types
typedef Basic2DVector< Scalar >	Point2D
typedef PixelRecoRange< float >	Range
typedef PixelRecoRange< Scalar >	RangeD
typedef TkRotation2D< Scalar >	Rotation
using	Scalar = double

Public Member Functions
void	init (const GlobalPoint &P1, const GlobalPoint &P2, Scalar ip, Scalar curv)
void	init (Scalar x1, Scalar y1, Scalar x2, Scalar y2, Scalar ip, Scalar curv)
Range	operator() (Scalar radius, int charge) const
Range	operator() (Scalar radius) const
Range	rangeRPhi (Scalar radius, int charge) const
Range	rangeRPhi (Scalar radius) const
	ThirdHitPredictionFromInvParabola ()
	ThirdHitPredictionFromInvParabola (Scalar x1, Scalar y1, Scalar x2, Scalar y2, Scalar ip, Scalar curv, Scalar tolerance)
	ThirdHitPredictionFromInvParabola (const GlobalPoint &P1, const GlobalPoint &P2, Scalar ip, Scalar curv, Scalar tolerance)

Private Member Functions
Scalar	coeffA (Scalar impactParameter) const
Scalar	coeffB (Scalar impactParameter) const
void	findPointAtCurve (Scalar radius, Scalar ip, Scalar &u, Scalar &v) const
Scalar	ipFromCurvature (Scalar curvature, bool pos) const
Scalar	predV (Scalar u, Scalar ip) const
Point2D	transform (Point2D const &p) const
Point2D	transformBack (Point2D const &p) const

Private Attributes
Scalar	dv
bool	emptyM
bool	emptyP
Scalar	overDu
Scalar	pv
Scalar	su
RangeD	theIpRangeMinus
RangeD	theIpRangePlus
Rotation	theRotation
Scalar	theTolerance
Scalar	u1u2

Detailed Description

Definition at line 27 of file ThirdHitPredictionFromInvParabola.h.

Member Typedef Documentation

typedef Basic2DVector<Scalar> ThirdHitPredictionFromInvParabola::Point2D

Definition at line 34 of file ThirdHitPredictionFromInvParabola.h.

typedef PixelRecoRange<float> ThirdHitPredictionFromInvParabola::Range

Definition at line 32 of file ThirdHitPredictionFromInvParabola.h.

typedef PixelRecoRange<Scalar> ThirdHitPredictionFromInvParabola::RangeD

Definition at line 33 of file ThirdHitPredictionFromInvParabola.h.

typedef TkRotation2D<Scalar> ThirdHitPredictionFromInvParabola::Rotation

Definition at line 31 of file ThirdHitPredictionFromInvParabola.h.

using ThirdHitPredictionFromInvParabola::Scalar = double

Definition at line 30 of file ThirdHitPredictionFromInvParabola.h.

Constructor & Destructor Documentation

ThirdHitPredictionFromInvParabola::ThirdHitPredictionFromInvParabola ( )

inline

Definition at line 36 of file ThirdHitPredictionFromInvParabola.h.

{}

ThirdHitPredictionFromInvParabola::ThirdHitPredictionFromInvParabola ( Scalar  x1, Scalar  y1, Scalar  x2, Scalar  y2, Scalar  ip, Scalar  curv, Scalar  tolerance  )

inline

Definition at line 37 of file ThirdHitPredictionFromInvParabola.h.

References funct::abs(), and init().

                                     : theTolerance(tolerance)
  {
    init(x1,y1,x2,y2,ip,std::abs(curv));
  }

ThirdHitPredictionFromInvParabola::ThirdHitPredictionFromInvParabola	(	const GlobalPoint &	P1,
		const GlobalPoint &	P2,
		Scalar	ip,
		Scalar	curv,
		Scalar	tolerance
	)

Definition at line 20 of file ThirdHitPredictionFromInvParabola.cc.

References funct::abs(), and init().

  : theTolerance(tolerance)
{
  init(P1,P2,ip,std::abs(curv));
}

Member Function Documentation

ThirdHitPredictionFromInvParabola::Scalar ThirdHitPredictionFromInvParabola::coeffA ( Scalar  impactParameter ) const

inlineprivate

Definition at line 98 of file ThirdHitPredictionFromInvParabola.h.

References EnergyCorrector::c, overDu, pv, and u1u2.

Referenced by findPointAtCurve(), and predV().

{
  auto c = -pv*overDu;
  return c - u1u2*impactParameter;
}

ThirdHitPredictionFromInvParabola::Scalar ThirdHitPredictionFromInvParabola::coeffB ( Scalar  impactParameter ) const

inlineprivate

Definition at line 105 of file ThirdHitPredictionFromInvParabola.h.

References EnergyCorrector::c, dv, overDu, and su.

Referenced by findPointAtCurve(), and predV().

{
  auto c = dv*overDu;
  return c - su*impactParameter;
}

void ThirdHitPredictionFromInvParabola::findPointAtCurve ( Scalar  radius, Scalar  ip, Scalar &  u, Scalar &  v  ) const

inlineprivate

Definition at line 128 of file ThirdHitPredictionFromInvParabola.h.

References funct::A, a, alpha, EnergyCorrector::c, coeffA(), coeffB(), delta, alignCSCRings::r, mathSSE::sqrt(), and findQualityFiles::v.

Referenced by rangeRPhi().

{
  //
  // assume u=(1-alpha^2/2)/r v=alpha/r
  // solve qudratic equation neglecting aplha^4 term
  //
  Scalar A = coeffA(ip);
  Scalar B = coeffB(ip);

  // Scalar overR = 1./r;
  Scalar ipOverR = ip/r; // *overR;

  Scalar a = 0.5*B+ipOverR;
  Scalar c = -B+A*r-ipOverR;

  Scalar delta = 1-4*a*c;
  // Scalar sqrtdelta = (delta > 0) ? std::sqrt(delta) : 0.;
  Scalar sqrtdelta = std::sqrt(delta);
  //  Scalar alpha = (-1+sqrtdelta)/(2*a);
  Scalar alpha = (-2*c)/(1+sqrtdelta);

  v = alpha;  // *overR
  Scalar d2 = 1. - v*v;  // overR*overR - v*v
  // u = (d2 > 0) ? std::sqrt(d2) : 0.;
  u = std::sqrt(d2);

  // u,v not rotated! not multiplied by 1/r
}

void ThirdHitPredictionFromInvParabola::init ( const GlobalPoint &  P1, const GlobalPoint &  P2, Scalar  ip, Scalar  curv  )

inline

Definition at line 59 of file ThirdHitPredictionFromInvParabola.h.

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

Referenced by ThirdHitPredictionFromInvParabola().

                                                                                      {
     init( P1.x(), P1.y(), P2.x(),P2.y(),ip,curv);
  }

void ThirdHitPredictionFromInvParabola::init	(	Scalar	x1,
		Scalar	y1,
		Scalar	x2,
		Scalar	y2,
		Scalar	ip,
		Scalar	curv
	)

Definition at line 28 of file ThirdHitPredictionFromInvParabola.cc.

References funct::abs(), dv, PixelRecoRange< T >::empty(), emptyM, emptyP, PixelRecoRange< T >::intersection(), ipFromCurvature(), bookConverter::max, min(), overDu, p1, p2, pv, su, theIpRangeMinus, theIpRangePlus, theRotation, transform(), u1u2, Basic2DVector< T >::x(), and Basic2DVector< T >::y().

                                                                                                               {
//  GlobalVector aX = GlobalVector( P2.x()-P1.x(), P2.y()-P1.y(), 0.).unit();
 
  Point2D p1(x1,y1);
  Point2D p2(x2,y2);
  theRotation = Rotation(p1);
  p1 = transform(p1);  // (1./P1.xy().mag(),0); 
  p2 = transform(p2);

 
  u1u2 = p1.x()*p2.x();
  overDu = 1./(p2.x() - p1.x());
  pv = p1.y()*p2.x() - p2.y()*p1.x();
  dv = p2.y() - p1.y();
  su = p2.x() + p1.x();

  ip = std::abs(ip);
  RangeD ipRange(-ip, ip); 

  
  Scalar ipIntyPlus = ipFromCurvature(0.,true);
  Scalar ipCurvPlus = ipFromCurvature(curv, true);
  Scalar ipCurvMinus = ipFromCurvature(curv, false);

  
  RangeD ipRangePlus(std::min(ipIntyPlus, ipCurvPlus),std::max(ipIntyPlus, ipCurvPlus)); 
  RangeD ipRangeMinus(std::min(-ipIntyPlus, ipCurvMinus),std::max(-ipIntyPlus, ipCurvMinus));

  theIpRangePlus  = ipRangePlus.intersection(ipRange);
  theIpRangeMinus = ipRangeMinus.intersection(ipRange);

  emptyP =  theIpRangePlus.empty();
  emptyM =  theIpRangeMinus.empty();

}

ThirdHitPredictionFromInvParabola::Scalar ThirdHitPredictionFromInvParabola::ipFromCurvature ( Scalar  curvature, bool  pos  ) const

inlineprivate

Definition at line 112 of file ThirdHitPredictionFromInvParabola.h.

References overDu, pv, and u1u2.

Referenced by init().

{
  Scalar overU1u2 = 1./u1u2;
  Scalar inInf = -pv*overDu*overU1u2;
  return (pos? inInf : -inInf) -curvature*overU1u2*0.5;
}

Range ThirdHitPredictionFromInvParabola::operator() ( Scalar  radius, int  charge  ) const

inline

Definition at line 47 of file ThirdHitPredictionFromInvParabola.h.

References rangeRPhi().

{ return rangeRPhi(radius,charge); } 

Range ThirdHitPredictionFromInvParabola::operator() ( Scalar  radius ) const

inline

Definition at line 49 of file ThirdHitPredictionFromInvParabola.h.

References rangeRPhi().

{ return rangeRPhi(radius); } 

ThirdHitPredictionFromInvParabola::Scalar ThirdHitPredictionFromInvParabola::predV ( Scalar  u, Scalar  ip  ) const

inlineprivate

Definition at line 121 of file ThirdHitPredictionFromInvParabola.h.

References EnergyCorrector::c, coeffA(), and coeffB().

{
  auto c =  -( coeffA(ip) - coeffB(ip*u) - ip*u*u);
  return c;
}

ThirdHitPredictionFromInvParabola::Range ThirdHitPredictionFromInvParabola::rangeRPhi	(	Scalar	radius,
		int	charge
	)		const

Definition at line 67 of file ThirdHitPredictionFromInvParabola.cc.

References Basic2DVector< T >::barePhi(), PixelRecoRange< T >::empty(), findPointAtCurve(), i, PixelRecoRange< T >::intersection(), PixelRecoRange< T >::max(), PixelRecoRange< T >::min(), diffTwoXMLs::r1, diffTwoXMLs::r2, TkRotation2D< T >::rotateBack(), std::swap(), theIpRangeMinus, theIpRangePlus, theRotation, theTolerance, and findQualityFiles::v.

Referenced by operator()().

{
  bool pos =  icharge>0;

  RangeD ip = (pos) ? theIpRangePlus : theIpRangeMinus;


  //  it will vectorize with gcc 4.7 (with -O3 -fno-math-errno)
  // change sign as intersect assume -ip for negative charge...
  Scalar ipv[2]={(pos)? ip.min() : -ip.min() ,(pos)? ip.max() : -ip.max()};
  Scalar u[2], v[2];
  for (int i=0; i!=2; ++i)
    findPointAtCurve(radius,ipv[i],u[i],v[i]);

  // 
  Scalar phi1 = theRotation.rotateBack(Point2D(u[0],v[0])).barePhi();
  Scalar phi2 = phi1+(v[1]-v[0]); 
  
  if (ip.empty()) {
    Range r1(phi1*radius-theTolerance, phi1*radius+theTolerance); 
    Range r2(phi2*radius-theTolerance, phi2*radius+theTolerance); 
    return r1.intersection(r2);
  }

  if (phi2<phi1) std::swap(phi1, phi2); 
  return Range(radius*phi1-theTolerance, radius*phi2+theTolerance);
  

}

ThirdHitPredictionFromInvParabola::Range ThirdHitPredictionFromInvParabola::rangeRPhi

(

Scalar

radius

)

const

Definition at line 148 of file ThirdHitPredictionFromInvParabola.cc.

References relativeConstraints::empty, emptyM, emptyP, findPointAtCurve(), helper::getRange(), i, PixelRecoRange< T >::intersection(), PixelRecoRange< T >::max(), bookConverter::max, PixelRecoRange< T >::min(), min(), diffTwoXMLs::r1, diffTwoXMLs::r2, theIpRangeMinus, theIpRangePlus, theRotation, theTolerance, findQualityFiles::v, TkRotation2D< T >::x(), and TkRotation2D< T >::y().

{

  auto getRange = [&](Scalar phi1, Scalar phi2, bool empty)->RangeD {
    
    if (empty) {
      RangeD r1(phi1*radius-theTolerance, phi1*radius+theTolerance); 
      RangeD r2(phi2*radius-theTolerance, phi2*radius+theTolerance); 
      return r1.intersection(r2);
    }
    
    return RangeD(radius*std::min(phi1,phi2)-theTolerance, radius*std::max(phi1,phi2)+theTolerance);
  };


  //  it will vectorize with gcc 4.7 (with -O3 -fno-math-errno)
  // change sign as intersect assume -ip for negative charge...
  Scalar ipv[4]={theIpRangePlus.min(), -theIpRangeMinus.min(), theIpRangePlus.max(),  -theIpRangeMinus.max()};
  Scalar u[4], v[4];
  for (int i=0; i<4; ++i)
    findPointAtCurve(radius,ipv[i],u[i],v[i]);

  // 
  auto xr = theRotation.x();
  auto yr = theRotation.y();

  Scalar phi1[2],phi2[2];
  for (int i=0; i<2; ++i) {
    auto x =  xr[0]*u[i] + yr[0]*v[i];
    auto y =  xr[1]*u[i] + yr[1]*v[i];
    phi1[i] = atan2_fast(y,x);
    phi2[i] = phi1[i]+(v[i+2]-v[i]); 
  }


  return getRange(phi1[1],phi2[1],emptyM).sum(getRange(phi1[0],phi2[0],emptyP));

  /*
  Scalar phi1P = theRotation.rotateBack(Point2D(u[0],v[0])).barePhi();
  Scalar phi2P= phi1P+(v[2]-v[0]); 

  Scalar phi1M = theRotation.rotateBack(Point2D(u[1],v[1])).barePhi();
  Scalar phi2M = phi1M+(v[3]-v[1]); 


  return getRange(phi1M,phi2M,emptyM).sum(getRange(phi1P,phi2P,emptyP));

  */

}

Point2D ThirdHitPredictionFromInvParabola::transform ( Point2D const &  p ) const

inlineprivate

Definition at line 73 of file ThirdHitPredictionFromInvParabola.h.

References Basic2DVector< T >::mag2(), TkRotation2D< T >::rotate(), and theRotation.

Referenced by init().

                                             {
    return theRotation.rotate(p)/p.mag2();
  }

Point2D ThirdHitPredictionFromInvParabola::transformBack ( Point2D const &  p ) const

inlineprivate

Definition at line 77 of file ThirdHitPredictionFromInvParabola.h.

References Basic2DVector< T >::mag2(), TkRotation2D< T >::rotateBack(), and theRotation.

                                                 {
    return theRotation.rotateBack(p)/p.mag2();
  }

Member Data Documentation

Scalar ThirdHitPredictionFromInvParabola::dv

private

Definition at line 84 of file ThirdHitPredictionFromInvParabola.h.

Referenced by coeffB(), and init().

bool ThirdHitPredictionFromInvParabola::emptyM

private

Definition at line 91 of file ThirdHitPredictionFromInvParabola.h.

Referenced by init(), and rangeRPhi().

bool ThirdHitPredictionFromInvParabola::emptyP

private

Definition at line 91 of file ThirdHitPredictionFromInvParabola.h.

Referenced by init(), and rangeRPhi().

Scalar ThirdHitPredictionFromInvParabola::overDu

private

Definition at line 84 of file ThirdHitPredictionFromInvParabola.h.

Referenced by coeffA(), coeffB(), init(), and ipFromCurvature().

Scalar ThirdHitPredictionFromInvParabola::pv

private

Definition at line 84 of file ThirdHitPredictionFromInvParabola.h.

Referenced by coeffA(), init(), and ipFromCurvature().

Scalar ThirdHitPredictionFromInvParabola::su

private

Definition at line 84 of file ThirdHitPredictionFromInvParabola.h.

Referenced by coeffB(), and init().

RangeD ThirdHitPredictionFromInvParabola::theIpRangeMinus

private

Definition at line 89 of file ThirdHitPredictionFromInvParabola.h.

Referenced by init(), and rangeRPhi().

RangeD ThirdHitPredictionFromInvParabola::theIpRangePlus

private

Definition at line 89 of file ThirdHitPredictionFromInvParabola.h.

Referenced by init(), and rangeRPhi().

Rotation ThirdHitPredictionFromInvParabola::theRotation

private

Definition at line 83 of file ThirdHitPredictionFromInvParabola.h.

Referenced by init(), rangeRPhi(), transform(), and transformBack().

Scalar ThirdHitPredictionFromInvParabola::theTolerance

private

Definition at line 90 of file ThirdHitPredictionFromInvParabola.h.

Referenced by rangeRPhi().

Scalar ThirdHitPredictionFromInvParabola::u1u2

private

Definition at line 84 of file ThirdHitPredictionFromInvParabola.h.

Referenced by coeffA(), init(), and ipFromCurvature().