jaxspec 0.0.2__py3-none-any.whl → 0.0.4__py3-none-any.whl

jaxspec/model/abc.py

@@ -1,5 +1,6 @@
 from __future__ import annotations
 import haiku as hk
+import jax
 import jax.numpy as jnp
 import networkx as nx
 from haiku._src import base
@@ -106,7 +107,7 @@ class SpectralModel:
     def params(self):
         return self.transformed_func_photon.init(None, jnp.ones(10), jnp.ones(10))
 
-    def photon_flux(self, *args, **kwargs):
+    def photon_flux(self, params, e_low, e_high, n_points=2):
         r"""
         Compute the expected counts between $E_\min$ and $E_\max$ by integrating the model.
 
@@ -114,16 +115,27 @@ class SpectralModel:
         \int _{E_\min}^{E_\max}\text{d}E ~ \mathcal{M}\left( E \right)
         \quad \left[\frac{\text{photons}}{\text{cm}^2\text{s}}\right]$$
 
+        Parameters:
+            params : The parameters of the model.
+            e_low : The lower bound of the energy bins.
+            e_high : The upper bound of the energy bins.
+            n_points : The number of points used to integrate the model in each bin.
+
         !!! info
             This method is internally used in the inference process and should not be used directly. See
-            [`photon_flux`](/references/results/#jaxspec.analysis.results.ChainResult.photon_flux) to compute
+            [`photon_flux`][jaxspec.analysis.results.ChainResult.photon_flux] to compute
             the photon flux associated with a set of fitted parameters in a
-            [`ChainResult`](/references/results/#jaxspec.analysis.results.ChainResult)
+            [`ChainResult`][jaxspec.analysis.results.ChainResult]
             instead.
         """
-        return self.transformed_func_photon.apply(*args, **kwargs)
 
-    def energy_flux(self, *args, **kwargs):
+        params = jax.tree_map(lambda x: jnp.asarray(x), params)
+        e_low = jnp.asarray(e_low)
+        e_high = jnp.asarray(e_high)
+
+        return self.transformed_func_photon.apply(params, e_low, e_high, n_points=n_points)
+
+    def energy_flux(self, params, e_low, e_high, n_points=2):
         r"""
         Compute the expected energy flux between $E_\min$ and $E_\max$ by integrating the model.
 
@@ -131,6 +143,12 @@ class SpectralModel:
         \int _{E_\min}^{E_\max}\text{d}E ~ E ~ \mathcal{M}\left( E \right)
         \quad \left[\frac{\text{keV}}{\text{cm}^2\text{s}}\right]$$
 
+        Parameters:
+            params : The parameters of the model.
+            e_low : The lower bound of the energy bins.
+            e_high : The upper bound of the energy bins.
+            n_points : The number of points used to integrate the model in each bin.
+
         !!! info
             This method is internally used in the inference process and should not be used directly. See
             [`energy_flux`](/references/results/#jaxspec.analysis.results.ChainResult.energy_flux) to compute
@@ -138,7 +156,12 @@ class SpectralModel:
             [`ChainResult`](/references/results/#jaxspec.analysis.results.ChainResult)
             instead.
         """
-        return self.transformed_func_energy.apply(*args, **kwargs)
+
+        params = jax.tree_map(lambda x: jnp.asarray(x), params)
+        e_low = jnp.asarray(e_low)
+        e_high = jnp.asarray(e_high)
+
+        return self.transformed_func_energy.apply(params, e_low, e_high, n_points=n_points)
 
     def build_namespace(self):
         """

jaxspec/model/additive.py

@@ -6,10 +6,9 @@ import haiku as hk
 import jax
 import jax.numpy as jnp
 import jax.scipy as jsp
-import numpy as np
-import importlib.resources
-from jax.lax import dynamic_slice_in_dim as jax_slice
-from functools import partial
+import astropy.units as u
+import astropy.constants
+
 from .abc import ModelComponent
 from haiku.initializers import Constant as HaikuConstant
 from ..util.integrate import integrate_interval
@@ -72,7 +71,7 @@ class Powerlaw(AdditiveComponent):
     ??? abstract "Parameters"
         * $\alpha$ : Photon index of the power law $\left[\text{dimensionless}\right]$
         * $E_0$ : Reference energy fixed at 1 keV $\left[ \mathrm{keV}\right]$
-        * $K$ : Normalization at 1 keV $\left[\frac{\text{photons}}{\text{cm}^2\text{s}}\right]$
+        * $K$ : Normalization at the reference energy (1 keV) $\left[\frac{\text{photons}}{\text{cm}^2\text{s}}\right]$
     """
 
     def continuum(self, energy):
@@ -135,8 +134,8 @@ class Logparabola(AdditiveComponent):
     ??? abstract "Parameters"
         * $a$ : Slope of the LogParabola at the pivot energy $\left[\text{dimensionless}\right]$
         * $b$ : Curve parameter of the LogParabola $\left[\text{dimensionless}\right]$
-        * $K$ : Normalization at the pivot energy $\left[\frac{\text{photons}}{\text{cm}^2\text{s}}\right]$
         * $E_{\text{Pivot}}$ : Pivot energy fixed at 1 keV $\left[ \mathrm{keV}\right]$
+        * $K$ : Normalization at the pivot energy (1keV) $\left[\frac{\text{photons}}{\text{cm}^2\text{s}}\right]$
     """
 
     # TODO : conform with xspec definition
@@ -188,16 +187,10 @@ class Blackbodyrad(AdditiveComponent):
 
 class Gauss(AdditiveComponent):
     r"""
-    A Gaussian line profile
-
-    $$\mathcal{M}\left( E \right) = \frac{K}{\sqrt{2\pi\sigma^2}}\exp\left(\frac{-(E-E_L)^2}{2\sigma^2}\right)$$
+    A Gaussian line profile. If the width is $$\leq 0$$ then it is treated as a delta function.
+    The `Zgauss` variant computes a redshifted Gaussian.
 
-    The primitive is defined as :
-
-    $$
-    \int \mathcal{M}\left( E \right) \text{d}E =
-    \frac{K}{2}\left( 1+\text{erf} \frac{(E-E_L)}{\sqrt{2}\sigma} \right)
-    $$
+    $$\mathcal{M}\left( E \right) = \frac{K}{\sigma \sqrt{2 \pi}}\exp\left(\frac{-(E-E_L)^2}{2\sigma^2}\right)$$
 
     ??? abstract "Parameters"
         * $E_L$ : Energy of the line $\left[\text{keV}\right]$
@@ -213,6 +206,7 @@ class Gauss(AdditiveComponent):
         return norm * jsp.stats.norm.pdf(energy, loc=line_energy, scale=sigma)
 
 
+"""
 class APEC(AdditiveComponent):
     def __init__(self, name="apec"):
         super(APEC, self).__init__(name=name)
@@ -307,6 +301,7 @@ class APEC(AdditiveComponent):
         )
 
         return (continuum + pseudo) * 1e14 * norm
+"""
 
 
 class Cutoffpl(AdditiveComponent):
@@ -319,7 +314,7 @@ class Cutoffpl(AdditiveComponent):
         * $\alpha$ : Photon index of the power law $\left[\text{dimensionless}\right]$
         * $\beta$ : Folding energy of the exponential cutoff $\left[\text{keV}\right]$
         * $E_0$ : Reference energy fixed at 1 keV $\left[ \mathrm{keV}\right]$
-        * $K$ : Normalization at 1 keV $\left[\frac{\text{photons}}{\text{cm}^2\text{s}}\right]$
+        * $K$ : Normalization at the reference energy (1 keV) $\left[\frac{\text{photons}}{\text{cm}^2\text{s}}\right]$
     """
 
     def continuum(self, energy):
@@ -369,23 +364,93 @@ class Diskbb(AdditiveComponent):
     `Diskpbb` with $p=0.75$
 
     ??? abstract "Parameters"
+        * $T_{\text{in}}$ : Temperature at inner disk radius $\left[ \mathrm{keV}\right]$
         * $\text{norm}$ : $\cos i(r_{\text{in}}/d)^{2}$,
         where $r_{\text{in}}$ is "an apparent" inner disk radius $\left[\text{km}\right]$,
-        $d$ the distance to the source in units of $10 \text{kpc}$,
-        $i$ the angle of the disk ($i=0$ is face-on)
-        * $T_{\text{in}}$ : Temperature at inner disk radius $\left[ \mathrm{keV}\right]$
+        $d$ the distance to the source in units of $10 \text{kpc}$, $i$ the angle of the disk ($i=0$ is face-on)
     """
 
     def continuum(self, energy):
-        norm = hk.get_parameter("norm", [], init=HaikuConstant(1))
         p = 0.75
-        tin = hk.get_parameter("Tin", [], init=HaikuConstant(1))
         tout = 0.0
+        tin = hk.get_parameter("Tin", [], init=HaikuConstant(1))
+        norm = hk.get_parameter("norm", [], init=HaikuConstant(1))
 
         # Tout is set to 0 as it is evaluated at R=infinity
         def integrand(kT, e, tin, p):
             return e**2 * (kT / tin) ** (-2 / p - 1) / (jnp.exp(e / kT) - 1)
 
         integral = integrate_interval(integrand)
+        return norm * 2.78e-3 * (0.75 / p) / tin * jnp.vectorize(lambda e: integral(tout, tin, e, tin, p))(energy)
+
+
+class Agauss(AdditiveComponent):
+    r"""
+    A simple Gaussian line profile in wavelength space.
+    If the width is $\leq 0$ then it is treated as a delta function.
+    The `Zagauss` variant computes a redshifted Gaussian.
+
+    $$\mathcal{M}\left( \lambda \right) =
+    \frac{K}{\sigma \sqrt{2 \pi}} \exp\left(\frac{-(\lambda - \lambda_L)^2}{2 \sigma^2}\right)$$
+
+    ??? abstract "Parameters"
+        * $\lambda_L$ : Wavelength of the line in Angström $\left[\unicode{x212B}\right]$
+        * $\sigma$ : Width of the line width in Angström  $\left[\unicode{x212B}\right]$
+        * $K$ : Normalization $\left[\frac{\text{photons}}{\text{cm}^2\text{s}}\right]$
+    """
+
+    def continuum(self, energy) -> (jax.Array, jax.Array):
+        hc = (astropy.constants.h * astropy.constants.c).to(u.angstrom * u.keV).value
+        line_wavelength = hk.get_parameter("Lambda_l", [], init=HaikuConstant(hc))
+        sigma = hk.get_parameter("sigma", [], init=HaikuConstant(0.001))
+        norm = hk.get_parameter("norm", [], init=HaikuConstant(1))
+
+        return norm * jsp.stats.norm.pdf(hc / energy, loc=line_wavelength, scale=sigma)
+
+
+class Zagauss(AdditiveComponent):
+    r"""
+    A redshifted Gaussian line profile in wavelength space.
+    If the width is $\leq 0$ then it is treated as a delta function.
+
+    $$\mathcal{M}\left( \lambda \right) =
+    \frac{K (1+z)}{\sigma \sqrt{2 \pi}} \exp\left(\frac{-(\lambda/(1+z) - \lambda_L)^2}{2 \sigma^2}\right)$$
+
+    ??? abstract "Parameters"
+        * $\lambda_L$ : Wavelength of the line in Angström $\left[\text{\AA}\right]$
+        * $\sigma$ : Width of the line width in Angström $\left[\text{\AA}\right]$
+        * $K$ : Normalization $\left[\frac{\text{photons}}{\text{cm}^2\text{s}}\right]$
+        * $z$ : Redshift [dimensionless]
+    """
+
+    def continuum(self, energy) -> (jax.Array, jax.Array):
+        hc = (astropy.constants.h * astropy.constants.c).to(u.angstrom * u.keV).value
+        line_wavelength = hk.get_parameter("Lambda_l", [], init=HaikuConstant(hc))
+        sigma = hk.get_parameter("sigma", [], init=HaikuConstant(0.001))
+        norm = hk.get_parameter("norm", [], init=HaikuConstant(1))
+        redshift = hk.get_parameter("redshift", [], init=HaikuConstant(0))
+
+        return norm * (1 + redshift) * jsp.stats.norm.pdf((hc / energy) / (1 + redshift), loc=line_wavelength, scale=sigma)
+
+
+class Zgauss(AdditiveComponent):
+    r"""
+    A redshifted Gaussian line profile. If the width is $\leq 0$ then it is treated as a delta function.
+
+    $$\mathcal{M}\left( E \right) =
+    \frac{K}{(1+z) \sigma \sqrt{2 \pi}}\exp\left(\frac{-(E(1+z)-E_L)^2}{2\sigma^2}\right)$$
+
+    ??? abstract "Parameters"
+        * $E_L$ : Energy of the line $\left[\text{keV}\right]$
+        * $\sigma$ : Width of the line $\left[\text{keV}\right]$
+        * $K$ : Normalization $\left[\frac{\text{photons}}{\text{cm}^2\text{s}}\right]$
+        * $z$ : Redshift [dimensionless]
+    """
+
+    def continuum(self, energy) -> (jax.Array, jax.Array):
+        line_energy = hk.get_parameter("E_l", [], init=HaikuConstant(1))
+        sigma = hk.get_parameter("sigma", [], init=HaikuConstant(1))
+        norm = hk.get_parameter("norm", [], init=HaikuConstant(1))
+        redshift = hk.get_parameter("redshift", [], init=HaikuConstant(0))
 
-        return norm * 2.78e-3 * (0.75 / p) / tin * jax.vmap(lambda e: integral(tout, tin, e, tin, p))(energy)
+        return (norm / (1 + redshift)) * jsp.stats.norm.pdf(energy * (1 + redshift), loc=line_energy, scale=sigma)

jaxspec/model/background.py

@@ -61,6 +61,7 @@ class BackgroundWithError(BackgroundModel):
         return countrate
 
 
+'''
 class ConjugateBackground(BackgroundModel):
     r"""
     This class fit an expected rate $\lambda$ in each bin of the background spectrum. Assuming a Gamma prior
@@ -87,15 +88,14 @@ class ConjugateBackground(BackgroundModel):
 
     def numpyro_model(self, energy, observed_counts, name: str = "bkg", observed=True):
         # Gamma in numpyro is parameterized by concentration and rate (alpha/beta)
-        alpha = observed_counts + 1
-        beta = 1
+        # alpha = observed_counts + 1
+        # beta = 1
 
         with numpyro.plate(f"{name}_plate", len(observed_counts)):
-            countrate = numpyro.sample(
-                f"{name}", dist.NegativeBinomialProbs(alpha, probs=beta / (beta + 1)), obs=observed_counts if observed else None
-            )
+            countrate = numpyro.sample(f"{name}", dist.Gamma(2 * observed_counts + 1, 2), obs=None)
 
         return countrate
+'''
 
 
 class GaussianProcessBackground(BackgroundModel):

jaxspec/model/multiplicative.py

@@ -29,9 +29,9 @@ class Expfac(MultiplicativeComponent):
     $$
 
     ??? abstract "Parameters"
-        * $A$ : amplitude of the modification $\left[\text{dimensionless}\right]$
-        * $f$ : exponential factor $\left[\text{keV}^{-1}\right]$
-        * $E_c$ : start energy of modification $\left[\text{keV}\right]$
+        * $A$ : Amplitude of the modification $\left[\text{dimensionless}\right]$
+        * $f$ : Exponential factor $\left[\text{keV}^{-1}\right]$
+        * $E_c$ : Start energy of modification $\left[\text{keV}\right]$
 
     """
 
@@ -45,12 +45,21 @@ class Expfac(MultiplicativeComponent):
 
 class Tbabs(MultiplicativeComponent):
     r"""
-    The Tuebingen-Boulder ISM absorption model.
+    The Tuebingen-Boulder ISM absorption model. This model calculates the cross section for X-ray absorption by the ISM
+    as the sum of the cross sections for X-ray absorption due to the gas-phase ISM, the grain-phase ISM,
+    and the molecules in the ISM.
+
+    $$
+        \mathcal{M}(E) = \exp^{-N_{\text{H}}\sigma(E)}
+    $$
 
     ??? abstract "Parameters"
-        * $N_H$ : equivalent hydrogen column density
+        * $N_{\text{H}}$ : Equivalent hydrogen column density
             $\left[\frac{\text{atoms}~10^{22}}{\text{cm}^2}\right]$
 
+    !!! note
+        Abundances and cross-sections $\sigma$ can be found in Wilms et al. (2000).
+
     """
 
     ref = importlib.resources.files("jaxspec") / "tables/xsect_tbabs_wilm.fits"
@@ -71,7 +80,7 @@ class Phabs(MultiplicativeComponent):
     A photoelectric absorption model.
 
     ??? abstract "Parameters"
-        * $N_H$ : equivalent hydrogen column density
+        * $N_{\text{H}}$ : Equivalent hydrogen column density
             $\left[\frac{\text{atoms}~10^{22}}{\text{cm}^2}\right]$
 
     """
@@ -94,7 +103,7 @@ class Wabs(MultiplicativeComponent):
     A photo-electric absorption using Wisconsin (Morrison & McCammon 1983) cross-sections.
 
     ??? abstract "Parameters"
-        * $N_H$ : equivalent hydrogen column density
+        * $N_{\text{H}}$ : Equivalent hydrogen column density
             $\left[\frac{\text{atoms}~10^{22}}{\text{cm}^2}\right]$
 
     """
@@ -122,9 +131,9 @@ class Gabs(MultiplicativeComponent):
     $$
 
     ??? abstract "Parameters"
-        * $\tau$ : absorption strength $\left[\text{dimensionless}\right]$
-        * $\sigma$ : absorption width $\left[\text{keV}\right]$
-        * $E_0$ : absorption center $\left[\text{keV}\right]$
+        * $\tau$ : Absorption strength $\left[\text{dimensionless}\right]$
+        * $\sigma$ : Absorption width $\left[\text{keV}\right]$
+        * $E_0$ : Absorption center $\left[\text{keV}\right]$
 
     !!! note
         The optical depth at line center is $\tau/(\sqrt{2 \pi} \sigma)$.
@@ -149,8 +158,8 @@ class Highecut(MultiplicativeComponent):
     $$
 
     ??? abstract "Parameters"
-        * $E_c$ : cutoff energy $\left[\text{keV}\right]$
-        * $E_f$ : e-folding energy $\left[\text{keV}\right]$
+        * $E_c$ : Cutoff energy $\left[\text{keV}\right]$
+        * $E_f$ : Folding energy $\left[\text{keV}\right]$
     """
 
     def continuum(self, energy):
@@ -170,9 +179,9 @@ class Zedge(MultiplicativeComponent):
     $$
 
     ??? abstract "Parameters"
-        * $E_c$ : threshold energy
-        * $E_f$ : absorption depth at the threshold
-        * $z$ : redshift
+        * $E_c$ : Threshold energy
+        * $E_f$ : Absorption depth at the threshold
+        * $z$ : Redshift [dimensionless]
     """
 
     def continuum(self, energy):
@@ -181,3 +190,35 @@ class Zedge(MultiplicativeComponent):
         z = hk.get_parameter("z", [], init=HaikuConstant(0))
 
         return jnp.where(energy <= E_c, 1.0, jnp.exp(-D * (energy * (1 + z) / E_c) ** 3))
+
+
+class Tbpcf(MultiplicativeComponent):
+    r"""
+    Partial covering model using the Tuebingen-Boulder ISM absorption model (for more details, see `Tbabs`).
+
+    $$
+        \mathcal{M}(E) = f \exp^{-N_{\text{H}}\sigma(E)} + (1-f)
+    $$
+
+    ??? abstract "Parameters"
+        * $N_{\text{H}}$ : Equivalent hydrogen column density
+            $\left[\frac{\text{atoms}~10^{22}}{\text{cm}^2}\right]$
+        * $f$ : Partial covering fraction, ranges between 0 and 1 [dimensionless]
+
+    !!! note
+        Abundances and cross-sections $\sigma$ can be found in Wilms et al. (2000).
+
+    """
+
+    ref = importlib.resources.files("jaxspec") / "tables/xsect_tbabs_wilm.fits"
+    with importlib.resources.as_file(ref) as path:
+        table = Table.read(path)
+    energy = jnp.asarray(np.array(table["ENERGY"]), dtype=np.float32)
+    sigma = jnp.asarray(np.array(table["SIGMA"]), dtype=np.float32)
+
+    def continuum(self, energy):
+        nh = hk.get_parameter("N_H", [], init=HaikuConstant(1))
+        f = hk.get_parameter("f", [], init=HaikuConstant(0.2))
+        sigma = jnp.interp(energy, self.energy, self.sigma, left=1e9, right=0.0)
+
+        return f * jnp.exp(-nh * sigma) + (1 - f)

{jaxspec-0.0.2.dist-info → jaxspec-0.0.4.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: jaxspec
-Version: 0.0.2
+Version: 0.0.4
 Summary: jaxspec is a bayesian spectral fitting library for X-ray astronomy.
 License: MIT
 Author: sdupourque
@@ -10,11 +10,11 @@ Classifier: License :: OSI Approved :: MIT License
 Classifier: Programming Language :: Python :: 3
 Classifier: Programming Language :: Python :: 3.10
 Classifier: Programming Language :: Python :: 3.11
-Requires-Dist: arviz (>=0.17.0,<0.18.0)
+Requires-Dist: arviz (>=0.17.1,<0.18.0)
 Requires-Dist: astropy (>=6.0.0,<7.0.0)
 Requires-Dist: chainconsumer (>=1.0.0,<2.0.0)
 Requires-Dist: cmasher (>=1.6.3,<2.0.0)
-Requires-Dist: dm-haiku (>=0.0.11,<0.0.12)
+Requires-Dist: dm-haiku (>=0.0.11,<0.0.13)
 Requires-Dist: gpjax (>=0.8.0,<0.9.0)
 Requires-Dist: jax (>=0.4.23,<0.5.0)
 Requires-Dist: jaxlib (>=0.4.23,<0.5.0)
@@ -23,15 +23,19 @@ Requires-Dist: matplotlib (>=3.8.0,<4.0.0)
 Requires-Dist: mkdocstrings (>=0.24.0,<0.25.0)
 Requires-Dist: networkx (>=3.1,<4.0)
 Requires-Dist: numpy (>=1.26.1,<2.0.0)
-Requires-Dist: numpyro (>=0.13.2,<0.14.0)
+Requires-Dist: numpyro (>=0.13.2,<0.15.0)
 Requires-Dist: pandas (>=2.2.0,<3.0.0)
+Requires-Dist: scipy (<1.13)
 Requires-Dist: seaborn (>=0.13.1,<0.14.0)
 Requires-Dist: simpleeval (>=0.9.13,<0.10.0)
+Requires-Dist: sparse (>=0.15.1,<0.16.0)
 Requires-Dist: tinygp (>=0.3.0,<0.4.0)
 Description-Content-Type: text/markdown
 
 # jaxspec
 
+[![PyPI - Version](https://img.shields.io/pypi/v/jaxspec?style=for-the-badge&logo=pypi&color=rgb(37%2C%20150%2C%20190))](https://pypi.org/project/jaxspec/)
+[![Python package](https://img.shields.io/pypi/pyversions/jaxspec?style=for-the-badge)](https://pypi.org/project/jaxspec/)
 [![Read the Docs](https://img.shields.io/readthedocs/jaxspec?style=for-the-badge)](https://jaxspec.readthedocs.io/en/latest/)
 [![Codecov](https://img.shields.io/codecov/c/github/renecotyfanboy/jaxspec?style=for-the-badge)](https://app.codecov.io/gh/renecotyfanboy/jaxspec)
 [![Slack](https://img.shields.io/badge/Slack-4A154B?style=for-the-badge&logo=slack&logoColor=white)](https://join.slack.com/t/jaxspec/shared_invite/zt-2cuxkdl2f-t0EEAKP~HBEHKvIUZJL2sg)

{jaxspec-0.0.2.dist-info → jaxspec-0.0.4.dist-info}/RECORD

@@ -1,8 +1,8 @@
 jaxspec/__init__.py,sha256=fKMzN6U8LmdsTrUn_-w4An0X-ta38xB7AdBjQWbcWZU,147
 jaxspec/analysis/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
 jaxspec/analysis/compare.py,sha256=PZ9WdRe1kdIlNyzoD3vkGvcawe44cToA6KB9FfUdzWg,726
-jaxspec/analysis/results.py,sha256=WIe4OJahQuaEISCyJQ1Cy8btk7Z4NUOuEyD-vTaIiG4,15783
-jaxspec/data/__init__.py,sha256=RGL3JE_gu25W9BWSjxu4-ngPEzEz-Clz-_3oNQF345c,263
+jaxspec/analysis/results.py,sha256=IjTj72rJaovG6wR4iuBB6KJQYBBsCTEhxD00R9rHWs8,15790
+jaxspec/data/__init__.py,sha256=5MmiOFyDBg_lBMN0piUY0ULN0gD_mW33rurVOYXebdA,497
 jaxspec/data/example_data/MOS1.arf,sha256=kBetpqOR1G-bVpuNAzj7q7YqU5fnwp6woW5OAtSGgVk,34560
 jaxspec/data/example_data/MOS1.pha,sha256=fh_2ZFRbq0_c4e-UdocVtNh6ObJSth4HDnFCflyGkqw,83520
 jaxspec/data/example_data/MOS1.rmf,sha256=LRE40iwxxTmKTu0RcLC4iwc1Ds_senrvT1UIcctDCa4,14846400
@@ -18,25 +18,26 @@ jaxspec/data/example_data/PN.pha,sha256=qYXgYHa_Bg06UzHyBBOvqSCKjXcfpsZx6JGKGGBX
 jaxspec/data/example_data/PN.rmf,sha256=kbqe-C2oufc-anmd_gl7h8aKcCCsbFqg3NQGe_nLQoc,3962880
 jaxspec/data/example_data/PN_spectrum_grp20.fits,sha256=qYXgYHa_Bg06UzHyBBOvqSCKjXcfpsZx6JGKGGBXfJA,138240
 jaxspec/data/example_data/PNbackground_spectrum.fits,sha256=VeAX4MGbMkJF_vBJ3_KnouSbmjkWZ8qcT2Z8T2g7H0k,120960
+jaxspec/data/example_data/fakeit.pha,sha256=IhkeWkE-b3ELECd_Uasjo9h3cXgcjCYH20wDpXJ8LMk,60480
 jaxspec/data/grouping.py,sha256=hhgBt-voiH0DDSyePacaIGsaMnrYbJM_-ZeU66keC7I,622
-jaxspec/data/instrument.py,sha256=1O3pWkm9Wb-lkjtQOEL0z84oKer3rACM5c9SwrgS-yI,3827
-jaxspec/data/obsconf.py,sha256=-qlzQlIBDDxlApVm2oaXxUXHPJo6sn6_AZMlCrkQXZ4,4698
-jaxspec/data/observation.py,sha256=Zy7AEGhiKFBYu63saesbscP3Cm0keMpaTusM5X63zyc,5298
-jaxspec/data/ogip.py,sha256=8DYaIOe1aDy7CFPWruX8SMMSz2H6vhABOLLRIQXpHtA,6674
-jaxspec/data/util.py,sha256=l7bqBToox4GeCTd9abZ7fmNby5v3keVlWTGIVk9WeF8,7823
-jaxspec/fit.py,sha256=5G4OdSy9jRbZl2_gMjhMhznJxpUXGAXCPNa28PtFIg8,8986
+jaxspec/data/instrument.py,sha256=0Ef3zhNT7ca-nHtRCGKXpbkyXpgVVpNtO8_XrMGWnyU,3984
+jaxspec/data/obsconf.py,sha256=tnXCXim6eBjZbvNbx2ViRJ3wzQhNacXRcCrLZlPTdY0,7378
+jaxspec/data/observation.py,sha256=1UnFu5ihZp9z-vP_I7tsFY8jhhIJunv46JyuE-acrg0,6394
+jaxspec/data/ogip.py,sha256=d0OEEkznA7s8xPkqPrnfIvedTdIRsQuBnRK4wHoQ17M,8793
+jaxspec/data/util.py,sha256=BsIVpmx2kAvAn-w6uSjC793OmPAfmDKMGGLDuk9yztY,7923
+jaxspec/fit.py,sha256=FhnedL7-_9eovErw2SESH5sbMhU-qJNdkqxkoVsLglo,9531
 jaxspec/model/__init__.py,sha256=AbpHGcgLb-kRsJGnwFEktk7uzpZOCcBY74-YBdrKVGs,1
-jaxspec/model/abc.py,sha256=mVqWo5BtXpEfnLU5X-auC4uNnpFhOSDR-DprOIUlSHk,16982
-jaxspec/model/additive.py,sha256=iBrwscqhNYKExH_VX9OSWzmbvsKZECl41R9gaSmjzDs,15325
-jaxspec/model/background.py,sha256=IcTy7vS_n_KlTSQ7vomWtT8p2FNuZCthmx9KU_B8T4k,6943
+jaxspec/model/abc.py,sha256=O5DwgcXWfuGSYxTk_TmZcZXak1lLuEvBzJXUTd3E9tE,17851
+jaxspec/model/additive.py,sha256=iu_n820tOdKbG5BKsFM-nVOa2_sKPWBSQ_TKS1bljSk,18656
+jaxspec/model/background.py,sha256=zej99rVfcRb75T85o3u4qeYQIgnFwGtxK8niZJ8S5mM,6872
 jaxspec/model/list.py,sha256=0RPAoscVz_zM1CWdx_Gd5wfrQWV5Nv4Kd4bSXu2ayUA,860
-jaxspec/model/multiplicative.py,sha256=5zpaexxVVgLIitGk1BBttzcGeKxeMEuy6WOPUtgGvng,5933
+jaxspec/model/multiplicative.py,sha256=sAKDkiplhdY7TsaPk7gwkR18dcXmU2nytgBCiHNBPMk,7537
 jaxspec/tables/xsect_phabs_aspl.fits,sha256=Pq_7oqYuOmEeCk4f9KVzQtfVdvAj17u2MnENx1uaUBk,86400
 jaxspec/tables/xsect_tbabs_wilm.fits,sha256=PPReRcnWccTE_BKDFLfDposw8Jbu3ms-sIv1UiSkSTU,86400
 jaxspec/tables/xsect_wabs_angr.fits,sha256=mzBzpHejC1LiB_LEv3mvxq4Zq7qPIHGQrExpcCT3QHM,86400
 jaxspec/util/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
 jaxspec/util/integrate.py,sha256=6oGYhRTevgOo47hbzdcomOCShC2cpKgF8AP7bkx-Zsw,4379
-jaxspec-0.0.2.dist-info/LICENSE.md,sha256=2q5XoWzddts5IqzIcgYYMOL21puU3MfO8gvT3Ype1eQ,1073
-jaxspec-0.0.2.dist-info/METADATA,sha256=vH1D1IeCcDSfrcVf_rg5mLfTi4wUSbFQ5rlgHaZWJ58,2845
-jaxspec-0.0.2.dist-info/WHEEL,sha256=d2fvjOD7sXsVzChCqf0Ty0JbHKBaLYwDbGQDwQTnJ50,88
-jaxspec-0.0.2.dist-info/RECORD,,
+jaxspec-0.0.4.dist-info/LICENSE.md,sha256=2q5XoWzddts5IqzIcgYYMOL21puU3MfO8gvT3Ype1eQ,1073
+jaxspec-0.0.4.dist-info/METADATA,sha256=JSU_vEIQjYwXmTQBXPrXL4ZC8Xfg3s60ktPDmBoB_3E,3194
+jaxspec-0.0.4.dist-info/WHEEL,sha256=d2fvjOD7sXsVzChCqf0Ty0JbHKBaLYwDbGQDwQTnJ50,88
+jaxspec-0.0.4.dist-info/RECORD,,