CRPropa / CRPropa3

CRPropa is a public astrophysical simulation framework for propagating extraterrestrial ultra-high energy particles. https://crpropa.github.io/CRPropa3/
https://crpropa.desy.de
GNU General Public License v3.0
68 stars 67 forks source link

New feature in CRPropa 3.2: Targeting #380

Closed wilson-2015 closed 2 years ago

wilson-2015 commented 2 years ago

The simulations of UHECRs are emitted isotropically from their sources and undergo deflections through IGMF before detection at the observer sphere. Obtaining enough statistics in this scenario is usually computationally intensive. Personally I have been solving this problem by calculating the cosmic ray emission angle at the source and injecting particles through a small opening angle.

A new approach has been described in CRPropa 3.2 known as targeting method that maximizes the number of cosmic ray hits, based on von Misses Fisher distribution. This method is also described in this paper POS(ICRC2019)447.

I have looked around and could not get a way of implementing this in my simulation script. Any further guideline on how to do this or a snippet code will be highly appreciated.

lukasmerten commented 2 years ago

The CRPropa 3.2 has not yet been published but will contain a more detailed description of the targeting feature. In addition, there is already an example ipython notebook included in the current master branch.

wilson-2015 commented 2 years ago

Thanks Lukas, this is exactly what I was looking for.