| 100 |
|
seleMan_.setSelectionSet(evaluator_.evaluate()); |
| 101 |
|
} |
| 102 |
|
|
| 103 |
+ |
|
| 104 |
+ |
int runningTot = 0; |
| 105 |
|
for (sd = seleMan_.beginSelected(i); sd != NULL; |
| 106 |
|
sd = seleMan_.nextSelected(i)) { |
| 107 |
< |
|
| 107 |
> |
|
| 108 |
|
Vector3d pos = sd->getPos(); |
| 109 |
|
|
| 110 |
|
Vector3d r1 = CenterOfMass - pos; |
| 136 |
|
atot += count_[i]; |
| 137 |
|
|
| 138 |
|
for(int i = 0; i < count_.size(); ++i) { |
| 139 |
< |
histogram_[i] = double(count_[i] / atot); |
| 139 |
> |
histogram_[i] = double(count_[i]) / double(atot); |
| 140 |
|
} |
| 141 |
|
} |
| 142 |
|
|
| 149 |
|
rdfStream << "#nFrames:\t" << nProcessed_ << "\n"; |
| 150 |
|
rdfStream << "#selection: (" << selectionScript_ << ")\n"; |
| 151 |
|
rdfStream << "#cos(theta)\tp(cos(theta))\n"; |
| 152 |
+ |
RealType dct = 2.0 / histogram_.size(); |
| 153 |
|
for (int i = 0; i < histogram_.size(); ++i) { |
| 154 |
< |
RealType ct = -1.0 + i / histogram_.size(); |
| 155 |
< |
rdfStream << ct << "\t" << histogram_[i] << "\n"; |
| 154 |
> |
RealType ct = -1.0 + (2.0 * i + 1) / (histogram_.size()); |
| 155 |
> |
rdfStream << ct << "\t" << histogram_[i]/dct << "\n"; |
| 156 |
|
} |
| 157 |
|
|
| 158 |
|
} else { |
