diff --git a/optika/sensors/_materials/_materials.py b/optika/sensors/_materials/_materials.py
index 4aaa14f..b15c11d 100644
--- a/optika/sensors/_materials/_materials.py
+++ b/optika/sensors/_materials/_materials.py
@@ -1,11 +1,13 @@
 import numpy as np
 import astropy.units as u
 import named_arrays as na
+import optika
 
 __all__ = [
     "energy_bandgap",
     "energy_electron_hole",
     "quantum_yield_ideal",
+    "quantum_efficiency_effective",
 ]
 
 energy_bandgap = 1.12 * u.eV
@@ -81,3 +83,301 @@ def quantum_yield_ideal(
     result = np.where(energy > energy_bandgap, result, 0)
 
     return result
+
+
+def quantum_efficiency_effective(
+    wavelength: u.Quantity | na.AbstractScalar,
+    direction: na.AbstractCartesian3dVectorArray,
+    thickness_oxide: u.Quantity | na.AbstractScalar,
+    thickness_implant: u.Quantity | na.AbstractScalar,
+    thickness_substrate: u.Quantity | na.AbstractScalar,
+    cce_backsurface: u.Quantity | na.AbstractScalar,
+    formula_oxide: str = "SiO2",
+    n_ambient: complex | na.AbstractScalar = 1,
+    n_substrate: None | complex | na.AbstractScalar = None,
+    normal: None | na.AbstractCartesian3dVectorArray = None,
+) -> na.AbstractScalar:
+    r"""
+    Calculate the effective quantum efficiency of a back-illuminated detector.
+
+    Parameters
+    ----------
+    wavelength
+        The wavelength of the incident light.
+    direction
+        The propagation direction of the incident light.
+    thickness_oxide
+        The thickness of the oxide layer on the back surface of the sensor.
+    thickness_implant
+        The thickness of the implant layer.
+    thickness_substrate
+        The thickness of the silicon substrate.
+    cce_backsurface
+        The differential charge collection efficiency on the back surface
+        of the sensor.
+    formula_oxide
+        The chemical formula of the oxide layer.
+        Default is silicon dioxide.
+    n_ambient
+        Optional complex refractive index of the ambient propagation medium.
+        Default is :math:`1`, the refractive index of vacuum.
+    n_substrate
+        Optional complex refractive index of the implant region and substrate.
+        If :obj:`None`, then the refractive index of silicon is used.
+    normal
+        The vector perpendicular to the surface of the sensor.
+        If :obj:`None`, then the normal is assumed to be :math:`-\hat{z}`
+
+    Examples
+    --------
+    Reproduce Figure 12 from :cite:t:`Stern1994`, the modeled quantum efficiency
+    of a Tektronix TK512CB :math:`512 \times 512` pixel backilluminated CCD.
+
+    .. jupyter-execute::
+
+        import matplotlib.pyplot as plt
+        import numpy as np
+        import astropy.units as u
+        import named_arrays as na
+        import optika
+
+        # Define an array of wavelengths with which to sample the EQE
+        wavelength = na.geomspace(10, 10000, axis="wavelength", num=1001) * u.AA
+
+        # Assume normal incidence
+        direction = na.Cartesian3dVectorArray(0, 0, 1)
+
+        # Store the fit parameters found in Stern 1994
+        thickness_oxide = 50 * u.AA
+        thickness_implant = 2317 * u.AA
+        thickness_substrate = 7 * u.um
+        cce_backsurface = 0.21
+
+        # Compute the effective quantum efficiency
+        eqe = optika.sensors.quantum_efficiency_effective(
+            wavelength=wavelength,
+            direction=direction,
+            thickness_oxide=thickness_oxide,
+            thickness_implant=thickness_implant,
+            thickness_substrate=thickness_substrate,
+            cce_backsurface=cce_backsurface,
+        )
+
+        # Compute the maximum theoretical quantum efficiency
+        eqe_max = optika.sensors.quantum_efficiency_effective(
+            wavelength=wavelength,
+            direction=direction,
+            thickness_oxide=thickness_oxide,
+            thickness_implant=thickness_implant,
+            thickness_substrate=thickness_substrate,
+            cce_backsurface=1,
+        )
+
+        # Plot the effective and maximum quantum efficiency
+        fig, ax = plt.subplots(constrained_layout=True)
+        na.plt.plot(
+            wavelength,
+            eqe,
+            ax=ax,
+            label="effective quantum efficiency",
+        );
+        na.plt.plot(
+            wavelength,
+            eqe_max,
+            ax=ax,
+            label="maximum quantum efficiency",
+        );
+        ax.set_xscale("log");
+        ax.set_xlabel(f"wavelength ({wavelength.unit:latex_inline})");
+        ax.set_ylabel("efficiency");
+        ax.legend();
+
+    Plot the EQE as a function of wavelength for normal and oblique incidence
+
+    .. jupyter-execute::
+
+        angle = na.linspace(0, 30, axis="angle", num=2) * u.deg
+        direction = na.Cartesian3dVectorArray(
+            x=np.sin(angle),
+            y=0,
+            z=np.cos(angle),
+        )
+
+        eqe = optika.sensors.quantum_efficiency_effective(
+            wavelength=wavelength,
+            direction=direction,
+            thickness_oxide=thickness_oxide,
+            thickness_implant=thickness_implant,
+            thickness_substrate=thickness_substrate,
+            cce_backsurface=cce_backsurface,
+        )
+
+        # Plot the results
+        fig, ax = plt.subplots(constrained_layout=True)
+        angle_str = angle.value.astype(str).astype(object)
+        na.plt.plot(
+            wavelength,
+            eqe,
+            ax=ax,
+            axis="wavelength",
+            label=r"$\theta$ = " + angle_str + f"{angle.unit:latex_inline}",
+        );
+        ax.set_xscale("log");
+        ax.set_xlabel(f"wavelength ({wavelength.unit:latex_inline})");
+        ax.set_ylabel("efficiency");
+        ax.legend();
+
+    Notes
+    -----
+    Our goal is to recover Equation 11 in :cite:t:`Stern1994`, along wtih the
+    correction for photons lost by transmission through the entire CCD substrate.
+    From inspecting Equations 6 and 9 in :cite:t:`Stern1994`,
+    we can see that the effective quantum efficiency is:
+
+    .. math::
+        :label: eqe-definition
+
+        \text{EQE} = T_\lambda \int_0^\infty \alpha \eta(x) e^{-\alpha x} \; dx
+
+    where :math:`T_\lambda` is the net transmission of photons through the backsurface
+    oxide layer (accounting for both absorption and reflections) calculated using
+    :func:`optika.materials.multilayer_efficiency`,
+    :math:`\alpha` is the absorption coefficient of silicon,
+    :math:`x` is the distance from the backsurface,
+    and :math:`\eta(x)` is the differential charge collection efficiency (CCE).
+
+    :cite:t:`Stern1994` assumes that the differential CCE takes the following
+    linear form,
+
+    .. math::
+        :label: differential-cce
+
+        \eta(x) = \begin{cases}
+            \eta_0 + (1 - \eta_0) x / W, & 0 < x < W \\
+            1, & W < x < D, \\
+            0, & D < x
+        \end{cases}
+
+    where :math:`\eta_0` is the differential CCE at the backsurface,
+    :math:`W` is the thickness of the implant region,
+    and :math:`D` is the total thickness of the silicon substrate.
+
+    Plugging Equation :eq:`differential-cce` into Equation :eq:`eqe-definition`
+    and integrating yields
+
+    .. math::
+        :label: eqe
+
+        \text{EQE} &= \alpha T_\lambda \left\{
+                        \int_0^W \left[ \eta_0 + \left( \frac{1 - \eta_0}{W} \right) x \right] e^{-\alpha x} \; dx
+                        + \int_W^D e^{-\alpha x} \; dx
+        \right\} \\
+        &= \alpha T_\lambda \left\{
+            \eta_0 \int_0^W e^{-\alpha x} \; dx
+            + \left( \frac{1 - \eta_0}{W} \right) \int_0^W x e^{-\alpha x} \; dx
+            + \int_W^D e^{-\alpha x} \; dx
+        \right\} \\
+        &= \alpha T_\lambda \left\{
+            -\left[ \frac{\eta_0}{\alpha} e^{-\alpha x} \right|_0^W
+            - \left( \frac{1 - \eta_0}{W} \right) \left[ \left( \frac{\alpha x + 1}{\alpha^2} \right) e^{-\alpha x} \right|_0^W
+            - \left[ \frac{1}{\alpha} e^{-\alpha x} \right|_W^D
+        \right\} \\
+        &= T_\lambda \left\{
+            - \left[ \eta_0 (e^{-\alpha W} - 1) \right]
+            - \left( \frac{1 - \eta_0}{\alpha W} \right) \left[ (\alpha W + 1) e^{-\alpha W} - 1 \right]
+            - \left[ e^{-\alpha D} - e^{-\alpha W} \right]
+        \right\} \\
+        &= T_\lambda \left\{
+            - \eta_0 e^{-\alpha W}
+            + \eta_0
+            - e^{-\alpha W}
+            + \eta_0 e^{-\alpha W}
+            + \left( \frac{1 - \eta_0}{\alpha W} \right) (1 - e^{-\alpha W})
+            - e^{-\alpha D}
+            + e^{-\alpha W}
+        \right\} \\
+        &= T_\lambda \left\{
+            \eta_0
+            + \left( \frac{1 - \eta_0}{\alpha W} \right) (1 - e^{-\alpha W})
+            - e^{-\alpha D}
+        \right\} \\
+
+    Equation :eq:`eqe` is equivalent to Equation 11 in :cite:t:`Stern1994`,
+    with the addition of an :math:`e^{-\alpha W}-e^{-\alpha D}` term which represents photons
+    that traveled all the way through the silicon substrate without interacting.
+
+    Equation :eq:`eqe` is only valid for normally-incident light.
+    We can generalize it to obliquely-incident light by making the substitution
+
+    .. math::
+        :label: x-oblique
+
+        x \rightarrow \frac{x}{\cos \theta}
+
+    where :math:`\theta` is the angle between the propagation direction
+    inside the silicon substrate and the normal vector.
+
+    Substituting :eq:`x-oblique` into Equation :eq:`eqe` and solving yields
+
+    .. math::
+        :label: eqe-oblique
+
+        \text{EQE} = T_\lambda \left\{
+            \eta_0
+            + \left( \frac{1 - \eta_0}{\alpha W \sec \theta} \right) (1 - e^{-\alpha W \sec \theta})
+            - e^{-\alpha D \sec \theta}
+        \right\} \\
+
+    """
+
+    if n_substrate is None:
+        substrate = optika.chemicals.Chemical("Si")
+        index_refraction_substrate = na.interp(
+            x=wavelength,
+            xp=substrate.index_refraction.inputs,
+            fp=substrate.index_refraction.outputs,
+        )
+        wavenumber_substrate = na.interp(
+            x=wavelength,
+            xp=substrate.wavenumber.inputs,
+            fp=substrate.wavenumber.outputs,
+        )
+        n_substrate = index_refraction_substrate + wavenumber_substrate * 1j
+
+    if normal is None:
+        normal = na.Cartesian3dVectorArray(0, 0, -1)
+
+    reflectivity, transmissivity = optika.materials.multilayer_efficiency(
+        material_layers=na.ScalarArray(np.array([formula_oxide]), axes="_layer"),
+        thickness_layers=na.stack([thickness_oxide], axis="_layer"),
+        axis_layers="_layer",
+        wavelength_ambient=wavelength,
+        direction_ambient=direction,
+        n_ambient=n_ambient,
+        n_substrate=n_substrate,
+        normal=normal,
+    )
+
+    wavenumber_substrate = np.imag(n_substrate)
+    absorption_substrate = 4 * np.pi * wavenumber_substrate / wavelength
+
+    direction_substrate = optika.materials.snells_law(
+        wavelength=wavelength,
+        direction=direction,
+        index_refraction=np.real(n_ambient),
+        index_refraction_new=np.real(n_substrate),
+        normal=normal,
+    )
+
+    cos_theta = -direction_substrate @ normal
+
+    z0 = absorption_substrate * thickness_implant / cos_theta
+    exp_z0 = np.exp(-z0)
+
+    term_1 = cce_backsurface
+    term_2 = ((1 - cce_backsurface) / z0) * (1 - exp_z0)
+    term_3 = -np.exp(-absorption_substrate * thickness_substrate / cos_theta)
+
+    result = transmissivity * (term_1 + term_2 + term_3)
+
+    return result
diff --git a/optika/sensors/_materials/_materials_test.py b/optika/sensors/_materials/_materials_test.py
index 13992f1..12fc2ba 100644
--- a/optika/sensors/_materials/_materials_test.py
+++ b/optika/sensors/_materials/_materials_test.py
@@ -18,3 +18,61 @@ def test_quantum_yield_ideal(
 ):
     result = optika.sensors.quantum_yield_ideal(wavelength)
     assert np.all(result == result_expected)
+
+
+@pytest.mark.parametrize(
+    argnames="wavelength",
+    argvalues=[
+        304 * u.AA,
+        na.linspace(100, 200, axis="wavelength", num=4) * u.AA,
+    ],
+)
+@pytest.mark.parametrize(
+    argnames="direction",
+    argvalues=[
+        na.Cartesian3dVectorArray(0, 0, 1),
+    ],
+)
+@pytest.mark.parametrize(
+    argnames="thickness_oxide",
+    argvalues=[
+        10 * u.AA,
+    ],
+)
+@pytest.mark.parametrize(
+    argnames="thickness_implant",
+    argvalues=[
+        1000 * u.AA,
+    ],
+)
+@pytest.mark.parametrize(
+    argnames="thickness_substrate",
+    argvalues=[
+        1 * u.um,
+    ],
+)
+@pytest.mark.parametrize(
+    argnames="cce_backsurface",
+    argvalues=[
+        0.2,
+        1,
+    ],
+)
+def test_quantum_efficiency_effective(
+    wavelength: u.Quantity | na.AbstractScalar,
+    direction: na.AbstractCartesian3dVectorArray,
+    thickness_oxide: u.Quantity | na.AbstractScalar,
+    thickness_implant: u.Quantity | na.AbstractScalar,
+    thickness_substrate: u.Quantity | na.AbstractScalar,
+    cce_backsurface: u.Quantity | na.AbstractScalar,
+):
+    result = optika.sensors.quantum_efficiency_effective(
+        wavelength=wavelength,
+        direction=direction,
+        thickness_oxide=thickness_oxide,
+        thickness_implant=thickness_implant,
+        thickness_substrate=thickness_substrate,
+        cce_backsurface=cce_backsurface,
+    )
+    assert np.all(result >= 0)
+    assert np.all(result <= 1)
diff --git a/optika/sensors/_sensors.py b/optika/sensors/_sensors.py
index b14802b..564bad0 100644
--- a/optika/sensors/_sensors.py
+++ b/optika/sensors/_sensors.py
@@ -3,13 +3,11 @@
 """
 from typing import TypeVar
 import abc
-import numpy as np
 import astropy.units as u
 import named_arrays as na
 import optika
 
 __all__ = [
-    "quantum_efficiency_effective",
     "AbstractImagingSensor",
     "AbstractCCD",
 ]
@@ -18,304 +16,6 @@
 MaterialT = TypeVar("MaterialT", bound=optika.materials.AbstractMaterial)
 
 
-def quantum_efficiency_effective(
-    wavelength: u.Quantity | na.AbstractScalar,
-    direction: na.AbstractCartesian3dVectorArray,
-    thickness_oxide: u.Quantity | na.AbstractScalar,
-    thickness_implant: u.Quantity | na.AbstractScalar,
-    thickness_substrate: u.Quantity | na.AbstractScalar,
-    cce_backsurface: u.Quantity | na.AbstractScalar,
-    formula_oxide: str = "SiO2",
-    n_ambient: complex | na.AbstractScalar = 1,
-    n_substrate: None | complex | na.AbstractScalar = None,
-    normal: None | na.AbstractCartesian3dVectorArray = None,
-) -> na.AbstractScalar:
-    r"""
-    Calculate the effective quantum efficiency of a back-illuminated detector.
-
-    Parameters
-    ----------
-    wavelength
-        The wavelength of the incident light.
-    direction
-        The propagation direction of the incident light.
-    thickness_oxide
-        The thickness of the oxide layer on the back surface of the sensor.
-    thickness_implant
-        The thickness of the implant layer.
-    thickness_substrate
-        The thickness of the silicon substrate.
-    cce_backsurface
-        The differential charge collection efficiency on the back surface
-        of the sensor.
-    formula_oxide
-        The chemical formula of the oxide layer.
-        Default is silicon dioxide.
-    n_ambient
-        Optional complex refractive index of the ambient propagation medium.
-        Default is :math:`1`, the refractive index of vacuum.
-    n_substrate
-        Optional complex refractive index of the implant region and substrate.
-        If :obj:`None`, then the refractive index of silicon is used.
-    normal
-        The vector perpendicular to the surface of the sensor.
-        If :obj:`None`, then the normal is assumed to be :math:`-\hat{z}`
-
-    Examples
-    --------
-    Reproduce Figure 12 from :cite:t:`Stern1994`, the modeled quantum efficiency
-    of a Tektronix TK512CB :math:`512 \times 512` pixel backilluminated CCD.
-
-    .. jupyter-execute::
-
-        import matplotlib.pyplot as plt
-        import numpy as np
-        import astropy.units as u
-        import named_arrays as na
-        import optika
-
-        # Define an array of wavelengths with which to sample the EQE
-        wavelength = na.geomspace(10, 10000, axis="wavelength", num=1001) * u.AA
-
-        # Assume normal incidence
-        direction = na.Cartesian3dVectorArray(0, 0, 1)
-
-        # Store the fit parameters found in Stern 1994
-        thickness_oxide = 50 * u.AA
-        thickness_implant = 2317 * u.AA
-        thickness_substrate = 7 * u.um
-        cce_backsurface = 0.21
-
-        # Compute the effective quantum efficiency
-        eqe = optika.sensors.quantum_efficiency_effective(
-            wavelength=wavelength,
-            direction=direction,
-            thickness_oxide=thickness_oxide,
-            thickness_implant=thickness_implant,
-            thickness_substrate=thickness_substrate,
-            cce_backsurface=cce_backsurface,
-        )
-
-        # Compute the maximum theoretical quantum efficiency
-        eqe_max = optika.sensors.quantum_efficiency_effective(
-            wavelength=wavelength,
-            direction=direction,
-            thickness_oxide=thickness_oxide,
-            thickness_implant=thickness_implant,
-            thickness_substrate=thickness_substrate,
-            cce_backsurface=1,
-        )
-
-        # Plot the effective and maximum quantum efficiency
-        fig, ax = plt.subplots(constrained_layout=True)
-        na.plt.plot(
-            wavelength,
-            eqe,
-            ax=ax,
-            label="effective quantum efficiency",
-        );
-        na.plt.plot(
-            wavelength,
-            eqe_max,
-            ax=ax,
-            label="maximum quantum efficiency",
-        );
-        ax.set_xscale("log");
-        ax.set_xlabel(f"wavelength ({wavelength.unit:latex_inline})");
-        ax.set_ylabel("efficiency");
-        ax.legend();
-
-    Plot the EQE as a function of wavelength for normal and oblique incidence
-
-    .. jupyter-execute::
-
-        angle = na.linspace(0, 30, axis="angle", num=2) * u.deg
-        direction = na.Cartesian3dVectorArray(
-            x=np.sin(angle),
-            y=0,
-            z=np.cos(angle),
-        )
-
-        eqe = optika.sensors.quantum_efficiency_effective(
-            wavelength=wavelength,
-            direction=direction,
-            thickness_oxide=thickness_oxide,
-            thickness_implant=thickness_implant,
-            thickness_substrate=thickness_substrate,
-            cce_backsurface=cce_backsurface,
-        )
-
-        # Plot the results
-        fig, ax = plt.subplots(constrained_layout=True)
-        angle_str = angle.value.astype(str).astype(object)
-        na.plt.plot(
-            wavelength,
-            eqe,
-            ax=ax,
-            axis="wavelength",
-            label=r"$\theta$ = " + angle_str + f"{angle.unit:latex_inline}",
-        );
-        ax.set_xscale("log");
-        ax.set_xlabel(f"wavelength ({wavelength.unit:latex_inline})");
-        ax.set_ylabel("efficiency");
-        ax.legend();
-
-    Notes
-    -----
-    Our goal is to recover Equation 11 in :cite:t:`Stern1994`, along wtih the
-    correction for photons lost by transmission through the entire CCD substrate.
-    From inspecting Equations 6 and 9 in :cite:t:`Stern1994`,
-    we can see that the effective quantum efficiency is:
-
-    .. math::
-        :label: eqe-definition
-
-        \text{EQE} = T_\lambda \int_0^\infty \alpha \eta(x) e^{-\alpha x} \; dx
-
-    where :math:`T_\lambda` is the net transmission of photons through the backsurface
-    oxide layer (accounting for both absorption and reflections) calculated using
-    :func:`optika.materials.multilayer_efficiency`,
-    :math:`\alpha` is the absorption coefficient of silicon,
-    :math:`x` is the distance from the backsurface,
-    and :math:`\eta(x)` is the differential charge collection efficiency (CCE).
-
-    :cite:t:`Stern1994` assumes that the differential CCE takes the following
-    linear form,
-
-    .. math::
-        :label: differential-cce
-
-        \eta(x) = \begin{cases}
-            \eta_0 + (1 - \eta_0) x / W, & 0 < x < W \\
-            1, & W < x < D, \\
-            0, & D < x
-        \end{cases}
-
-    where :math:`\eta_0` is the differential CCE at the backsurface,
-    :math:`W` is the thickness of the implant region,
-    and :math:`D` is the total thickness of the silicon substrate.
-
-    Plugging Equation :eq:`differential-cce` into Equation :eq:`eqe-definition`
-    and integrating yields
-
-    .. math::
-        :label: eqe
-
-        \text{EQE} &= \alpha T_\lambda \left\{
-                        \int_0^W \left[ \eta_0 + \left( \frac{1 - \eta_0}{W} \right) x \right] e^{-\alpha x} \; dx
-                        + \int_W^D e^{-\alpha x} \; dx
-        \right\} \\
-        &= \alpha T_\lambda \left\{
-            \eta_0 \int_0^W e^{-\alpha x} \; dx
-            + \left( \frac{1 - \eta_0}{W} \right) \int_0^W x e^{-\alpha x} \; dx
-            + \int_W^D e^{-\alpha x} \; dx
-        \right\} \\
-        &= \alpha T_\lambda \left\{
-            -\left[ \frac{\eta_0}{\alpha} e^{-\alpha x} \right|_0^W
-            - \left( \frac{1 - \eta_0}{W} \right) \left[ \left( \frac{\alpha x + 1}{\alpha^2} \right) e^{-\alpha x} \right|_0^W
-            - \left[ \frac{1}{\alpha} e^{-\alpha x} \right|_W^D
-        \right\} \\
-        &= T_\lambda \left\{
-            - \left[ \eta_0 (e^{-\alpha W} - 1) \right]
-            - \left( \frac{1 - \eta_0}{\alpha W} \right) \left[ (\alpha W + 1) e^{-\alpha W} - 1 \right]
-            - \left[ e^{-\alpha D} - e^{-\alpha W} \right]
-        \right\} \\
-        &= T_\lambda \left\{
-            - \eta_0 e^{-\alpha W}
-            + \eta_0
-            - e^{-\alpha W}
-            + \eta_0 e^{-\alpha W}
-            + \left( \frac{1 - \eta_0}{\alpha W} \right) (1 - e^{-\alpha W})
-            - e^{-\alpha D}
-            + e^{-\alpha W}
-        \right\} \\
-        &= T_\lambda \left\{
-            \eta_0
-            + \left( \frac{1 - \eta_0}{\alpha W} \right) (1 - e^{-\alpha W})
-            - e^{-\alpha D}
-        \right\} \\
-
-    Equation :eq:`eqe` is equivalent to Equation 11 in :cite:t:`Stern1994`,
-    with the addition of an :math:`e^{-\alpha W}-e^{-\alpha D}` term which represents photons
-    that traveled all the way through the silicon substrate without interacting.
-
-    Equation :eq:`eqe` is only valid for normally-incident light.
-    We can generalize it to obliquely-incident light by making the substitution
-
-    .. math::
-        :label: x-oblique
-
-        x \rightarrow \frac{x}{\cos \theta}
-
-    where :math:`\theta` is the angle between the propagation direction
-    inside the silicon substrate and the normal vector.
-
-    Substituting :eq:`x-oblique` into Equation :eq:`eqe` and solving yields
-
-    .. math::
-        :label: eqe-oblique
-
-        \text{EQE} = T_\lambda \left\{
-            \eta_0
-            + \left( \frac{1 - \eta_0}{\alpha W \sec \theta} \right) (1 - e^{-\alpha W \sec \theta})
-            - e^{-\alpha D \sec \theta}
-        \right\} \\
-
-    """
-
-    if n_substrate is None:
-        substrate = optika.chemicals.Chemical("Si")
-        index_refraction_substrate = na.interp(
-            x=wavelength,
-            xp=substrate.index_refraction.inputs,
-            fp=substrate.index_refraction.outputs,
-        )
-        wavenumber_substrate = na.interp(
-            x=wavelength,
-            xp=substrate.wavenumber.inputs,
-            fp=substrate.wavenumber.outputs,
-        )
-        n_substrate = index_refraction_substrate + wavenumber_substrate * 1j
-
-    if normal is None:
-        normal = na.Cartesian3dVectorArray(0, 0, -1)
-
-    reflectivity, transmissivity = optika.materials.multilayer_efficiency(
-        material_layers=na.ScalarArray(np.array([formula_oxide]), axes="_layer"),
-        thickness_layers=na.stack([thickness_oxide], axis="_layer"),
-        axis_layers="_layer",
-        wavelength_ambient=wavelength,
-        direction_ambient=direction,
-        n_ambient=n_ambient,
-        n_substrate=n_substrate,
-        normal=normal,
-    )
-
-    wavenumber_substrate = np.imag(n_substrate)
-    absorption_substrate = 4 * np.pi * wavenumber_substrate / wavelength
-
-    direction_substrate = optika.materials.snells_law(
-        wavelength=wavelength,
-        direction=direction,
-        index_refraction=np.real(n_ambient),
-        index_refraction_new=np.real(n_substrate),
-        normal=normal,
-    )
-
-    cos_theta = -direction_substrate @ normal
-
-    z0 = absorption_substrate * thickness_implant / cos_theta
-    exp_z0 = np.exp(-z0)
-
-    term_1 = cce_backsurface
-    term_2 = ((1 - cce_backsurface) / z0) * (1 - exp_z0)
-    term_3 = -np.exp(-absorption_substrate * thickness_substrate / cos_theta)
-
-    result = transmissivity * (term_1 + term_2 + term_3)
-
-    return result
-
-
 class AbstractImagingSensor(
     optika.surfaces.AbstractSurface[
         None,
diff --git a/optika/sensors/_tests/test_sensors.py b/optika/sensors/_tests/test_sensors.py
index d418c1b..771a183 100644
--- a/optika/sensors/_tests/test_sensors.py
+++ b/optika/sensors/_tests/test_sensors.py
@@ -5,59 +5,4 @@
 import optika
 
 
-@pytest.mark.parametrize(
-    argnames="wavelength",
-    argvalues=[
-        304 * u.AA,
-        na.linspace(100, 200, axis="wavelength", num=4) * u.AA,
-    ],
-)
-@pytest.mark.parametrize(
-    argnames="direction",
-    argvalues=[
-        na.Cartesian3dVectorArray(0, 0, 1),
-    ],
-)
-@pytest.mark.parametrize(
-    argnames="thickness_oxide",
-    argvalues=[
-        10 * u.AA,
-    ],
-)
-@pytest.mark.parametrize(
-    argnames="thickness_implant",
-    argvalues=[
-        1000 * u.AA,
-    ],
-)
-@pytest.mark.parametrize(
-    argnames="thickness_substrate",
-    argvalues=[
-        1 * u.um,
-    ],
-)
-@pytest.mark.parametrize(
-    argnames="cce_backsurface",
-    argvalues=[
-        0.2,
-        1,
-    ],
-)
-def test_quantum_efficiency_effective(
-    wavelength: u.Quantity | na.AbstractScalar,
-    direction: na.AbstractCartesian3dVectorArray,
-    thickness_oxide: u.Quantity | na.AbstractScalar,
-    thickness_implant: u.Quantity | na.AbstractScalar,
-    thickness_substrate: u.Quantity | na.AbstractScalar,
-    cce_backsurface: u.Quantity | na.AbstractScalar,
-):
-    result = optika.sensors.quantum_efficiency_effective(
-        wavelength=wavelength,
-        direction=direction,
-        thickness_oxide=thickness_oxide,
-        thickness_implant=thickness_implant,
-        thickness_substrate=thickness_substrate,
-        cce_backsurface=cce_backsurface,
-    )
-    assert np.all(result >= 0)
-    assert np.all(result <= 1)
+