The correlation function measured in ultrarelativistic nuclear collisions is non-Gaussian. By making use of models we discuss and assess how much various effects can influence its shape. In particular, we focus on the parametrisations expressed with the help of Lévy-stable distributions. We show that the Lévy index may deviate substantially from 2 due to non-critical effects such as non-spherical shape, resonance decays, event-by-event fluctuations and functional dependence on Qinv or similar.