Changeset 23df833 in sasmodels for example

Timestamp:

Oct 30, 2018 6:19:38 PM (6 years ago)

Author:

Paul Kienzle <pkienzle@…>

Branches:

master, core_shell_microgels, magnetic_model, ticket-1257-vesicle-product, ticket_1156, ticket_1265_superball, ticket_822_more_unit_tests

Children:

5024a56

Parents:

aa8c6e0

Message:

Document direct call to Fq in user guide

File:

• example/cylinder_eval.py (modified) (1 diff)

example/cylinder_eval.py

@@ -3 +3 @@
 """
 
-from numpy import logspace
+from numpy import logspace, sqrt
 from matplotlib import pyplot as plt
 from sasmodels.core import load_model
-from sasmodels.direct_model import call_kernel
+from sasmodels.direct_model import call_kernel, call_Fq
 
 model = load_model('cylinder')
 q = logspace(-3, -1, 200)
 kernel = model.make_kernel([q])
-Iq = call_kernel(kernel, dict(radius=200.))
-plt.loglog(q, Iq)
+pars = {'radius': 200, 'radius_pd': 0.1, 'scale': 2}
+Iq = call_kernel(kernel, pars)
+F, Fsq, Reff, V, Vratio = call_Fq(kernel, pars)
+plt.loglog(q, Iq, label='2 I(q)')
+plt.loglog(q, F**2/V, label='<F(q)>^2/V')
+plt.loglog(q, Fsq/V, label='<F^2(q)>/V')
 plt.xlabel('q (1/A)')
-plt.ylabel('I(q)')
+plt.ylabel('I(q) (1/cm)')
 plt.title('Cylinder with radius 200.')
+plt.legend()
 plt.show()