CGEM BOSS 6.6.5.i
BESIII Offline Software System
Loading...
Searching...
No Matches
KalmanFit::Lpav Class Reference

#include <Lpav.h>

+ Inheritance diagram for KalmanFit::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 Lpavoperator= (const Lpav &)
 
const Lpavoperator= (const Lpar &)
 
const Lpavoperator+= (const Lpav &)
 
- Public Member Functions inherited from KalmanFit::Lpar
 Lpar ()
 
virtual ~Lpar ()
 
const Lparoperator= (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)
 

Detailed Description

Constructor & Destructor Documentation

◆ Lpav()

◆ ~Lpav()

Lpav::~Lpav ( )
virtual

Definition at line 62 of file KalFitAlg/KalFitAlg-00-15-20/src/lpav/Lpav.cxx.

63{
64}

Member Function Documentation

◆ add_point()

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

◆ add_point_frac()

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

Definition at line 584 of file KalFitAlg/KalFitAlg-00-15-20/src/lpav/Lpav.cxx.

584 {
585 register double wi = w * a;
586 m_wsum += wi;
587 m_xsum += wi * xi;
588 m_ysum += wi * yi;
589 m_xxsum += wi * xi * xi;
590 m_yysum += wi * yi * yi;
591 m_xysum += wi * xi * yi;
592 register double rri = ( xi * xi + yi * yi );
593 register double wrri = wi * rri;
594 m_xrrsum += wrri * xi;
595 m_yrrsum += wrri * yi;
596 m_rrrrsum += wrri * rri;
597 m_nc += a;
598}

◆ calculate_average() [1/2]

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

Definition at line 87 of file KalFitAlg/KalFitAlg-00-15-20/src/lpav/Lpav.cxx.

87 {
88 if(m_wsum<=0) return;
89 m_wsum_temp = m_wsum + wi;
90 double rri(xi * xi + yi * yi);
91 double wrri(wi * rri);
92 double wsum_inv(1/m_wsum_temp);
93 m_xav = (m_xsum + wi * xi) * wsum_inv;
94 m_yav = (m_ysum + wi * yi) * wsum_inv;
95
96 double xxav((m_xxsum + wi * xi * xi) * wsum_inv);
97 double yyav((m_yysum + wi * yi * yi) * wsum_inv);
98 double xyav((m_xysum + wi * xi * yi) * wsum_inv);
99 double xrrav((m_xrrsum + xi * wrri) * wsum_inv);
100 double yrrav((m_yrrsum + yi * wrri) * wsum_inv);
101 double rrrrav((m_rrrrsum + wrri * rri) * wsum_inv);
102
103 calculate_average_n(xxav, yyav, xyav, xrrav, yrrav, rrrrav);
104
105}

◆ calculate_average() [2/2]

void Lpav::calculate_average ( void )

Definition at line 107 of file KalFitAlg/KalFitAlg-00-15-20/src/lpav/Lpav.cxx.

107 {
108 if(m_wsum<=0) return;
109 m_wsum_temp = m_wsum;
110 double wsum_inv(1/m_wsum_temp);
111 m_xav = m_xsum * wsum_inv;
112 m_yav = m_ysum * wsum_inv;
113
114 double xxav(m_xxsum * wsum_inv);
115 double yyav(m_yysum * wsum_inv);
116 double xyav(m_xysum * wsum_inv);
117 double xrrav(m_xrrsum * wsum_inv);
118 double yrrav(m_yrrsum * wsum_inv);
119 double rrrrav(m_rrrrsum * wsum_inv);
120
121 calculate_average_n(xxav, yyav, xyav, xrrav, yrrav, rrrrav);
122}

Referenced by fit(), and fit().

◆ calculate_average3() [1/2]

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

Definition at line 193 of file KalFitAlg/KalFitAlg-00-15-20/src/lpav/Lpav.cxx.

193 {
194 if(m_wsum<=0) return;
195 m_wsum_temp = m_wsum + wi;
196 double wsum_inv(1/m_wsum_temp);
197 double rri(xi * xi + yi * yi);
198 m_xav = (m_xsum + wi * xi) * wsum_inv;
199 m_yav = (m_ysum + wi * yi) * wsum_inv;
200
201 m_rscale = 1;
202 m_cosrot = 1;
203 m_sinrot = 0;
204 m_xxavp = (m_xxsum + wi * xi * xi) * wsum_inv;
205 m_xyavp = (m_xysum + wi * xi * yi) * wsum_inv;
206 m_yyavp = (m_yysum + wi * yi * yi) * wsum_inv;
207 double wrri(wi * rri);
208 m_xrravp = (m_xrrsum + xi * wrri) * wsum_inv;
209 m_yrravp = (m_yrrsum + yi * wrri) * wsum_inv;
210 m_rrrravp = (m_rrrrsum + rri * wrri) * wsum_inv;
211}

◆ calculate_average3() [2/2]

void Lpav::calculate_average3 ( void )

Definition at line 213 of file KalFitAlg/KalFitAlg-00-15-20/src/lpav/Lpav.cxx.

213 {
214 if(m_wsum<=0) return;
215 m_wsum_temp = m_wsum;
216 double wsum_inv(1/m_wsum_temp);
217 m_xav = m_xsum * wsum_inv;
218 m_yav = m_ysum * wsum_inv;
219
220 m_rscale = 1;
221 m_cosrot = 1;
222 m_sinrot = 0;
223 m_xxavp = m_xxsum * wsum_inv;
224 m_xyavp = m_xysum * wsum_inv;
225 m_yyavp = m_yysum * wsum_inv;
226 m_xrravp = m_xrrsum * wsum_inv;
227 m_yrravp = m_yrrsum * wsum_inv;
228 m_rrrravp = m_rrrrsum * wsum_inv;
229}

Referenced by fit(), and fit().

◆ calculate_lpar()

double Lpav::calculate_lpar ( void )

Definition at line 316 of file KalFitAlg/KalFitAlg-00-15-20/src/lpav/Lpav.cxx.

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

Referenced by fit(), and fit().

◆ calculate_lpar3()

double Lpav::calculate_lpar3 ( void )

Definition at line 369 of file KalFitAlg/KalFitAlg-00-15-20/src/lpav/Lpav.cxx.

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

Referenced by fit(), and fit().

◆ chi_deg()

double Lpav::chi_deg ( ) const

Definition at line 660 of file KalFitAlg/KalFitAlg-00-15-20/src/lpav/Lpav.cxx.

660 {
661 if (m_nc<=3) return -1;
662 else return m_chisq/(m_nc-3);
663}

◆ chisq()

double KalmanFit::Lpav::chisq ( ) const
inline

Definition at line 66 of file KalFitAlg/KalFitAlg-00-15-20/KalFitAlg/lpav/Lpav.h.

66{ return m_chisq; }

◆ clear()

void Lpav::clear ( )
inline

Definition at line 201 of file KalFitAlg/KalFitAlg-00-15-20/KalFitAlg/lpav/Lpav.h.

201 {
202 m_wsum = m_xsum = m_ysum = m_xxsum = m_yysum = m_xysum
203 = m_xrrsum = m_yrrsum = m_rrrrsum = m_rscale = m_nc = 0;
204 m_chisq = -1;
205}

Referenced by Lpav().

◆ cov()

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

Definition at line 444 of file KalFitAlg/KalFitAlg-00-15-20/src/lpav/Lpav.cxx.

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

Referenced by extrapolate(), and similarity().

◆ cov_c()

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

Definition at line 476 of file KalFitAlg/KalFitAlg-00-15-20/src/lpav/Lpav.cxx.

479{
480#else
481{
482 HepSymMatrix vret(3);
483#endif
484#ifdef HAVE_EXCEPTION
485 try {
486#endif
487 vret = cov(1).similarity(dldc());
488#ifdef HAVE_EXCEPTION
489 }
490 catch (Lpav::Singular) {
491 THROW(Lpav::cov_c1,Singular_c);
492 }
493#endif
494 if(inv==0) {
495// int i = vret.Inv();
496 int i;
497 vret.invert(i);
498 if (i!=0) {
499 std::cerr << "Lpav::cov_c:could not invert " << vret;
500#ifdef HAVE_EXCEPTION
501 THROW(Lpav::cov_c2,Singular_c);
502#endif
503 }
504 }
505 return vret;
506}

◆ delta_chisq()

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

Definition at line 665 of file KalFitAlg/KalFitAlg-00-15-20/src/lpav/Lpav.cxx.

665 {
666 double sim = similarity(x,y);
667 if(sim<0) return -1;
668 double d = d0(x,y);
669 double delta = std::sqrt(d) * w / (1 + sim * w);
670 return delta;
671}

◆ extrapolate()

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

Definition at line 508 of file KalFitAlg/KalFitAlg-00-15-20/src/lpav/Lpav.cxx.

508 {
509 double x, y;
510 if (m_chisq<0) return -1;
511 if (xy(r, x, y)!=0) return -1;
512 phi = std::atan2(y,x);
513 if (phi<0) phi += (2*M_PI);
514 HepVector v(4);
515 v(1) = x;
516 v(2) = y;
517 v(3) = 1;
518 v(4) = r * r;
519// HepSymMatrix l = cov().similarityT(v);
520#ifdef HAVE_EXCEPTION
521 try {
522#endif
523// HepSymMatrix l = cov().similarity(v.T());
524// // cout << "delta d^2=" << l(1,1);
525// if (l(1,1)>0) {
526 double l = cov().similarity(v);
527 if(l>0) {
528 double ls = std::sqrt(l);
529 dphi = ls / r;
530 // cout << " delta d=" << ls << " dphi=" << dphi;
531 }
532#ifdef HAVE_EXCEPTION
533 }
534 catch (Lpav::Singular) {
535 return -1;
536 }
537#endif
538// cout << endl;
539 return 0;
540}
Double_t x[10]
**********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
#define M_PI
Definition TConstant.h:4
double phi(double r, int dir=0) const

◆ fit() [1/2]

double Lpav::fit ( void )

Definition at line 428 of file KalFitAlg/KalFitAlg-00-15-20/src/lpav/Lpav.cxx.

428 {
429 if (m_nc<=3) return -1;
430 m_chisq = -1;
431 double q;
432 if (m_nc<4) {
434 q = calculate_lpar3();
435 if (q>0) m_chisq = q * m_wsum_temp * m_rscale * m_rscale;
436 } else {
438 q = calculate_lpar();
439 if (q>0) m_chisq = q * m_wsum_temp * m_rscale * m_rscale;
440 }
441 return m_chisq;
442}
****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

◆ fit() [2/2]

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

Definition at line 412 of file KalFitAlg/KalFitAlg-00-15-20/src/lpav/Lpav.cxx.

412 {
413 if (m_nc<=3) return -1;
414 m_chisq = -1;
415 double q;
416 if (m_nc<4) {
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 } else {
421 calculate_average(x,y,w);
422 q = calculate_lpar();
423 if (q>0) m_chisq = q * m_wsum_temp * m_rscale * m_rscale;
424 }
425 return m_chisq;
426}

◆ nc()

double KalmanFit::Lpav::nc ( ) const
inline

Definition at line 60 of file KalFitAlg/KalFitAlg-00-15-20/KalFitAlg/lpav/Lpav.h.

60{ return m_nc; }

◆ operator+=()

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

Definition at line 616 of file KalFitAlg/KalFitAlg-00-15-20/src/lpav/Lpav.cxx.

616 {
617 m_wsum += la1.m_wsum;
618 m_xsum += la1.m_xsum;
619 m_ysum += la1.m_ysum;
620 m_xxsum += la1.m_xxsum;
621 m_yysum += la1.m_yysum;
622 m_xysum += la1.m_xysum;
623 m_xrrsum += la1.m_xrrsum;
624 m_yrrsum += la1.m_yrrsum;
625 m_rrrrsum += la1.m_rrrrsum;
626 m_nc += la1.m_nc;
627 return *this;
628}

◆ operator=() [1/2]

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

◆ operator=() [2/2]

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

Definition at line 140 of file KalFitAlg/KalFitAlg-00-15-20/KalFitAlg/lpav/Lpav.h.

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

◆ prob()

double Lpav::prob ( ) const

Definition at line 651 of file KalFitAlg/KalFitAlg-00-15-20/src/lpav/Lpav.cxx.

651 {
652 if (m_nc<=3) return 0;
653 if (m_chisq<0) return 0;
654 float c = m_chisq;
655 int nci = (int)m_nc - 3;
656 double p = (double) prob_(&c, &nci);
657 return p;
658}
float prob_(float *, int *)

◆ similarity()

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

Definition at line 542 of file KalFitAlg/KalFitAlg-00-15-20/src/lpav/Lpav.cxx.

542 {
543 if (m_nc<=3) return -1;
544 HepVector v(4);
545 v(1) = x;
546 v(2) = y;
547 v(3) = 1;
548 v(4) = x * x + y * y;
549 double l;
550#ifdef HAVE_EXCEPTION
551 try {
552#endif
553 l = cov().similarity(v);
554#ifdef HAVE_EXCEPTION
555 }
556 catch (Lpav::Singular) {
557 return -1;
558 }
559#endif
560 return l;
561}

Referenced by delta_chisq().

Friends And Related Symbol Documentation

◆ operator+

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

Definition at line 630 of file KalFitAlg/KalFitAlg-00-15-20/src/lpav/Lpav.cxx.

633{
634#else
635{
636 Lpav la;
637#endif
638 la.m_wsum = la1.m_wsum + la2.m_wsum;
639 la.m_xsum = la1.m_xsum + la2.m_xsum;
640 la.m_ysum = la1.m_ysum + la2.m_ysum;
641 la.m_xxsum = la1.m_xxsum + la2.m_xxsum;
642 la.m_yysum = la1.m_yysum + la2.m_yysum;
643 la.m_xysum = la1.m_xysum + la2.m_xysum;
644 la.m_xrrsum = la1.m_xrrsum + la2.m_xrrsum;
645 la.m_yrrsum = la1.m_yrrsum + la2.m_yrrsum;
646 la.m_rrrrsum = la1.m_rrrrsum + la2.m_rrrrsum;
647 la.m_nc = la1.m_nc + la2.m_nc;
648 return la;
649}

◆ operator<<

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

Definition at line 241 of file KalFitAlg/KalFitAlg-00-15-20/src/lpav/Lpav.cxx.

241 {
242// o << "wsum=" << a.m_wsum << " xsum=" << a.m_xsum << " ysum=" << a.m_ysum
243// << " xxsum=" << a.m_xxsum << " xysum=" << a.m_xysum
244// << " yysum=" << a.m_yysum
245// << " xrrsum=" << a.m_xrrsum << " yrrsum=" << a.m_yrrsum
246// << " rrrrsum=" << a.m_rrrrsum;
247// o << " rscale=" << a.m_rscale
248// << " xxavp=" << a.m_xxavp << " yyavp=" << a.m_yyavp
249// << " xrravp=" << a.m_xrravp << " yrravp=" << a.m_yrravp
250// << " rrrravp=" << a.m_rrrravp << " cosrot=" << a.m_cosrot
251// << " sinrot=" << a.m_sinrot
252// << endl;
253 o << " nc=" << a.m_nc << " chisq=" << a.m_chisq << " " << (Lpar&) a;
254 return o;
255}

The documentation for this class was generated from the following files: