Afterglow is an important phenomenon in luminescent materials and can be desired (e.g., persistent phosphors) or undesired (e.g., scintillators). Understanding and predicting afterglow is often based on analysis of thermally stimulated luminescence (TSL) glow curves, assuming the presence of one or more discrete trap states. Here we present a new approach for the description of the time-dependent afterglow from TSL glow curves using a model with a distribution of trap depths. The method is based on the deconvolution of the energy dependent density of occupied traps derived from TSL glow curves using Tikhonov regularization. To test the validity of this new approach, the procedure is applied to experimental TSL and afterglow data for Lu₁Gd₂Ga₃Al₂O₁₂:Ce ceramics codoped with 40 ppm of Yb³⁺ or Eu³⁺ traps. The experimentally measured afterglow curves are compared with simulations based on models with and without the continuous trap depth distribution. The analysis clearly demonstrates the presence of a distribution of trap depths and shows that the new approach gives a more accurate description of the experimentally observed afterglow. The new method will be especially useful in understanding and reducing undesired afterglow in scintillators.