Light Curves

Setup to make the output clean for the docs:

%%capture
from threeML import silence_logs
import warnings
warnings.filterwarnings("ignore")
silence_logs()
import matplotlib.pyplot as plt
%matplotlib inline
from jupyterthemes import jtplot
jtplot.style(context="talk", fscale=1, ticks=True, grid=False)

Gamma-Ray Bursts are transient sources with a typical duration between milliseconds and a few tens of seconds. Therefore they are nicely visible in light curves. In the following we will see how we can get the light curve of a real GRB as seen by an INTEGRAL/SPI detector.

First we have to define the rough time of the GRB.

from astropy.time import Time
grbtime = Time("2012-07-11T02:44:53", format='isot', scale='utc')

Next we need to define the bounds of the energy bins we want to use.

import numpy as np
ebounds = np.geomspace(20,8000,100)

Now we can construct the time series.

from pyspi.utils.data_builder.time_series_builder import TimeSeriesBuilderSPI
det = 0
tsb = TimeSeriesBuilderSPI.from_spi_grb(f"SPIDet{det}",
                                        det,
                                        grbtime,
                                        ebounds=ebounds,
                                        sgl_type="both",
                                        )

We can now plot the light curves for visualization, in which we can clearly see a transient source in this case.

fig = tsb.view_lightcurve(-50,250)