13 #include "msdevstudio/MSconfig.h"
30 using namespace hippodraw;
50 namespace dp = hippodraw::DataPoint2DTuple;
65 for(
int i = 0; i < num_bins; i++ )
68 const vector < double > & row =
m_ntuple -> getRow ( i );
71 double x = row [ ix ];
72 double hwidth = 0.5 * row [ ie ];
73 double edgeL = x - hwidth;
74 double edgeR = x + hwidth;
77 double theta =
m_function-> integrate( edgeL, edgeR );
80 double n = 2.0 * hwidth * row [ iy ];
104 sum1 += n * log ( theta );
109 result = sum1 - sum2;
113 return ( -2 * result );