#include <TangentCircle.h>

Public Member Functions
int	charge (float magz)

double	curvatureError ()

GlobalVector	directionAtVertex ()
	Return the direction at the vertex. More...

GlobalPoint	innerPoint () const

GlobalPoint	outerPoint () const

double	rho () const

	TangentCircle ()

	TangentCircle (const GlobalVector &direction, const GlobalPoint &innerPoint, const GlobalPoint &outerPoint)

	TangentCircle (const GlobalPoint &outerPoint, const GlobalPoint &innerPoint, const GlobalPoint &vertexPoint)
	Copy of FastCircle. More...

	TangentCircle (const TangentCircle &primCircle, const GlobalPoint &outerPoint, const GlobalPoint &innerPoint)
	Calculate the parameters of a circle which pass by 2 points (innerPoint and outerPoint) and which is tangent to primCircle. More...

double	vertexError () const

GlobalPoint	vertexPoint () const

double	x0 () const

double	y0 () const

Private Member Functions
int	chargeLocally (float magz, GlobalVector v) const

GlobalVector	direction (const GlobalPoint &point) const

GlobalPoint	getPosition (const TangentCircle &circle, const GlobalPoint &initalPosition, double theta, int direction) const

double	isTangent (const TangentCircle &primCircle, const TangentCircle &secCircle) const

Private Attributes
int	theCharge

GlobalVector	theDirectionAtVertex

GlobalPoint	theInnerPoint

GlobalPoint	theOuterPoint

double	theRho

double	theVertexError

GlobalPoint	theVertexPoint

double	theX0

double	theY0

bool	valid

Detailed Description

Definition at line 7 of file TangentCircle.h.

Constructor & Destructor Documentation

◆ TangentCircle() [1/4]

TangentCircle::TangentCircle ( )

inline

Definition at line 13 of file TangentCircle.h.

       : theInnerPoint(),
         theOuterPoint(),
         theVertexPoint(),
         theDirectionAtVertex(),
         theX0(0),
         theY0(0),
         theRho(0),
         theVertexError(0),
         valid(false),
         theCharge(0) {}

◆ TangentCircle() [2/4]

TangentCircle::TangentCircle	(	const GlobalVector &	direction,
		const GlobalPoint &	innerPoint,
		const GlobalPoint &	outerPoint
	)

Calculate the circle from 2 points on the circle (the vertex=innerPoint and the outerPoint) and the tangent direction at the inner point

Definition at line 7 of file TangentCircle.cc.

References funct::cos(), bTagMiniDQMDeepCSV::denominator, direction(), SurfaceOrientation::inner, PV3DBase< T, PVType, FrameType >::mag(), mag(), SurfaceOrientation::outer, PV3DBase< T, PVType, FrameType >::perp2(), funct::sin(), theCharge, theDirectionAtVertex, theInnerPoint, theOuterPoint, theRho, theVertexError, theX0, theY0, valid, PV3DBase< T, PVType, FrameType >::x(), testProducerWithPsetDescEmpty_cfi::x1, testProducerWithPsetDescEmpty_cfi::x2, PV3DBase< T, PVType, FrameType >::y(), testProducerWithPsetDescEmpty_cfi::y1, and testProducerWithPsetDescEmpty_cfi::y2.

     : theInnerPoint(inner), theOuterPoint(outer), theVertexPoint(inner) {
   if (theInnerPoint.perp2() > theOuterPoint.perp2()) {
     valid = false;
   } else
     valid = true;
 
   double x1 = inner.x();
   double y1 = inner.y();
   double x2 = outer.x();
   double y2 = outer.y();
   double alpha1 = (direction.y() != 0) ? atan(-direction.x() / direction.y()) : PI / 2;
   double denominator = 2 * ((x1 - x2) * cos(alpha1) + (y1 - y2) * sin(alpha1));
   theRho = (denominator != 0) ? ((x1 - x2) * (x1 - x2) + (y1 - y2) * (y1 - y2)) / denominator : 1E12;
 
   // TODO : variable not yet calculated look in nucl.C
   theX0 = 1E10;
   theY0 = 1E10;
 
   theDirectionAtVertex = direction;
   theDirectionAtVertex /= theDirectionAtVertex.mag();
 
   //theCharge = (theRho>0) ? -1 : 1;
 
   theCharge = 0;
   theRho = fabs(theRho);
 
   theVertexError = (theInnerPoint - theOuterPoint).mag();
 }

◆ TangentCircle() [3/4]

TangentCircle::TangentCircle	(	const GlobalPoint &	outerPoint,
		const GlobalPoint &	innerPoint,
		const GlobalPoint &	vertexPoint
	)

Copy of FastCircle.

Definition at line 37 of file TangentCircle.cc.

References innerPoint(), FastCircle::isValid(), outerPoint(), PV3DBase< T, PVType, FrameType >::perp2(), FastCircle::rho(), theCharge, theDirectionAtVertex, theInnerPoint, theOuterPoint, theRho, theVertexError, theX0, theY0, valid, vertexPoint(), FastCircle::x0(), and FastCircle::y0().

     : theInnerPoint(innerPoint), theOuterPoint(outerPoint), theVertexPoint(vertexPoint) {
   FastCircle circle(outerPoint, innerPoint, vertexPoint);
   theX0 = circle.x0();
   theY0 = circle.y0();
   theRho = circle.rho();
   theVertexError = 0;
   theCharge = 0;
   theDirectionAtVertex = GlobalVector(1000, 1000, 1000);
   if (theInnerPoint.perp2() > theOuterPoint.perp2() || !circle.isValid()) {
     valid = false;
   } else
     valid = true;
 }

◆ TangentCircle() [4/4]

TangentCircle::TangentCircle	(	const TangentCircle &	primCircle,
		const GlobalPoint &	outerPoint,
		const GlobalPoint &	innerPoint
	)

Calculate the parameters of a circle which pass by 2 points (innerPoint and outerPoint) and which is tangent to primCircle.

Definition at line 54 of file TangentCircle.cc.

References DeadROC_duringRun::dir, getPosition(), mps_fire::i, innerPoint(), isTangent(), LogDebug, mag(), outerPoint(), PV3DBase< T, PVType, FrameType >::perp2(), rho(), alignCSCRings::s, theCharge, theDirectionAtVertex, theInnerPoint, theOuterPoint, theRho, tauSpinnerTable_cfi::theta, theVertexError, theVertexPoint, theX0, theY0, valid, bphysicsOniaDQM_cfi::vertex, vertexPoint(), PV3DBase< T, PVType, FrameType >::x(), x0(), PV3DBase< T, PVType, FrameType >::y(), and y0().

                                                             {
   if (theInnerPoint.perp2() > theOuterPoint.perp2()) {
     valid = false;
   } else
     valid = true;
 
   int NITER = 10;
 
   // Initial vertex used = outerPoint of the primary circle (should be the first estimation of the nuclear interaction position)
   GlobalPoint InitialVertex(primCircle.outerPoint().x(), primCircle.outerPoint().y(), 0);
   GlobalPoint SecInnerPoint(innerPoint.x(), innerPoint.y(), 0);
   GlobalPoint SecOuterPoint(outerPoint.x(), outerPoint.y(), 0);
 
   // distance between the initial vertex and the inner point of the secondary circle
   double s = (SecInnerPoint - InitialVertex).mag();
   double deltaTheta = s / primCircle.rho();
 
   double minTangentCondition = 1E12;
   TangentCircle theCorrectSecCircle;
   GlobalPoint vertex = InitialVertex;
   int dir = 1;
   double theta = deltaTheta / (NITER - 1);
 
   for (int i = 0; i < NITER; i++) {
     // get the circle which pass through outerPoint, innerPoint and the vertex
     TangentCircle secCircle(SecOuterPoint, SecInnerPoint, vertex);
 
     // get a value relative to the tangentness of the 2 circles
     double minCond = isTangent(primCircle, secCircle);
 
     // double dirDiff = (primCircle.direction(vertex) - secCircle.direction(vertex)).mag();
     // if( dirDiff > 1) dirDiff = 2-dirDiff;
 
     if (minCond < minTangentCondition) {
       minTangentCondition = minCond;
       theCorrectSecCircle = secCircle;
       vertex = getPosition(primCircle, secCircle.vertexPoint(), theta, dir);
       if (i == 0 && ((vertex - SecInnerPoint).mag() > (InitialVertex - SecInnerPoint).mag())) {
         dir = -1;
         vertex = getPosition(primCircle, InitialVertex, theta, dir);
         LogDebug("NuclearSeedGenerator") << "Change direction to look for vertex"
                                          << "\n";
       }
     } else
       break;
   }
   theInnerPoint = theCorrectSecCircle.innerPoint();
   theOuterPoint = theCorrectSecCircle.outerPoint();
   theVertexPoint = theCorrectSecCircle.vertexPoint();
   theX0 = theCorrectSecCircle.x0();
   theY0 = theCorrectSecCircle.y0();
   theRho = theCorrectSecCircle.rho();
   theCharge = 0;
   theDirectionAtVertex = GlobalVector(1000, 1000, 1000);
 
   theVertexError = s / NITER;
 }

Member Function Documentation

◆ charge()

int TangentCircle::charge ( float magz )

Definition at line 205 of file TangentCircle.cc.

References chargeLocally(), directionAtVertex(), F(), LogDebug, theCharge, theVertexPoint, theX0, theY0, findQualityFiles::v, PV3DBase< T, PVType, FrameType >::x(), and PV3DBase< T, PVType, FrameType >::y().

Referenced by TangentHelix::charge().

                                     {
   if (theCharge != 0)
     return theCharge;
 
   if (theX0 > 1E9 || theY0 > 1E9)
     theCharge = chargeLocally(magz, directionAtVertex());
   else {
     GlobalPoint center(theX0, theY0, 0);
     GlobalVector u = center - theVertexPoint;
     GlobalVector v = directionAtVertex();
 
     // F = force vector
     GlobalVector F(v.y() * magz, -v.x() * magz, 0);
     if (u.x() * F.x() + u.y() * F.y() > 0)
       theCharge = -1;
     else
       theCharge = 1;
 
     if (theCharge != chargeLocally(magz, v)) {
       LogDebug("NuclearSeedGenerator") << "Inconsistency in calculation of the charge"
                                        << "\n";
     }
   }
   return theCharge;
 }

◆ chargeLocally()

int TangentCircle::chargeLocally	(	float	magz,
		GlobalVector	v
	)		const

private

Definition at line 231 of file TangentCircle.cc.

References theOuterPoint, theVertexPoint, findQualityFiles::v, PV3DBase< T, PVType, FrameType >::x(), and PV3DBase< T, PVType, FrameType >::y().

Referenced by charge().

                                                                  {
   GlobalVector u = theOuterPoint - theVertexPoint;
   double tz = v.x() * u.y() - v.y() * u.x();
 
   if (tz * magz > 0)
     return 1;
   else
     return -1;
 }

◆ curvatureError()

double TangentCircle::curvatureError ( )

Definition at line 193 of file TangentCircle.cc.

References directionAtVertex(), mag(), rho(), theInnerPoint, theOuterPoint, theVertexError, and theVertexPoint.

Referenced by TangentHelix::curvatureError().

                                      {
   if ((theInnerPoint - theVertexPoint).mag() < theVertexError) {
     TangentCircle circle1(directionAtVertex(), theVertexPoint - theVertexError * directionAtVertex(), theOuterPoint);
     TangentCircle circle2(directionAtVertex(), theVertexPoint + theVertexError * directionAtVertex(), theOuterPoint);
     return fabs(1 / circle1.rho() - 1 / circle2.rho());
   } else {
     TangentCircle circle1(theOuterPoint, theInnerPoint, theVertexPoint - theVertexError * directionAtVertex());
     TangentCircle circle2(theOuterPoint, theInnerPoint, theVertexPoint + theVertexError * directionAtVertex());
     return fabs(1 / circle1.rho() - 1 / circle2.rho());
   }
 }

◆ direction()

GlobalVector TangentCircle::direction ( const GlobalPoint & point ) const

private

Definition at line 125 of file TangentCircle.cc.

References change_name::diff, DeadROC_duringRun::dir, LogDebug, mag(), point, theInnerPoint, theOuterPoint, theX0, and theY0.

Referenced by directionAtVertex(), and TangentCircle().

                                                                     {
   if (theY0 > 1E9 || theX0 > 1E9) {
     LogDebug("NuclearSeedGenerator") << "Center of TangentCircle not calculated but used !!!"
                                      << "\n";
   }
 
   // calculate the direction perpendicular to the vector v = point - center_of_circle
   GlobalVector dir(point.y() - theY0, theX0 - point.x(), 0);
 
   dir /= dir.mag();
 
   // Check the sign :
   GlobalVector fastDir = theOuterPoint - theInnerPoint;
   double diff = (dir - fastDir).mag();
   double sum = (dir + fastDir).mag();
 
   if (sum < diff)
     dir = (-1) * dir;
 
   return dir;
 }

◆ directionAtVertex()

GlobalVector TangentCircle::directionAtVertex ( )

Return the direction at the vertex.

Definition at line 147 of file TangentCircle.cc.

References direction(), theDirectionAtVertex, theVertexPoint, and PV3DBase< T, PVType, FrameType >::x().

Referenced by charge(), curvatureError(), and TangentHelix::directionAtVertex().

                                               {
   if (theDirectionAtVertex.x() > 999)
     theDirectionAtVertex = direction(theVertexPoint);
   return theDirectionAtVertex;
 }

◆ getPosition()

GlobalPoint TangentCircle::getPosition	(	const TangentCircle &	circle,
		const GlobalPoint &	initalPosition,
		double	theta,
		int	direction
	)		const

private

Definition at line 153 of file TangentCircle.cc.

References isotrackApplyRegressor::alpha, HLT_2024v14_cff::beta, funct::cos(), DeadROC_duringRun::dir, CustomPhysics_cfi::gamma, MainPageGenerator::l, rho(), Validation_hcalonly_cfi::sign, funct::sin(), tauSpinnerTable_cfi::theta, PV3DBase< T, PVType, FrameType >::x(), x0(), testProducerWithPsetDescEmpty_cfi::x2, PV3DBase< T, PVType, FrameType >::y(), y0(), and testProducerWithPsetDescEmpty_cfi::y2.

Referenced by TangentCircle().

                                                       {
   int sign[3];
   double x2 = initalPosition.x();
   double y2 = initalPosition.y();
 
   if ((x2 > circle.x0()) && dir > 0) {
     sign[0] = 1;
     sign[1] = -1;
     sign[2] = -1;
   }
   if ((x2 > circle.x0()) && dir < 0) {
     sign[0] = 1;
     sign[1] = 1;
     sign[2] = 1;
   }
   if ((x2 < circle.x0()) && dir > 0) {
     sign[0] = -1;
     sign[1] = 1;
     sign[2] = -1;
   }
   if ((x2 < circle.x0()) && dir < 0) {
     sign[0] = -1;
     sign[1] = -1;
     sign[2] = 1;
   }
 
   double l = 2 * circle.rho() * sin(theta / 2);
   double alpha = atan((y2 - circle.y0()) / (x2 - circle.x0()));
   double beta = PI / 2 - theta / 2;
   double gamma = PI + sign[2] * alpha - beta;
 
   double xnew = x2 + sign[0] * l * cos(gamma);
   double ynew = y2 + sign[1] * l * sin(gamma);
 
   return GlobalPoint(xnew, ynew, 0);
 }

◆ innerPoint()

GlobalPoint TangentCircle::innerPoint ( ) const

inline

Definition at line 46 of file TangentCircle.h.

References theInnerPoint.

Referenced by TangentCircle().

46 { return theInnerPoint; }

TangentCircle::theInnerPoint

GlobalPoint theInnerPoint

Definition: TangentCircle.h:57

◆ isTangent()

double TangentCircle::isTangent	(	const TangentCircle &	primCircle,
		const TangentCircle &	secCircle
	)		const

private

Definition at line 114 of file TangentCircle.cc.

References SiStripPI::min, rho(), x0(), and y0().

Referenced by TangentCircle().

                                                                                                      {
   // return a value that should be equal to 0 if primCircle and secCircle are tangent
 
   double distanceBetweenCircle = (primCircle.x0() - secCircle.x0()) * (primCircle.x0() - secCircle.x0()) +
                                  (primCircle.y0() - secCircle.y0()) * (primCircle.y0() - secCircle.y0());
   double RadiusSum = (primCircle.rho() + secCircle.rho()) * (primCircle.rho() + secCircle.rho());
   double RadiusDifference = (primCircle.rho() - secCircle.rho()) * (primCircle.rho() - secCircle.rho());
 
   return std::min(fabs(RadiusSum - distanceBetweenCircle), fabs(RadiusDifference - distanceBetweenCircle));
 }