source: sasview/sansview/perspectives/fitting/models.py @ 6aa47df

import wx
import imp
import os,sys,math
import os.path

(ModelEvent, EVT_MODEL) = wx.lib.newevent.NewEvent()
class ModelInfo(object):
    def __init__(self,model,description=None):
        self.model=model
        self.description=description
    def set_description(self, descrition):
        self.description =str(description)
    def get_description(self):
        return self.description
def log(message):
    print message
    out = open("plugins.log", 'a')
    out.write("%10g:  %s\n" % (time.clock(), message))
    out.close()

def findModels():
    print "looking for models"
    if os.path.isdir('plugins'):
        return _findModels('plugins')
    return []
    
def _findModels(dir):
    # List of plugin objects
    plugins = []
    # Go through files in plug-in directory
    try:
        list = os.listdir(dir)
        for item in list:
            toks = os.path.splitext(os.path.basename(item))
            if toks[1]=='.py' and not toks[0]=='__init__':
                name = toks[0]
            
                path = [os.path.abspath(dir)]
                file = None
                try:
                    (file, path, info) = imp.find_module(name, path)
                    module = imp.load_module( name, file, item, info )
                    if hasattr(module, "Model"):
                        try:
                            plugins.append(module.Model)
                        except:
                            log("Error accessing Model in %s\n  %s" % (name, sys.exc_value))
                except:
                    log("Error accessing Model in %s\n  %s" % (name, sys.exc_value))
                finally:
                    if not file==None:
                        file.close()
    except:
        pass
    return plugins
class ModelManager:
    
    ## Dictionary of models
    model_list = {}
    indep_model_list = {}
    model_list_box = {}
    custom_models={}
    plugins=[]
    indep_model=[]
    ## Event owner
    event_owner = None
    
    def _getModelList(self):
        """
            List of models we want to make available by default
            for this application
            
            @param id: first event ID to register the menu events with
            @return: the next free event ID following the new menu events
        """
        self.model_list = {}
        self.model_list_box = {}
        
        model_info="shape-based models"
        
        from sans.models.SphereModel import SphereModel
        self.model_list[str(wx.NewId())] =  ModelInfo(SphereModel , model_info)
        
        from sans.models.CylinderModel import CylinderModel
        self.model_list[str(wx.NewId())] = ModelInfo(CylinderModel , model_info)
      
        from sans.models.CoreShellModel import CoreShellModel
        self.model_list[str(wx.NewId())] = ModelInfo(CoreShellModel , model_info)
        
        from sans.models.CoreShellCylinderModel import CoreShellCylinderModel
        self.model_list[str(wx.NewId())] = ModelInfo(CoreShellCylinderModel , model_info)
        
        from sans.models.EllipticalCylinderModel import EllipticalCylinderModel
        self.model_list[str(wx.NewId())] = ModelInfo(EllipticalCylinderModel , model_info)
        
        from sans.models.EllipsoidModel import EllipsoidModel
        self.model_list[str(wx.NewId())] = ModelInfo(EllipsoidModel , model_info)
        
        from sans.guitools.LineModel import LineModel
        self.model_list[str(wx.NewId())]  = ModelInfo(LineModel , model_info)
        
        
        model_info="shape-independent models"
        
        from sans.models.BEPolyelectrolyte import BEPolyelectrolyte
        self.indep_model.append( ModelInfo( BEPolyelectrolyte , model_info) )
        
        from sans.models.DABModel import DABModel
        self.indep_model.append( ModelInfo(DABModel , model_info+
        "\n Evaluates F(x) = scale/( 1 + (x*L)^2 )^(2) +bkd") )
        
        from sans.models.DebyeModel import DebyeModel
        self.indep_model.append( ModelInfo(DebyeModel , model_info+
        "\n Evaluates F(x) = 2( exp(-x)+x -1 )/x**2") )
        
        from sans.models.FractalModel import FractalModel
        class FractalAbsModel(FractalModel):
            def _Fractal(self, x):
                return FractalModel._Fractal(self, math.fabs(x))
        self.indep_model.append( ModelInfo(FractalAbsModel , model_info) )
        
        from sans.models.LorentzModel import LorentzModel
        self.indep_model.append( ModelInfo(LorentzModel , model_info+
                "\n Evaluates F(x)= scale/( 1 + (x*L)^2 ) + bkd ") ) 
            
        from sans.models.PowerLawModel import PowerLawModel
        class PowerLawAbsModel(PowerLawModel):
            def _PowerLaw(self, x):
                try:
                    return PowerLawModel._PowerLaw(self, math.fabs(x))
                except:
                    print sys.exc_value  
        self.indep_model.append( ModelInfo(PowerLawAbsModel , model_info+
        "\n Evaluates abs(F(x)) \n with F(x) = scale* (x)^(m) + bkd ") )
                    
       
        from sans.models.TeubnerStreyModel import TeubnerStreyModel
        self.indep_model.append( ModelInfo(TeubnerStreyModel , model_info) )
        
        
          
        
        
      
        
        #Looking for plugins
        self.plugins = findModels()
       
        return 0

    
    def populate_menu(self, modelmenu, event_owner):
        """
            Populate a menu with our models
            
            @param id: first menu event ID to use when binding the menu events
            @param modelmenu: wx.Menu object to populate
            @param event_owner: wx object to bind the menu events to
            @return: the next free event ID following the new menu events
        """
        self._getModelList()
        self.event_owner = event_owner
        
        shape_submenu= wx.Menu() 
        indep_submenu = wx.Menu()
        added_models = wx.Menu()
        
        for id_str,value in self.model_list.iteritems():
            item = self.model_list[id_str]
            name = item.model.__name__
            if hasattr(item, "name"):
                name = item.model.name
            self.model_list_box[name] =value.model
            shape_submenu.Append(int(id_str), name, name)
            wx.EVT_MENU(event_owner, int(id_str), self._on_model)
        modelmenu.AppendMenu(wx.NewId(), "Shapes...", shape_submenu, "List of shape-based models")
        
        id = wx.NewId()
        if len(self.indep_model_list) == 0:
            for item in self.indep_model:
                if item not in self.indep_model_list.values():
                    self.indep_model_list[str(id)] = item
                    self.model_list[str(id)]=item
                    if hasattr(item, "name"):
                        name = item.model.name
                    else:
                        name = item.model.__name__
                    indep_submenu.Append(id,name, name)
                    self.model_list_box[name] =item.model
                    wx.EVT_MENU(event_owner, int(id), self._on_model)
                    id = wx.NewId()         
        modelmenu.AppendMenu(wx.NewId(), "Shape-independent...", indep_submenu, "List of shape-independent models")
        
        
        
        model_info="additional models"
        id = wx.NewId()
        if len(self.custom_models) == 0:
            for item in self.plugins:
                if item not in self.custom_models.values():
                    self.custom_models[str(id)] = item
                    
                    self.model_list[str(id)]=ModelInfo(item,model_info)
                    if hasattr(item, "name"):
                        name = item.name
                    else:
                        name = item.__name__
                    added_models.Append(id, name, name)
                    self.model_list_box[name] =item
                    wx.EVT_MENU(event_owner, int(id), self._on_model)
                    id = wx.NewId()
        modelmenu.AppendMenu(wx.NewId(),"Added models...", added_models, "List of additional models")
        return 0
    
    def _on_model(self, evt):
        """
            React to a model menu event
            @param event: wx menu event
        """
        if str(evt.GetId()) in self.model_list.keys():
            # Notify the application manager that a new model has been set
            #self.app_manager.set_model(self.model_list[str(evt.GetId())]())
            
            #TODO: post a model event to update all panels that need
            #evt = ModelEvent(model=self.model_list[str(evt.GetId())]())
           
            model = self.model_list[str(evt.GetId())]
            evt = ModelEvent(modelinfo=model)
            wx.PostEvent(self.event_owner, evt)
        
    def get_model_list(self):    
        """ @ return dictionary of models for fitpanel use """
        return self.model_list_box