We show that the analysis of post-transit photocurrent i(t) in a multi-trapping context to determine the density of trapping states g(E) is capable of resolving features less than kT in width. A commonly used method uses a Laplace inversion of i(t) data giving the well-known result g(E) ∼ t i(t) but employs a delta function approximation for trap release times, which results in loss of energy resolution. We show that it is possible to retain the exponential distribution function for trap release time and solve the multi-trapping rate equations directly, giving significantly improved resolution. The analysis is performed on computer generated posttransit data for distributed and discrete traps, and compared with the earlier method and other related Fourier transform methods for determining g(E). In addition, the versatility of the new method in handling cases with either distributed traps or with discrete traps means that it can be applied to disordered materials or to crystalline materials with well-defined defect levels.