| 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_; |