An essential part in the modeling of kilonova lightcurves is the computation of ejecta properties based on the binary parameters. Below are some simple code snippets that allow for such a computation.

 import matplotlib.pyplot as plt
 from matplotlib.colors import BoundaryNorm
 import numpy as np
 from gwemlightcurves.EjectaFits import KaKy2016
 cmap = plt.get_cmap('Blues')
 #
 Q    = np.linspace(1.0, 15.0, num=25)
 C    = np.linspace(0.1, 0.22, num=25)
 x, y = np.meshgrid(Q, C)
 mej  = KaKy2016.calc_meje(x,0,y,1.2,1.0)
 cf   = plt.contourf(x,y, mej,cmap=cmap)
 # Add title and axis names
 plt.title('Ejecta mass estimate for BHNS based on Kawaguchi et al.')
 plt.xlabel('Mass ratio')
 plt.ylabel('Compactness')
 cbar = plt.colorbar(cf)
 

Within this code, we compute the ejecta mass based on the work of Kawaguchi el al. 2016 for various mass ratios and neutron star compactnesses, as shown in the plot below.