Observing the Galactic Underworld: Code supplement to predicting photometry and astrometry from compact remnant microlensing events

Code and data used/generated throughout the research for the Observing the Galactic Underworld paper.

This project contains the following underlying code and data:

• calculate_incidence.py; code which loads in a population of objects, calculates their path on the sky, queries Gaia for objects near that path and then calculates microlensing event.  
• plot_microlensing_results.py; collection of plotting functions used to visualise the distribution of microlensing events. Most plots in the paper are generated by this script.
• Plot Preparation.pynb; Jupyter notebook for exploring the data and generating plots. Any plots in the paper which are not gerneated by plot_microlensing_results.py are generated by this notebook.
• 10000_year_lensing_results_all_GUW_paper_mag.ecsv; microlensing events from due to black holes and neutron stars over a 10,000 year period. The original population was undersampled by 1,000 so dataset is roughly equivalent to 10 years of microlensing data.
• 100000_year_lensing_results_milkyway_all_MW_paper_mag.ecsv; microlensing events from due to stars over a 100,000 year period. The original population was undersampled by 1,000,000 so dataset is roughly equivalent to 0.1 years of microlensing data
• 32_year_Gaia_x101_lensing_results_all_mag.ecsv; microlensing events due to black holes and neutron stars over a 12 year observation period, a 15 year gap and a 5 year further observation period. 101 such periods were combined together to form this file and, to reduce the file size, non events were not recorded. The original population was undersampled by a factor of 1,000 so this is data is undersampled by roughly a factor of 10.
• kicked_remnants_igoshev_young_7.8_DC_integrated_final_ECSN.csv; population of black holes and neutron stars, as outputted by Sweeney et al. (2020). This population is undersampled by a factor of 10,000.
• milkyway_f1e-6.ebf; population of stars, as generated by GALAXIA. This population is undersampled by a factor of 1,000,000.

Data file containing microlensing events from due to black holes, where black holes have a mass of 20 solar masses, over a 10,000 year period is available from the author upon reasonable request. Attempts were made to upload it here but the upload would not complete.