Matrices

[1]:

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)

scattering_matrices = nx.Matrices(experiment = exp,
                                  energy = 14400)

print(scattering_matrices())

[[[0.71751273-0.58897687j 0.        +0.j        ]
  [0.        +0.j         0.71751273-0.58897687j]]

 [[0.66100408+0.40190631j 0.        +0.j        ]
  [0.        +0.j         0.66100408+0.40190631j]]]