get_albedo_matrix#

sunkit_spex.models.physical.albedo.get_albedo_matrix( energy_edges: Annotated[Quantity, Unit('keV')], theta: Annotated[Quantity, Unit('deg')], anisotropy=1, )[source]#

Get albedo correction matrix.

Matrix used to correct a photon spectrum for the component reflected by the solar atmosphere following interpolated to given angle and energy indices.

Parameters:

energy_edges – Energy edges associated with the spectrum
theta – Angle between Sun-observer line and X-ray source
anisotropy – Ratio of the flux in observer direction to the flux downwards, 1 for an isotropic source

Example

>>> import astropy.units as u
>>> import numpy as np
>>> from sunkit_spex.models.physical.albedo import get_albedo_matrix
>>> e = np.linspace(5,  500, 5)*u.keV
>>> albedo_matrix = get_albedo_matrix(e,theta=45*u.deg)
>>> albedo_matrix
array([[7.64944936e-03, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00],
       [7.17787454e-01, 1.54795970e-10, 0.00000000e+00, 0.00000000e+00],
       [5.22059171e-01, 3.02951100e-01, 1.46291699e-13, 0.00000000e+00],
       [4.52582540e-01, 3.69821128e-01, 1.13435321e-01, 5.95953019e-15]])