Added optika.sensors.TektronixTK512CBMaterial class.

optika/sensors/_materials/__init__.py

@@ -1,3 +1,4 @@
 from ._materials import *
+from ._tektronix_tk512cb import *
 from ._e2v_ccd97 import *
 from ._e2v_ccd_aia import *

optika/sensors/_materials/_tektronix_tk512cb/__init__.py

@@ -0,0 +1 @@
+from ._tektronix_tk512cb import *

optika/sensors/_materials/_tektronix_tk512cb/_tektronix_tk512cb.py

@@ -0,0 +1,137 @@
+import astropy.units as u
+import named_arrays as na
+from .._materials import AbstractStern1994BackilluminatedCCDMaterial
+
+__all__ = [
+    "TektronixTK512CBMaterial",
+]
+
+
+class TektronixTK512CBMaterial(
+    AbstractStern1994BackilluminatedCCDMaterial,
+):
+    """
+    A model of the light-sensitive material of a Tektronix TK512CB sensor based on
+    measurements by :cite:t:`Stern1994`.
+
+    Examples
+    --------
+
+    Plot the measured TK512CB quantum efficiency vs the fitted
+    quantum efficiency calculated using the method of :cite:t:`Stern1994`.
+
+    .. jupyter-execute::
+
+        import matplotlib.pyplot as plt
+        import astropy.units as u
+        import astropy.visualization
+        import named_arrays as na
+        import optika
+
+        # Create a new instance of the e2v CCD97 light-sensitive material
+        material_ccd = optika.sensors.TektronixTK512CBMaterial()
+
+        # Store the wavelengths at which the QE was measured
+        wavelength_measured = material_ccd.quantum_efficiency_measured.inputs
+
+        # Store the QE measurements
+        qe_measured = material_ccd.quantum_efficiency_measured.outputs
+
+        # Define a grid of wavelengths with which to evaluate the fitted QE
+        wavelength_fit = na.geomspace(10, 10000, axis="wavelength", num=1001) * u.AA
+
+        # Evaluate the fitted QE using the given wavelengths
+        qe_fit = material_ccd.quantum_efficiency_effective(
+            rays=optika.rays.RayVectorArray(
+                wavelength=wavelength_fit,
+                direction=na.Cartesian3dVectorArray(0, 0, 1),
+            ),
+            normal=na.Cartesian3dVectorArray(0, 0, -1),
+        )
+
+        # Plot the measured QE vs the fitted QE
+        with astropy.visualization.quantity_support():
+            fig, ax = plt.subplots(constrained_layout=True)
+            na.plt.scatter(
+                wavelength_measured,
+                qe_measured,
+                label="measured",
+            )
+            na.plt.plot(
+                wavelength_fit,
+                qe_fit,
+                label="fit",
+            )
+            ax.set_xscale("log")
+            ax.set_xlabel(f"wavelength ({wavelength_fit.unit:latex_inline})")
+            ax.set_ylabel("quantum efficiency")
+            ax.legend()
+
+    The thickness of the oxide layer found by the fit is
+
+    .. jupyter-execute::
+
+        material_ccd.thickness_oxide
+
+    The thickness of the implant layer found by the fit is
+
+    .. jupyter-execute::
+
+        material_ccd.thickness_implant
+
+    The thickness of the substrate found by the fit is
+
+    .. jupyter-execute::
+
+        material_ccd.thickness_substrate
+
+    And the differential charge collection efficiency at the backsurface
+    found by the fit is
+
+    .. jupyter-execute::
+
+        material_ccd.cce_backsurface
+    """
+
+    @property
+    def quantum_efficiency_measured(self) -> na.FunctionArray:
+        wavelength = [
+            13.3,
+            23.6,
+            44.7,
+            67.6,
+            114.0,
+            135.5,
+            171.4,
+            256.0,
+            303.8,
+            461.0,
+            584.0,
+            736.0,
+            1215.5,
+            2537.0,
+            3650.0,
+            4050.0,
+        ] * u.AA
+        qe = [
+            0.91,
+            0.80,
+            0.48,
+            0.32,
+            0.42,
+            0.86,
+            0.82,
+            0.60,
+            0.58,
+            0.53,
+            0.30,
+            0.085,
+            0.055,
+            0.06,
+            0.09,
+            0.29,
+        ] * u.dimensionless_unscaled
+        return na.FunctionArray(
+            inputs=na.ScalarArray(wavelength, axes="wavelength"),
+            outputs=na.ScalarArray(qe, axes="wavelength"),
+        )

optika/sensors/_materials/_tektronix_tk512cb/_tektronix_tk512cb_test.py

@@ -0,0 +1,15 @@
+import pytest
+import optika
+from .._materials_test import AbstractTestAbstractStern1994BackilluminatedCCDMaterial
+
+
+@pytest.mark.parametrize(
+    argnames="a",
+    argvalues=[
+        optika.sensors.TektronixTK512CBMaterial(),
+    ],
+)
+class TestTektronixTK512CBMaterial(
+    AbstractTestAbstractStern1994BackilluminatedCCDMaterial,
+):
+    pass