template<class HistT>
| double Impala::Core::Histogram::FitWeibullMarginalOld< HistT >::AndersonDarling | ( | ) | const [inline, private] |
Definition at line 184 of file FitWeibullMarginalOld.h.
References Impala::Core::Histogram::FitWeibullMarginalOld< HistT >::A2(), Impala::Core::Histogram::FitWeibullMarginalOld< HistT >::mBeta, Impala::Core::Histogram::FitWeibullMarginalOld< HistT >::mBins, Impala::Core::Histogram::FitWeibullMarginalOld< HistT >::mDx, Impala::Core::Histogram::FitWeibullMarginalOld< HistT >::mGamma, Impala::Core::Histogram::FitWeibullMarginalOld< HistT >::mHist, Impala::Core::Histogram::FitWeibullMarginalOld< HistT >::mNorm, and Impala::Core::Histogram::FitWeibullMarginalOld< HistT >::mPrecision.
Referenced by Impala::Core::Histogram::FitWeibullMarginalOld< HistT >::A2().
00185 { 00186 double A2 = 0; 00187 double x; 00188 int i, bins = 0; 00189 00190 double Fi = 0.0; 00191 00192 for (i=0, x = mDx; i<mBins; i++, x += mDx) 00193 { 00194 if (mHist[i]>mPrecision) 00195 { 00196 double xt = pow(x/mBeta, mGamma); 00197 double F = exp(-xt); 00198 double dF = (mGamma/mBeta)*pow(x/mBeta, mGamma-1.0)*F*mDx; 00199 00200 Fi += mHist[i]/mNorm; 00201 F = 1.0-F; 00202 if ((F>0.0) && (F<1.0)) 00203 A2 += (pow(Fi-F,2.0)/(F*(1-F)))*dF; 00204 bins++; 00205 } 00206 } 00207 00208 return bins*A2; 00209 }
Here is the call graph for this function:
