from bumps.names import *

from sasmodels.data import load_data
from sasmodels.core import load_model
from sasmodels.bumps_model import Model, Experiment

# Spherical particle data, not ellipsoids
data = load_data('../../sasview/sasview/test/sesans_data/sphere2micron.ses')
data.oriented = True

kernel = load_model("ellipsoid")

model = Model(
    kernel,
    scale=0.08, background=0,
    sld=.291, sld_solvent=7.105,
    radius_polar=10000, radius_polar_pd=0.222296, radius_polar_pd_n=0,
    radius_equatorial=2600, radius_equatorial_pd=0.28, radius_equatorial_pd_n=0,
    theta=60, theta_pd=0, theta_pd_n=0,
    phi=60, phi_pd=0, phi_pd_n=0,
    )

# SET THE FITTING PARAMETERS
model.radius_polar.range(1000, 100000)
model.radius_equatorial.range(1000, 100000)
model.theta.range(0, 90)
model.phi.range(0, 180)
model.background.range(0,1000)
model.scale.range(0, 10)

M = Experiment(data=data, model=model)
problem = FitProblem(M)