BOSS 8.0.0
BESIII Offline Software System
Loading...
Searching...
No Matches
Lpar Class Reference

#include <Lpar.h>

Inheritance diagram for Lpar:

Public Member Functions

 Lpar ()
virtual ~Lpar ()
const Lpar & operator= (const Lpar &)
void neg ()
void circle (double x1, double y1, double x2, double y2, double x3, double y3)
double kappa () const
double radius () const
HepVector center () const
double s (double x, double y) const
double d (double x, double y) const
double dr (double x, double y) const
double s (double r, int dir=0) const
double phi (double r, int dir=0) const
int sd (double r, double x, double y, double limit, double &s, double &d) const
HepVector Hpar (const HepPoint3D &pivot) const
 Lpar ()
virtual ~Lpar ()
const Lpar & operator= (const Lpar &)
void neg ()
void circle (double x1, double y1, double x2, double y2, double x3, double y3)
double kappa () const
double radius () const
HepVector center () const
double s (double x, double y) const
double d (double x, double y) const
double dr (double x, double y) const
double s (double r, int dir=0) const
double phi (double r, int dir=0) const
int sd (double r, double x, double y, double limit, double &s, double &d) const
HepVector Hpar (const HepPoint3D &pivot) const
 Lpar ()
virtual ~Lpar ()
const Lpar & operator= (const Lpar &)
void neg ()
void circle (double x1, double y1, double x2, double y2, double x3, double y3)
double kappa () const
double radius () const
HepVector center () const
double s (double x, double y) const
double d (double x, double y) const
double dr (double x, double y) const
double s (double r, int dir=0) const
double phi (double r, int dir=0) const
int sd (double r, double x, double y, double limit, double &s, double &d) const
HepVector Hpar (const HepPoint3D &pivot) const

Friends

class Lpav
class Lpar::Cpar
std::ostream & operator<< (std::ostream &o, Lpar &s)
int intersect (const Lpar &lp1, const Lpar &lp2, HepVector &v1, HepVector &v2)
std::ostream & operator<< (std::ostream &o, Lpar &s)
int intersect (const Lpar &lp1, const Lpar &lp2, HepVector &v1, HepVector &v2)
std::ostream & operator<< (std::ostream &o, Lpar &s)
int intersect (const Lpar &lp1, const Lpar &lp2, HepVector &v1, HepVector &v2)

Detailed Description

Definition at line 47 of file InstallArea/x86_64-el9-gcc13-dbg/include/TrackUtil/Lpar.h.

Constructor & Destructor Documentation

◆ Lpar() [1/3]

Lpar::Lpar ( )
inline

Definition at line 162 of file InstallArea/x86_64-el9-gcc13-dbg/include/TrackUtil/Lpar.h.

162 {
163 m_alpha = 0;
164 m_beta = 1;
165 m_gamma = 0;
166 m_kappa = 0;
167 StatusCode scmgn = Gaudi::svcLocator()->service( "MagneticFieldSvc", m_pmgnIMF );
168 if ( scmgn != StatusCode::SUCCESS )
169 { std::cout << "Unable to open Magnetic field service" << std::endl; }
170}

Referenced by Hpar(), intersect, Lpar::Cpar, Lpav, operator<<, Lpav::operator<<, operator=(), and Lpav::operator=().

◆ ~Lpar() [1/3]

Lpar::~Lpar ( )
virtual

Definition at line 52 of file TrackUtil/src/Lpar.cxx.

52{}

◆ Lpar() [2/3]

Lpar::Lpar ( )

◆ ~Lpar() [2/3]

virtual Lpar::~Lpar ( )
virtual

◆ Lpar() [3/3]

Lpar::Lpar ( )

◆ ~Lpar() [3/3]

virtual Lpar::~Lpar ( )
virtual

Member Function Documentation

◆ center() [1/3]

HepVector Lpar::center ( ) const

Definition at line 241 of file TrackUtil/src/Lpar.cxx.

244{
245#else
246{
247 HepVector v( 3 );
248#endif
249 v( 1 ) = xc();
250 v( 2 ) = yc();
251 v( 3 ) = 0;
252 return ( v );
253}
v
**********Class see also m_nmax DOUBLE PRECISION m_amel DOUBLE PRECISION m_x2 DOUBLE PRECISION m_alfinv DOUBLE PRECISION m_Xenph INTEGER m_KeyWtm INTEGER m_idyfs DOUBLE PRECISION m_zini DOUBLE PRECISION m_q2 DOUBLE PRECISION m_Wt_KF DOUBLE PRECISION m_WtCut INTEGER m_KFfin *COMMON c_KarLud $ !Input CMS energy[GeV] $ !CMS energy after beam spread beam strahlung[GeV] $ !Beam energy spread[GeV] $ !z boost due to beam spread $ !electron beam mass *ff pair spectrum $ !minimum v
Definition KarLud.h:35

Referenced by TCircleFitter::fit(), intersect, and FTWire::z().

◆ center() [2/3]

HepVector Lpar::center ( ) const

◆ center() [3/3]

HepVector Lpar::center ( ) const

◆ circle() [1/3]

void Lpar::circle ( double x1,
double y1,
double x2,
double y2,
double x3,
double y3 )

Definition at line 75 of file TrackUtil/src/Lpar.cxx.

75 {
76 double a;
77 double b;
78 double c;
79 double delta = ( x1 - x2 ) * ( y1 - y3 ) - ( y1 - y2 ) * ( x1 - x3 );
80 if ( delta == 0 )
81 {
82 //
83 // three points are on a line.
84 //
85 m_kappa = 0;
86 double r12sq = ( x1 - x2 ) * ( x1 - x2 ) + ( y1 - y2 ) * ( y1 - y2 );
87 if ( r12sq > 0 )
88 {
89 double r12 = sqrt( r12sq );
90 m_beta = -( x1 - x2 ) / r12;
91 m_alpha = ( y1 - y2 ) / r12;
92 m_gamma = -( m_alpha * x1 + m_beta * y1 );
93 }
94 else
95 {
96 double r13sq = ( x1 - x3 ) * ( x1 - x3 ) + ( y1 - y3 ) * ( y1 - y3 );
97 if ( r13sq > 0 )
98 {
99 double r13 = sqrt( r13sq );
100 m_beta = -( x1 - x3 ) / r13;
101 m_alpha = ( y1 - y3 ) / r13;
102 m_gamma = -( m_alpha * x3 + m_beta * y3 );
103 }
104 else
105 {
106 double r23sq = ( x2 - x3 ) * ( x2 - x3 ) + ( y2 - y3 ) * ( y2 - y3 );
107 if ( r23sq > 0 )
108 {
109 double r23 = sqrt( r23sq );
110 m_beta = -( x2 - x3 ) / r23;
111 m_alpha = ( y2 - y3 ) / r23;
112 m_gamma = -( m_alpha * x3 + m_beta * y3 );
113 }
114 else
115 {
116 m_alpha = 1;
117 m_beta = 0;
118 m_gamma = 0;
119 }
120 }
121 }
122 }
123 else
124 {
125 double r1sq = x1 * x1 + y1 * y1;
126 double r2sq = x2 * x2 + y2 * y2;
127 double r3sq = x3 * x3 + y3 * y3;
128 a = 0.5 * ( ( y1 - y3 ) * ( r1sq - r2sq ) - ( y1 - y2 ) * ( r1sq - r3sq ) ) / delta;
129 b = 0.5 * ( -( x1 - x3 ) * ( r1sq - r2sq ) + ( x1 - x2 ) * ( r1sq - r3sq ) ) / delta;
130 double csq = ( x1 - a ) * ( x1 - a ) + ( y1 - b ) * ( y1 - b );
131 c = sqrt( csq );
132 double csq2 = ( x2 - a ) * ( x2 - a ) + ( y2 - b ) * ( y2 - b );
133 double csq3 = ( x3 - a ) * ( x3 - a ) + ( y3 - b ) * ( y3 - b );
134 m_kappa = 1 / ( 2 * c );
135 m_alpha = -2 * a * m_kappa;
136 m_beta = -2 * b * m_kappa;
137 m_gamma = ( a * a + b * b - c * c ) * m_kappa;
138 }
139}

◆ circle() [2/3]

void Lpar::circle ( double x1,
double y1,
double x2,
double y2,
double x3,
double y3 )

◆ circle() [3/3]

void Lpar::circle ( double x1,
double y1,
double x2,
double y2,
double x3,
double y3 )

◆ d() [1/3]

double Lpar::d ( double x,
double y ) const
inline

Definition at line 218 of file InstallArea/x86_64-el9-gcc13-dbg/include/TrackUtil/Lpar.h.

218 {
219 double dd = d0( x, y );
220 const double approx_limit = 0.2;
221 if ( std::fabs( m_kappa * dd ) > approx_limit ) return -1;
222 return dd * ( 1 - m_kappa * dd );
223}

Referenced by Lpav::delta_chisq(), intersect, and sd().

◆ d() [2/3]

double Lpar::d ( double x,
double y ) const
inline

◆ d() [3/3]

double Lpar::d ( double x,
double y ) const
inline

◆ dr() [1/3]

double Lpar::dr ( double x,
double y ) const
inline

Definition at line 225 of file InstallArea/x86_64-el9-gcc13-dbg/include/TrackUtil/Lpar.h.

225 {
226 double dx = xc() - x;
227 double dy = yc() - y;
228 double r = 0.5 / std::fabs( m_kappa );
229 return std::fabs( std::sqrt( dx * dx + dy * dy ) - r );
230}

◆ dr() [2/3]

double Lpar::dr ( double x,
double y ) const
inline

◆ dr() [3/3]

double Lpar::dr ( double x,
double y ) const
inline

◆ Hpar() [1/3]

HepVector Lpar::Hpar ( const HepPoint3D & pivot) const
inline

Definition at line 280 of file InstallArea/x86_64-el9-gcc13-dbg/include/TrackUtil/Lpar.h.

280 {
281 HepVector a( 5 );
282 double m_bField = -10000 * ( m_pmgnIMF->getReferField() );
283 double m_alpha = 10000. / 2.99792458 / m_bField;
284 double dd = d0( pivot.x(), pivot.y() );
285 a( 1 ) = dd * ( m_kappa * dd - 1 );
286 a( 2 ) = ( m_kappa > 0 )
287 ? std::atan2( yc() - pivot.y(), xc() - pivot.x() ) + M_PI
288 : std::atan2( pivot.y() - yc(), pivot.x() - xc() ) - M_PI + 2 * M_PI;
289 // a(3) = -2.0*BELLE_ALPHA*m_kappa;
290 a( 3 ) = -2.0 * m_alpha * m_kappa;
291 a( 4 ) = 0;
292 a( 5 ) = 0;
293 return a;
294}
M_PI
#define M_PI
Definition TConstant.h:4

◆ Hpar() [2/3]

HepVector Lpar::Hpar ( const HepPoint3D & pivot) const
inline

◆ Hpar() [3/3]

HepVector Lpar::Hpar ( const HepPoint3D & pivot) const
inline

◆ kappa() [1/3]

double Lpar::kappa ( ) const
inline

Definition at line 70 of file InstallArea/x86_64-el9-gcc13-dbg/include/TrackUtil/Lpar.h.

70{ return m_kappa; }

Referenced by TCircleFitter::fit(), intersect, s(), and s().

◆ kappa() [2/3]

double Lpar::kappa ( ) const
inline

Definition at line 70 of file InstallArea/x86_64-el9-gcc13-opt/include/TrackUtil/Lpar.h.

70{ return m_kappa; }

◆ kappa() [3/3]

double Lpar::kappa ( ) const
inline

Definition at line 70 of file Reconstruction/TrackUtil/include/TrackUtil/Lpar.h.

70{ return m_kappa; }

◆ neg() [1/3]

void Lpar::neg ( )
inline

Definition at line 207 of file InstallArea/x86_64-el9-gcc13-dbg/include/TrackUtil/Lpar.h.

207 {
208 m_alpha = -m_alpha;
209 m_beta = -m_beta;
210 m_gamma = -m_gamma;
211 m_kappa = -m_kappa;
212}

Referenced by Lpav::calculate_lpar(), and Lpav::calculate_lpar3().

◆ neg() [2/3]

void Lpar::neg ( )
inline

◆ neg() [3/3]

void Lpar::neg ( )
inline

◆ operator=() [1/3]

const Lpar & Lpar::operator= ( const Lpar & l)
inline

Definition at line 179 of file InstallArea/x86_64-el9-gcc13-dbg/include/TrackUtil/Lpar.h.

179 {
180 if ( this != &l )
181 {
182 m_alpha = l.m_alpha;
183 m_beta = l.m_beta;
184 m_gamma = l.m_gamma;
185 m_kappa = l.m_kappa;
186 }
187 return *this;
188}

Referenced by Lpav::operator=().

◆ operator=() [2/3]

const Lpar & Lpar::operator= ( const Lpar & )
inline

◆ operator=() [3/3]

const Lpar & Lpar::operator= ( const Lpar & )
inline

◆ phi() [1/3]

double Lpar::phi ( double r,
int dir = 0 ) const

Definition at line 199 of file TrackUtil/src/Lpar.cxx.

199 {
200 double x, y;
201 if ( !xy( r, x, y, dir ) ) return -1;
202 double p = atan2( y, x );
203 if ( p < 0 ) p += ( 2 * M_PI );
204 return p;
205}

Referenced by Lpav::extrapolate(), and sd().

◆ phi() [2/3]

double Lpar::phi ( double r,
int dir = 0 ) const

◆ phi() [3/3]

double Lpar::phi ( double r,
int dir = 0 ) const

◆ radius() [1/3]

double Lpar::radius ( ) const
inline

Definition at line 71 of file InstallArea/x86_64-el9-gcc13-dbg/include/TrackUtil/Lpar.h.

71{ return 0.5 / std::fabs( m_kappa ); }

Referenced by TCircleFitter::fit(), and FTWire::z().

◆ radius() [2/3]

double Lpar::radius ( ) const
inline

Definition at line 71 of file InstallArea/x86_64-el9-gcc13-opt/include/TrackUtil/Lpar.h.

71{ return 0.5 / std::fabs( m_kappa ); }

◆ radius() [3/3]

double Lpar::radius ( ) const
inline

Definition at line 71 of file Reconstruction/TrackUtil/include/TrackUtil/Lpar.h.

71{ return 0.5 / std::fabs( m_kappa ); }

◆ s() [1/6]

double Lpar::s ( double r,
int dir = 0 ) const

Definition at line 232 of file TrackUtil/src/Lpar.cxx.

232 {
233 double d0 = da();
234 if ( fabs( r ) < fabs( d0 ) ) return -1;
235 double b = fabs( kappa() ) * sqrt( ( r * r - d0 * d0 ) / ( 1 + 2 * kappa() * d0 ) );
236 if ( fabs( b ) > 1 ) return -1;
237 if ( dir == 0 ) return asin( b ) / fabs( kappa() );
238 return ( M_PI - asin( b ) ) / fabs( kappa() );
239}

◆ s() [2/6]

double Lpar::s ( double r,
int dir = 0 ) const

◆ s() [3/6]

double Lpar::s ( double r,
int dir = 0 ) const

◆ s() [4/6]

double Lpar::s ( double x,
double y ) const

Definition at line 220 of file TrackUtil/src/Lpar.cxx.

220 {
221 double xh, yh, xx, yy;
222 xhyh( x, y, xh, yh );
223 double fk = fabs( kappa() );
224 if ( fk == 0 ) return 0;
225 yy = 2 * fk * ( alpha() * yh - beta() * xh );
226 xx = 2 * kappa() * ( alpha() * xh + beta() * yh ) + xi2();
227 double sp = atan2( yy, xx );
228 if ( sp < 0 ) sp += ( 2 * M_PI );
229 return sp / 2 / fk;
230}
alpha
double alpha
Definition EvtPhokharaDef.hh:17

Referenced by intersect, Lpar::Cpar, Lpav::operator+=(), operator<<, Lpav::operator=(), and sd().

◆ s() [5/6]

double Lpar::s ( double x,
double y ) const

◆ s() [6/6]

double Lpar::s ( double x,
double y ) const

◆ sd() [1/3]

int Lpar::sd ( double r,
double x,
double y,
double limit,
double & s,
double & d ) const
inline

Definition at line 262 of file InstallArea/x86_64-el9-gcc13-dbg/include/TrackUtil/Lpar.h.

262 {
263 if ( ( x * yc() - y * xc() ) * m_kappa < 0 ) return 0;
264 double dd = d0( x, y );
265 d = dd * ( 1 - m_kappa * dd );
266 double d_cross_limit = d * limit;
267 if ( d_cross_limit < 0 || d_cross_limit > limit * limit ) return 0;
268 double rc = std::sqrt( m_alpha * m_alpha + m_beta * m_beta ) / ( 2 * m_kappa );
269 double rho = 1. / ( -2 * m_kappa );
270 double cosPhi = ( rc * rc + rho * rho - r * r ) / ( -2 * rc * rho );
271 double phi = std::acos( cosPhi );
272 s = std::fabs( rho ) * phi;
273 d *= r / ( std::fabs( rc ) * std::sin( phi ) );
274 if ( abs( d ) > abs( limit ) ) return 0;
275 d_cross_limit = d * limit;
276 if ( d_cross_limit > limit * limit ) return 0;
277 return 1;
278}
Lpar::s
double s(double x, double y) const
Definition TrackUtil/src/Lpar.cxx:220
Lpar::phi
double phi(double r, int dir=0) const
Definition TrackUtil/src/Lpar.cxx:199

Referenced by FTSegment::update3D().

◆ sd() [2/3]

int Lpar::sd ( double r,
double x,
double y,
double limit,
double & s,
double & d ) const
inline

◆ sd() [3/3]

int Lpar::sd ( double r,
double x,
double y,
double limit,
double & s,
double & d ) const
inline

◆ intersect [1/3]

int intersect ( const Lpar & lp1,
const Lpar & lp2,
HepVector & v1,
HepVector & v2 )
friend

Definition at line 255 of file TrackUtil/src/Lpar.cxx.

255 {
256 HepVector cen1( lp1.center() );
257 HepVector cen2( lp2.center() );
258 double dx = cen1( 1 ) - cen2( 1 );
259 double dy = cen1( 2 ) - cen2( 2 );
260 double dc = sqrt( dx * dx + dy * dy );
261 if ( dc < fabs( 0.5 / lp1.kappa() ) + fabs( 0.5 / lp2.kappa() ) )
262 {
263 double a1 = std::sqrt( lp1.alpha() ) + std::sqrt( lp1.beta() );
264 double a2 = std::sqrt( lp2.alpha() ) + std::sqrt( lp2.beta() );
265 double a3 = lp1.alpha() * lp2.alpha() + lp1.beta() * lp2.beta();
266 double det = lp1.alpha() * lp2.beta() - lp1.beta() * lp2.alpha();
267 if ( fabs( det ) > 1e-12 )
268 {
269 double c1 = a2 * std::sqrt( lp1.kappa() ) + a1 * std::sqrt( lp2.kappa() ) -
270 2.0 * a3 * lp1.kappa() * lp2.kappa();
271 if ( c1 != 0 )
272 {
273 double cinv = 1.0 / c1;
274 double c2 = std::sqrt( a3 ) - 0.5 * ( a1 + a2 ) -
275 2.0 * a3 * ( lp1.gamma() * lp2.kappa() + lp2.gamma() * lp1.kappa() );
276 double c3 = a2 * std::sqrt( lp1.gamma() ) + a1 * std::sqrt( lp2.gamma() ) -
277 2.0 * a3 * lp1.gamma() * lp2.gamma();
278 double root = std::sqrt( c2 ) - 4.0 * c1 * c3;
279 if ( root >= 0 )
280 {
281 root = sqrt( root );
282 double rad2[2];
283 rad2[0] = 0.5 * cinv * ( -c2 - root );
284 rad2[1] = 0.5 * cinv * ( -c2 + root );
285 double ab1 = -( lp2.beta() * lp1.gamma() - lp1.beta() * lp2.gamma() );
286 double ab2 = ( lp2.alpha() * lp1.gamma() - lp1.alpha() * lp2.gamma() );
287 double ac1 = -( lp2.beta() * lp1.kappa() - lp1.beta() * lp2.kappa() );
288 double ac2 = ( lp2.alpha() * lp1.kappa() - lp1.alpha() * lp2.kappa() );
289 double dinv = 1.0 / det;
290 v1( 1 ) = dinv * ( ab1 + ac1 * rad2[0] );
291 v1( 2 ) = dinv * ( ab2 + ac2 * rad2[0] );
292 v1( 3 ) = 0;
293 v2( 1 ) = dinv * ( ab1 + ac1 * rad2[1] );
294 v2( 2 ) = dinv * ( ab2 + ac2 * rad2[1] );
295 v2( 3 ) = 0;
296 double d1 = lp1.d( v1( 1 ), v1( 2 ) );
297 double d2 = lp2.d( v1( 1 ), v1( 2 ) );
298 double d3 = lp1.d( v2( 1 ), v2( 2 ) );
299 double d4 = lp2.d( v2( 1 ), v2( 2 ) );
300 double r = sqrt( rad2[0] );
301 Lpar::Cpar cp1( lp1 );
302 Lpar::Cpar cp2( lp2 );
303 for ( int j = 0; j < 2; j++ )
304 {
305 double s1, s2;
306 if ( j == 0 )
307 {
308 s1 = lp1.s( v1( 1 ), v1( 2 ) );
309 s2 = lp2.s( v1( 1 ), v1( 2 ) );
310 }
311 else
312 {
313 s1 = lp1.s( v2( 1 ), v2( 2 ) );
314 s2 = lp2.s( v2( 1 ), v2( 2 ) );
315 }
316 double phi1 = cp1.fi() + 2 * cp1.cu() * s1;
317 double phi2 = cp2.fi() + 2 * cp2.cu() * s2;
318 double f = ( 1 + 2 * cp1.cu() * cp1.da() ) * ( 1 + 2 * cp2.cu() * cp2.da() ) *
319 cos( cp1.fi() - cp2.fi() );
320 f -= 2 * ( lp1.gamma() * lp2.kappa() + lp2.gamma() * lp1.kappa() );
321 double cosphi12 = f;
322 }
323 return 2;
324 }
325 }
326 }
327 }
328 return 0;
329}
root
std::string root
Definition CalibModel.cxx:40
f
TFile f("ana_bhabha660a_dqa_mcPat_zy_old.root")
phi2
Double_t phi2
Definition DataBase/tau_mode.c:8
phi1
Double_t phi1
Definition DataBase/tau_mode.c:8
dc
int dc[18]
Definition EvtPycont.cc:63
cos
double cos(const BesAngle a)
Definition InstallArea/x86_64-el9-gcc13-dbg/include/MdcGeom/BesAngle.h:187
Lpar::center
HepVector center() const
Definition TrackUtil/src/Lpar.cxx:241

◆ intersect [2/3]

int intersect ( const Lpar & lp1,
const Lpar & lp2,
HepVector & v1,
HepVector & v2 )
friend

Definition at line 255 of file TrackUtil/src/Lpar.cxx.

255 {
256 HepVector cen1( lp1.center() );
257 HepVector cen2( lp2.center() );
258 double dx = cen1( 1 ) - cen2( 1 );
259 double dy = cen1( 2 ) - cen2( 2 );
260 double dc = sqrt( dx * dx + dy * dy );
261 if ( dc < fabs( 0.5 / lp1.kappa() ) + fabs( 0.5 / lp2.kappa() ) )
262 {
263 double a1 = std::sqrt( lp1.alpha() ) + std::sqrt( lp1.beta() );
264 double a2 = std::sqrt( lp2.alpha() ) + std::sqrt( lp2.beta() );
265 double a3 = lp1.alpha() * lp2.alpha() + lp1.beta() * lp2.beta();
266 double det = lp1.alpha() * lp2.beta() - lp1.beta() * lp2.alpha();
267 if ( fabs( det ) > 1e-12 )
268 {
269 double c1 = a2 * std::sqrt( lp1.kappa() ) + a1 * std::sqrt( lp2.kappa() ) -
270 2.0 * a3 * lp1.kappa() * lp2.kappa();
271 if ( c1 != 0 )
272 {
273 double cinv = 1.0 / c1;
274 double c2 = std::sqrt( a3 ) - 0.5 * ( a1 + a2 ) -
275 2.0 * a3 * ( lp1.gamma() * lp2.kappa() + lp2.gamma() * lp1.kappa() );
276 double c3 = a2 * std::sqrt( lp1.gamma() ) + a1 * std::sqrt( lp2.gamma() ) -
277 2.0 * a3 * lp1.gamma() * lp2.gamma();
278 double root = std::sqrt( c2 ) - 4.0 * c1 * c3;
279 if ( root >= 0 )
280 {
281 root = sqrt( root );
282 double rad2[2];
283 rad2[0] = 0.5 * cinv * ( -c2 - root );
284 rad2[1] = 0.5 * cinv * ( -c2 + root );
285 double ab1 = -( lp2.beta() * lp1.gamma() - lp1.beta() * lp2.gamma() );
286 double ab2 = ( lp2.alpha() * lp1.gamma() - lp1.alpha() * lp2.gamma() );
287 double ac1 = -( lp2.beta() * lp1.kappa() - lp1.beta() * lp2.kappa() );
288 double ac2 = ( lp2.alpha() * lp1.kappa() - lp1.alpha() * lp2.kappa() );
289 double dinv = 1.0 / det;
290 v1( 1 ) = dinv * ( ab1 + ac1 * rad2[0] );
291 v1( 2 ) = dinv * ( ab2 + ac2 * rad2[0] );
292 v1( 3 ) = 0;
293 v2( 1 ) = dinv * ( ab1 + ac1 * rad2[1] );
294 v2( 2 ) = dinv * ( ab2 + ac2 * rad2[1] );
295 v2( 3 ) = 0;
296 double d1 = lp1.d( v1( 1 ), v1( 2 ) );
297 double d2 = lp2.d( v1( 1 ), v1( 2 ) );
298 double d3 = lp1.d( v2( 1 ), v2( 2 ) );
299 double d4 = lp2.d( v2( 1 ), v2( 2 ) );
300 double r = sqrt( rad2[0] );
301 Lpar::Cpar cp1( lp1 );
302 Lpar::Cpar cp2( lp2 );
303 for ( int j = 0; j < 2; j++ )
304 {
305 double s1, s2;
306 if ( j == 0 )
307 {
308 s1 = lp1.s( v1( 1 ), v1( 2 ) );
309 s2 = lp2.s( v1( 1 ), v1( 2 ) );
310 }
311 else
312 {
313 s1 = lp1.s( v2( 1 ), v2( 2 ) );
314 s2 = lp2.s( v2( 1 ), v2( 2 ) );
315 }
316 double phi1 = cp1.fi() + 2 * cp1.cu() * s1;
317 double phi2 = cp2.fi() + 2 * cp2.cu() * s2;
318 double f = ( 1 + 2 * cp1.cu() * cp1.da() ) * ( 1 + 2 * cp2.cu() * cp2.da() ) *
319 cos( cp1.fi() - cp2.fi() );
320 f -= 2 * ( lp1.gamma() * lp2.kappa() + lp2.gamma() * lp1.kappa() );
321 double cosphi12 = f;
322 }
323 return 2;
324 }
325 }
326 }
327 }
328 return 0;
329}

◆ intersect [3/3]

int intersect ( const Lpar & lp1,
const Lpar & lp2,
HepVector & v1,
HepVector & v2 )
friend

Definition at line 255 of file TrackUtil/src/Lpar.cxx.

255 {
256 HepVector cen1( lp1.center() );
257 HepVector cen2( lp2.center() );
258 double dx = cen1( 1 ) - cen2( 1 );
259 double dy = cen1( 2 ) - cen2( 2 );
260 double dc = sqrt( dx * dx + dy * dy );
261 if ( dc < fabs( 0.5 / lp1.kappa() ) + fabs( 0.5 / lp2.kappa() ) )
262 {
263 double a1 = std::sqrt( lp1.alpha() ) + std::sqrt( lp1.beta() );
264 double a2 = std::sqrt( lp2.alpha() ) + std::sqrt( lp2.beta() );
265 double a3 = lp1.alpha() * lp2.alpha() + lp1.beta() * lp2.beta();
266 double det = lp1.alpha() * lp2.beta() - lp1.beta() * lp2.alpha();
267 if ( fabs( det ) > 1e-12 )
268 {
269 double c1 = a2 * std::sqrt( lp1.kappa() ) + a1 * std::sqrt( lp2.kappa() ) -
270 2.0 * a3 * lp1.kappa() * lp2.kappa();
271 if ( c1 != 0 )
272 {
273 double cinv = 1.0 / c1;
274 double c2 = std::sqrt( a3 ) - 0.5 * ( a1 + a2 ) -
275 2.0 * a3 * ( lp1.gamma() * lp2.kappa() + lp2.gamma() * lp1.kappa() );
276 double c3 = a2 * std::sqrt( lp1.gamma() ) + a1 * std::sqrt( lp2.gamma() ) -
277 2.0 * a3 * lp1.gamma() * lp2.gamma();
278 double root = std::sqrt( c2 ) - 4.0 * c1 * c3;
279 if ( root >= 0 )
280 {
281 root = sqrt( root );
282 double rad2[2];
283 rad2[0] = 0.5 * cinv * ( -c2 - root );
284 rad2[1] = 0.5 * cinv * ( -c2 + root );
285 double ab1 = -( lp2.beta() * lp1.gamma() - lp1.beta() * lp2.gamma() );
286 double ab2 = ( lp2.alpha() * lp1.gamma() - lp1.alpha() * lp2.gamma() );
287 double ac1 = -( lp2.beta() * lp1.kappa() - lp1.beta() * lp2.kappa() );
288 double ac2 = ( lp2.alpha() * lp1.kappa() - lp1.alpha() * lp2.kappa() );
289 double dinv = 1.0 / det;
290 v1( 1 ) = dinv * ( ab1 + ac1 * rad2[0] );
291 v1( 2 ) = dinv * ( ab2 + ac2 * rad2[0] );
292 v1( 3 ) = 0;
293 v2( 1 ) = dinv * ( ab1 + ac1 * rad2[1] );
294 v2( 2 ) = dinv * ( ab2 + ac2 * rad2[1] );
295 v2( 3 ) = 0;
296 double d1 = lp1.d( v1( 1 ), v1( 2 ) );
297 double d2 = lp2.d( v1( 1 ), v1( 2 ) );
298 double d3 = lp1.d( v2( 1 ), v2( 2 ) );
299 double d4 = lp2.d( v2( 1 ), v2( 2 ) );
300 double r = sqrt( rad2[0] );
301 Lpar::Cpar cp1( lp1 );
302 Lpar::Cpar cp2( lp2 );
303 for ( int j = 0; j < 2; j++ )
304 {
305 double s1, s2;
306 if ( j == 0 )
307 {
308 s1 = lp1.s( v1( 1 ), v1( 2 ) );
309 s2 = lp2.s( v1( 1 ), v1( 2 ) );
310 }
311 else
312 {
313 s1 = lp1.s( v2( 1 ), v2( 2 ) );
314 s2 = lp2.s( v2( 1 ), v2( 2 ) );
315 }
316 double phi1 = cp1.fi() + 2 * cp1.cu() * s1;
317 double phi2 = cp2.fi() + 2 * cp2.cu() * s2;
318 double f = ( 1 + 2 * cp1.cu() * cp1.da() ) * ( 1 + 2 * cp2.cu() * cp2.da() ) *
319 cos( cp1.fi() - cp2.fi() );
320 f -= 2 * ( lp1.gamma() * lp2.kappa() + lp2.gamma() * lp1.kappa() );
321 double cosphi12 = f;
322 }
323 return 2;
324 }
325 }
326 }
327 }
328 return 0;
329}

◆ Lpar::Cpar

Lpar::Cpar
friend

Definition at line 110 of file InstallArea/x86_64-el9-gcc13-dbg/include/TrackUtil/Lpar.h.

◆ Lpav

Lpav
friend

Definition at line 83 of file InstallArea/x86_64-el9-gcc13-dbg/include/TrackUtil/Lpar.h.

Referenced by Hpar(), and Lpav.

◆ operator<< [1/3]

std::ostream & operator<< ( std::ostream & o,
Lpar & s )
friend

Definition at line 339 of file TrackUtil/src/Lpar.cxx.

339 {
340 return o << " al=" << s.m_alpha << " be=" << s.m_beta << " ka=" << s.m_kappa
341 << " ga=" << s.m_gamma;
342}

◆ operator<< [2/3]

std::ostream & operator<< ( std::ostream & o,
Lpar & s )
friend

Definition at line 339 of file TrackUtil/src/Lpar.cxx.

339 {
340 return o << " al=" << s.m_alpha << " be=" << s.m_beta << " ka=" << s.m_kappa
341 << " ga=" << s.m_gamma;
342}

◆ operator<< [3/3]

std::ostream & operator<< ( std::ostream & o,
Lpar & s )
friend

Definition at line 339 of file TrackUtil/src/Lpar.cxx.

339 {
340 return o << " al=" << s.m_alpha << " be=" << s.m_beta << " ka=" << s.m_kappa
341 << " ga=" << s.m_gamma;
342}
The documentation for this class was generated from the following files: