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