# Draft minutes for Diagnostic Meeting 20210927

def _spectral_desnity_max(...):
    ...

def _spectral_density_kappa(...):
    ...

class Thomson:
    _spectal_types = {
        "maxwellian": _spectral_density_max,
        "kappa": _spectral_density_kappa
    }
    
    def __init__(*, distribution_type="maxwellian", **kwargs):
        self.parameters = kwargs
        self._spectral_density_func = self._spectral_types[distribution_type]
    
    def spectral_density(self):
        return self._spectral_densithy_func(**self.parameters)
    
    def power(self):
        ...

def thomson_fitting(wavelength, signal, type="maxwellian"):
    ...
    return Thomson()

def spectral_density(type="maxwellian", **kwargs):
    ...

def spectral_power(type="maxwellain", **kwargs):
    sd = spectral_density(type=type, **kwargs)
    return modify*sd