| 122 |
|
// z is possibly symmetric around 0 |
| 123 |
|
z = pos.z(); |
| 124 |
|
|
| 125 |
< |
std::size_t whichRBin = int(r / dr); |
| 126 |
< |
std::size_t whichZBin = int( (zLen/2.0 + z) / dz); |
| 125 |
> |
int whichRBin = int(r / dr); |
| 126 |
> |
int whichZBin = int( (zLen/2.0 + z) / dz); |
| 127 |
|
|
| 128 |
< |
if ((whichRBin < nRBins_) && (whichZBin >= 0) && (whichZBin < nZBins_)) |
| 129 |
< |
histo[whichRBin][whichZBin] += sd->getMass(); |
| 128 |
> |
if ((whichRBin < nRBins_) && (whichZBin >= 0) && (whichZBin < nZBins_)) { |
| 129 |
> |
std::size_t i = static_cast<std::size_t>(whichRBin); |
| 130 |
> |
std::size_t j = static_cast<std::size_t>(whichZBin); |
| 131 |
> |
histo[i][j] += sd->getMass(); |
| 132 |
> |
} |
| 133 |
|
|
| 134 |
|
} |
| 135 |
|
|