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

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.

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

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

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.

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

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

Member Function Documentation

coeffA()

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

inlineprivate

Definition at line 92 of file ThirdHitPredictionFromInvParabola.h.

References HltBtagPostValidation_cff::c, overDu, and u1u2.

Referenced by findPointAtCurve(), and predV().

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

coeffB()

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

inlineprivate

Definition at line 97 of file ThirdHitPredictionFromInvParabola.h.

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

Referenced by findPointAtCurve(), and predV().

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

findPointAtCurve()

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

inlineprivate

Definition at line 114 of file ThirdHitPredictionFromInvParabola.h.

References A, a, simBeamSpotPI::alpha, B, HltBtagPostValidation_cff::c, coeffA(), coeffB(), dumpMFGeometry_cfg::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
}

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.

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);
  }

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.

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

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

ipFromCurvature()

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

inlineprivate

Definition at line 102 of file ThirdHitPredictionFromInvParabola.h.

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

Referenced by init().

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

operator()() [1/2]

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

inline

Definition at line 55 of file ThirdHitPredictionFromInvParabola.h.

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

{ return rangeRPhi(radius, charge); }

operator()() [2/2]

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

inline

Definition at line 57 of file ThirdHitPredictionFromInvParabola.h.

References CosmicsPD_Skims::radius, and rangeRPhi().

{ return rangeRPhi(radius); }

predV()

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

inlineprivate

Definition at line 109 of file ThirdHitPredictionFromInvParabola.h.

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

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

rangeRPhi() [1/2]

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

Definition at line 70 of file ThirdHitPredictionFromInvParabola.cc.

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()().

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

rangeRPhi() [2/2]

ThirdHitPredictionFromInvParabola::Range ThirdHitPredictionFromInvParabola::rangeRPhi

(

Scalar

radius

)

const

Definition at line 99 of file ThirdHitPredictionFromInvParabola.cc.

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(), mkfit::Config::xr, y, and TkRotation2D< T >::y().

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

transform()

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

inlineprivate

Definition at line 76 of file ThirdHitPredictionFromInvParabola.h.

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

Referenced by init().

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

transformBack()

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

inlineprivate

Definition at line 78 of file ThirdHitPredictionFromInvParabola.h.

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

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

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.

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().