Changes in example/sesansfit.py [9c117a2:0ac3db5] in sasmodels


Ignore:
File:
1 edited

Legend:

Unmodified
Added
Removed
  • example/sesansfit.py

    r9c117a2 r0ac3db5  
    44from sasmodels import bumps_model as sas 
    55kernel = sas.load_model('sphere', dtype='single') 
    6 #kernel = sas.load_model('triaxial_ellipsoid', dtype='single') 
    76 
    87 
     
    1211    loader=Loader() 
    1312    data=loader.load('testsasview1.ses') 
    14     data.x /=10 
    1513     
    1614#    data = load_sesans('mydatfile.pz') 
     
    2220    err_data = np.ones_like(SElength)*0.03 
    2321 
    24     class Sample: 
    25         zacceptance = 0.1 # [A^-1] 
    26         thickness = 0.2 # [cm] 
    27          
    2822    class SESANSData1D: 
    2923        #q_zmax = 0.23 # [A^-1] 
    30         lam = 0.2 # [nm] 
     24        zacceptance = 0.1 # [A^-1] 
     25        lam = 2e-10 # [m] 
     26        thickness = 0.2 # [cm] 
    3127        x = SElength 
    3228        y = data 
    3329        dy = err_data 
    34         sample = Sample() 
    35     data = SESANSData1D() 
     30    data = SesansData() 
     31print dir(data) 
    3632 
    3733radius = 1000 
    3834data.Rmax = 3*radius # [A] 
    3935 
    40 ##  Sphere parameters 
    41  
    4236phi = Parameter(0.1, name="phi") 
    4337model = sas.BumpsModel(data, kernel, 
    4438    scale=phi*(1-phi), sld=7.0, solvent_sld=1.0, radius=radius) 
    45 phi.range(0.001,0.90) 
    46 #model.radius.pmp(40) 
    47 model.radius.range(100,10000) 
    48 #model.sld.pmp(5) 
    49 #model.background 
    50 #model.radius_pd=0 
    51 #model.radius_pd_n=0 
     39phi.pmp(10) 
     40model.radius.pmp(40) 
     41model.sld.pm(2) 
     42model.background.range(0,5) 
     43 
    5244 
    5345if False: # have sans data 
     
    5749    problem = FitProblem(model) 
    5850 
    59  
    60 ### Tri-Axial Ellipsoid 
    61 # 
    62 #phi = Parameter(0.1, name='phi') 
    63 #model = sas.BumpsModel(data, kernel, 
    64 #    scale=phi*(1-phi), sld=7.0, solvent_sld=1.0, radius=radius) 
    65 #phi.range(0.001,0.90) 
    66 ##model.radius.pmp(40) 
    67 #model.radius.range(100,10000) 
    68 ##model.sld.pmp(5) 
    69 ##model.background 
    70 ##model.radius_pd=0 
    71 ##model.radius_pd_n=0 
    72 # 
    73 #if False: # have sans data 
    74 #    sansmodel = sas.BumpsModel(sans_data, kernel, **model.parameters()) 
    75 #    problem = FitProblem([model, sansmodel]) 
    76 #else: 
    77 #    problem = FitProblem(model) 
Note: See TracChangeset for help on using the changeset viewer.