An analytical model that approximates Bragg curves due to protons in water is extended to lithium fluoride (LiF), both in bulk and lower-density homogeneous thin film form. To that purpose, Monte Carlo simulations of proton range for energies from 5 to 65 MeV are best fitted so that optimal values of the model parameters are deduced for LiF. Tests are performed by comparing analytically computed Bragg curves with corresponding ones simulated by Monte Carlo software in LiF, bulk and lower density, at two example proton energies and energy-straggling widths. The model is then applied to the reconstruction of the energy spectra of proton beams by best fitting experimental Bragg profiles due to visible photoluminescence of optically-active F2 and F3+ color centers created in LiF crystals by the ionizing interaction of protons with the material.