ruling_measurement

esis.flights.f1.optics.gratings.rulings.ruling_measurement(num_distribution=11)

Load a model of the rulings based on the total efficiency of the gratings.

The total efficiency of the gratings is given by efficiency_vs_wavelength(), and the efficiency of the multilayer coating is given by multilayer_fit(). This function takes the ratio of those two functions to estimate the efficiency of the rulings.

Parameters:

num_distribution (int) – The number of Monte Carlo samples to draw when computing uncertainties.

Examples

Compare the as-designed efficiency of the rulings to the measured efficiency of the rulings.

import numpy as np
import matplotlib.pyplot as plt
import astropy.units as u
import astropy.visualization
import named_arrays as na
import optika
from esis.flights.f1.optics import gratings

# Define an array of wavelengths with which to sample the efficiency
wavelength = na.geomspace(500, 700, axis="wavelength", num=1001) * u.AA

# Define the incidence angle to be the same as the Horiba technical proposal
angle = 1.3 * u.deg

# Define the incident rays from the wavelength array
rays = optika.rays.RayVectorArray(
    wavelength=wavelength,
    position=0 * u.mm,
    direction=na.Cartesian3dVectorArray(
        x=np.sin(angle),
        y=0,
        z=np.cos(angle),
    ),
)

# Define the surface normal
normal = na.Cartesian3dVectorArray(0, 0, -1)

# Initialize the ESIS diffraction grating ruling model
ruling_design = gratings.rulings.ruling_design(num_distribution=0)
ruling_measurement = gratings.rulings.ruling_measurement(num_distribution=0)

# Compute the efficiency of the grating rulings
efficiency_design = ruling_design.efficiency(
    rays=rays,
    normal=normal,
)
efficiency_measurement = ruling_measurement.efficiency(
    rays=rays,
    normal=normal,
)

# Plot the efficiency vs wavelength
with astropy.visualization.quantity_support():
    fig, ax = plt.subplots(constrained_layout=True)
    na.plt.plot(
        wavelength,
        efficiency_design,
        ax=ax,
        color="tab:blue",
        label="design",
    );
    na.plt.plot(
        wavelength,
        efficiency_measurement,
        ax=ax,
        color="tab:orange",
        label="measurement",
    );
    ax.set_xlabel(f"wavelength ({ax.get_xlabel()})");
    ax.set_ylabel("efficiency");
    ax.legend();

References to esis.flights.f1.optics.gratings.rulings.ruling_measurement
ruling_measurement