source: sasview/sansview/perspectives/fitting/models.py @ 7381f34

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 7381f34 was a1b2471, checked in by Jae Cho <jhjcho@…>, 14 years ago

added sld plot for onion model and etc…

  • Property mode set to 100644
File size: 18.1 KB
RevLine 
[33afff7]1
[d89f09b]2import wx
[bb18ef1]3import wx.lib.newevent
[b30f001]4import imp
[442895f]5import os,sys,math
[d89f09b]6import os.path
[8d78399]7
[d89f09b]8(ModelEvent, EVT_MODEL) = wx.lib.newevent.NewEvent()
[e7b1ccf]9from sans.guicomm.events import StatusEvent 
[33afff7]10# Time is needed by the log method
11import time
[2dbb681]12
[33afff7]13# Explicitly import from the pluginmodel module so that py2exe
14# places it in the distribution. The Model1DPlugin class is used
15# as the base class of plug-in models.
16from sans.models.pluginmodel import Model1DPlugin
[00561739]17   
[b30f001]18def log(message):
[5062bbf]19    """
20    """
[b30f001]21    out = open("plugins.log", 'a')
22    out.write("%10g%s\n" % (time.clock(), message))
23    out.close()
24
[aa92772]25def findModels():
[5062bbf]26    """
27    """
[33afff7]28    log("looking for models in: %s/plugins" % os.getcwd())
[b30f001]29    if os.path.isdir('plugins'):
30        return _findModels('plugins')
31    return []
[aa92772]32   
[bfe4644]33def _check_plugin(model, name):
34    """
[5062bbf]35    Do some checking before model adding plugins in the list
36   
37    :param model: class model to add into the plugin list
38    :param name:name of the module plugin
39   
40    :return model: model if valid model or None if not valid
41   
[bfe4644]42    """
43    #Check is the plugin is of type Model1DPlugin
44    if not issubclass(model, Model1DPlugin):
45        msg= "Plugin %s must be of type Model1DPlugin \n"%str(name)
46        log(msg)
[8d78399]47        return None
[bfe4644]48    if model.__name__!="Model":
49        msg= "Plugin %s class name must be Model \n"%str(name)
50        log(msg)
[8d78399]51        return None
[bfe4644]52    try:
53        new_instance= model()
54    except:
55        msg="Plugin %s error in __init__ \n\t: %s %s\n"%(str(name),
56                                    str(sys.exc_type),sys.exc_value)
57        log(msg)
[8d78399]58        return None
[bfe4644]59   
60    new_instance= model() 
61    if hasattr(new_instance,"function"):
62        try:
63           value=new_instance.function()
64        except:
65           msg="Plugin %s: error writing function \n\t :%s %s\n "%(str(name),
66                                    str(sys.exc_type),sys.exc_value)
67           log(msg)
[8d78399]68           return None
[bfe4644]69    else:
70       msg="Plugin  %s needs a method called function \n"%str(name)
71       log(msg)
[8d78399]72       return None
[bfe4644]73    return model
74 
75 
[1c66bc5]76def _findModels(dir):
[5062bbf]77    """
78    """
[1c66bc5]79    # List of plugin objects
80    plugins = []
81    # Go through files in plug-in directory
82    try:
83        list = os.listdir(dir)
84        for item in list:
85            toks = os.path.splitext(os.path.basename(item))
86            if toks[1]=='.py' and not toks[0]=='__init__':
87                name = toks[0]
88           
89                path = [os.path.abspath(dir)]
90                file = None
91                try:
92                    (file, path, info) = imp.find_module(name, path)
93                    module = imp.load_module( name, file, item, info )
94                    if hasattr(module, "Model"):
95                        try:
[bfe4644]96                            if _check_plugin(module.Model, name)!=None:
97                                plugins.append(module.Model)
[1c66bc5]98                        except:
[bfe4644]99                            msg="Error accessing Model"
100                            msg+="in %s\n  %s %s\n" % (name,
101                                    str(sys.exc_type), sys.exc_value)
102                            log(msg)
[1c66bc5]103                except:
[bfe4644]104                    msg="Error accessing Model"
105                    msg +=" in %s\n  %s %s \n" %(name,
106                                    str(sys.exc_type), sys.exc_value)
107                    log(msg)
[1c66bc5]108                finally:
[a92d51b]109             
[1c66bc5]110                    if not file==None:
111                        file.close()
112    except:
[33afff7]113        # Don't deal with bad plug-in imports. Just skip.
[1c66bc5]114        pass
115    return plugins
[bb18ef1]116
117class ModelList(object):
118    """
[5062bbf]119    Contains dictionary of model and their type
[bb18ef1]120    """
121    def __init__(self):
[5062bbf]122        """
123        """
124        self.mydict = {}
[bb18ef1]125       
126    def set_list(self, name, mylist):
127        """
[5062bbf]128        :param name: the type of the list
129        :param mylist: the list to add
130       
[bb18ef1]131        """
132        if name not in self.mydict.keys():
133            self.mydict[name] = mylist
134           
135           
136    def get_list(self):
137        """
[5062bbf]138        return all the list stored in a dictionary object
[bb18ef1]139        """
140        return self.mydict
141       
[d89f09b]142class ModelManager:
[5062bbf]143    """
144    """
[bb18ef1]145    ## external dict for models
146    model_combobox = ModelList()
147    ## Dictionary of form models
148    form_factor_dict = {}
149    ## dictionary of other
150    struct_factor_dict = {}
151    ##list of form factors
152    shape_list =[]
153    ## independent shape model list
154    shape_indep_list = []
155    ##list of structure factors
156    struct_list= []
[376916c]157    ##list of model allowing multiplication
158    multiplication_factor=[]
[e87f9fc]159    ##list of multifunctional shapes
160    multi_func_list =[]
[bb18ef1]161    ## list of added models
[b30f001]162    plugins=[]
[bb18ef1]163    ## Event owner (guiframe)
[d89f09b]164    event_owner = None
165   
166    def _getModelList(self):
167        """
[5062bbf]168        List of models we want to make available by default
169        for this application
170   
171        :return: the next free event ID following the new menu events
[e7b1ccf]172       
[d89f09b]173        """
[bb18ef1]174        ## form factor
[442895f]175        from sans.models.SphereModel import SphereModel
[bb18ef1]176        self.shape_list.append(SphereModel)
[376916c]177        self.multiplication_factor.append(SphereModel)
[442895f]178       
[e87f9fc]179        from sans.models.SphereExpShellModel import SphereExpShellModel
180        self.shape_list.append(SphereExpShellModel)
181        self.multiplication_factor.append(SphereExpShellModel)
182        self.multi_func_list.append(SphereExpShellModel)
183       
[ce07fa8]184        from sans.models.FuzzySphereModel import FuzzySphereModel
185        self.shape_list.append(FuzzySphereModel)
186        self.multiplication_factor.append(FuzzySphereModel)
187       
[442895f]188        from sans.models.CoreShellModel import CoreShellModel
[bb18ef1]189        self.shape_list.append(CoreShellModel)
[5eb9154]190        self.multiplication_factor.append(CoreShellModel)
[4523b68]191       
192        from sans.models.CoreMultiShellModel import CoreMultiShellModel
193        self.shape_list.append(CoreMultiShellModel)
194        self.multiplication_factor.append(CoreMultiShellModel)
[a1b2471]195        self.multi_func_list.append(CoreMultiShellModel)
[cee6867]196       
[eddff027]197        from sans.models.VesicleModel import VesicleModel
198        self.shape_list.append(VesicleModel)
[5eb9154]199        self.multiplication_factor.append(VesicleModel)
200       
201        from sans.models.MultiShellModel import MultiShellModel
202        self.shape_list.append(MultiShellModel)
203        self.multiplication_factor.append(MultiShellModel)
[eddff027]204       
[cee6867]205        from sans.models.BinaryHSModel import BinaryHSModel
206        self.shape_list.append(BinaryHSModel)
207       
[eddff027]208        from sans.models.CylinderModel import CylinderModel
209        self.shape_list.append(CylinderModel)
210        self.multiplication_factor.append(CylinderModel)
211       
212        from sans.models.CoreShellCylinderModel import CoreShellCylinderModel
213        self.shape_list.append(CoreShellCylinderModel)
[5eb9154]214        self.multiplication_factor.append(CoreShellCylinderModel)
[cee6867]215       
216        from sans.models.HollowCylinderModel import HollowCylinderModel
217        self.shape_list.append(HollowCylinderModel)
[5eb9154]218        self.multiplication_factor.append(HollowCylinderModel)
[eddff027]219             
220        from sans.models.FlexibleCylinderModel import FlexibleCylinderModel
221        self.shape_list.append(FlexibleCylinderModel)
[72f719b]222
[ce07fa8]223        from sans.models.FlexCylEllipXModel import FlexCylEllipXModel
224        self.shape_list.append(FlexCylEllipXModel)
[eddff027]225       
226        from sans.models.StackedDisksModel import StackedDisksModel
227        self.shape_list.append(StackedDisksModel)
[5eb9154]228        self.multiplication_factor.append(StackedDisksModel)
[eddff027]229       
230        from sans.models.ParallelepipedModel import ParallelepipedModel
231        self.shape_list.append(ParallelepipedModel)
[72f719b]232        self.multiplication_factor.append(ParallelepipedModel)
[cee6867]233       
[442895f]234        from sans.models.EllipticalCylinderModel import EllipticalCylinderModel
[bb18ef1]235        self.shape_list.append(EllipticalCylinderModel)
[72f719b]236        self.multiplication_factor.append(EllipticalCylinderModel)
237               
[442895f]238        from sans.models.EllipsoidModel import EllipsoidModel
[bb18ef1]239        self.shape_list.append(EllipsoidModel)
[376916c]240        self.multiplication_factor.append(EllipsoidModel)
[eddff027]241     
242        from sans.models.CoreShellEllipsoidModel import CoreShellEllipsoidModel
243        self.shape_list.append(CoreShellEllipsoidModel)
[5eb9154]244        self.multiplication_factor.append(CoreShellEllipsoidModel)
[bb18ef1]245         
[e65050e]246        from sans.models.TriaxialEllipsoidModel import TriaxialEllipsoidModel
247        self.shape_list.append(TriaxialEllipsoidModel)
[9002927]248        self.multiplication_factor.append(TriaxialEllipsoidModel)
[e65050e]249       
250        from sans.models.LamellarModel import LamellarModel
251        self.shape_list.append(LamellarModel)
252       
253        from sans.models.LamellarFFHGModel import LamellarFFHGModel
254        self.shape_list.append(LamellarFFHGModel)
255       
256        from sans.models.LamellarPSModel import LamellarPSModel
257        self.shape_list.append(LamellarPSModel)
[7a69683]258     
[e65050e]259        from sans.models.LamellarPSHGModel import LamellarPSHGModel
260        self.shape_list.append(LamellarPSHGModel)
[7a69683]261     
[376916c]262        ## Structure factor
[8346667]263        from sans.models.SquareWellStructure import SquareWellStructure
[bb18ef1]264        self.struct_list.append(SquareWellStructure)
[8346667]265       
266        from sans.models.HardsphereStructure import HardsphereStructure
[bb18ef1]267        self.struct_list.append(HardsphereStructure)
268         
[8346667]269        from sans.models.StickyHSStructure import StickyHSStructure
[bb18ef1]270        self.struct_list.append(StickyHSStructure)
[8346667]271       
272        from sans.models.HayterMSAStructure import HayterMSAStructure
[bb18ef1]273        self.struct_list.append(HayterMSAStructure)
[5062bbf]274   
275        ##shape-independent models 
[ce07fa8]276        from sans.models.PowerLawAbsModel import PowerLawAbsModel
277        self.shape_indep_list.append( PowerLawAbsModel )
278       
[442895f]279        from sans.models.BEPolyelectrolyte import BEPolyelectrolyte
[bb18ef1]280        self.shape_indep_list.append(BEPolyelectrolyte )
281        self.form_factor_dict[str(wx.NewId())] =  [SphereModel]
[ce07fa8]282
[442895f]283        from sans.models.DABModel import DABModel
[bb18ef1]284        self.shape_indep_list.append(DABModel )
[442895f]285       
[ce07fa8]286        from sans.models.DebyeModel import DebyeModel
287        self.shape_indep_list.append(DebyeModel )
288       
[f39511b]289        from sans.models.GuinierModel import GuinierModel
[bb18ef1]290        self.shape_indep_list.append(GuinierModel )
[f39511b]291       
[ce07fa8]292        from sans.models.FractalModel import FractalModel
293        self.shape_indep_list.append(FractalModel )
[bb18ef1]294       
[ce07fa8]295        from sans.models.LorentzModel import LorentzModel
296        self.shape_indep_list.append( LorentzModel) 
[442895f]297       
[cee6867]298        from sans.models.PeakGaussModel import PeakGaussModel
299        self.shape_indep_list.append(PeakGaussModel)
300       
301        from sans.models.PeakLorentzModel import PeakLorentzModel
302        self.shape_indep_list.append(PeakLorentzModel)
303       
[ce07fa8]304        from sans.models.Poly_GaussCoil import Poly_GaussCoil
305        self.shape_indep_list.append(Poly_GaussCoil)
306                 
307        from sans.models.PorodModel import PorodModel
308        self.shape_indep_list.append(PorodModel )
[442895f]309       
[ce07fa8]310        #FractalModel (a c-model)will be used.
311        #from sans.models.FractalAbsModel import FractalAbsModel
312        #self.shape_indep_list.append(FractalAbsModel)
[81bece4]313       
[442895f]314        from sans.models.TeubnerStreyModel import TeubnerStreyModel
[bb18ef1]315        self.shape_indep_list.append(TeubnerStreyModel )
[eddff027]316       
317        from sans.models.LineModel import LineModel
318        self.shape_indep_list.append(LineModel)
[5062bbf]319       
[49b7efa]320        #Looking for plugins
321        self.plugins = findModels()
[d89f09b]322        return 0
323
324   
325    def populate_menu(self, modelmenu, event_owner):
326        """
[5062bbf]327        Populate a menu with our models
328       
329        :param id: first menu event ID to use when binding the menu events
330        :param modelmenu: wx.Menu object to populate
331        :param event_owner: wx object to bind the menu events to
332       
333        :return: the next free event ID following the new menu events
334       
[d89f09b]335        """
[bb18ef1]336        ## Fill model lists
[d89f09b]337        self._getModelList()
[bb18ef1]338        ## store reference to model menu of guiframe
339        self.modelmenu = modelmenu
340        ## guiframe reference
[d89f09b]341        self.event_owner = event_owner
[bb18ef1]342       
343        shape_submenu = wx.Menu()
344        shape_indep_submenu = wx.Menu()
345        structure_factor = wx.Menu()
[b30f001]346        added_models = wx.Menu()
[376916c]347        multip_models = wx.Menu()
[bb18ef1]348        ## create menu with shape
[5062bbf]349        self._fill_simple_menu(menuinfo=["Shapes",shape_submenu," simple shape"],
350                         list1=self.shape_list)
[376916c]351       
[5062bbf]352        self._fill_simple_menu(menuinfo=["Shape-Independent",shape_indep_submenu,
[bb18ef1]353                                    "List of shape-independent models"],
[5062bbf]354                         list1=self.shape_indep_list )
[bb18ef1]355       
[5062bbf]356        self._fill_simple_menu(menuinfo=["Structure Factors",structure_factor,
[bb18ef1]357                                          "List of Structure factors models" ],
[5062bbf]358                                list1=self.struct_list)
[bb18ef1]359       
[5062bbf]360        self._fill_plugin_menu(menuinfo=["Customized Models", added_models,
[376916c]361                                            "List of additional models"],
[5062bbf]362                                 list1=self.plugins)
[376916c]363       
364        self._fill_menu(menuinfo=["P(Q)*S(Q)",multip_models,
365                                  "mulplication of 2 models"],
[5062bbf]366                                   list1=self.multiplication_factor ,
367                                   list2= self.struct_list)
[d89f09b]368        return 0
369   
[5062bbf]370    def _fill_plugin_menu(self, menuinfo, list1):
[bfe4644]371        """
[5062bbf]372        fill the plugin menu with costumized models
[bfe4644]373        """
374        if len(list1)==0:
375            id = wx.NewId() 
376            msg= "No model available check plugins.log for errors to fix problem"
377            menuinfo[1].Append(int(id),"Empty",msg)
378        self._fill_simple_menu( menuinfo,list1)
379       
[5062bbf]380    def _fill_simple_menu(self, menuinfo, list1):
[bb18ef1]381        """
[5062bbf]382        Fill the menu with list item
383       
384        :param modelmenu: the menu to fill
385        :param menuinfo: submenu item for the first column of this modelmenu
386                         with info.Should be a list :
387                         [name(string) , menu(wx.menu), help(string)]
388        :param list1: contains item (form factor )to fill modelmenu second column
389       
[bb18ef1]390        """
391        if len(list1)>0:
392            self.model_combobox.set_list(menuinfo[0],list1)
[e7b1ccf]393           
[bb18ef1]394            for item in list1:
[e7b1ccf]395                try:
396                    id = wx.NewId() 
397                    struct_factor=item()
398                    struct_name = struct_factor.__class__.__name__
399                    if hasattr(struct_factor, "name"):
400                        struct_name = struct_factor.name
401                       
402                    menuinfo[1].Append(int(id),struct_name,struct_name)
403                    if not  item in self.struct_factor_dict.itervalues():
404                        self.struct_factor_dict[str(id)]= item
405                    wx.EVT_MENU(self.event_owner, int(id), self._on_model)
406                except:
407                    msg= "Error Occured: %s"%sys.exc_value
408                    wx.PostEvent(self.event_owner, StatusEvent(status=msg))
[bb18ef1]409               
410        id = wx.NewId()         
411        self.modelmenu.AppendMenu(id, menuinfo[0],menuinfo[1],menuinfo[2])
412       
[5062bbf]413    def _fill_menu(self, menuinfo, list1, list2):
[bb18ef1]414        """
[5062bbf]415        Fill the menu with list item
416       
417        :param menuinfo: submenu item for the first column of this modelmenu
418                         with info.Should be a list :
419                         [name(string) , menu(wx.menu), help(string)]
420        :param list1: contains item (form factor )to fill modelmenu second column
421        :param list2: contains item (Structure factor )to fill modelmenu
422                third column
423               
[bb18ef1]424        """
425        if len(list1)>0:
426            self.model_combobox.set_list(menuinfo[0],list1)
427           
428            for item in list1:   
429                form_factor= item()
430                form_name = form_factor.__class__.__name__
431                if hasattr(form_factor, "name"):
432                    form_name = form_factor.name
433                ### store form factor to return to other users   
434                newmenu= wx.Menu()
435                if len(list2)>0:
436                    for model  in list2:
437                        id = wx.NewId()
438                        struct_factor = model()
439                        name = struct_factor.__class__.__name__
440                        if hasattr(struct_factor, "name"):
441                            name = struct_factor.name
442                        newmenu.Append(id,name, name)
443                        wx.EVT_MENU(self.event_owner, int(id), self._on_model)
444                        ## save form_fact and struct_fact
445                        self.form_factor_dict[int(id)] = [form_factor,struct_factor]
446                       
447                form_id= wx.NewId()   
448                menuinfo[1].AppendMenu(int(form_id), form_name,newmenu,menuinfo[2])
449        id=wx.NewId()
450        self.modelmenu.AppendMenu(id,menuinfo[0],menuinfo[1], menuinfo[2])
451       
[d89f09b]452    def _on_model(self, evt):
453        """
[5062bbf]454        React to a model menu event
455       
456        :param event: wx menu event
457       
[d89f09b]458        """
[bb18ef1]459        if int(evt.GetId()) in self.form_factor_dict.keys():
460            from sans.models.MultiplicationModel import MultiplicationModel
461            model1, model2 = self.form_factor_dict[int(evt.GetId())]
[5062bbf]462            model = MultiplicationModel(model1, model2)   
[bb18ef1]463        else:
464            model= self.struct_factor_dict[str(evt.GetId())]()
465        evt = ModelEvent( model= model )
466        wx.PostEvent(self.event_owner, evt)
[d89f09b]467       
468    def get_model_list(self):   
[5062bbf]469        """
470        return dictionary of models for fitpanel use
471       
472        """
[376916c]473        self.model_combobox.set_list("multiplication", self.multiplication_factor)
[e87f9fc]474        self.model_combobox.set_list("Multi-Functions", self.multi_func_list)
[bb18ef1]475        return self.model_combobox
[d89f09b]476   
[376916c]477 
[bb18ef1]478       
479   
[d89f09b]480   
[bb18ef1]481 
Note: See TracBrowser for help on using the repository browser.