Changeset 8cdb9f1 in sasmodels


Ignore:
Timestamp:
Aug 4, 2014 5:25:51 PM (11 years ago)
Author:
HMP1 <helen.park@…>
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:
1726b21
Parents:
a42fec0
Message:

organizing fit

File:
1 edited

Legend:

Unmodified
Added
Removed

  • fit.py

    ra42fec0 r8cdb9f1  
    2020#set_half(data, 'right') 
    2121 
    22 """ 
    2322 
    24 model = SasModel(data, OneDGpuCylinder, 
    25 scale=0.0013, 
    26 radius=105, 
    27 length=1000, 
    28 background=21, 
    29 sldCyl=.291e-6,sldSolv=7.105e-6, 
    30 radius_pd=0.1,radius_pd_n=10,radius_pd_nsigma=0, 
    31 length_pd=0.1,length_pd_n=5,length_pd_nsigma=0, 
    32 bolim=0.0, 
    33 uplim=90) #bottom limit, upper limit of angle integral 
    34 """ 
    35 """ 
    36 model = SasModel(data, GpuEllipse, 
    37 scale=0.08, 
    38 radius_a=15, radius_b=800, 
    39 sldEll=.291e-6, sldSolv=7.105e-6, 
    40 background=0, 
    41 axis_theta=90, axis_phi=0, 
    42 axis_theta_pd=15, axis_theta_pd_n=40, axis_theta_pd_nsigma=3, 
    43  
    44 radius_a_pd=0.222296, radius_a_pd_n=1, radius_a_pd_nsigma=0, 
    45 radius_b_pd=.000128, radius_b_pd_n=1, radius_b_pd_nsigma=0, 
    46 axis_phi_pd=2.63698e-05, axis_phi_pd_n=20, axis_phi_pd_nsigma=0, 
    47 dtype='float') 
     23if 0: 
     24    model = SasModel(data, OneDGpuCylinder, 
     25    scale=0.0013, radius=105, length=1000, 
     26    background=21, sldCyl=.291e-6,sldSolv=7.105e-6, 
     27    radius_pd=0.1,radius_pd_n=10,radius_pd_nsigma=0, 
     28    length_pd=0.1,length_pd_n=5,length_pd_nsigma=0, 
     29    bolim=0.0, uplim=90) #bottom limit, upper limit of angle integral 
    4830 
    4931 
    50 # SET THE FITTING PARAMETERS 
    51 model.radius_a.range(15, 1000) 
    52 model.radius_b.range(15, 1000) 
    53 #model.axis_theta_pd.range(0, 360) 
    54 #model.background.range(0,1000) 
    55 #model.scale.range(0, 1) 
    56 """ 
     32if 0: 
     33    model = SasModel(data, GpuEllipse, 
     34    scale=0.08, 
     35    radius_a=15, radius_b=800, 
     36    sldEll=.291e-6, sldSolv=7.105e-6, 
     37    background=0, 
     38    axis_theta=90, axis_phi=0, 
     39    axis_theta_pd=15, axis_theta_pd_n=40, axis_theta_pd_nsigma=3, 
     40    radius_a_pd=0.222296, radius_a_pd_n=1, radius_a_pd_nsigma=0, 
     41    radius_b_pd=.000128, radius_b_pd_n=1, radius_b_pd_nsigma=0, 
     42    axis_phi_pd=2.63698e-05, axis_phi_pd_n=20, axis_phi_pd_nsigma=0, 
     43    dtype='float') 
    5744 
    58 """ 
    59 model = SasModel(data, GpuLamellar, 
    60 scale=0.08, 
    61 bi_thick=19.2946, 
    62 sld_bi=5.38e-6,sld_sol=7.105e-6, 
    63 background=0.003, 
    64 bi_thick_pd= 0.37765, bi_thick_pd_n=10, bi_thick_pd_nsigma=3, 
    65 dtype='float') 
    6645 
    67 # SET THE FITTING PARAMETERS 
    68 #model.bi_thick.range(0, 1000) 
    69 #model.scale.range(0, 1) 
    70 #model.bi_thick_pd.range(0, 1000) 
    71 #model.background.range(0, 1000) 
    72 model.sld_bi.range(0, 1) 
     46    # SET THE FITTING PARAMETERS 
     47    model.radius_a.range(15, 1000) 
     48    model.radius_b.range(15, 1000) 
     49    #model.axis_theta_pd.range(0, 360) 
     50    #model.background.range(0,1000) 
     51    #model.scale.range(0, 1) 
    7352 
    74 """ 
     53 
     54if 0: 
     55    model = SasModel(data, GpuLamellar, 
     56    scale=0.08, 
     57    bi_thick=19.2946, 
     58    sld_bi=5.38e-6,sld_sol=7.105e-6, 
     59    background=0.003, 
     60    bi_thick_pd= 0.37765, bi_thick_pd_n=10, bi_thick_pd_nsigma=3, 
     61    dtype='float') 
     62 
     63    # SET THE FITTING PARAMETERS 
     64    #model.bi_thick.range(0, 1000) 
     65    #model.scale.range(0, 1) 
     66    #model.bi_thick_pd.range(0, 1000) 
     67    #model.background.range(0, 1000) 
     68    model.sld_bi.range(0, 1) 
     69 
     70 
    7571if 1: 
    7672    model = SasModel(data, GpuCylinder, 
    77     scale=0.0104, 
    78     radius=92.5, 
    79     length=798.3, 
    80     sldCyl=.29e-6, 
    81     sldSolv=7.105e-6, 
    82     background=5, 
    83     cyl_theta=0, 
    84     cyl_phi=0, 
    85     cyl_theta_pd=22.11, 
    86     cyl_theta_pd_n=20, 
    87     cyl_theta_pd_nsigma=3, 
    88     radius_pd=.0084, 
    89     radius_pd_n=10, 
    90     radius_pd_nsigma=3, 
    91     length_pd=0.493, 
    92     length_pd_n=10, 
    93     length_pd_nsigma=3, 
    94     cyl_phi_pd=0, 
    95     cyl_phi_pd_n=1, 
    96     cyl_phi_pd_nsigma=3, 
     73    scale=0.0104, radius=92.5, length=798.3, 
     74    sldCyl=.29e-6, sldSolv=7.105e-6, background=5, 
     75    cyl_theta=0, cyl_phi=0, 
     76    cyl_theta_pd=22.11, cyl_theta_pd_n=20, cyl_theta_pd_nsigma=3, 
     77    radius_pd=.0084, radius_pd_n=10, radius_pd_nsigma=3, 
     78    length_pd=0.493, length_pd_n=10, length_pd_nsigma=3, 
     79    cyl_phi_pd=0, cyl_phi_pd_n=1, cyl_phi_pd_nsigma=3, 
    9780    dtype='float') 
    9881 
    … …  
    11497if 0: 
    11598    model = SasModel(data, GpuCoreShellCylinder, 
    116                      scale= .00031, radius=19.5, thickness=30, length=22, 
    117                      core_sld=7.105e-6, shell_sld=.291e-6, solvent_sld=7.105e-6, 
    118                      background=0.2, axis_theta=0, axis_phi=0, 
     99    scale= .00031, radius=19.5, thickness=30, length=22, 
     100    core_sld=7.105e-6, shell_sld=.291e-6, solvent_sld=7.105e-6, 
     101    background=0.2, axis_theta=0, axis_phi=0, 
    119102 
    120                      radius_pd=0.26, radius_pd_n=10, radius_pd_nsigma=3, 
    121                      length_pd=0.26, length_pd_n=10, length_pd_nsigma=3, 
    122                      thickness_pd=1, thickness_pd_n=1, thickness_pd_nsigma=1, 
    123                      axis_theta_pd=1, axis_theta_pd_n=10, axis_theta_pd_nsigma=3, 
    124                      axis_phi_pd=0.1, axis_phi_pd_n=1, axis_phi_pd_nsigma=1, 
    125                      dtype='float') 
     103    radius_pd=0.26, radius_pd_n=10, radius_pd_nsigma=3, 
     104    length_pd=0.26, length_pd_n=10, length_pd_nsigma=3, 
     105    thickness_pd=1, thickness_pd_n=1, thickness_pd_nsigma=1, 
     106    axis_theta_pd=1, axis_theta_pd_n=10, axis_theta_pd_nsigma=3, 
     107    axis_phi_pd=0.1, axis_phi_pd_n=1, axis_phi_pd_nsigma=1, 
     108    dtype='float') 
    126109 
    127110    # SET THE FITTING PARAMETERS 
    … …  
    139122 
    140123 
    141 """ 
    142 model = SasModel(data, GpuCapCylinder, 
    143                  scale=1, rad_cyl=20, rad_cap=40, length=400, 
    144                  sld_capcyl=1e-6, sld_solv=6.3e-6, 
    145                  background=0, theta=0, phi=0, 
    146                  rad_cyl_pd=.1, rad_cyl_pd_n=1, rad_cyl_pd_nsigma=0, 
    147                  rad_cap_pd=.1, rad_cap_pd_n=1, rad_cap_pd_nsigma=0, 
    148                  length_pd=.1, length_pd_n=1, length_pd_nsigma=0, 
    149                  theta_pd=.1, theta_pd_n=1, theta_pd_nsigma=0, 
    150                  phi_pd=.1, phi_pd_n=1, phi_pd_nsigma=0, 
    151                  dtype='float') 
     124if 0: 
     125    model = SasModel(data, GpuCapCylinder, 
     126    scale=1, rad_cyl=20, rad_cap=40, length=400, 
     127    sld_capcyl=1e-6, sld_solv=6.3e-6, 
     128    background=0, theta=0, phi=0, 
     129    rad_cyl_pd=.1, rad_cyl_pd_n=1, rad_cyl_pd_nsigma=0, 
     130    rad_cap_pd=.1, rad_cap_pd_n=1, rad_cap_pd_nsigma=0, 
     131    length_pd=.1, length_pd_n=1, length_pd_nsigma=0, 
     132    theta_pd=.1, theta_pd_n=1, theta_pd_nsigma=0, 
     133    phi_pd=.1, phi_pd_n=1, phi_pd_nsigma=0, 
     134    dtype='float') 
    152135 
    153 model.scale.range(0, 1) 
     136    model.scale.range(0, 1) 
    154137 
    155 """ 
    156 """ 
    157 model = SasModel(data, GpuTriEllipse, 
    158                  scale=0.08, axisA=15, axisB=20, axisC=500, 
    159                  sldEll=7.105e-6, sldSolv=.291e-6, 
    160                  background=5, theta=0, phi=0, psi=0, 
    161                  theta_pd=20, theta_pd_n=40, theta_pd_nsigma=3, 
    162                  phi_pd=.1, phi_pd_n=1, phi_pd_nsigma=0, 
    163                  psi_pd=30, psi_pd_n=1, psi_pd_nsigma=0, 
    164                  axisA_pd=.1, axisA_pd_n=1, axisA_pd_nsigma=0, 
    165                  axisB_pd=.1, axisB_pd_n=1, axisB_pd_nsigma=0, 
    166                  axisC_pd=.1, axisC_pd_n=1, axisC_pd_nsigma=0, dtype='float') 
    167138 
    168 # SET THE FITTING PARAMETERS 
    169 model.axisA.range(15, 1000) 
    170 model.axisB.range(15, 1000) 
    171 #model.axisC.range(15, 1000) 
    172 #model.background.range(0,1000) 
    173 #model.theta_pd.range(0, 360) 
    174 #model.phi_pd.range(0, 360) 
    175 #model.psi_pd.range(0, 360) 
     139if 0: 
     140    model = SasModel(data, GpuTriEllipse, 
     141    scale=0.08, axisA=15, axisB=20, axisC=500, 
     142    sldEll=7.105e-6, sldSolv=.291e-6, 
     143    background=5, theta=0, phi=0, psi=0, 
     144    theta_pd=20, theta_pd_n=40, theta_pd_nsigma=3, 
     145    phi_pd=.1, phi_pd_n=1, phi_pd_nsigma=0, 
     146    psi_pd=30, psi_pd_n=1, psi_pd_nsigma=0, 
     147    axisA_pd=.1, axisA_pd_n=1, axisA_pd_nsigma=0, 
     148    axisB_pd=.1, axisB_pd_n=1, axisB_pd_nsigma=0, 
     149    axisC_pd=.1, axisC_pd_n=1, axisC_pd_nsigma=0, dtype='float') 
    176150 
    177 """ 
     151    # SET THE FITTING PARAMETERS 
     152    model.axisA.range(15, 1000) 
     153    model.axisB.range(15, 1000) 
     154    #model.axisC.range(15, 1000) 
     155    #model.background.range(0,1000) 
     156    #model.theta_pd.range(0, 360) 
     157    #model.phi_pd.range(0, 360) 
     158    #model.psi_pd.range(0, 360) 
     159 
     160 
    178161 
    179162 
Note: See TracChangeset for help on using the changeset viewer.