===== Examples for the computation of lightcurves =====

This code snippet should allow the computation of the lightcurves using the simple, analytical model derived by Brian Metzger: 

   import matplotlib
   import matplotlib.pyplot as plt
   matplotlib.use('TkAgg')
   %matplotlib inline
   from gwemlightcurves.sampler.model import Me2017_model
   t_Me, lbol_Me, mag_Me = Me2017_model(mej=1e-2,vej=0.1,beta=3,kappa_r=25)
   for i in range(1,8):
       plt.plot(t_Me, mag_Me[i])
   plt.xlim(1, 14)
   plt.ylim(-10, -15)
   plt.show();


{{:lightcurve_1.png?400|}}


Similarly you could also call other models, e.g., here the one derived by [[https://arxiv.org/pdf/1612.03665.pdf|Dietrich and Ujevic in 2017]]:

   from gwemlightcurves.sampler.model import DiUj2017_model_ejecta
   import numpy as np
   t, lbol, mag = DiUj2017_model_ejecta(mej=1e-1,vej=0.1,th=0.3,ph=90)
   print(np.size(mag[1]))
   for i in range(1,8):
       plt.plot(t, mag[i])
   plt.xlim(4, 20)
   plt.ylim(-10, -18)
   plt.show();

{{:lightcurve_2.png?400|}}