template<class HistT>
| double Impala::Core::Histogram::FitWeibullMarginal< HistT >::AndersonDarling | ( | ) | const [inline, private] |
Definition at line 252 of file FitWeibullMarginal.h.
References Impala::Core::Histogram::FitWeibullMarginal< HistT >::A2(), Impala::Core::Histogram::FitWeibullMarginal< HistT >::mBeta, Impala::Core::Histogram::FitWeibullMarginal< HistT >::mBins, Impala::Core::Histogram::FitWeibullMarginal< HistT >::mDx, Impala::Core::Histogram::FitWeibullMarginal< HistT >::mGamma, Impala::Core::Histogram::FitWeibullMarginal< HistT >::mHist, Impala::Core::Histogram::FitWeibullMarginal< HistT >::mNorm, and Impala::Core::Histogram::FitWeibullMarginal< HistT >::mPrecision.
Referenced by Impala::Core::Histogram::FitWeibullMarginal< HistT >::A2().
00253 { 00254 double A2 = 0; 00255 double x; 00256 int i, n=0; 00257 00258 double Fi = 0.0; 00259 00260 for (i=0, x = mDx; i<mBins; i++, x += mDx) 00261 { 00262 if (mHist[i] > mPrecision) 00263 { 00264 double xt = pow(x/mBeta, mGamma); 00265 double F = exp(-xt); 00266 double dF = (mGamma/mBeta)*pow(x/mBeta, mGamma-1.0)*F*mDx; 00267 00268 Fi += mHist[i]/mNorm; 00269 F = 1.0-F; 00270 if ((F>0.0) && (F<1.0)) 00271 { 00272 A2 += ((Fi-F)*(Fi-F)/(F*(1-F)))*dF; 00273 n++; 00274 } 00275 } 00276 } 00277 00278 return n*A2; 00279 }
Here is the call graph for this function:
