We present a detailed examination of thermalization after inflation for perturbative inflaton decay. Different interactions among particles in the plasma of inflaton decay products are considered and it will be shown that 2 -> 2 scatterings and particle decay are the important ones. We show that thermalization occurs after decays dominate scatterings, and that depending on the typical mass scale of inflaton decay products, different situations may arise. In particular, thermalization may be delayed until late times, in which case the bounds from thermal gravitino production on supersymmetric models of inflation are considerably relaxed. We will also consider the case where the observable sector consists only of the MSSM matter content, and point out that flat directions with large vevs may result in earlier thermalization of the plasma and push the reheat temperature towards its upper limit.