Impala Documentation

template<class HistT>

double Impala::Core::Histogram::FitWeibull< HistT >::AndersonDarling ( ) const [inline, private]

Definition at line 167 of file FitWeibull.h.

References Impala::Core::Histogram::FitWeibull< HistT >::A2(), Impala::Core::Histogram::FitWeibull< HistT >::mBeta, Impala::Core::Histogram::FitWeibull< HistT >::mBins, Impala::Core::Histogram::FitWeibull< HistT >::mDx, Impala::Core::Histogram::FitWeibull< HistT >::mGamma, Impala::Core::Histogram::FitWeibull< HistT >::mHist, Impala::Core::Histogram::FitWeibull< HistT >::mNorm, and Impala::Core::Histogram::FitWeibull< HistT >::mPrecision.

Referenced by Impala::Core::Histogram::FitWeibull< HistT >::A2().

00168     { 
00169         double A2 = 0;
00170         double x;
00171         int i, n=0;
00172 
00173         double Fi = 0.0;
00174 
00175         for (i=0, x = mDx; i<mBins; i++, x += mDx)
00176         {
00177             if (mHist[i] > mPrecision)
00178             {
00179                 double xt = pow(x/mBeta, mGamma);
00180                 double F = exp(-xt);
00181                 double dF = (mGamma/mBeta)*pow(x/mBeta, mGamma-1.0)*F*mDx;
00182 
00183                 Fi += mHist[i]/mNorm;
00184                 F = 1.0-F;
00185                 if ((F>0.0) && (F<1.0)) {
00186                     A2 += ( (Fi-F)*(Fi-F)/(F*(1-F)))*dF;
00187                     n++;
00188                 }
00189             }
00190         }
00191 
00192         return n*A2;
00193     }

Here is the call graph for this function: