Purpose: In proton therapy, complex density heterogeneities within the beam path constitute a challenge to dose calculation algorithms. This might question the reliability of dose distributions predicted by treatment planning systems based on analytical dose calculation. For cases in which substantial dose errors are expected, resorting to Monte Carlo dose calculations might be essential to ensure a successful treatment outcome and therefore the benefit is worth a presumably long computation time. The aim of this study was to define an indicator for the accuracy of dose delivery based on analytical dose calculations in treatment planning systems for small proton therapy fields to identify those patients for which Monte Carlo dose calculation is warranted. Methods: Fourteen patients treated at our facility with small passively scattered proton beams (apertures diameters below 7 cm) were selected. Plans were generated in the XiO treatment planning system in combination with a pencil beam algorithm developed at the Massachusetts General Hospital and compared to Monte Carlo dose calculations. Differences in the dose to the 50% of the gross tumor volume (D50, GTV) were assessed in a field-by-field basis. A simple and fast methodology was developed to quantify the inhomogeneity of the tissue traversed by a single small proton beam using a heterogeneity index (HI)-a concept presented by Plugfelder et al. [Med. Phys. 34, 1506-1513 (2007)] for scanned proton beams. Finally, the potential correlation between the error made by the pencil beam based treatment planning algorithm for each field and the level of tissue heterogeneity traversed by the proton beam given by the HI was evaluated. Results: Discrepancies up to 5.4% were found in D50 for single fields, although dose differences were within clinical tolerance levels (<3%) when combining all of the fields involved in the treatment. The discrepancies found for each field exhibited a strong correlation to their associated HI-values (Spearman's ρ = 0.8, p < 0.0001); the higher the level of tissue inhomogeneities for a particular field, the larger the error by the analytical algorithm. With the established correlation a threshold for HI can be set by choosing a tolerance level of 2-3%-commonly accepted in radiotherapy. Conclusions: The HI is a good indicator for the accuracy of proton field delivery in terms of GTV prescription dose coverage when small fields are delivered. Each HI-value was obtained from the CT image in less than 3 min on a computer with 2 GHz CPU allowing implementation of this methodology in clinical routine. For HI-values exceeding the threshold, either a change in beam direction (if feasible) or a recalculation of the dose with Monte Carlo would be highly recommended.