source: sasview/sansmodels/src/sans/models/DABModel.py @ 1bbf2ac

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 1bbf2ac was 988130c6, checked in by Gervaise Alina <gervyh@…>, 16 years ago

self.fixed field added in models but need to change wrapper generator

  • Property mode set to 100644
File size: 2.8 KB
Line 
1#!/usr/bin/env python
2"""
3    Provide F(x) = scale/( 1 + (x*L)^2 )^(2) + background
4    DAB (Debye Anderson Brumberger) function as a BaseComponent model
5"""
6
7from sans.models.BaseComponent import BaseComponent
8import math
9
10class DABModel(BaseComponent):
11   
12    """
13        Class that evaluates a DAB model.
14       
15        F(x) = scale/( 1 + (x*L)^2 )^(2) + background
16       
17        The model has three parameters:
18            L             =  Correlation Length
19            scale         =  scale factor
20            background    =  incoherent background
21    """
22       
23    def __init__(self):
24        """ Initialization """
25       
26        # Initialize BaseComponent first, then sphere
27        BaseComponent.__init__(self)
28       
29        ## Name of the model
30        self.name = "DAB_Model"
31        self.description="""F(x) = scale/( 1 + (x*L)^2 )^(2) + background
32       
33        The model has three parameters:
34        L             =  Correlation Length
35        scale         =  scale factor
36        background    =  incoherent background"""
37        ## Define parameters
38        self.params = {}
39        self.params['length']             = 50.0
40        self.params['scale']              = 1.0
41        self.params['background']         = 0.0
42
43        ## Parameter details [units, min, max]
44        self.details = {}
45        self.details['length']            = ['', None, None]
46        self.details['scale']             = ['', None, None]
47        self.details['background']        = ['', None, None]
48        #list of parameter that cannot be fitted
49        self.fixed= []     
50    def _DAB(self, x):
51        """
52            Evaluate  F(x) = scale/( 1 + (x*L)^2 )^(2) + background
53           
54        """
55        return self.params['scale']/math.pow(( 1 + math.pow(x * self.params['length'],2)),2) \
56         + self.params['background']
57       
58   
59    def run(self, x = 0.0):
60        """ Evaluate the model
61            @param x: input q-value (float or [float, float] as [r, theta])
62            @return: (DAB value)
63        """
64        if x.__class__.__name__ == 'list':
65            return self._DAB(x[0])
66        elif x.__class__.__name__ == 'tuple':
67            raise ValueError, "Tuples are not allowed as input to BaseComponent models"
68        else:
69            return self._DAB(x)
70   
71    def runXY(self, x = 0.0):
72        """ Evaluate the model
73            @param x: input q-value (float or [float, float] as [qx, qy])
74            @return: DAB value
75        """
76        if x.__class__.__name__ == 'list':
77            q = math.sqrt(x[0]**2 + x[1]**2)
78            return self._DAB(q)
79        elif x.__class__.__name__ == 'tuple':
80            raise ValueError, "Tuples are not allowed as input to BaseComponent models"
81        else:
82            return self._DAB(x)
Note: See TracBrowser for help on using the repository browser.