Simple Sample

import nexus as nx
import numpy as np
import matplotlib.pyplot as plt

site = nx.Hyperfine(magnetic_field = 33,
                        isotropic = True)

sample = nx.SimpleSample(thickness=3000,
                         composition = [["Fe", 1]],
                         density = 7.874,
                         isotope = nx.lib.moessbauer.Fe57,
                         abundance = 0.02119,
                         lamb_moessbauer = 0.796,
                         hyperfine_sites = [site])

beam = nx.Beam(polarization = 0)

exp = nx.Experiment(beam = beam,
                    objects = [sample],
                    isotope = nx.lib.moessbauer.Fe57)

velocities = np.linspace(-10, 10, 512)

moessbauer_spectrum = nx.MoessbauerSpectrum(experiment = exp,
                                            velocity = velocities,
                                            id  = "my first spectrum")

moessbauer_spectrum.Calculate()

moessbauer_spectrum.Plot(velocity = True)