===== Examples for the computation of ejecta properties =====

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. 

{{:ejecta_1.png?400|}}