Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Lightmediator #13

Merged
merged 62 commits into from
Apr 8, 2024
Merged

Lightmediator #13

merged 62 commits into from
Apr 8, 2024

Conversation

temken
Copy link
Owner

@temken temken commented Apr 8, 2024

No description provided.

Timon Emken and others added 30 commits August 10, 2022 08:54
…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.
… the computation of the scattering rate tabulation
…he relative speed. Use the quicker simulation function if it doesn't.
…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)
…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
…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.
… used in the simulations anymore, might delete later.)
… (if possible, i.e. if the nucleus and electron rate scale identically).
… section (e or N) is zero.

In that case a rescaling is also possible and should be preferred over a new tabulation.
…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.
@temken temken merged commit e768346 into dev Apr 8, 2024
1 check passed
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

1 participant