Photothermal therapy using nanoparticles is a promising new approach for the treatment of cancer. The principle is to utilise plasmonic nanoparticle light interaction for efficient heat conversion. However, there are many hurdles to overcome before it can be accepted in clinical practice. One issue is a current poor characterization of the thermal dose that is distributed over the tumour region and the surrounding normal tissue. Here, we use Monte Carlo simulations of photon radiative transfer through tissue and subsequent heat diffusion calculations, to model the spatial thermal dose in a skin cancer model. We validate our heat rise simulations against experimental data from the literature and estimate the concentration of nanorods in the tumor that are associated with the heat rise. We use the cumulative equivalent minutes at 43 °C (CEM43) metric to analyse the percentage cell kill across the tumour and the surrounding normal tissue. Overall, we show that computer simulations of photothermal therapy are an invaluable tool to fully characterize thermal dose within tumour and normal tissue.