| 1 | from bumps.names import * |
|---|
| 2 | |
|---|
| 3 | from sasmodels.data import load_data |
|---|
| 4 | from sasmodels.core import load_model |
|---|
| 5 | from sasmodels.bumps_model import Model, Experiment |
|---|
| 6 | |
|---|
| 7 | # Spherical particle data, not ellipsoids |
|---|
| 8 | data = load_data('../../sasview/sasview/test/sesans_data/sphere2micron.ses') |
|---|
| 9 | data.oriented = True |
|---|
| 10 | |
|---|
| 11 | kernel = load_model("ellipsoid") |
|---|
| 12 | |
|---|
| 13 | model = Model( |
|---|
| 14 | kernel, |
|---|
| 15 | scale=0.08, background=0, |
|---|
| 16 | sld=.291, sld_solvent=7.105, |
|---|
| 17 | radius_polar=10000, radius_polar_pd=0.222296, radius_polar_pd_n=0, |
|---|
| 18 | radius_equatorial=2600, radius_equatorial_pd=0.28, radius_equatorial_pd_n=0, |
|---|
| 19 | theta=60, theta_pd=0, theta_pd_n=0, |
|---|
| 20 | phi=60, phi_pd=0, phi_pd_n=0, |
|---|
| 21 | ) |
|---|
| 22 | |
|---|
| 23 | # SET THE FITTING PARAMETERS |
|---|
| 24 | model.radius_polar.range(1000, 100000) |
|---|
| 25 | model.radius_equatorial.range(1000, 100000) |
|---|
| 26 | model.theta.range(0, 90) |
|---|
| 27 | model.phi.range(0, 180) |
|---|
| 28 | model.background.range(0,1000) |
|---|
| 29 | model.scale.range(0, 10) |
|---|
| 30 | |
|---|
| 31 | M = Experiment(data=data, model=model) |
|---|
| 32 | problem = FitProblem(M) |
|---|
