…to include long range interactions

…he scattering rates

This slows down the interpolation of the scattering rates considerably, but it is necessary for long range interactions. This should be optimized later.

…ass for light mediators

… the computation of the scattering rate tabulation

…he relative speed. Use the quicker simulation function if it doesn't.

…function including unit tests

…ge x.

…q) from previously sampling the scattering angle and the target velocity.

NOT COMPLETE YET

Medium function still accounts for electron density only (and wrongly at that for nuclei, needs to be fixed later)

…not in trajectory simulator class)

…for vDM = 0.0 (which happens during interpolation).

…transform sampling function for q sampling

NOTE: Medium function not included in scattering rates at this point

…Min and qMax are now their own arguments

…s can act as a target argument for the scattering rates and reduce the number of arguments of these functions a lot!

… is found dynamically.

…functions. I want to give both std::vector<Isotope> and std::vector<Solar_Isotope> as argument. Not too happy with this solution, but not sure about a better alternative.

…s to the sampling of q.

For a positive cos_theta, the pdf of q is 0 everywhere, and this gave an error. But cos_theta was sampled from a wrong PDF (due to a faulty integration using Gauss-Legendre instead of Gauss-Kronrod).

… the simulation code. Update the main function for testing of sampling.

…ith progress bar

… used in the simulations anymore, might delete later.)

… (if possible, i.e. if the nucleus and electron rate scale identically).

… zero when computing the total rate.

… section (e or N) is zero.

In that case a rescaling is also possible and should be preferred over a new tabulation.

…nfig file

…ts to be shown in the terminal

…for the fact when the number of scatterings reaches the maximum value allowed.

Particles that scatter more often are considered as "captured" in the statistics.