#include <zav.h>
Public Member Functions | |
double | a () const |
double | a () const |
double | a () const |
double | a () const |
void | add (double, double, double) |
void | add (double, double, double) |
void | add (double, double, double) |
void | add (double, double, double) |
double | b () const |
double | b () const |
double | b () const |
double | b () const |
double | calculate () |
double | calculate () |
double | calculate () |
double | calculate () |
double | chisq () const |
double | chisq () const |
double | chisq () const |
double | chisq () const |
void | clear (void) |
void | clear (void) |
void | clear (void) |
void | clear (void) |
double | d (double s, double z) const |
double | d (double s, double z) const |
double | d (double s, double z) const |
double | d (double s, double z) const |
int | nc () const |
int | nc () const |
int | nc () const |
int | nc () const |
void | set (const zav *) |
void | set (const zav *) |
void | set (const zav *) |
void | set (const zav *) |
double | z (double s) const |
double | z (double s) const |
double | z (double s) const |
double | z (double s) const |
zav (int) | |
zav () | |
zav (int) | |
zav () | |
zav (int) | |
zav () | |
zav (int) | |
zav () | |
Private Attributes | |
double | _a |
double | _b |
double | _c11 |
double | _c21 |
double | _c22 |
double | _chisq |
int | _nc |
double | _sav |
double | _sig_inv |
double | _ssav |
double | _szav |
double | _w |
double | _zav |
double | _zzav |
Friends | |
std::ostream & | operator<< (std::ostream &, const zav &) |
std::ostream & | operator<< (std::ostream &, const zav &) |
std::ostream & | operator<< (std::ostream &, const zav &) |
std::ostream & | operator<< (std::ostream &, const zav &) |
|
00047 { 00048 _a = _b = _w = _sav = _ssav = _zav = _szav = _zzav = 0; 00049 _chisq = -1; 00050 _c22 = _c21 = _c11 = _sig_inv = 0; 00051 _nc = 0; 00052 }
|
|
00015 {} // dummy for one of the chain constructor
|
|
|
|
00014 {} // dummy for one of the chain constructor
|
|
|
|
00015 {} // dummy for one of the chain constructor
|
|
|
|
00014 {} // dummy for one of the chain constructor
|
|
00020 { return _a; }
|
|
00021 { return _a; }
|
|
00020 { return _a; }
|
|
00021 { return _a; }
|
|
|
|
|
|
|
|
00083 { 00084 _w += w; 00085 double sw = s * w; 00086 _sav += sw; 00087 _ssav += sw * s; 00088 double zw = z * w; 00089 _zav += zw; 00090 _szav += zw * s; 00091 _zzav += zw * z; 00092 _chisq = -1; 00093 _nc++; 00094 }
|
|
00021 { return _b; }
|
|
00022 { return _b; }
|
|
00021 { return _b; }
|
|
00022 { return _b; }
|
|
|
|
|
|
|
|
00115 { 00116 double sig = _ssav * _w - _sav * _sav; 00117 if (sig!=0) { 00118 _sig_inv = 1/sig; 00119 _a = ( _szav * _w - _sav * _zav ) * _sig_inv; 00120 _b = ( _ssav * _zav - _sav * _szav ) * _sig_inv; 00121 _chisq = _zzav - 2 * _a * _szav - 2 * _b * _zav + _a * _a * _ssav 00122 + _b * _b * _w + 2 * _a * _b * _sav; 00123 _c11 = _w * _sig_inv; 00124 _c21 = - _sav * _sig_inv; 00125 _c22 = _ssav * _sig_inv; 00126 } else { 00127 _sig_inv = 0; 00128 _c11 = _c21 = _c22 = 0; 00129 _chisq = -1; 00130 } 00131 if(_nc==2) { 00132 _chisq = 0; 00133 } 00134 return _chisq; 00135 }
|
|
00016 { return _chisq; }
|
|
00017 { return _chisq; }
|
|
00016 { return _chisq; }
|
|
00017 { return _chisq; }
|
|
|
|
|
|
|
|
00054 { 00055 _w = _sav = _ssav = _zav = _szav = _zzav = 0; 00056 _chisq = -1; 00057 _c22 = _c21 = _c11 = _sig_inv = 0; 00058 _nc = 0; 00059 }
|
|
|
|
|
|
|
|
|
|
00024 { return _nc; }
|
|
00025 { return _nc; }
|
|
00024 { return _nc; }
|
|
00025 { return _nc; }
|
|
|
|
|
|
|
|
00061 { 00062 if(c) { 00063 _w = c->_w; 00064 _sav = c->_sav; 00065 _ssav = c->_ssav; 00066 _zav = c->_zav; 00067 _szav = c->_szav; 00068 _zzav = c->_zzav; 00069 _sig_inv = c->_sig_inv; 00070 _c11 = c->_c11; 00071 _c21 = c->_c21; 00072 _c22 = c->_c22; 00073 _nc = c->_nc; 00074 } else { 00075 _w = _sav = _ssav = _zav = _szav = _zzav = 00076 _sig_inv = _c11 = _c21 = _c22 = 0; 00077 _nc = 0; 00078 } 00079 _a = _b = 0; 00080 _chisq = -1; 00081 }
|
|
|
|
|
|
|
|
|
|
00137 { 00138 o << " zav::w=" << z._w << " sav=" << z._sav << " zav=" << z._zav 00139 << " nc=" << z._nc << " chisq=" << z._chisq << " a=" << z._a 00140 << " b=" << z._b << " c11=" << z._c11 << " c21=" << z._c21 00141 << " c22=" << z._c22 << " sig_inv=" << z._sig_inv << std::endl; 00142 return o; 00143 }
|
|
00137 { 00138 o << " zav::w=" << z._w << " sav=" << z._sav << " zav=" << z._zav 00139 << " nc=" << z._nc << " chisq=" << z._chisq << " a=" << z._a 00140 << " b=" << z._b << " c11=" << z._c11 << " c21=" << z._c21 00141 << " c22=" << z._c22 << " sig_inv=" << z._sig_inv << std::endl; 00142 return o; 00143 }
|
|
00137 { 00138 o << " zav::w=" << z._w << " sav=" << z._sav << " zav=" << z._zav 00139 << " nc=" << z._nc << " chisq=" << z._chisq << " a=" << z._a 00140 << " b=" << z._b << " c11=" << z._c11 << " c21=" << z._c21 00141 << " c22=" << z._c22 << " sig_inv=" << z._sig_inv << std::endl; 00142 return o; 00143 }
|
|
00137 { 00138 o << " zav::w=" << z._w << " sav=" << z._sav << " zav=" << z._zav 00139 << " nc=" << z._nc << " chisq=" << z._chisq << " a=" << z._a 00140 << " b=" << z._b << " c11=" << z._c11 << " c21=" << z._c21 00141 << " c22=" << z._c22 << " sig_inv=" << z._sig_inv << std::endl; 00142 return o; 00143 }
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|