\end{cases} """ # %% # Example plot # ------------ # Here is an example plot of the model: from astropy import units as u import matplotlib.pyplot as plt from gammapy.modeling.models import BrokenPowerLawSpectralModel, Models, SkyModel energy_range = [0.1, 100] * u.TeV model = BrokenPowerLawSpectralModel( index1=1.5, index2=2.5, amplitude="1e-12 TeV-1 cm-2 s-1", ebreak="1 TeV", ) model.plot(energy_range) plt.grid(which="both") # %% # YAML representation # ------------------- # Here is an example YAML file using the model: model = SkyModel(spectral_model=model, name="broken-power-law-model") models = Models([model]) print(models.to_yaml())
index=1.8 * u.Unit(""), amplitude=4 / u.cm ** 2 / u.s / u.TeV, reference=1 * u.TeV, lambda_=0.1 / u.TeV, alpha=0.8, ), val_at_2TeV=u.Quantity(0.871694294554192, "cm-2 s-1 TeV-1"), integral_1_10TeV=u.Quantity(3.026342, "cm-2 s-1"), eflux_1_10TeV=u.Quantity(7.38652453, "TeV cm-2 s-1"), e_peak=1.7677669529663684 * u.TeV, ), dict( name="bpl", model=BrokenPowerLawSpectralModel( index1=1.5 * u.Unit(""), index2=2.5 * u.Unit(""), amplitude=4 / u.cm ** 2 / u.s / u.TeV, ebreak=0.5 * u.TeV, ), val_at_2TeV=u.Quantity(0.125, "cm-2 s-1 TeV-1"), integral_1_10TeV=u.Quantity(0.45649740094103286, "cm-2 s-1"), eflux_1_10TeV=u.Quantity(0.9669999668731384, "TeV cm-2 s-1"), ), dict( name="sbpl", model=SmoothBrokenPowerLawSpectralModel( index1=1.5 * u.Unit(""), index2=2.5 * u.Unit(""), amplitude=4 / u.cm ** 2 / u.s / u.TeV, ebreak=0.5 * u.TeV, reference=1 * u.TeV, beta=1,