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_; |
123 |
|
rdfStream << "selection2: (" << selectionScript2_ << ")\n"; |
124 |
|
rdfStream << "#nAngleBins =" << nAngleBins_ << "deltaCosAngle = " << deltaCosAngle_ << "\n"; |
125 |
|
for (unsigned int i = 0; i < avgGofr_.size(); ++i) { |
126 |
< |
RealType cosAngle1 = -1.0 + (i + 0.5)*deltaCosAngle_; |
127 |
< |
|
126 |
> |
// RealType cosAngle1 = -1.0 + (i + 0.5)*deltaCosAngle_; |
127 |
> |
|
128 |
|
for(unsigned int j = 0; j < avgGofr_[i].size(); ++j) { |
129 |
< |
RealType cosAngle2 = -1.0 + (j + 0.5)*deltaCosAngle_; |
129 |
> |
// RealType cosAngle2 = -1.0 + (j + 0.5)*deltaCosAngle_; |
130 |
|
rdfStream <<avgGofr_[i][j]/nProcessed_ << "\t"; |
131 |
|
} |
132 |
|
|