The electron amplification and transport within a photomultiplier tube (PMT) has been investigated by developing an in-house Monte Carlo simulation code. The secondary electron emission in the dynodes is implemented via an effective electron model and the Modified Vaughan’s model, whereas the transport is computed with the Boris leapfrog algorithm. The PMT gain, rise time and transit time have been studied as a function of supply voltage and external magnetostatic field. A good agreement with experimental measurements using a Hamamatsu R13408-100 PMT was obtained. The simulations have been conducted following different treatments of the underlying geometry: three-dimensional, two-dimensional and intermediate (2.5D). The validity of these approaches is compared. The developed framework will help in understanding the behavior of PMTs under highly intense and irregular illumination or varying external magnetic fields, as in the case of prompt gamma-ray measurements during pencil-beam proton therapy; and aid in optimizing the design of voltage dividers with behavioral circuit models.