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

#include <Lpav.h>

Inheritance diagram for Lpav:

Classes

class  Singular
class  Singular_c

Public Member Functions

 Lpav ()
virtual ~Lpav ()
void calculate_average (void)
void calculate_average3 (void)
void calculate_average (double x, double y, double w=1)
void calculate_average3 (double x, double y, double w=1)
double calculate_lpar (void)
double calculate_lpar3 (void)
double fit ()
double fit (double x, double y, double w=1)
void clear ()
void add_point (double x, double y, double w=1)
void add_point_frac (double x, double y, double w, double f)
double nc () const
HepSymMatrix cov (int=0) const
HepSymMatrix cov_c (int=0) const
int extrapolate (double, double &, double &) const
double similarity (double, double) const
double delta_chisq (double x, double y, double w=1) const
double chisq () const
double prob () const
double chi_deg () const
const Lpav & operator= (const Lpav &)
const Lpav & operator= (const Lpar &)
const Lpav & operator+= (const Lpav &)
 Lpav ()
virtual ~Lpav ()
void calculate_average (void)
void calculate_average3 (void)
void calculate_average (double x, double y, double w=1)
void calculate_average3 (double x, double y, double w=1)
double calculate_lpar (void)
double calculate_lpar3 (void)
double fit ()
double fit (double x, double y, double w=1)
void clear ()
void add_point (double x, double y, double w=1)
void add_point_frac (double x, double y, double w, double f)
double nc () const
HepSymMatrix cov (int=0) const
HepSymMatrix cov_c (int=0) const
int extrapolate (double, double &, double &) const
double similarity (double, double) const
double delta_chisq (double x, double y, double w=1) const
double chisq () const
double prob () const
double chi_deg () const
const Lpav & operator= (const Lpav &)
const Lpav & operator= (const Lpar &)
const Lpav & operator+= (const Lpav &)
 Lpav ()
virtual ~Lpav ()
void calculate_average (void)
void calculate_average3 (void)
void calculate_average (double x, double y, double w=1)
void calculate_average3 (double x, double y, double w=1)
double calculate_lpar (void)
double calculate_lpar3 (void)
double fit ()
double fit (double x, double y, double w=1)
void clear ()
void add_point (double x, double y, double w=1)
void add_point_frac (double x, double y, double w, double f)
double nc () const
HepSymMatrix cov (int=0) const
HepSymMatrix cov_c (int=0) const
int extrapolate (double, double &, double &) const
double similarity (double, double) const
double delta_chisq (double x, double y, double w=1) const
double chisq () const
double prob () const
double chi_deg () const
const Lpav & operator= (const Lpav &)
const Lpav & operator= (const Lpar &)
const Lpav & operator+= (const Lpav &)
Public Member Functions inherited from Lpar
 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

std::ostream & operator<< (std::ostream &o, const Lpav &s)
Lpav operator+ (const Lpav &la1, const Lpav &la2)
std::ostream & operator<< (std::ostream &o, const Lpav &s)
Lpav operator+ (const Lpav &la1, const Lpav &la2)
std::ostream & operator<< (std::ostream &o, const Lpav &s)
Lpav operator+ (const Lpav &la1, const Lpav &la2)

Detailed Description

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

Constructor & Destructor Documentation

◆ Lpav() [1/3]

Lpav::Lpav ( )

Definition at line 49 of file TrackUtil/src/Lpav.cxx.

49{ clear(); }

Referenced by operator+, operator+=(), operator<<, operator=(), and operator=().

◆ ~Lpav() [1/3]

Lpav::~Lpav ( )
virtual

Definition at line 55 of file TrackUtil/src/Lpav.cxx.

55{}

◆ Lpav() [2/3]

Lpav::Lpav ( )

◆ ~Lpav() [2/3]

virtual Lpav::~Lpav ( )
virtual

◆ Lpav() [3/3]

Lpav::Lpav ( )

◆ ~Lpav() [3/3]

virtual Lpav::~Lpav ( )
virtual

Member Function Documentation

◆ add_point() [1/3]

void Lpav::add_point ( double x,
double y,
double w = 1 )

Referenced by TCircleFitter::fit().

◆ add_point() [2/3]

void Lpav::add_point ( double x,
double y,
double w = 1 )

◆ add_point() [3/3]

void Lpav::add_point ( double x,
double y,
double w = 1 )

◆ add_point_frac() [1/3]

void Lpav::add_point_frac ( double x,
double y,
double w,
double f )

Definition at line 590 of file TrackUtil/src/Lpav.cxx.

590 {
591 register double wi = w * a;
592 m_wsum += wi;
593 m_xsum += wi * xi;
594 m_ysum += wi * yi;
595 m_xxsum += wi * xi * xi;
596 m_yysum += wi * yi * yi;
597 m_xysum += wi * xi * yi;
598 register double rri = ( xi * xi + yi * yi );
599 register double wrri = wi * rri;
600 m_xrrsum += wrri * xi;
601 m_yrrsum += wrri * yi;
602 m_rrrrsum += wrri * rri;
603 m_nc += a;
604}
w
double w
Definition EvtPhokharaDef.hh:28

◆ add_point_frac() [2/3]

void Lpav::add_point_frac ( double x,
double y,
double w,
double f )

◆ add_point_frac() [3/3]

void Lpav::add_point_frac ( double x,
double y,
double w,
double f )

◆ calculate_average() [1/6]

void Lpav::calculate_average ( double x,
double y,
double w = 1 )

Definition at line 78 of file TrackUtil/src/Lpav.cxx.

78 {
79 if ( m_wsum <= 0 ) return;
80 m_wsum_temp = m_wsum + wi;
81 double rri( xi * xi + yi * yi );
82 double wrri( wi * rri );
83 double wsum_inv( 1 / m_wsum_temp );
84 m_xav = ( m_xsum + wi * xi ) * wsum_inv;
85 m_yav = ( m_ysum + wi * yi ) * wsum_inv;
86
87 double xxav( ( m_xxsum + wi * xi * xi ) * wsum_inv );
88 double yyav( ( m_yysum + wi * yi * yi ) * wsum_inv );
89 double xyav( ( m_xysum + wi * xi * yi ) * wsum_inv );
90 double xrrav( ( m_xrrsum + xi * wrri ) * wsum_inv );
91 double yrrav( ( m_yrrsum + yi * wrri ) * wsum_inv );
92 double rrrrav( ( m_rrrrsum + wrri * rri ) * wsum_inv );
93
94 calculate_average_n( xxav, yyav, xyav, xrrav, yrrav, rrrrav );
95}

◆ calculate_average() [2/6]

void Lpav::calculate_average ( double x,
double y,
double w = 1 )

◆ calculate_average() [3/6]

void Lpav::calculate_average ( double x,
double y,
double w = 1 )

◆ calculate_average() [4/6]

void Lpav::calculate_average ( void )

Definition at line 97 of file TrackUtil/src/Lpav.cxx.

97 {
98 if ( m_wsum <= 0 ) return;
99 m_wsum_temp = m_wsum;
100 double wsum_inv( 1 / m_wsum_temp );
101 m_xav = m_xsum * wsum_inv;
102 m_yav = m_ysum * wsum_inv;
103
104 double xxav( m_xxsum * wsum_inv );
105 double yyav( m_yysum * wsum_inv );
106 double xyav( m_xysum * wsum_inv );
107 double xrrav( m_xrrsum * wsum_inv );
108 double yrrav( m_yrrsum * wsum_inv );
109 double rrrrav( m_rrrrsum * wsum_inv );
110
111 calculate_average_n( xxav, yyav, xyav, xrrav, yrrav, rrrrav );
112}

Referenced by fit(), and fit().

◆ calculate_average() [5/6]

void Lpav::calculate_average ( void )

◆ calculate_average() [6/6]

void Lpav::calculate_average ( void )

◆ calculate_average3() [1/6]

void Lpav::calculate_average3 ( double x,
double y,
double w = 1 )

Definition at line 186 of file TrackUtil/src/Lpav.cxx.

186 {
187 if ( m_wsum <= 0 ) return;
188 m_wsum_temp = m_wsum + wi;
189 double wsum_inv( 1 / m_wsum_temp );
190 double rri( xi * xi + yi * yi );
191 m_xav = ( m_xsum + wi * xi ) * wsum_inv;
192 m_yav = ( m_ysum + wi * yi ) * wsum_inv;
193
194 m_rscale = 1;
195 m_cosrot = 1;
196 m_sinrot = 0;
197 m_xxavp = ( m_xxsum + wi * xi * xi ) * wsum_inv;
198 m_xyavp = ( m_xysum + wi * xi * yi ) * wsum_inv;
199 m_yyavp = ( m_yysum + wi * yi * yi ) * wsum_inv;
200 double wrri( wi * rri );
201 m_xrravp = ( m_xrrsum + xi * wrri ) * wsum_inv;
202 m_yrravp = ( m_yrrsum + yi * wrri ) * wsum_inv;
203 m_rrrravp = ( m_rrrrsum + rri * wrri ) * wsum_inv;
204}

◆ calculate_average3() [2/6]

void Lpav::calculate_average3 ( double x,
double y,
double w = 1 )

◆ calculate_average3() [3/6]

void Lpav::calculate_average3 ( double x,
double y,
double w = 1 )

◆ calculate_average3() [4/6]

void Lpav::calculate_average3 ( void )

Definition at line 206 of file TrackUtil/src/Lpav.cxx.

206 {
207 if ( m_wsum <= 0 ) return;
208 m_wsum_temp = m_wsum;
209 double wsum_inv( 1 / m_wsum_temp );
210 m_xav = m_xsum * wsum_inv;
211 m_yav = m_ysum * wsum_inv;
212
213 m_rscale = 1;
214 m_cosrot = 1;
215 m_sinrot = 0;
216 m_xxavp = m_xxsum * wsum_inv;
217 m_xyavp = m_xysum * wsum_inv;
218 m_yyavp = m_yysum * wsum_inv;
219 m_xrravp = m_xrrsum * wsum_inv;
220 m_yrravp = m_yrrsum * wsum_inv;
221 m_rrrravp = m_rrrrsum * wsum_inv;
222}

Referenced by fit(), and fit().

◆ calculate_average3() [5/6]

void Lpav::calculate_average3 ( void )

◆ calculate_average3() [6/6]

void Lpav::calculate_average3 ( void )

◆ calculate_lpar() [1/3]

double Lpav::calculate_lpar ( void )

Definition at line 305 of file TrackUtil/src/Lpav.cxx.

305 {
306 double lambda = solve_lambda();
307 // changed on Oct-13-93
308 // if (lambda<=0) return -1;
309 if ( lambda < 0 ) return -1;
310 double h11 = m_xxavp - lambda;
311 double h22 = m_yyavp - lambda;
312 if ( h11 == 0.0 ) return -1;
313 double h14 = m_xrravp;
314 double h24 = m_yrravp;
315 double h34 = 1 + 2 * lambda;
316 double rootsq = ( h14 * h14 / h11 / h11 ) + 4 * h34;
317 if ( std::fabs( h22 ) > std::fabs( h24 ) )
318 {
319 if ( h22 == 0.0 ) return -1;
320 double ratio = h24 / h22;
321 rootsq += ratio * ratio;
322 m_kappa = 1 / std::sqrt( rootsq );
323 m_beta = -ratio * m_kappa;
324 }
325 else
326 {
327 if ( h24 == 0.0 ) return -1;
328 double ratio = h22 / h24;
329 rootsq = 1 + ratio * ratio * rootsq;
330 m_beta = 1 / std::sqrt( rootsq );
331 m_beta = h24 > 0 ? -m_beta : m_beta;
332 m_kappa = -ratio * m_beta;
333 }
334 m_alpha = -( h14 / h11 ) * m_kappa;
335 m_gamma = -h34 * m_kappa;
336 // if (lambda<0.0001) {
337 // cout << " lambda=" << lambda << " h34=" << h34
338 // << " rootsq=" << rootsq << " h22=" << h22
339 // << " h11=" << h11 << " h14=" << h14 << " h24=" << h24 <<
340 // " " << *this << endl;
341 // }
342 //
343 // C TRANSFORM THESE INTO THE LAB COORDINATE SYSTEM
344 //
345 // C FIRST GET KAPPA AND GAMMA BACK TO REAL DIMENSIONS
346 //
347 scale( m_rscale );
348 //
349 // C NEXT ROTATE ALPHA AND BETA
350 //
351 rotate( m_cosrot, -m_sinrot );
352 //
353 // C THEN TRANSLATE BY (XAV,YAV)
354 //
355 move( -m_xav, -m_yav );
356 if ( m_yrravp < 0 ) neg();
357 if ( lambda >= 0 ) m_chisq = lambda * m_wsum_temp * m_rscale * m_rscale;
358 return lambda;
359}

Referenced by fit(), and fit().

◆ calculate_lpar() [2/3]

double Lpav::calculate_lpar ( void )

◆ calculate_lpar() [3/3]

double Lpav::calculate_lpar ( void )

◆ calculate_lpar3() [1/3]

double Lpav::calculate_lpar3 ( void )

Definition at line 361 of file TrackUtil/src/Lpav.cxx.

361 {
362 double lambda = solve_lambda3();
363 // changed on Oct-13-93
364 // if (lambda<=0) return -1;
365 if ( lambda < 0 ) return -1;
366 double h11 = m_xxavp - lambda;
367 double h22 = m_yyavp - lambda;
368 double h14 = m_xrravp;
369 double h24 = m_yrravp;
370 m_gamma = 0;
371 double h12 = m_xyavp;
372 double det = h11 * h22 - h12 * h12;
373 if ( det != 0 )
374 {
375 double r1 = ( h14 * h22 - h24 * h12 ) / ( det );
376 double r2 = ( h24 * h11 - h14 * h12 ) / ( det );
377 double kinvsq = r1 * r1 + r2 * r2;
378 m_kappa = std::sqrt( 1 / kinvsq );
379 if ( h11 != 0 ) m_alpha = -m_kappa * r1;
380 else m_alpha = 1;
381 if ( h22 != 0 ) m_beta = -m_kappa * r2;
382 else m_beta = 1;
383 }
384 else
385 {
386 m_kappa = 0;
387 if ( h11 != 0 && h22 != 0 )
388 {
389 m_beta = 1 / std::sqrt( 1 + h12 * h12 / h11 / h11 );
390 m_alpha = std::sqrt( 1 - m_beta * m_beta );
391 }
392 else if ( h11 != 0 )
393 {
394 m_beta = 1;
395 m_alpha = 0;
396 }
397 else
398 {
399 m_beta = 0;
400 m_alpha = 1;
401 }
402 }
403 if ( ( m_alpha * m_xav + m_beta * m_yav ) * ( m_beta * m_xav - m_alpha * m_yav ) < 0 ) neg();
404 // if (std::fabs(m_alpha)<0.01 && std::fabs(m_beta)<0.01) {
405 // cout << " lambda=" << lambda << " " << *this << endl;
406 // }
407 if ( lambda >= 0 ) m_chisq = lambda * m_wsum_temp * m_rscale * m_rscale;
408 return lambda;
409}

Referenced by fit(), and fit().

◆ calculate_lpar3() [2/3]

double Lpav::calculate_lpar3 ( void )

◆ calculate_lpar3() [3/3]

double Lpav::calculate_lpar3 ( void )

◆ chi_deg() [1/3]

double Lpav::chi_deg ( ) const

Definition at line 666 of file TrackUtil/src/Lpav.cxx.

666 {
667 if ( m_nc <= 3 ) return -1;
668 else return m_chisq / ( m_nc - 3 );
669}

◆ chi_deg() [2/3]

double Lpav::chi_deg ( ) const

◆ chi_deg() [3/3]

double Lpav::chi_deg ( ) const

◆ chisq() [1/3]

double Lpav::chisq ( ) const
inline

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

58{ return m_chisq; }

◆ chisq() [2/3]

double Lpav::chisq ( ) const
inline

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

58{ return m_chisq; }

◆ chisq() [3/3]

double Lpav::chisq ( ) const
inline

Definition at line 58 of file Reconstruction/TrackUtil/include/TrackUtil/Lpav.h.

58{ return m_chisq; }

◆ clear() [1/3]

void Lpav::clear ( )
inline

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

198 {
199 m_wsum = m_xsum = m_ysum = m_xxsum = m_yysum = m_xysum = m_xrrsum = m_yrrsum = m_rrrrsum =
200 m_rscale = m_nc = 0;
201 m_chisq = -1;
202}

Referenced by Lpav().

◆ clear() [2/3]

void Lpav::clear ( )
inline

◆ clear() [3/3]

void Lpav::clear ( )
inline

◆ cov() [1/3]

HepSymMatrix Lpav::cov ( int inv = 0) const

Definition at line 449 of file TrackUtil/src/Lpav.cxx.

452{
453#else
454{
455 HepSymMatrix vret( 4 );
456#endif
457 vret( 1, 1 ) = m_xxsum;
458 vret( 2, 1 ) = m_xysum;
459 vret( 2, 2 ) = m_yysum;
460 vret( 3, 1 ) = m_xsum;
461 vret( 3, 2 ) = m_ysum;
462 vret( 3, 3 ) = m_wsum;
463 vret( 4, 1 ) = m_xrrsum;
464 vret( 4, 2 ) = m_yrrsum;
465 vret( 4, 3 ) = m_xxsum + m_yysum;
466 vret( 4, 4 ) = m_rrrrsum;
467 if ( inv == 0 )
468 {
469 // int i=vret.Inv();
470 int i;
471 vret.invert( i );
472 if ( i != 0 )
473 {
474 std::cerr << "Lpav::cov:could not invert nc=" << m_nc << vret;
475#ifdef HAVE_EXCEPTION
476 THROW( Lpav::cov, Singular );
477#endif
478 }
479 }
480 return vret;
481}
Lpav::cov
HepSymMatrix cov(int=0) const
Definition TrackUtil/src/Lpav.cxx:449

Referenced by cov(), cov_c(), extrapolate(), and similarity().

◆ cov() [2/3]

HepSymMatrix Lpav::cov ( int = 0) const

◆ cov() [3/3]

HepSymMatrix Lpav::cov ( int = 0) const

◆ cov_c() [1/3]

HepSymMatrix Lpav::cov_c ( int inv = 0) const

Definition at line 483 of file TrackUtil/src/Lpav.cxx.

486{
487#else
488{
489 HepSymMatrix vret( 3 );
490#endif
491#ifdef HAVE_EXCEPTION
492 try
493 {
494#endif
495 vret = cov( 1 ).similarity( dldc() );
496#ifdef HAVE_EXCEPTION
497 } catch ( Lpav::Singular )
498 { THROW( Lpav::cov_c1, Singular_c ); }
499#endif
500 if ( inv == 0 )
501 {
502 // int i = vret.Inv();
503 int i;
504 vret.invert( i );
505 if ( i != 0 )
506 {
507 std::cerr << "Lpav::cov_c:could not invert " << vret;
508#ifdef HAVE_EXCEPTION
509 THROW( Lpav::cov_c2, Singular_c );
510#endif
511 }
512 }
513 return vret;
514}
Lpav::Singular_c
Definition InstallArea/x86_64-el9-gcc13-dbg/include/TrackUtil/Lpav.h:72

◆ cov_c() [2/3]

HepSymMatrix Lpav::cov_c ( int = 0) const

◆ cov_c() [3/3]

HepSymMatrix Lpav::cov_c ( int = 0) const

◆ delta_chisq() [1/3]

double Lpav::delta_chisq ( double x,
double y,
double w = 1 ) const

Definition at line 671 of file TrackUtil/src/Lpav.cxx.

671 {
672 double sim = similarity( x, y );
673 if ( sim < 0 ) return -1;
674 double d = d0( x, y );
675 double delta = std::sqrt( d ) * w / ( 1 + sim * w );
676 return delta;
677}
Lpav::similarity
double similarity(double, double) const
Definition TrackUtil/src/Lpav.cxx:550

◆ delta_chisq() [2/3]

double Lpav::delta_chisq ( double x,
double y,
double w = 1 ) const

◆ delta_chisq() [3/3]

double Lpav::delta_chisq ( double x,
double y,
double w = 1 ) const

◆ extrapolate() [1/3]

int Lpav::extrapolate ( double r,
double & phi,
double & dphi ) const

Definition at line 516 of file TrackUtil/src/Lpav.cxx.

516 {
517 double x, y;
518 if ( m_chisq < 0 ) return -1;
519 if ( xy( r, x, y ) != 0 ) return -1;
520 phi = std::atan2( y, x );
521 if ( phi < 0 ) phi += ( 2 * M_PI );
522 HepVector v( 4 );
523 v( 1 ) = x;
524 v( 2 ) = y;
525 v( 3 ) = 1;
526 v( 4 ) = r * r;
527// HepSymMatrix l = cov().similarityT(v);
528#ifdef HAVE_EXCEPTION
529 try
530 {
531#endif
532 // HepSymMatrix l = cov().similarity(v.T());
533 // // cout << "delta d^2=" << l(1,1);
534 // if (l(1,1)>0) {
535 double l = cov().similarity( v );
536 if ( l > 0 )
537 {
538 double ls = std::sqrt( l );
539 dphi = ls / r;
540 // cout << " delta d=" << ls << " dphi=" << dphi;
541 }
542#ifdef HAVE_EXCEPTION
543 } catch ( Lpav::Singular )
544 { return -1; }
545#endif
546 // cout << endl;
547 return 0;
548}
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
M_PI
#define M_PI
Definition TConstant.h:4
Lpar::phi
double phi(double r, int dir=0) const
Definition TrackUtil/src/Lpar.cxx:199

◆ extrapolate() [2/3]

int Lpav::extrapolate ( double ,
double & ,
double &  ) const

◆ extrapolate() [3/3]

int Lpav::extrapolate ( double ,
double & ,
double &  ) const

◆ fit() [1/6]

double Lpav::fit ( void )

Definition at line 430 of file TrackUtil/src/Lpav.cxx.

430 {
431 if ( m_nc <= 3 ) return -1;
432 m_chisq = -1;
433 double q;
434 if ( m_nc < 4 )
435 {
436 calculate_average3();
437 q = calculate_lpar3();
438 if ( q > 0 ) m_chisq = q * m_wsum_temp * m_rscale * m_rscale;
439 }
440 else
441 {
442 calculate_average();
443 q = calculate_lpar();
444 if ( q > 0 ) m_chisq = q * m_wsum_temp * m_rscale * m_rscale;
445 }
446 return m_chisq;
447}
q
****INTEGER imax DOUBLE PRECISION m_pi *DOUBLE PRECISION m_amfin DOUBLE PRECISION m_Chfin DOUBLE PRECISION m_Xenph DOUBLE PRECISION m_sinw2 DOUBLE PRECISION m_GFermi DOUBLE PRECISION m_MfinMin DOUBLE PRECISION m_ta2 INTEGER m_out INTEGER m_KeyFSR INTEGER m_KeyQCD *COMMON c_Semalib $ !copy of input $ !CMS energy $ !beam mass $ !final mass $ !beam charge $ !final charge $ !smallest final mass $ !Z mass $ !Z width $ !EW mixing angle $ !Gmu Fermi $ alphaQED at q
Definition KKsem.h:33
Lpav::calculate_lpar3
double calculate_lpar3(void)
Definition TrackUtil/src/Lpav.cxx:361
Lpav::calculate_average3
void calculate_average3(void)
Definition TrackUtil/src/Lpav.cxx:206
Lpav::calculate_average
void calculate_average(void)
Definition TrackUtil/src/Lpav.cxx:97
Lpav::calculate_lpar
double calculate_lpar(void)
Definition TrackUtil/src/Lpav.cxx:305

Referenced by TCircleFitter::fit().

◆ fit() [2/6]

double Lpav::fit ( )

◆ fit() [3/6]

double Lpav::fit ( )

◆ fit() [4/6]

double Lpav::fit ( double x,
double y,
double w = 1 )

Definition at line 411 of file TrackUtil/src/Lpav.cxx.

411 {
412 if ( m_nc <= 3 ) return -1;
413 m_chisq = -1;
414 double q;
415 if ( m_nc < 4 )
416 {
417 calculate_average3( x, y, w );
418 double q = calculate_lpar3();
419 if ( q > 0 ) m_chisq = q * m_wsum_temp * m_rscale * m_rscale;
420 }
421 else
422 {
423 calculate_average( x, y, w );
424 q = calculate_lpar();
425 if ( q > 0 ) m_chisq = q * m_wsum_temp * m_rscale * m_rscale;
426 }
427 return m_chisq;
428}

◆ fit() [5/6]

double Lpav::fit ( double x,
double y,
double w = 1 )

◆ fit() [6/6]

double Lpav::fit ( double x,
double y,
double w = 1 )

◆ nc() [1/3]

double Lpav::nc ( ) const
inline

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

52{ return m_nc; }

◆ nc() [2/3]

double Lpav::nc ( ) const
inline

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

52{ return m_nc; }

◆ nc() [3/3]

double Lpav::nc ( ) const
inline

Definition at line 52 of file Reconstruction/TrackUtil/include/TrackUtil/Lpav.h.

52{ return m_nc; }

◆ operator+=() [1/3]

const Lpav & Lpav::operator+= ( const Lpav & la1)

Definition at line 622 of file TrackUtil/src/Lpav.cxx.

622 {
623 m_wsum += la1.m_wsum;
624 m_xsum += la1.m_xsum;
625 m_ysum += la1.m_ysum;
626 m_xxsum += la1.m_xxsum;
627 m_yysum += la1.m_yysum;
628 m_xysum += la1.m_xysum;
629 m_xrrsum += la1.m_xrrsum;
630 m_yrrsum += la1.m_yrrsum;
631 m_rrrrsum += la1.m_rrrrsum;
632 m_nc += la1.m_nc;
633 return *this;
634}

◆ operator+=() [2/3]

const Lpav & Lpav::operator+= ( const Lpav & )

◆ operator+=() [3/3]

const Lpav & Lpav::operator+= ( const Lpav & )

◆ operator=() [1/6]

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

◆ operator=() [2/6]

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

◆ operator=() [3/6]

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

◆ operator=() [4/6]

const Lpav & Lpav::operator= ( const Lpav & lp)
inline

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

137 {
138 Lpar::operator=( lp );
139 m_wsum = lp.m_wsum;
140 m_xsum = lp.m_xsum;
141 m_ysum = lp.m_ysum;
142 m_xxsum = lp.m_xxsum;
143 m_yysum = lp.m_yysum;
144 m_xysum = lp.m_xysum;
145 m_xrrsum = lp.m_xrrsum;
146 m_yrrsum = lp.m_yrrsum;
147 m_rrrrsum = lp.m_rrrrsum;
148
149 m_wsum_temp = lp.m_wsum_temp;
150 m_xav = lp.m_xav;
151 m_yav = lp.m_yav;
152 m_xyavp = lp.m_xyavp;
153
154 m_rscale = lp.m_rscale;
155 m_xxavp = lp.m_xxavp;
156 m_yyavp = lp.m_yyavp;
157 m_xrravp = lp.m_xrravp;
158 m_yrravp = lp.m_yrravp;
159 m_rrrravp = lp.m_rrrravp;
160 m_sinrot = lp.m_sinrot;
161 m_cosrot = lp.m_cosrot;
162
163 m_nc = lp.m_nc;
164 m_chisq = lp.m_chisq;
165 return *this;
166}

◆ operator=() [5/6]

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

◆ operator=() [6/6]

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

◆ prob() [1/3]

double Lpav::prob ( ) const

Definition at line 657 of file TrackUtil/src/Lpav.cxx.

657 {
658 if ( m_nc <= 3 ) return 0;
659 if ( m_chisq < 0 ) return 0;
660 float c = m_chisq;
661 int nci = (int)m_nc - 3;
662 double p = (double)prob_( &c, &nci );
663 return p;
664}
prob_
float prob_(float *, int *)

◆ prob() [2/3]

double Lpav::prob ( ) const

◆ prob() [3/3]

double Lpav::prob ( ) const

◆ similarity() [1/3]

double Lpav::similarity ( double x,
double y ) const

Definition at line 550 of file TrackUtil/src/Lpav.cxx.

550 {
551 if ( m_nc <= 3 ) return -1;
552 HepVector v( 4 );
553 v( 1 ) = x;
554 v( 2 ) = y;
555 v( 3 ) = 1;
556 v( 4 ) = x * x + y * y;
557 double l;
558#ifdef HAVE_EXCEPTION
559 try
560 {
561#endif
562 l = cov().similarity( v );
563#ifdef HAVE_EXCEPTION
564 } catch ( Lpav::Singular )
565 { return -1; }
566#endif
567 return l;
568}

Referenced by delta_chisq().

◆ similarity() [2/3]

double Lpav::similarity ( double ,
double  ) const

◆ similarity() [3/3]

double Lpav::similarity ( double ,
double  ) const

◆ operator+ [1/3]

Lpav operator+ ( const Lpav & la1,
const Lpav & la2 )
friend

Definition at line 636 of file TrackUtil/src/Lpav.cxx.

639{
640#else
641{
642 Lpav la;
643#endif
644 la.m_wsum = la1.m_wsum + la2.m_wsum;
645 la.m_xsum = la1.m_xsum + la2.m_xsum;
646 la.m_ysum = la1.m_ysum + la2.m_ysum;
647 la.m_xxsum = la1.m_xxsum + la2.m_xxsum;
648 la.m_yysum = la1.m_yysum + la2.m_yysum;
649 la.m_xysum = la1.m_xysum + la2.m_xysum;
650 la.m_xrrsum = la1.m_xrrsum + la2.m_xrrsum;
651 la.m_yrrsum = la1.m_yrrsum + la2.m_yrrsum;
652 la.m_rrrrsum = la1.m_rrrrsum + la2.m_rrrrsum;
653 la.m_nc = la1.m_nc + la2.m_nc;
654 return la;
655}

◆ operator+ [2/3]

Lpav operator+ ( const Lpav & la1,
const Lpav & la2 )
friend

Definition at line 636 of file TrackUtil/src/Lpav.cxx.

639{
640#else
641{
642 Lpav la;
643#endif
644 la.m_wsum = la1.m_wsum + la2.m_wsum;
645 la.m_xsum = la1.m_xsum + la2.m_xsum;
646 la.m_ysum = la1.m_ysum + la2.m_ysum;
647 la.m_xxsum = la1.m_xxsum + la2.m_xxsum;
648 la.m_yysum = la1.m_yysum + la2.m_yysum;
649 la.m_xysum = la1.m_xysum + la2.m_xysum;
650 la.m_xrrsum = la1.m_xrrsum + la2.m_xrrsum;
651 la.m_yrrsum = la1.m_yrrsum + la2.m_yrrsum;
652 la.m_rrrrsum = la1.m_rrrrsum + la2.m_rrrrsum;
653 la.m_nc = la1.m_nc + la2.m_nc;
654 return la;
655}

◆ operator+ [3/3]

Lpav operator+ ( const Lpav & la1,
const Lpav & la2 )
friend

Definition at line 636 of file TrackUtil/src/Lpav.cxx.

639{
640#else
641{
642 Lpav la;
643#endif
644 la.m_wsum = la1.m_wsum + la2.m_wsum;
645 la.m_xsum = la1.m_xsum + la2.m_xsum;
646 la.m_ysum = la1.m_ysum + la2.m_ysum;
647 la.m_xxsum = la1.m_xxsum + la2.m_xxsum;
648 la.m_yysum = la1.m_yysum + la2.m_yysum;
649 la.m_xysum = la1.m_xysum + la2.m_xysum;
650 la.m_xrrsum = la1.m_xrrsum + la2.m_xrrsum;
651 la.m_yrrsum = la1.m_yrrsum + la2.m_yrrsum;
652 la.m_rrrrsum = la1.m_rrrrsum + la2.m_rrrrsum;
653 la.m_nc = la1.m_nc + la2.m_nc;
654 return la;
655}

◆ operator<< [1/3]

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

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

232 {
233 // o << "wsum=" << a.m_wsum << " xsum=" << a.m_xsum << " ysum=" << a.m_ysum
234 // << " xxsum=" << a.m_xxsum << " xysum=" << a.m_xysum
235 // << " yysum=" << a.m_yysum
236 // << " xrrsum=" << a.m_xrrsum << " yrrsum=" << a.m_yrrsum
237 // << " rrrrsum=" << a.m_rrrrsum;
238 // o << " rscale=" << a.m_rscale
239 // << " xxavp=" << a.m_xxavp << " yyavp=" << a.m_yyavp
240 // << " xrravp=" << a.m_xrravp << " yrravp=" << a.m_yrravp
241 // << " rrrravp=" << a.m_rrrravp << " cosrot=" << a.m_cosrot
242 // << " sinrot=" << a.m_sinrot
243 // << endl;
244 o << " nc=" << a.m_nc << " chisq=" << a.m_chisq << " " << (Lpar&)a;
245 return o;
246}

◆ operator<< [2/3]

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

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

232 {
233 // o << "wsum=" << a.m_wsum << " xsum=" << a.m_xsum << " ysum=" << a.m_ysum
234 // << " xxsum=" << a.m_xxsum << " xysum=" << a.m_xysum
235 // << " yysum=" << a.m_yysum
236 // << " xrrsum=" << a.m_xrrsum << " yrrsum=" << a.m_yrrsum
237 // << " rrrrsum=" << a.m_rrrrsum;
238 // o << " rscale=" << a.m_rscale
239 // << " xxavp=" << a.m_xxavp << " yyavp=" << a.m_yyavp
240 // << " xrravp=" << a.m_xrravp << " yrravp=" << a.m_yrravp
241 // << " rrrravp=" << a.m_rrrravp << " cosrot=" << a.m_cosrot
242 // << " sinrot=" << a.m_sinrot
243 // << endl;
244 o << " nc=" << a.m_nc << " chisq=" << a.m_chisq << " " << (Lpar&)a;
245 return o;
246}

◆ operator<< [3/3]

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

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

232 {
233 // o << "wsum=" << a.m_wsum << " xsum=" << a.m_xsum << " ysum=" << a.m_ysum
234 // << " xxsum=" << a.m_xxsum << " xysum=" << a.m_xysum
235 // << " yysum=" << a.m_yysum
236 // << " xrrsum=" << a.m_xrrsum << " yrrsum=" << a.m_yrrsum
237 // << " rrrrsum=" << a.m_rrrrsum;
238 // o << " rscale=" << a.m_rscale
239 // << " xxavp=" << a.m_xxavp << " yyavp=" << a.m_yyavp
240 // << " xrravp=" << a.m_xrravp << " yrravp=" << a.m_yrravp
241 // << " rrrravp=" << a.m_rrrravp << " cosrot=" << a.m_cosrot
242 // << " sinrot=" << a.m_sinrot
243 // << endl;
244 o << " nc=" << a.m_nc << " chisq=" << a.m_chisq << " " << (Lpar&)a;
245 return o;
246}
The documentation for this class was generated from the following files: