74 |
|
} |
75 |
|
} |
76 |
|
|
77 |
< |
void GofAngle2::initalizeHistogram() { |
77 |
> |
void GofAngle2::initializeHistogram() { |
78 |
|
npairs_ = 0; |
79 |
|
for (unsigned int i = 0; i < histogram_.size(); ++i) |
80 |
|
std::fill(histogram_[i].begin(), histogram_[i].end(), 0); |
124 |
|
RealType cosAngle2 = dot(dipole1, dipole2); |
125 |
|
|
126 |
|
RealType halfBin = (nAngleBins_ - 1) * 0.5; |
127 |
< |
int angleBin1 = halfBin * (cosAngle1 + 1.0); |
128 |
< |
int angleBin2 = halfBin * (cosAngle2 + 1.0); |
127 |
> |
int angleBin1 = int(halfBin * (cosAngle1 + 1.0)); |
128 |
> |
int angleBin2 = int(halfBin * (cosAngle2 + 1.0)); |
129 |
|
|
130 |
|
++histogram_[angleBin1][angleBin2]; |
131 |
|
++npairs_; |