An energy sampling method for Monte Carlo integration of the adjoint photon transport equation is described. The method is based on a biasing scheme, making the statistical weight associated with the adjoint transport simulating particle depend only on current energy, without past history memory. Events considered are Compton scattering and absorption. The efficiency of the algorithm presented is investigated through numerical comparisons with other known sampling procedures.