| 134 |
|
void pAngle::processHistogram() { |
| 135 |
|
|
| 136 |
|
int atot = 0; |
| 137 |
< |
for(int i = 0; i < count_.size(); ++i) |
| 137 |
> |
for(unsigned int i = 0; i < count_.size(); ++i) |
| 138 |
|
atot += count_[i]; |
| 139 |
|
|
| 140 |
< |
for(int i = 0; i < count_.size(); ++i) { |
| 140 |
> |
for(unsigned int i = 0; i < count_.size(); ++i) { |
| 141 |
|
histogram_[i] = double(count_[i] / double(atot)); |
| 142 |
|
} |
| 143 |
|
} |
| 152 |
|
rdfStream << "#selection: (" << selectionScript_ << ")\n"; |
| 153 |
|
rdfStream << "#cos(theta)\tp(cos(theta))\n"; |
| 154 |
|
RealType dct = 2.0 / histogram_.size(); |
| 155 |
< |
for (int i = 0; i < histogram_.size(); ++i) { |
| 155 |
> |
for (unsigned int i = 0; i < histogram_.size(); ++i) { |
| 156 |
|
RealType ct = -1.0 + (2.0 * i + 1) / (histogram_.size()); |
| 157 |
|
rdfStream << ct << "\t" << histogram_[i]/dct << "\n"; |
| 158 |
|
} |