polarization example

import nexus as nx
import numpy as np

layer_Fe = nx.Layer(id = "Fe",
                    material = nx.Material.Template(nx.lib.material.Fe),
                    thickness = 3000,
                    roughness = 100
                    )

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

layer_Fe.material.hyperfine_sites = [site]

sample = nx.Sample(layers = layer_Fe)

beam  = nx.Beam()
beam.LinearSigma()

analyzer = nx.Analyzer()
analyzer.LinearSigma()

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



energy_spectrum = nx.EnergySpectrum(experiment = exp,
                                    detuning = np.linspace(-100, 100, 2001)
                                   )

energy_spectrum.Calculate()

energy_spectrum.Plot()



analyzer.LinearPi()

energy_spectrum.Calculate()

energy_spectrum.Plot()