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) const
Range	operator() (Scalar radius, int charge) const
Range	rangeRPhi (Scalar radius) const
Range	rangeRPhi (Scalar radius, int charge) const
	ThirdHitPredictionFromInvParabola ()
	ThirdHitPredictionFromInvParabola (const GlobalPoint &P1, const GlobalPoint &P2, Scalar ip, Scalar curv, Scalar tolerance)
	ThirdHitPredictionFromInvParabola (Scalar x1, Scalar y1, Scalar x2, Scalar y2, 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

Friends
int	test::PixelTriplets_InvPrbl_prec::test ()
int	test::PixelTriplets_InvPrbl_t::test ()

Detailed Description

Definition at line 33 of file ThirdHitPredictionFromInvParabola.h.

Member Typedef Documentation

Point2D

typedef Basic2DVector<Scalar> ThirdHitPredictionFromInvParabola::Point2D

Definition at line 43 of file ThirdHitPredictionFromInvParabola.h.

Range

typedef PixelRecoRange<float> ThirdHitPredictionFromInvParabola::Range

Definition at line 41 of file ThirdHitPredictionFromInvParabola.h.

RangeD

typedef PixelRecoRange<Scalar> ThirdHitPredictionFromInvParabola::RangeD

Definition at line 42 of file ThirdHitPredictionFromInvParabola.h.

Rotation

typedef TkRotation2D<Scalar> ThirdHitPredictionFromInvParabola::Rotation

Definition at line 40 of file ThirdHitPredictionFromInvParabola.h.

Scalar

using ThirdHitPredictionFromInvParabola::Scalar = double

Definition at line 39 of file ThirdHitPredictionFromInvParabola.h.

Constructor & Destructor Documentation

ThirdHitPredictionFromInvParabola() [1/3]

ThirdHitPredictionFromInvParabola::ThirdHitPredictionFromInvParabola ( )

inline

Definition at line 45 of file ThirdHitPredictionFromInvParabola.h.

{}

ThirdHitPredictionFromInvParabola() [2/3]

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

inline

Definition at line 46 of file ThirdHitPredictionFromInvParabola.h.

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

References funct::abs(), init(), testProducerWithPsetDescEmpty_cfi::x1, testProducerWithPsetDescEmpty_cfi::x2, testProducerWithPsetDescEmpty_cfi::y1, and testProducerWithPsetDescEmpty_cfi::y2.

ThirdHitPredictionFromInvParabola() [3/3]

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

Definition at line 32 of file ThirdHitPredictionFromInvParabola.cc.

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

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

Member Function Documentation

coeffA()

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

inlineprivate

Definition at line 92 of file ThirdHitPredictionFromInvParabola.h.

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

References c, overDu, pv, and u1u2.

Referenced by findPointAtCurve(), and predV().

coeffB()

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

inlineprivate

Definition at line 97 of file ThirdHitPredictionFromInvParabola.h.

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

References c, dv, overDu, and su.

Referenced by findPointAtCurve(), and predV().

findPointAtCurve()

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

inlineprivate

Definition at line 114 of file ThirdHitPredictionFromInvParabola.h.

                                                                                                        {
  //
  // 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
}

References a, alpha, TtFullHadDaughter::B, c, coeffA(), coeffB(), dumpMFGeometry_cfg::delta, alignCSCRings::r, mathSSE::sqrt(), and findQualityFiles::v.

Referenced by rangeRPhi().

init() [1/2]

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

inline

Definition at line 65 of file ThirdHitPredictionFromInvParabola.h.

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

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

Referenced by ThirdHitPredictionFromInvParabola().

init() [2/2]

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

Definition at line 38 of file ThirdHitPredictionFromInvParabola.cc.

                                                                                                               {
  //  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();
}

References funct::abs(), dv, PixelRecoRange< T >::empty(), emptyM, emptyP, PixelRecoRange< T >::intersection(), ipFromCurvature(), SiStripPI::max, min(), overDu, p1, p2, pv, su, theIpRangeMinus, theIpRangePlus, theRotation, transform(), u1u2, testProducerWithPsetDescEmpty_cfi::x1, testProducerWithPsetDescEmpty_cfi::x2, testProducerWithPsetDescEmpty_cfi::y1, and testProducerWithPsetDescEmpty_cfi::y2.

ipFromCurvature()

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

inlineprivate

Definition at line 102 of file ThirdHitPredictionFromInvParabola.h.

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

References PixelRecoUtilities::curvature(), overDu, pv, and u1u2.

Referenced by init().

operator()() [1/2]

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

inline

Definition at line 57 of file ThirdHitPredictionFromInvParabola.h.

{ return rangeRPhi(radius); }

References CosmicsPD_Skims::radius, and rangeRPhi().

operator()() [2/2]

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

inline

Definition at line 55 of file ThirdHitPredictionFromInvParabola.h.

{ return rangeRPhi(radius, charge); }

References ALCARECOTkAlJpsiMuMu_cff::charge, CosmicsPD_Skims::radius, and rangeRPhi().

predV()

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

inlineprivate

Definition at line 109 of file ThirdHitPredictionFromInvParabola.h.

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

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

rangeRPhi() [1/2]

ThirdHitPredictionFromInvParabola::Range ThirdHitPredictionFromInvParabola::rangeRPhi

(

Scalar

radius

)

const

Definition at line 99 of file ThirdHitPredictionFromInvParabola.cc.

                                                                                                       {
  auto getRange = [&](Scalar phi1, Scalar phi2, bool empty) -> RangeD {
    if (phi2 < phi1)
      std::swap(phi1, phi2);
    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 * phi1 - theTolerance, radius * 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] = f_atan2f(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));
}

References relativeConstraints::empty, emptyM, emptyP, findPointAtCurve(), mps_fire::i, PixelRecoRange< T >::max(), PixelRecoRange< T >::min(), diffTwoXMLs::r1, diffTwoXMLs::r2, CosmicsPD_Skims::radius, std::swap(), theIpRangeMinus, theIpRangePlus, theRotation, theTolerance, findQualityFiles::v, x, TkRotation2D< T >::x(), y, and TkRotation2D< T >::y().

rangeRPhi() [2/2]

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

Definition at line 70 of file ThirdHitPredictionFromInvParabola.cc.

                                                                                                         {
  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 = f_phi(theRotation.rotateBack(Point2D(u[0], v[0])));
  Scalar phi2 = phi1 + (v[1] - v[0]);
 
  if (phi2 < phi1)
    std::swap(phi1, phi2);
 
  if (ip.empty()) {
    Range r1(phi1 * radius - theTolerance, phi1 * radius + theTolerance);
    Range r2(phi2 * radius - theTolerance, phi2 * radius + theTolerance);
    return r1.intersection(r2);  // this range can be empty
  }
 
  return Range(radius * phi1 - theTolerance, radius * phi2 + theTolerance);
}

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

Referenced by operator()().

transform()

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

inlineprivate

Definition at line 76 of file ThirdHitPredictionFromInvParabola.h.

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

References AlCaHLTBitMon_ParallelJobs::p, TkRotation2D< T >::rotate(), and theRotation.

Referenced by init().

transformBack()

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

inlineprivate

Definition at line 78 of file ThirdHitPredictionFromInvParabola.h.

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

References AlCaHLTBitMon_ParallelJobs::p, TkRotation2D< T >::rotateBack(), and theRotation.

Friends And Related Function Documentation

test::PixelTriplets_InvPrbl_prec::test

int test::PixelTriplets_InvPrbl_prec::test ( )

friend

test::PixelTriplets_InvPrbl_t::test

int test::PixelTriplets_InvPrbl_t::test ( )

friend

Member Data Documentation

dv

Scalar ThirdHitPredictionFromInvParabola::dv

private

Definition at line 82 of file ThirdHitPredictionFromInvParabola.h.

Referenced by coeffB(), and init().

emptyM

bool ThirdHitPredictionFromInvParabola::emptyM

private

Definition at line 89 of file ThirdHitPredictionFromInvParabola.h.

Referenced by init(), and rangeRPhi().

emptyP

bool ThirdHitPredictionFromInvParabola::emptyP

private

Definition at line 89 of file ThirdHitPredictionFromInvParabola.h.

Referenced by init(), and rangeRPhi().

overDu

Scalar ThirdHitPredictionFromInvParabola::overDu

private

Definition at line 82 of file ThirdHitPredictionFromInvParabola.h.

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

pv

Scalar ThirdHitPredictionFromInvParabola::pv

private

Definition at line 82 of file ThirdHitPredictionFromInvParabola.h.

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

su

Scalar ThirdHitPredictionFromInvParabola::su

private

Definition at line 82 of file ThirdHitPredictionFromInvParabola.h.

Referenced by coeffB(), and init().

theIpRangeMinus

RangeD ThirdHitPredictionFromInvParabola::theIpRangeMinus

private

Definition at line 87 of file ThirdHitPredictionFromInvParabola.h.

Referenced by init(), and rangeRPhi().

theIpRangePlus

RangeD ThirdHitPredictionFromInvParabola::theIpRangePlus

private

Definition at line 87 of file ThirdHitPredictionFromInvParabola.h.

Referenced by init(), and rangeRPhi().

theRotation

Rotation ThirdHitPredictionFromInvParabola::theRotation

private

Definition at line 81 of file ThirdHitPredictionFromInvParabola.h.

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

theTolerance

Scalar ThirdHitPredictionFromInvParabola::theTolerance

private

Definition at line 88 of file ThirdHitPredictionFromInvParabola.h.

Referenced by rangeRPhi().

u1u2

Scalar ThirdHitPredictionFromInvParabola::u1u2

private

Definition at line 82 of file ThirdHitPredictionFromInvParabola.h.

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