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