X-Ray Absorption in Cu-W

import roentgen
import roentgen.absorption
import astropy.units as u
import numpy as np
import matplotlib.pyplot as plt
from astropy.visualization import quantity_support

alloy = roentgen.absorption.Material(
    {'Cu': 0.25, 'W': 0.75},
    thickness=350*u.micron,
    density=14.5*u.g/(u.cm**3),
)

energy = np.linspace(1,400,1000) * u.keV

with quantity_support():
    plt.plot(energy, alloy.absorption(energy))