from refl1d.names import * from refl1d.stajconvert import load_mlayer # Load neutron model and data from staj file M = load_mlayer("De2_VATR.staj") # Set thickness/roughness fitting parameters to +/- 20 % # Set SLD to +/- 5% for all but the incident medium and the substrate. for L in M.sample[1:-1]: L.thickness.pmp(20) L.interface.pmp(20) L.material.rho.pmp(5) # Let the substrate SLD vary by 2% M.sample[0].material.rho.pmp(2) M.sample[0].interface.range(0,20) M.sample[1].interface.range(0,20) problem = FitProblem(M) problem.name = "Desorption 2"
from refl1d.names import * from refl1d.stajconvert import load_mlayer # Load neutron model and data from staj file # Layer names are ordered from substrate to surface, and defaults to # the names in the original staj file. # Model name defaults to the data file name layers=["sappire","MgO","MgHx1","MgHx2","Pd","air"] M = load_mlayer("De2_VATR.staj",layers=layers,name="n6hd2") # Set thickness/roughness fitting parameters to +/- 20 % # Set SLD to +/- 5% for all but the incident medium and the substrate. for L in M.sample[1:-1]: L.thickness.pmp(20) L.interface.pmp(20) L.material.rho.pmp(5) # Let the substrate SLD vary by 2% M.sample[0].material.rho.pmp(2) M.sample[0].interface.range(0,20) M.sample[1].interface.range(0,20) problem = FitProblem(M) problem.name = "Desorption 2"
from refl1d.names import * from refl1d.stajconvert import load_mlayer # Load neutron model and data from staj file # Layer names are ordered from substrate to surface, and defaults to # the names in the original staj file. # Model name defaults to the data file name layers = ["sappire", "MgO", "MgHx1", "MgHx2", "Pd", "air"] M = load_mlayer("De2_VATR.staj", layers=layers, name="n6hd2") # Set thickness/roughness fitting parameters to +/- 20 % # Set SLD to +/- 5% for all but the incident medium and the substrate. for L in M.sample[1:-1]: L.thickness.pmp(20) L.interface.pmp(20) L.material.rho.pmp(5) # Let the substrate SLD vary by 2% M.sample[0].material.rho.pmp(2) M.sample[0].interface.range(0, 20) M.sample[1].interface.range(0, 20) problem = FitProblem(M) problem.name = "Desorption 2"