| 138 |
|
} |
| 139 |
|
} |
| 140 |
|
|
| 141 |
< |
for (unsigned int i = 0; i < histo.size(); ++i) { |
| 142 |
< |
RealType ther = dr * (i + 0.5); |
| 143 |
< |
for(unsigned int j = 0; j < histo[i].size(); ++j) { |
| 144 |
< |
if (histo[i][j] <= threshDens_) { |
| 145 |
< |
RealType thez = dz * (j + 0.5); |
| 146 |
< |
cerr << ther << "\t" << thez << "\n"; |
| 147 |
< |
break; |
| 141 |
> |
for (unsigned int j = 0; j < nZBins_; ++j) { |
| 142 |
> |
RealType thez = dz * (j + 0.5); |
| 143 |
> |
bool aboveThresh = false; |
| 144 |
> |
for (unsigned int i = 0; i < nRBins_; ++i) { |
| 145 |
> |
RealType ther = dr * (i + 0.5); |
| 146 |
> |
if (histo[i][j] >= threshDens_) aboveThresh = true; |
| 147 |
> |
|
| 148 |
> |
if (aboveThresh && (histo[i][j] <= threshDens_)) { |
| 149 |
> |
cerr << thez << "\t" << ther << "\n"; |
| 150 |
> |
break; |
| 151 |
|
} |
| 152 |
|
} |
| 153 |
|
} |
| 154 |
< |
|
| 154 |
> |
|
| 155 |
|
// values_.push_back( acos(maxct)*(180.0/M_PI) ); |
| 156 |
|
|
| 157 |
|
} |
