Impact of using the UHECR arrival energies to constrain source associations

Code used in arXiv:1811.06464. The directories here are organised to follow the paper and each section is documented using jupyter notebooks. The main part of the project is a joint hierarchical model for the UHECR energies and arrival directions which in implemented in Stan. The model includes UHECR production, propagation and detection effects.

If you use this code in your work please cite the paper and software using the above DOI.

Dependencies

• Stan and PyStan
• CRPropa 3 - for fits to simulations run using CRPropa 3
• basemap - for sky plots
• fancy (v3.0 release) - a python wrapper to make things (a bit) cleaner
• Jupyter - for notebooks

Questions?

If you have any questions feel free to open an issue or get in touch ([email protected]).