import nexus as nx
import numpy as np
import matplotlib.pyplot as plt
iron = nx.Material.Template(nx.lib.material.Fe)
layer_Fe = nx.Layer(id = "Fe layer",
material = iron,
thickness = 3000)
sample1 = nx.Sample(layers = [layer_Fe])
layer_Pt = nx.Layer(id = "Pt layer",
material = nx.Material.Template(nx.lib.material.Pt),
thickness = 1000)
sample2 = nx.Sample(layers=[layer_Pt])
beam = nx.Beam()
beam.LinearSigma()
exp = nx.Experiment(beam = beam,
objects = [sample1, sample2],
isotope = nx.lib.moessbauer.Fe57)
intensities = nx.Intensities(experiment = exp,
energy = 14400)
print(intensities())