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

no need for fct qr

parent ea573f32
No related branches found
No related tags found
1 merge request!1953cleanup fit examples
...@@ -60,16 +60,6 @@ def get_simulation(q_axis, fitParams): ...@@ -60,16 +60,6 @@ def get_simulation(q_axis, fitParams):
return simulation return simulation
def qr(result):
"""
Return q and reflectivity arrays from simulation result.
"""
q = np.array(result.axis(0).binCenters())
r = np.array(result.npArray())
return q, r
#################################################################### ####################################################################
# Experimental data # Experimental data
#################################################################### ####################################################################
...@@ -162,8 +152,8 @@ if __name__ == '__main__': ...@@ -162,8 +152,8 @@ if __name__ == '__main__':
# Plot data with initial model # Plot data with initial model
q, r = qr(get_simulation(qzs, PInitial).simulate()) r = get_simulation(qzs, PInitial).simulate().npArray()
plot(q, r, data, "PtLayerFit_initial.pdf", PInitial | fixedP) plot(qzs, r, data, "PtLayerFit_initial.pdf", PInitial | fixedP)
# Fit # Fit
...@@ -189,7 +179,7 @@ if __name__ == '__main__': ...@@ -189,7 +179,7 @@ if __name__ == '__main__':
print("Fit Result:") print("Fit Result:")
print(fitResult) print(fitResult)
q, r = qr(get_simulation(qzs, fitParams=fitResult).simulate()) r = get_simulation(qzs, fitParams=fitResult).simulate().npArray()
plot(q, r, data, "PtLayerFit_fit.pdf", fitResult | fixedP) plot(qzs, r, data, "PtLayerFit_fit.pdf", fitResult | fixedP)
plt.show() plt.show()
...@@ -60,16 +60,6 @@ def get_simulation(q_axis, fitParams): ...@@ -60,16 +60,6 @@ def get_simulation(q_axis, fitParams):
return simulation return simulation
def qr(result):
"""
Return q and reflectivity arrays from simulation result.
"""
q = np.array(result.axis(0).binCenters())
r = np.array(result.npArray())
return q, r
#################################################################### ####################################################################
# Experimental data # Experimental data
#################################################################### ####################################################################
...@@ -162,8 +152,8 @@ if __name__ == '__main__': ...@@ -162,8 +152,8 @@ if __name__ == '__main__':
# Plot data with initial model # Plot data with initial model
q, r = qr(get_simulation(qzs, PInitial).simulate()) r = get_simulation(qzs, PInitial).simulate().npArray()
plot(q, r, data, "PtLayerFit_initial.pdf", PInitial | fixedP) plot(qzs, r, data, "PtLayerFit_initial.pdf", PInitial | fixedP)
# Fit # Fit
...@@ -189,7 +179,7 @@ if __name__ == '__main__': ...@@ -189,7 +179,7 @@ if __name__ == '__main__':
print("Fit Result:") print("Fit Result:")
print(fitResult) print(fitResult)
q, r = qr(get_simulation(qzs, fitParams=fitResult).simulate()) r = get_simulation(qzs, fitParams=fitResult).simulate().npArray()
plot(q, r, data, "PtLayerFit_fit.pdf", fitResult | fixedP) plot(qzs, r, data, "PtLayerFit_fit.pdf", fitResult | fixedP)
plt.show() plt.show()
...@@ -60,16 +60,6 @@ def get_simulation(q_axis, fitParams): ...@@ -60,16 +60,6 @@ def get_simulation(q_axis, fitParams):
return simulation return simulation
def qr(result):
"""
Return q and reflectivity arrays from simulation result.
"""
q = np.array(result.axis(0).binCenters())
r = np.array(result.npArray())
return q, r
#################################################################### ####################################################################
# Experimental data # Experimental data
#################################################################### ####################################################################
...@@ -162,8 +152,8 @@ if __name__ == '__main__': ...@@ -162,8 +152,8 @@ if __name__ == '__main__':
# Plot data with initial model # Plot data with initial model
q, r = qr(get_simulation(qzs, PInitial).simulate()) r = get_simulation(qzs, PInitial).simulate().npArray()
plot(q, r, data, "PtLayerFit_initial.pdf", PInitial | fixedP) plot(qzs, r, data, "PtLayerFit_initial.pdf", PInitial | fixedP)
# Fit # Fit
...@@ -189,7 +179,7 @@ if __name__ == '__main__': ...@@ -189,7 +179,7 @@ if __name__ == '__main__':
print("Fit Result:") print("Fit Result:")
print(fitResult) print(fitResult)
q, r = qr(get_simulation(qzs, fitParams=fitResult).simulate()) r = get_simulation(qzs, fitParams=fitResult).simulate().npArray()
plot(q, r, data, "PtLayerFit_fit.pdf", fitResult | fixedP) plot(qzs, r, data, "PtLayerFit_fit.pdf", fitResult | fixedP)
plt.show() plt.show()
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment