computeSelfFields(_cycl). Needed to have individual cell counters
for the three directions, plus last cell was missing.
Added line to take z position of beam into account when doing
fieldsolve with internal/external beam pipe. Here, the transversal
offset should be taken into account as well...TODO!