Skip to content
Snippets Groups Projects
Commit 7f37e61c authored by Wuttke, Joachim's avatar Wuttke, Joachim
Browse files

update & fix Examples/fit/specular/Honeycomb_fit.py

parent 3fb622ad
No related branches found
No related tags found
1 merge request!1201fix several fit examples
Pipeline #83788 passed
......@@ -75,18 +75,16 @@ def get_sample(parameters, sign, ms150=1):
return sample
def get_simulation(q_axis, fitParams, sign, ms150=False):
def get_simulation(q_axis, parameters, sign, ms150=False):
n_points = 25
n_sig = 3.0
parameters = dict(fitParams, **fixedParams)
q_distr = ba.RangedDistributionGaussian(n_points, n_sig)
dq = parameters["dq"]*q_axis
scan = ba.QzScan(q_axis)
scan.setAbsoluteQResolution(q_distr, dq)
scan.setAbsoluteQVectorResolution(q_distr, dq)
if ms150:
......@@ -97,15 +95,18 @@ def get_simulation(q_axis, fitParams, sign, ms150=False):
sample = get_sample(parameters=parameters, sign=sign, ms150=1)
simulation = ba.SpecularSimulation(scan, sample)
simulation.beam().setIntensity(parameters["intensity"])
simulation.setBackground(ba.ConstantBackground(5e-7))
return simulation
def run_simulation(q_axis, fitParams, *, sign, ms150=False):
simulation = get_simulation(q_axis, fitParams, sign, ms150)
return simulation.simulate()
parameters = dict(fitParams, **fixedParams)
simulation = get_simulation(q_axis, parameters, sign, ms150)
result = simulation.simulate()
result.data_field().scale(parameters["intensity"])
return result
def qr(result):
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment