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

manually normalize examples, repair persistence tests

parent 875fcc69
No related branches found
No related tags found
No related merge requests found
......@@ -34,13 +34,14 @@ def get_sample():
return multi_layer
def get_simulation(scan_size=500):
def get_simulation(sample, scan_size=500):
"""
Defines and returns a specular simulation.
"""
simulation = ba.SpecularSimulation()
scan = ba.AngularSpecScan(1.54*angstrom, scan_size, 0.0*deg, 2.0*deg)
simulation.setScan(scan)
simulation.setSample(sample)
return simulation
......@@ -48,13 +49,10 @@ def run_simulation():
"""
Runs simulation and returns its result.
"""
sample = get_sample()
simulation = get_simulation()
simulation.setSample(sample)
simulation = get_simulation(get_sample())
simulation.runSimulation()
return simulation.result()
if __name__ == '__main__':
results = run_simulation()
ba.plot_simulation_result(results, c='k')
ba.run_and_plot(get_simulation(get_sample()), c='k')
......@@ -7,6 +7,7 @@ with BornAgain.
import numpy as np
import bornagain as ba
from os import path
from matplotlib import pyplot as plt
# input parameters
wavelength = 1.54*ba.angstrom
......@@ -68,7 +69,7 @@ def create_real_data():
return ax_values, real_data
def get_simulation(scan_size=500):
def get_simulation(sample, scan_size=500):
"""
Returns a specular simulation with beam and detector defined.
"""
......@@ -81,39 +82,14 @@ def get_simulation(scan_size=500):
simulation = ba.SpecularSimulation()
simulation.setScan(scan)
return simulation
def run_simulation():
"""
Runs simulation and returns it.
"""
sample = get_sample()
simulation = get_simulation()
simulation.setSample(sample)
simulation.runSimulation()
return simulation.result()
def plot(results):
"""
:param results:
:return:
"""
from matplotlib import pyplot as plt
return simulation
ba.plot_simulation_result(results, noshow=True)
if __name__ == '__main__':
ba.run_and_plot(get_simulation(get_sample()), noshow=True)
genx_axis, genx_values = create_real_data()
plt.semilogy(genx_axis, genx_values, 'ko', markevery=300)
plt.legend(['BornAgain', 'GenX'], loc='upper right')
plt.show()
if __name__ == '__main__':
results = run_simulation()
plot(results)
......@@ -50,7 +50,7 @@ def get_sample():
return multi_layer
def get_simulation(scan_size=500):
def get_simulation(sample, scan_size=500):
"""
Returns a specular simulation with beam and detector defined.
"""
......@@ -64,21 +64,10 @@ def get_simulation(scan_size=500):
simulation = ba.SpecularSimulation()
simulation.setScan(scan)
simulation.setSample(sample)
return simulation
def run_simulation():
"""
Runs simulation and returns it.
"""
sample = get_sample()
simulation = get_simulation()
simulation.setSample(sample)
simulation.runSimulation()
return simulation.result()
if __name__ == '__main__':
results = run_simulation()
ba.plot_simulation_result(results)
ba.run_and_plot(get_simulation(get_sample()))
......@@ -42,7 +42,7 @@ def get_sample():
return multi_layer
def get_simulation(scan_size=500):
def get_simulation(sample, scan_size=500):
"""
Defines and returns specular simulation
with a qz-defined beam
......@@ -59,21 +59,10 @@ def get_simulation(scan_size=500):
simulation = ba.SpecularSimulation()
simulation.setScan(scan)
simulation.setSample(sample)
return simulation
def run_simulation():
"""
Runs simulation and returns its result.
"""
sample = get_sample()
simulation = get_simulation()
simulation.setSample(sample)
simulation.runSimulation()
return simulation.result()
if __name__ == '__main__':
result = run_simulation()
ba.plot_simulation_result(result)
ba.run_and_plot(get_simulation(get_sample()))
......@@ -40,7 +40,7 @@ def get_sample():
return multi_layer
def get_simulation(scan_size=500):
def get_simulation(sample, scan_size=500):
"""
Defines and returns specular simulation
with a qz-defined beam
......@@ -49,20 +49,9 @@ def get_simulation(scan_size=500):
scan = ba.QSpecScan(qzs)
simulation = ba.SpecularSimulation()
simulation.setScan(scan)
return simulation
def run_simulation():
"""
Runs simulation and returns its result.
"""
sample = get_sample()
simulation = get_simulation()
simulation.setSample(sample)
simulation.runSimulation()
return simulation.result()
return simulation
if __name__ == '__main__':
result = run_simulation()
ba.plot_simulation_result(result)
ba.run_and_plot(get_simulation(get_sample()))
......@@ -116,7 +116,8 @@ def plot(results):
plt.subplot(2, 2, nplot + 1)
ba.plot_colormap(hist, zlabel="")
plt.tight_layout()
plt.show()
if not "NOSHOW" in os.environ:
plt.show()
if __name__ == '__main__':
......
......@@ -62,4 +62,6 @@ test_example(varia/CustomFormFactor 2e-10)
test_example(varia/DepthProbe 2e-10)
test_example(varia/Interference1DLattice 1.5)
test_example(varia/OffSpecularSimulation 2e-10)
test_example(varia/SimulationParameters 2e-10)
# TODO replace after revision of parameter machinery
# test_example(varia/SimulationParameters 2e-10)
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