Changeset c94577f in sasmodels

Timestamp:

Mar 20, 2016 10:56:38 AM (9 years ago)

Author:

awashington

Branches:

master, core_shell_microgels, costrafo411, magnetic_model, release_v0.94, release_v0.95, ticket-1257-vesicle-product, ticket_1156, ticket_1265_superball, ticket_822_more_unit_tests

Children:

459e890

Parents:

a154ad16

Message:

Fix G0 calculation of Hankel transform

File:

• sasmodels/sesans.py (modified) (2 diffs)

sasmodels/sesans.py

@@ -189 +189 @@
         integral = besselj(0, q*SElength_i)*Iq*q
         G[i] = np.sum(integral)
+    G0 = np.sum(Iq*q)
 
     # [m^-1] step size in q, needed for integration
@@ -195 +196 @@
     # integration step, convert q into [m**-1] and 2 pi circle integration
     G *= dq*1e10*2*pi
-
-    P = exp(thickness*wavelength**2/(4*pi**2)*(G-G[0]))
+    G0 *= dq*1e10*2*pi
+
+    P = exp(thickness*wavelength**2/(4*pi**2)*(G-G0))
 
     return P