/******************************************************************************* * * Component: SasView_model * * %I * Written by: Jakob Garde, Torben Nielsen, Peter Willendrup * Date: 03.02.2016 * Origin: SasView, DTU, European Spallation Source ERIC * * This SANS sample exposes SasView's scattering kernels to McStas. In this way SasView's monodisperse scattering kernels can be call from McStas. * * %D * Sample for use in SANS instruments. The models describe mono disperse particles in thin solution. The sample geometry may have the shape: * * Shape: - A filled box with dimensions xwidth, yheight and zdepth. * - A cylinder with dimensions radius and yheight. * - A filled sphere given by radius. * * These parameters are mutually exclusive. * * Example using spheres in thin solution, with radius=200 AA and a delta_sld=0.6 fm/AA^3: * SasView_model(model_index=47, model_scale=1.0, model_pars={1, 7, 200}, model_abs=0.0, xwidth=0.01, yheight=0.01, zdepth=0.005,) * * * The algorithm of this component requires use of the ISO C standard c99, standard from McStas 2.7 and 3.0. * * The list of scattering models in SasView is called sasmodels. The list of McStas available SasView sasmodels are found in the table below. * * A few models may require manual documentation lookup using the above link to the SasView site. * * McStas does currently not support multiplication of formfactor models with structure factor models. * * The 2D scattering scattering kernels are denoted by modelname_xy. I.e. to evalulate scattering from aligned cylinders use model_index=10 to use cylinder_xy. * * MDOC * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
Model no.SasView nameParameters
0NoneNone
1barbell(sld, solvent_sld, bell_radius, radius, length)
2barbell_xy(sld, solvent_sld, bell_radius, radius, length, theta, phi)
3bcc_paracrystal(dnn, d_factor, radius, sld, solvent_sld)
4bcc_paracrystal_xy(dnn, d_factor, radius, sld, solvent_sld, theta, phi, psi)
5capped_cylinder(sld, solvent_sld, radius, cap_radius, length)
6capped_cylinder_xy(sld, solvent_sld, radius, cap_radius, length, theta, phi)
7core_shell_cylinder(core_sld, shell_sld, solvent_sld, radius, thickness, length)
8core_shell_cylinder_xy(core_sld, shell_sld, solvent_sld, radius, thickness, length, theta, phi)
9cylinder(sld, solvent_sld, radius, length)
10cylinder_xy(sld, solvent_sld, radius, length, theta, phi)
11dab(length)
12dab_xy(length)
13ellipsoid(sld, solvent_sld, rpolar, requatorial)
14ellipsoid_xy(sld, solvent_sld, rpolar, requatorial, theta, phi)
15fcc_paracrystal(dnn, d_factor, radius, sld, solvent_sld)
16fcc_paracrystal_xy(dnn, d_factor, radius, sld, solvent_sld, theta, phi, psi)
17flexible_cylinder_ex(length, kuhn_length, radius, axis_ratio, sld, solvent_sld)
18flexible_cylinder_ex_xy(length, kuhn_length, radius, axis_ratio, sld, solvent_sld)
19gaussian_peak(q0, sigma)
20gaussian_peak_xy(q0, sigma)
21guinier(rg)
22guinier_xy(rg)
23hardsphere(effect_radius, volfraction)
24hardsphere_xy(effect_radius, volfraction)
25HayterMSAsq(effect_radius, zz, VolFrac, Temp, csalt, dialec)
26HayterMSAsq_xy(effect_radius, charge, volfraction, temperature, saltconc, dielectconst)
27hollow_cylinder(radius, core_radius, length, sld, solvent_sld)
28hollow_cylinder_xy(radius, core_radius, length, sld, solvent_sld, theta, phi)
29lamellar(sld, solvent_sld, thickness)
30lamellar_xy(sld, solvent_sld, thickness)
31lamellar_FFHG(tail_length, head_length, sld, head_sld, solvent_sld)
32lamellar_FFHG_xy(tail_length, head_length, sld, head_sld, solvent_sld)
33lamellarCailleHG(tail_length, head_length, Nlayers, dd, Cp, tail_sld, head_sld, solvent_sld)
34lamellarCailleHG_xy(tail_length, head_length, Nlayers, spacing, Caille_parameter, sld, head_sld, solvent_sld)
35lamellarPC(th, Nlayers, davg, pd, sld, solvent_sld)
36lamellarPC_xy(thickness, Nlayers, spacing, spacing_polydisp, sld, solvent_sld)
37lamellarPS(del, Nlayers, dd, Cp, sld, solvent_sld)
38lamellarPS_xy(thickness, Nlayers, spacing, Caille_parameter, sld, solvent_sld)
39linear_pearls(radius, edge_sep, num_pearls, pearl_sld, solvent_sld)
40linear_pearls_xy(radius, edge_sep, num_pearls, pearl_sld, solvent_sld)
41lorentz(cor_length)
42lorentz_xy(cor_length)
43mass_fractal(radius, mass_dim, cutoff_length)
44mass_fractal_xy(radius, mass_dim, cutoff_length)
45mass_surface_fractal(mass_dim, surface_dim, cluster_rg, primary_rg)
46mass_surface_fractal_xy(mass_dim, surface_dim, cluster_rg, primary_rg)
47parallelepiped(sld, solvent_sld, a_side, b_side, c_side)
48parallelepiped_xy(sld, solvent_sld, a_side, b_side, c_side, theta, phi, psi)
49pearl_necklace(radius, edge_separation, string_thickness, number_of_pearls, sld, string_sld, solvent_sld)
50pearl_necklace_xy(radius, edge_separation, string_thickness, number_of_pearls, sld, string_sld, solvent_sld)
51sphere(sld, solvent_sld, radius)
52sphere_xy(sld, solvent_sld, radius)
53star_polymer(radius2, arms)
54star_polymer_xy(radius2, arms)
55stickyhardsphere(effect_radius, volfraction, perturb, stickiness)
56stickyhardsphere_xy(effect_radius, volfraction, perturb, stickiness)
57triaxial_ellipsoid(sld, solvent_sld, req_minor, req_major, rpolar)
58triaxial_ellipsoid_xy(sld, solvent_sld, req_minor, req_major, rpolar, theta, phi, psi)
* * * *

Available SasView sasmodels - Extended parameters description

* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
sasmodels nameParameters UnitsRange
barbell { 4, 1, 40, 20, 400 }  
barbell_xy { 4, 1, 40, 20, 400, 60, 60 }  
    
 Parameters:  
 Barbell scattering length density4e-6/Ang^2[-inf, inf]
 Solvent scattering length density1e-6/Ang^2[-inf, inf]
 Spherical bell radiusAng[0, inf]
 Cylindrical bar radiusAng[0, inf]
 Cylinder bar lengthAng[0, inf]
 In plane angledegrees[-inf, inf]
 Out of plane angledegrees[-inf, inf]
    
bcc { 220, 0.06, 40, 4, 1 }  
bcc_xy { 220, 0.06, 40, 4, 1, 60, 60, 60 }  
    
 Parameters:  
 Nearest neighbour distanceAng[-inf, inf]
 Paracrystal distortion factor [-inf, inf]
 Particle radiusAng[0, inf]
 Particle scattering length density1e-6/Ang^2[-inf, inf]
 Solvent scattering length density1e-6/Ang^2[-inf, inf]
 In plane angledegrees[-inf, inf]
 Out of plane angledegrees[-inf, inf]
 Out of plane angledegrees[-inf, inf]
    
capped_cylinder { 4, 1, 20, 20, 400 }  
capped_cylinder_xy { 4, 1, 20, 20, 400, 60, 60 }  
    
 Parameters:  
 Cylinder scattering length density1e-6/Ang^2[-inf, inf]
 Solvent scattering length density1e-6/Ang^2[-inf, inf]
 Cylinder radiusAng[0, inf]
 Cap radiusAng[0, inf]
 Cylinder lengthAng[0, inf]
 In plane angledegrees[-inf, inf]
 Out of plane angledegrees[-inf, inf]
    
core_shell_cylinder { 4, 4, 1, 20, 20, 400 }  
core_shell_cylinder_xy { 4, 4, 1, 20, 20, 400, 60, 60 }  
    
 Parameters:  
 Cylinder core scattering length density1e-6/Ang^2[-inf, inf]
 Cylinder shell scattering length density1e-6/Ang^2[-inf, inf]
 Solvent scattering length density1e-6/Ang^2[-inf, inf]
 Cylinder core radiusAng[0, inf]
 Cylinder shell thicknessAng[0, inf]
 Cylinder lengthAng[0, inf]
 In plane angledegrees[-inf, inf]
 Out of plane angledegrees[-inf, inf]
    
cylinder { 4, 1, 20, 400 }  
cylinder_xy { 4, 1, 20, 400, 60, 60 }  
    
 Parameters:  
 Cylinder scattering length density4e-6/Ang^2[-inf, inf]
 Solvent scattering length density1e-6/Ang^2[-inf, inf]
 Cylinder radiusAng[0, inf]
 Cylinder lengthAng[0, inf]
 In plane angledegrees[-inf, inf]
 Out of plane angledegrees[-inf, inf]
    
dab { 50.0 }  
dab_xy { 50.0 }  
    
 Parameters:  
 correlation lengthAng[0, inf]
    
ellipsoid { 4, 1, 20, 400 }  
ellipsoid_xy { 4, 1, 20, 400, 60, 60 }  
    
 Parameters:  
 Ellipsoid scattering length density1e-6/Ang^2[-inf, inf]
 Solvent scattering length density1e-6/Ang^2[-inf, inf]
 Polar radiusAng[0, inf]
 Equatorial radiusAng[0, inf]
 In plane angledegrees[-inf, inf]
 Out of plane angledegrees[-inf, inf]
    
fcc { 220, 0.06, 40, 4, 1 }  
fcc_xy { 220, 0.06, 40, 4, 1, 60, 60, 60 }  
    
 Parameters:  
 Nearest neighbour distanceAng[-inf, inf]
 Paracrystal distortion factor [-inf, inf]
 Particle radiusAng[0, inf]
 Particle scattering length density1e-6/Ang^2[-inf, inf]
 Solvent scattering length density1e-6/Ang^2[-inf, inf]
 In plane angledegrees[-inf, inf]
 Out of plane angledegrees[-inf, inf]
 Out of plane angledegrees[-inf, inf]
    
flexible_cylinder_ex { 1000.0, 100.0, 20.0, 1.5, 1.0, 6.3 }  
flexible_cylinder_ex_xy { 1000.0, 100.0, 20.0, 1.5, 1.0, 6.3 }  
    
 Parameters:  
 Length of the flexible cylinderAng[0, inf]
 Kuhn length of the flexible cylinderAng[0, inf]
 Radius of the flexible cylinderAng[0, inf]
 Axis_ratio (major_radius/radius [0, inf]
 Cylinder scattering length density1e-6/Ang^2[-inf, inf]
 Solvent scattering length density1e-6/Ang^2[-inf, inf]
    
gaussian_peak { 0.05, 0.005 }  
gaussian_peak_xy { 0.05, 0.005 }  
    
 Parameters:  
 Peak position1/Ang[-inf, inf]
 Peak width (standard deviation)1/Ang[0, inf]
    
guinier { 60.0 }  
guinier_xy { 60.0 }  
    
 Parameters:  
 Radius of GyrationAng[0, inf]
    
hardsphere { 50.0, 0.2 }  
hardsphere_xy { 50.0, 0.2 }  
    
 Parameters:  
 effective radius of hard sphereAng[0, inf]
 volume fraction of hard spheres [0, 0.74]
    
HayterMSAsq { 20.75, 19.0, 0.0192, 318.16, 0.0, 71.08 }  
HayterMSAsq_xy { 20.75, 19.0, 0.0192, 318.16, 0.0, 71.08 }  
    
 Parameters:  
 effective radius of hard sphereAng[0, inf]
 charge on sphere (in electrons)e[0, inf]
 volume fraction of spheres [0, 0.74]
 temperature, in Kelvin, for Debye length calculationK[0, inf]
 conc of salt, 1:1 electolyte, for Debye lengthM[-inf, inf]
 dielectric constant of solvent (default water), for Debye length [-inf, inf]
    
hollow_cylinder { 30.0, 20.0, 400.0, 6.3, 1 }  
hollow_cylinder_xy { 30.0, 20.0, 400.0, 6.3, 1, 90, 0 }  
    
 Parameters:  
 Cylinder radiusAng[0, inf]
 Hollow core radiusAng[0, inf]
 Cylinder lengthAng[0, inf]
 Cylinder sld1/Ang^2[-inf, inf]
 Solvent sld1/Ang^2[-inf, inf]
 Theta angledegrees[-360, 360]
 Phi angledegrees[-360, 360]
    
lamellar { 1, 6, 50 }  
lamellar_xy { 1, 6, 50 }  
    
 Parameters:  
 Layer scattering length density1e-6/Ang^2[-inf, inf]
 Solvent scattering length density1e-6/Ang^2[-inf, inf]
 Bilayer thicknessAng[0, inf]
    
lamellarCaille { 30.0, 20, 400.0, 0.1, 6.3, 1.0 }  
lamellarCaille_xy { 30.0, 20, 400.0, 0.1, 6.3, 1.0 }  
    
 Parameters:  
 sheet thicknessAng[0, inf]
 Number of layers [0, inf]
 d-spacing of Caille S(Q)Ang[0.0, inf]
 Caille parameter1/Ang^2[0.0, 0.8]
 layer scattering length density1e-6/Ang^2[-inf, inf]
 Solvent scattering length density1e-6/Ang^2[-inf, inf]
    
lamellarCailleHG { 10, 2, 30, 40.0, 0.001, 0.4, 2.0, 6 }  
lamellarCailleHG_xy { 10, 2, 30, 40.0, 0.001, 0.4, 2.0, 6 }  
    
 Parameters:  
 Tail thicknessAng[0, inf]
 head thicknessAng[0, inf]
 Number of layers [0, inf]
 d-spacing of Caille S(Q)Ang[0.0, inf]
 Caille parameter [0.0, 0.8]
 Tail scattering length density1e-6/Ang^2[-inf, inf]
 Head scattering length density1e-6/Ang^2[-inf, inf]
 Solvent scattering length density1e-6/Ang^2[-inf, inf]
    
lamellarFFHG { 15, 10, 0.4, 3.0, 6 }  
lamellarFFHG_xy { 15, 10, 0.4, 3.0, 6 }  
    
 Parameters:  
 Tail thicknessAng[0, inf]
 head thicknessAng[0, inf]
 Tail scattering length density1e-6/Ang^2[-inf, inf]
 Head scattering length density1e-6/Ang^2[-inf, inf]
 Solvent scattering length density1e-6/Ang^2[-inf, inf]
    
lamellarPC { 33.0, 20, 250.0, 0.0, 1.0, 6.34 }  
lamellarPC_xy { 33.0, 20, 250.0, 0.0, 1.0, 6.34 }  
    
 Parameters:  
 sheet thicknessAng[0, inf]
 Number of layers [0, inf]
 d-spacing of paracrystal stackAng[0.0, inf]
 d-spacing polydispersityAng[0.0, inf]
 layer scattering length density1e-6/Ang^2[-inf, inf]
 Solvent scattering length density1e-6/Ang^2[-inf, inf]
    
linear_pearls { 80.0, 350.0, 3.0, 1.0, 6.3 }  
linear_pearls_xy { 80.0, 350.0, 3.0, 1.0, 6.3 }  
    
 Parameters:  
 Radius of the pearlsAng[0, inf]
 Length of the string segment - surface to surfaceAng[0, inf]
 Number of the pearls [0, inf]
 SLD of the pearl spheres1e-6/Ang^2[-inf, inf]
 SLD of the solvent1e-6/Ang^2[-inf, inf]
    
lorentz { 50.0 }  
lorentz_xy { 50.0 }  
    
 Parameters:  
 Screening lengthAng[0, inf]
    
mass_fractal { 10.0, 1.9, 100.0 }  
mass_fractal_xy { 10.0, 1.9, 100.0 }  
    
 Parameters:  
 Particle radiusAng[0.0, inf]
 Mass fractal dimension [1.0, 6.0]
 Cut-off lengthAng[0.0, inf]
    
mass_surface_fractal { 1.8, 2.3, 86.7, 4000.0 }  
mass_surface_fractal_xy { 1.8, 2.3, 86.7, 4000.0 }  
    
 Parameters:  
 Mass fractal dimension [1e-16, 6.0]
 Surface fractal dimension [1e-16, 6.0]
 Cluster radius of gyrationAng[0.0, inf]
 Primary particle radius of gyrationAng[0.0, inf]
    
parallelepiped { 4, 1, 35, 75, 400 }  
parallelepiped_xy { 4, 1, 35, 75, 400, 60, 60, 60 }  
    
 Parameters:  
 Parallelepiped scattering length density1e-6/Ang^2[-inf, inf]
 Solvent scattering length density1e-6/Ang^2[-inf, inf]
 Shorter side of the parallelepipedAng[0, inf]
 Second side of the parallelepipedAng[0, inf]
 Larger side of the parallelepipedAng[0, inf]
 In plane angledegrees[-inf, inf]
 Out of plane angledegrees[-inf, inf]
 Rotation angle around its own c axis against q planedegrees[-inf, inf]
    
pearl_necklace { 80.0, 350.0, 2.5, 3, 1.0, 1.0, 6.3 }  
pearl_necklace_xy { 80.0, 350.0, 2.5, 3, 1.0, 1.0, 6.3 }  
    
 Parameters:  
 Mean radius of the chained spheresAngstrom[0, inf]
 Mean separation of chained particlesAngstrom[0, inf]
 Thickness of the chain linkageAngstrom[0, inf]
 Mean number of pearls in each necklacenone[0, inf]
 Scattering length density of the chained spheresAngstrom^2[-inf, inf]
 Scattering length density of the chain linkageAngstrom^2[-inf, inf]
 Scattering length density of the solventAngstrom^2[-inf, inf]
    
sphere { 1, 6, 50 }  
sphere_xy { 1, 6, 50 }  
    
 Parameters:  
 Layer scattering length density1e-6/Ang^2[-inf, inf]
 Solvent scattering length density1e-6/Ang^2[-inf, inf]
 Sphere radiusAng[0, inf]
    
star_polymer { 100.0, 3 }  
star_polymer_xy { 100.0, 3 }  
    
 Parameters:  
 Ensemble radius of gyration squared of an armAng[0.0, inf]
 Number of arms in the model [1.0, 6.0]
    
stickyhardsphere { 50.0, 0.2, 0.05, 0.2 }  
stickyhardsphere_xy { 50.0, 0.2, 0.05, 0.2 }  
    
 Parameters:  
 effective radius of hard sphereAng[0, inf]
 volume fraction of hard spheres [0, 0.74]
 perturbation parameter, epsilon [0.01, 0.1]
 stickiness, tau [-inf, inf]
    
triaxial_ellipsoid { 4, 1, 20, 400, 10 }  
triaxial_ellipsoid_xy { 4, 1, 20, 400, 10, 60, 60, 60 }  
    
 Parameters:  
 Ellipsoid scattering length density1e-6/Ang^2[-inf, inf]
 Solvent scattering length density1e-6/Ang^2[-inf, inf]
 Minor equitorial radiusAng[0, inf]
 Major equatorial radiusAng[0, inf]
 Polar radiusAng[0, inf]
 In plane angledegrees[-inf, inf]
 Out of plane angledegrees[-inf, inf]
 Out of plane angledegrees[-inf, inf]
    
*

*

*

Name convention: SasView - sasmodels

* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
SasView namesasmodels name P(Q) S(Q)MultiplyMultiplicity
BarBellbarbell - - --
BCCrystalbcc - - --
CappedCylindercapped_cylinder - - --
CoreShellCylindercore_shell_cylinder - - --
Cylindercylinder - - --
DABdab - - --
Ellipsoidellipsoid - - --
FCCrystalfcc - - --
FlexCylEllipXflexible_cylinder_ex - - --
PeakGaussgaussian_peak - - --
Guinierguinier - - --
HardsphereStructurehardsphere - Y --
HayterMSAStructureHayterMSAsq - Y --
HollowCylinderhollow_cylinder - - --
Lamellarlamellar - - --
LamellarPSlamellarCaille - - --
LamellarPSHGlamellarCailleHG - - --
LamellarFFHGlamellarFFHG - - --
LamellarPCrystallamellarPC - - --
LinearPearlslinear_pearls - - --
Lorentzlorentz - - --
MassFractalmass_fractal - - --
MassSurfaceFractalmass_surface_fractal - - --
Parallelepipedparallelepiped - - --
PearlNecklacepearl_necklace - - --
Spheresphere - - --
StarPolymerstar_polymer - - --
StickyHSStructurestickyhardsphere - Y --
TriaxialEllipsoidtriaxial_ellipsoid - - --
* * * MDOC_END * * %P * Definition parameters: * * model_index: [] Index of the applied sasview model. Recompile instrument for changes to take effect. * model_scale: [] Global scale factor for scattering kernel. For systems without inter-particle interference, the form factors can be related to the scattering intensity by the particle volume fraction. * model_pars: [] Model parameters are given as a set of comma-separated values enclosed by {}, e.g. {60} for the model index 21 (the guinier model). Consult the sasview docs for further info. * model_abs: [1/m] Absorption cross section density at 2200 m/s * * Input parameters: * * radius: [m] Outer radius of sample in (x,z) plane for cylinder/sphere * xwidth: [m] horiz. dimension of sample, as a width * yheight: [m] vert . dimension of sample, as a height for cylinder/box * zdepth: [m] depth of sample * target_index: [1] Relative index of component to focus at, e.g. next is +1 * focus_xw: [m] horiz. dimension of a rectangular area * focus_yh: [m] vert. dimension of a rectangular area * focus_aw: [deg] horiz. angular dimension of a rectangular area * focus_ah: [deg] vert. angular dimension of a rectangular area * focus_r: [m] Detector (disk-shaped) radius * * Optional parameters: * * target_x: [m] relative focus target position * target_y: [m] relative focus target position * target_z: [m] relative focus target position * * Variables calculated in the component: * * my_s: Attenuation factor due to scattering [m^-1] * my_a: Attenuation factor due to absorbtion [m^-1] * * %Link * http://www.sasview.org/sasview/user/models/model_functions.html * * %E *******************************************************************************/ DEFINE COMPONENT SasView_model SETTING PARAMETERS ( model_index=21, model_scale=1.0, vector model_pars={60,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0}, model_abs=0.5, xwidth=0, yheight=0, zdepth=0, radius=0, target_x=0, target_y=0, target_z=6, int target_index=0, focus_xw=0, focus_yh=0, focus_aw=0, focus_ah=0, focus_r=0) /* Neutron parameters: (x,y,z,vx,vy,vz,t,sx,sy,sz,p) */ SHARE %{ %include "sasview_proxy.c" %} DECLARE %{ double shape; double my_a_v; double modelpars[15]; %} INITIALIZE %{ int chosen = model_index; #ifndef SASmodel_index exit (fprintf (stderr, "SasView_model: %s: Your instrument was compiled without the define -DSASmodel_index \n" "ERROR Please recompile your instrument with -DSASmodel_index=%i to enable model_index=%i.\n", NAME_CURRENT_COMP, chosen, chosen)); #else if (!(SASmodel_index == model_index)) { exit (fprintf (stderr, "SasView_model: %s: Your instrument was compiled with -DSASmodel_index=%i \n" "ERROR Please recompile your instrument with -DSASmodel_index=%i to enable model_index=%i\n", NAME_CURRENT_COMP, SASmodel_index, chosen, chosen)); // int SASmodel_index=model_index; } #endif shape = -1; /* -1:no shape, 0:cyl, 1:box, 2:sphere */ if (xwidth && yheight && zdepth) shape = 1; else if (radius > 0 && yheight) shape = 0; else if (radius > 0 && !yheight) shape = 2; if (shape < 0) exit (fprintf (stderr, "SasView_model: %s: sample has invalid dimensions.\n" "ERROR Please check parameter values.\n", NAME_CURRENT_COMP)); /* now compute target coords if a component index is supplied */ if (!target_index && !target_x && !target_y && !target_z) target_index = 1; if (target_index) { Coords ToTarget; ToTarget = coords_sub (POS_A_COMP_INDEX (INDEX_CURRENT_COMP + target_index), POS_A_CURRENT_COMP); ToTarget = rot_apply (ROT_A_CURRENT_COMP, ToTarget); coords_get (ToTarget, &target_x, &target_y, &target_z); } if (!(target_x || target_y || target_z)) { printf ("SasView_model: %s: The target is not defined. Using direct beam (Z-axis).\n", NAME_CURRENT_COMP); target_z = 1; } my_a_v = model_abs * 2200 * 100; /* Is not yet divided by v. 100: Convert barns -> fm^2 */ int j; for (j = 0; j < 15; j++) { modelpars[j] = model_pars[j]; } %} TRACE %{ double t0, t1, v, l_full, l, l_1, dt, d_phi, theta, my_s; double aim_x = 0, aim_y = 0, aim_z = 1, axis_x, axis_y, axis_z; double arg, tmp_vx, tmp_vy, tmp_vz, vout_x, vout_y, vout_z; double f, solid_angle, vx_i, vy_i, vz_i, q, qx, qy, qz; char intersect = 0; /* Intersection neutron trajectory / sample (sample surface) */ if (shape == 0) intersect = cylinder_intersect (&t0, &t1, x, y, z, vx, vy, vz, radius, yheight); else if (shape == 1) intersect = box_intersect (&t0, &t1, x, y, z, vx, vy, vz, xwidth, yheight, zdepth); else if (shape == 2) intersect = sphere_intersect (&t0, &t1, x, y, z, vx, vy, vz, radius); if (intersect) { if (t0 < 0) ABSORB; /* Neutron enters at t=t0. */ v = sqrt (vx * vx + vy * vy + vz * vz); l_full = v * (t1 - t0); /* Length of full path through sample */ dt = rand01 () * (t1 - t0) + t0; /* Time of scattering */ PROP_DT (dt); /* Point of scattering */ l = v * (dt - t0); /* Penetration in sample */ vx_i = vx; vy_i = vy; vz_i = vz; if ((target_x || target_y || target_z)) { aim_x = target_x - x; /* Vector pointing at target (anal./det.) */ aim_y = target_y - y; aim_z = target_z - z; } if (focus_aw && focus_ah) { randvec_target_rect_angular (&vx, &vy, &vz, &solid_angle, aim_x, aim_y, aim_z, focus_aw, focus_ah, ROT_A_CURRENT_COMP); } else if (focus_xw && focus_yh) { randvec_target_rect (&vx, &vy, &vz, &solid_angle, aim_x, aim_y, aim_z, focus_xw, focus_yh, ROT_A_CURRENT_COMP); } else { randvec_target_circle (&vx, &vy, &vz, &solid_angle, aim_x, aim_y, aim_z, focus_r); } NORM (vx, vy, vz); vx *= v; vy *= v; vz *= v; qx = V2K * (vx_i - vx); qy = V2K * (vy_i - vy); qz = V2K * (vz_i - vz); q = sqrt (qx * qx + qy * qy + qz * qz); float Iq_out; Iq_out = getIq (q, qx, qy, modelpars); float vol; vol = getFormVol (modelpars); // Scale by 1.0E2 [SasView: 1/cm -> McStas: 1/m] Iq_out = model_scale * Iq_out / vol * 1.0E2; l_1 = v * t1; p *= l_full * solid_angle / (4 * PI) * Iq_out * exp (-my_a_v * (l + l_1) / v); SCATTER; } %} MCDISPLAY %{ if (shape == 0) { /* cylinder */ circle ("xz", 0, yheight / 2.0, 0, radius); circle ("xz", 0, -yheight / 2.0, 0, radius); line (-radius, -yheight / 2.0, 0, -radius, +yheight / 2.0, 0); line (+radius, -yheight / 2.0, 0, +radius, +yheight / 2.0, 0); line (0, -yheight / 2.0, -radius, 0, +yheight / 2.0, -radius); line (0, -yheight / 2.0, +radius, 0, +yheight / 2.0, +radius); } else if (shape == 1) { /* box */ double xmin = -0.5 * xwidth; double xmax = 0.5 * xwidth; double ymin = -0.5 * yheight; double ymax = 0.5 * yheight; double zmin = -0.5 * zdepth; double zmax = 0.5 * zdepth; multiline (5, xmin, ymin, zmin, xmax, ymin, zmin, xmax, ymax, zmin, xmin, ymax, zmin, xmin, ymin, zmin); multiline (5, xmin, ymin, zmax, xmax, ymin, zmax, xmax, ymax, zmax, xmin, ymax, zmax, xmin, ymin, zmax); line (xmin, ymin, zmin, xmin, ymin, zmax); line (xmax, ymin, zmin, xmax, ymin, zmax); line (xmin, ymax, zmin, xmin, ymax, zmax); line (xmax, ymax, zmin, xmax, ymax, zmax); } else if (shape == 2) { /* sphere */ circle ("xy", 0, 0.0, 0, radius); circle ("xz", 0, 0.0, 0, radius); circle ("yz", 0, 0.0, 0, radius); } %} END