from refl1d.names import *

# Turn off resolution bars in plots.  Only do this after you have plotted the
# data with resolution bars so you know it looks reasonable, and you are not
# fitting the sample_broadening parameter in the probe.
Probe.show_resolution = False

nickel = Material("Ni")
sample = silicon(0, 10) | nickel(125, 10) | air

sample[1].thickness.pm(50)

sample[0].interface.range(3, 12)
sample[1].interface.range(0, 20)

instrument = SNS.Liquids()
files = ["nifilm-tof-%d.dat" % d for d in (1, 2, 3, 4)]
probe = ProbeSet(instrument.load(f) for f in files)

M = Experiment(probe=probe, sample=sample)

problem = FitProblem(M)