In this work, we addressed the inhomogeneity problem in gamma spectrometry caused by hot particles, which are dispersed into environment from large nuclear reactor accidents such as at Chernobyl and Fukushima. Using Monte Carlo simulation, we have determined the response of a gamma spectrometer to individual and grouped hot particles randomly distributed in a soil matrix of 1-L and 0.6-L sample containers. By exploring the fact that the peak-to-total ratio of efficiencies in gamma spectrometry is an empirical parameter, we derived and verified a power-law relationship between the peak efficiency and peak-to-total ratio. This enabled creation of a novel calibration model which was demonstrated to reduce the bias range and bias standard deviation, caused by measuring hot particles, by several times, as compared with the homogeneous calibration. The new model is independent of the number, location, and distribution of hot particles in the samples. In this work, we demonstrated successful performance of the model for a single-peak 137 Cs radionuclide. An extension to multi-peak radionuclide was also derived.