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: