fig, ax = plt.subplots(1,1,figsize=(10, 7))
ax.plot(times, hemiLuminosity)
ax.set_xlabel("Time [s]")
ax.set_ylabel("Luminsoity [erg s$^{-1}$]")