72 |
|
} |
73 |
|
} |
74 |
|
|
75 |
< |
void GofAngle2::initalizeHistogram() { |
75 |
> |
void GofAngle2::initializeHistogram() { |
76 |
|
npairs_ = 0; |
77 |
|
for (unsigned int i = 0; i < histogram_.size(); ++i) |
78 |
|
std::fill(histogram_[i].begin(), histogram_[i].end(), 0); |
108 |
|
RealType cosAngle2 = dot(dipole1, dipole2); |
109 |
|
|
110 |
|
RealType halfBin = (nAngleBins_ - 1) * 0.5; |
111 |
< |
int angleBin1 = halfBin * (cosAngle1 + 1.0); |
112 |
< |
int angleBin2 = halfBin * (cosAngle2 + 1.0); |
111 |
> |
int angleBin1 = int(halfBin * (cosAngle1 + 1.0)); |
112 |
> |
int angleBin2 = int(halfBin * (cosAngle2 + 1.0)); |
113 |
|
|
114 |
|
++histogram_[angleBin1][angleBin2]; |
115 |
|
++npairs_; |