compute_scat

gdt.missions.fermi.gbm.localization.dol.legacy_functions.compute_scat(npoints, rgrid, cenergies, geodir, nai_az, nai_zen, nai_unit_vec, back_unit_vec, front_only=False)[source]

Function to compute response of NaI detectors over a selected set of energy channels. NOTE: Output of this function must be multiplied by a spectral shape in order to provide expected detector rates.

Parameters:

  • npoints (int32) – Number of points in the sky grid

  • rgrid (np.array((2, npoints), float32)) – Grid of (az, zen) positions on the sky

  • cenergies (np.ndarray(2, float32)) – Array with first and last index of chosen energy channel range

  • geodir (np.ndarray(3, np.float64)) – Vector pointing from spacecraft to center of the Earth

  • nai_az (np.ndarray(12, float32)) – Azimuth position of each NaI detector on spacecraft in degrees

  • nai_zen (np.ndarray(12, float32)) – Zenith position of each NaI detector on spacecraft in degrees

  • nai_unit_vec (np.ndarray(12, float32)) – Unit vectors along the direction of each NaI detector on spacecraft

  • back_unit_vec (np.ndarray(12, float32)) – Unit vector along the anti-direction of each NaI detector on spacecraft

  • front_only (bool) – Only compute forward scattering geometry when True

Returns:

  • front_scattered_rates (np.ndarray((12, npoints, 16), float32)) – Array with energy response to forward scattering computed separately for each location on the sky for each NaI detector

  • back_scattered_rates (np.ndarray((12, npoints, 16), float32)) – Array with energy response to backward scattering computed separately for each location on the sky for each NaI detector