Note
Click here to download the full example code or to run this example in your browser via Binder
Light curve temporal model#
This model parametrises a LightCurve time model.
from astropy.time import Time
from gammapy.modeling.models import (
LightCurveTemplateTemporalModel,
Models,
PowerLawSpectralModel,
SkyModel,
)
time_range = [Time("59100", format="mjd"), Time("59365", format="mjd")]
path = "$GAMMAPY_DATA/tests/models/light_curve/lightcrv_PKSB1222+216.fits"
light_curve_model = LightCurveTemplateTemporalModel.read(path)
light_curve_model.plot(time_range)
<AxesSubplot:xlabel='Time [iso]', ylabel='Norm / A.U.'>
YAML representation#
Here is an example YAML file using the model:
model = SkyModel(
spectral_model=PowerLawSpectralModel(),
temporal_model=light_curve_model,
name="light_curve_model",
)
models = Models([model])
print(models.to_yaml())
components:
- name: light_curve_model
type: SkyModel
spectral:
type: PowerLawSpectralModel
parameters:
- name: index
value: 2.0
- name: amplitude
value: 1.0e-12
unit: cm-2 s-1 TeV-1
- name: reference
value: 1.0
unit: TeV
temporal:
type: LightCurveTemplateTemporalModel
parameters:
- name: t_ref
value: 59000.5
unit: d
filename: /Users/terrier/Code/gammapy-dev/gammapy-data/tests/models/light_curve/lightcrv_PKSB1222+216.fits
format: table
unit: ''