The low energy region (< 50 eV) of the electron energy loss spectrum (EELS) can contain a great deal of spectral detail associated with excitations of the valence electrons. Calculation of the spectra from first principles can assist with interpretation and the most widely used method is the random phase approximation (RPA), usually neglecting local field effects (LFE). For KBr this approach is insufficient due to the importance of quasiparticle and excitonic effects. Calculations including these multi-electron effects are shown to give much improved agreement with the experimental spectra, and the inclusion of spin-orbit coupling (SOC) reproduces the excitonic doublet just above band-edge onset. A review of the complex theory behind these methods is given along with practical guidance on performing these calculations.