source: sasview/sansmodels/include/StickyHS.h @ 3fdb68a

ESS_GUIESS_GUI_DocsESS_GUI_batch_fittingESS_GUI_bumps_abstractionESS_GUI_iss1116ESS_GUI_iss879ESS_GUI_iss959ESS_GUI_openclESS_GUI_orderingESS_GUI_sync_sascalccostrafo411magnetic_scattrelease-4.1.1release-4.1.2release-4.2.2release_4.0.1ticket-1009ticket-1094-headlessticket-1242-2d-resolutionticket-1243ticket-1249ticket885unittest-saveload
Last change on this file since 3fdb68a was 6319646, checked in by Jae Cho <jhjcho@…>, 13 years ago

new includes for scale fix for P*S

  • Property mode set to 100644
File size: 1.6 KB
Line 
1#if !defined(StickyHS_h)
2#define StickyHS_h
3#include "parameters.hh"
4
5/**
6 * Structure definition for StickyHS_struct (struncture factor) parameters
7 */
8//[PYTHONCLASS] = StickyHSStructure
9//[DISP_PARAMS] = effect_radius
10//[DESCRIPTION] =<text> Structure Factor for interacting particles:                               .
11//
12//  The interaction potential is
13//
14//                     U(r)= inf , r < 2R
15//                         = -Uo  , 2R < r < 2R + w
16//                         = 0   , r >= 2R +w
17//
18//                                              R: effective radius of the hardsphere
19//                     stickiness = [exp(Uo/kT)]/(12*perturb)
20//                     perturb = w/(w+ 2R) , 0.01 =< w <= 0.1
21//                     w: The width of the square well ,w > 0
22//                                              v: The volume fraction , v > 0
23//
24//                     Ref: Menon, S. V. G.,et.al., J. Chem.
25//                      Phys., 1991, 95(12), 9186-9190.
26//                              </text>
27//[FIXED]= effect_radius.width
28
29class StickyHSStructure{
30public:
31  // Model parameters
32  /// effect radius of hardsphere [A]
33  //  [DEFAULT]=effect_radius=50.0 [A]
34  Parameter effect_radius;
35
36  /// Volume fraction
37  //  [DEFAULT]=volfraction= 0.1
38  Parameter volfraction;
39
40  /// Perturbation Parameter ( between 0.01 - 0.1)
41  //  [DEFAULT]=perturb=0.05
42  Parameter perturb;
43
44  /// Stickiness
45  //  [DEFAULT]=stickiness= 0.2
46  Parameter stickiness;
47
48  // Constructor
49  StickyHSStructure();
50
51  // Operators to get I(Q)
52  double operator()(double q);
53  double operator()(double qx, double qy);
54  double calculate_ER();
55  double calculate_VR();
56  double evaluate_rphi(double q, double phi);
57};
58
59#endif
Note: See TracBrowser for help on using the repository browser.