data-manipulation-utilities 0.1.7__tar.gz → 0.1.9__tar.gz

{data_manipulation_utilities-0.1.7 → data_manipulation_utilities-0.1.9}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.2
 Name: data_manipulation_utilities
-Version: 0.1.7
+Version: 0.1.9
 Description-Content-Type: text/markdown
 Requires-Dist: logzero
 Requires-Dist: PyYAML
@@ -121,6 +121,24 @@ samples:
 
 ## PDFs
 
+### Model building
+
+In order to do complex fits, one often needs PDFs with many parameters, which need to be added.
+In these PDFs certain parameters (e.g. $\mu$ or $\sigma$) need to be shared. This project provides
+`ModelFactory`, which can do this as shown below:
+
+```python
+from dmu.stats.model_factory import ModelFactory
+
+l_pdf = ['cbr'] + 2 * ['cbl']
+l_shr = ['mu', 'sg']
+mod   = ModelFactory(obs = Data.obs, l_pdf = l_pdf, l_shared=l_shr)
+pdf   = mod.get_pdf()
+```
+
+where the model is a sum of three `CrystallBall` PDFs, one with a right tail and two with a left tail.
+The `mu` and `sg` parameters are shared.
+
 ### Printing PDFs
 
 One can print a zfit PDF by doing:
@@ -231,6 +249,44 @@ likelihood :
     nbins : 100 #If specified, will do binned likelihood fit instead of unbinned
 ```
 
+## Minimizers
+
+These are alternative implementations of the minimizers in zfit meant to be used for special types of fits.
+
+### Anealing minimizer
+
+This minimizer is meant to be used for fits to models with many parameters, where multiple minima are expected in the
+likelihood. The minimizer use is illustrated in:
+
+```python
+from dmu.stats.minimizers  import AnealingMinimizer
+
+nll       = _get_nll()
+minimizer = AnealingMinimizer(ntries=10, pvalue=0.05)
+res       = minimizer.minimize(nll)
+```
+
+this will:
+
+- Take the `NLL` object.
+- Try fitting at most 10 times
+- After each fit, calculate the goodness of fit (in this case the p-value)
+- Stop when the number of tries has been exhausted or the p-value reached is higher than `0.05`
+- If the fit has not succeeded because of convergence, validity or goodness of fit issues, 
+randomize the parameters and try again.
+- If the desired goodness of fit has not been achieved, pick the best result.
+- Return the `FitResult` object and set the PDF to the final fit result.
+
+The $\chi^2/Ndof$ can also be used as in:
+
+```python
+from dmu.stats.minimizers  import AnealingMinimizer
+
+nll       = _get_nll()
+minimizer = AnealingMinimizer(ntries=10, chi2ndof=1.00)
+res       = minimizer.minimize(nll)
+```
+
 ## Fit plotting
 
 The class `ZFitPlotter` can be used to plot fits done with zfit. For a complete set of examples of how to use

{data_manipulation_utilities-0.1.7 → data_manipulation_utilities-0.1.9}/README.md

@@ -101,6 +101,24 @@ samples:
 
 ## PDFs
 
+### Model building
+
+In order to do complex fits, one often needs PDFs with many parameters, which need to be added.
+In these PDFs certain parameters (e.g. $\mu$ or $\sigma$) need to be shared. This project provides
+`ModelFactory`, which can do this as shown below:
+
+```python
+from dmu.stats.model_factory import ModelFactory
+
+l_pdf = ['cbr'] + 2 * ['cbl']
+l_shr = ['mu', 'sg']
+mod   = ModelFactory(obs = Data.obs, l_pdf = l_pdf, l_shared=l_shr)
+pdf   = mod.get_pdf()
+```
+
+where the model is a sum of three `CrystallBall` PDFs, one with a right tail and two with a left tail.
+The `mu` and `sg` parameters are shared.
+
 ### Printing PDFs
 
 One can print a zfit PDF by doing:
@@ -211,6 +229,44 @@ likelihood :
     nbins : 100 #If specified, will do binned likelihood fit instead of unbinned
 ```
 
+## Minimizers
+
+These are alternative implementations of the minimizers in zfit meant to be used for special types of fits.
+
+### Anealing minimizer
+
+This minimizer is meant to be used for fits to models with many parameters, where multiple minima are expected in the
+likelihood. The minimizer use is illustrated in:
+
+```python
+from dmu.stats.minimizers  import AnealingMinimizer
+
+nll       = _get_nll()
+minimizer = AnealingMinimizer(ntries=10, pvalue=0.05)
+res       = minimizer.minimize(nll)
+```
+
+this will:
+
+- Take the `NLL` object.
+- Try fitting at most 10 times
+- After each fit, calculate the goodness of fit (in this case the p-value)
+- Stop when the number of tries has been exhausted or the p-value reached is higher than `0.05`
+- If the fit has not succeeded because of convergence, validity or goodness of fit issues, 
+randomize the parameters and try again.
+- If the desired goodness of fit has not been achieved, pick the best result.
+- Return the `FitResult` object and set the PDF to the final fit result.
+
+The $\chi^2/Ndof$ can also be used as in:
+
+```python
+from dmu.stats.minimizers  import AnealingMinimizer
+
+nll       = _get_nll()
+minimizer = AnealingMinimizer(ntries=10, chi2ndof=1.00)
+res       = minimizer.minimize(nll)
+```
+
 ## Fit plotting
 
 The class `ZFitPlotter` can be used to plot fits done with zfit. For a complete set of examples of how to use

{data_manipulation_utilities-0.1.7 → data_manipulation_utilities-0.1.9}/pyproject.toml

@@ -1,6 +1,6 @@
 [project]
 name        = 'data_manipulation_utilities'
-version     = '0.1.7'
+version     = '0.1.9'
 readme      = 'README.md'
 dependencies= [
 'logzero',

{data_manipulation_utilities-0.1.7 → data_manipulation_utilities-0.1.9}/src/data_manipulation_utilities.egg-info/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.2
 Name: data_manipulation_utilities
-Version: 0.1.7
+Version: 0.1.9
 Description-Content-Type: text/markdown
 Requires-Dist: logzero
 Requires-Dist: PyYAML
@@ -121,6 +121,24 @@ samples:
 
 ## PDFs
 
+### Model building
+
+In order to do complex fits, one often needs PDFs with many parameters, which need to be added.
+In these PDFs certain parameters (e.g. $\mu$ or $\sigma$) need to be shared. This project provides
+`ModelFactory`, which can do this as shown below:
+
+```python
+from dmu.stats.model_factory import ModelFactory
+
+l_pdf = ['cbr'] + 2 * ['cbl']
+l_shr = ['mu', 'sg']
+mod   = ModelFactory(obs = Data.obs, l_pdf = l_pdf, l_shared=l_shr)
+pdf   = mod.get_pdf()
+```
+
+where the model is a sum of three `CrystallBall` PDFs, one with a right tail and two with a left tail.
+The `mu` and `sg` parameters are shared.
+
 ### Printing PDFs
 
 One can print a zfit PDF by doing:
@@ -231,6 +249,44 @@ likelihood :
     nbins : 100 #If specified, will do binned likelihood fit instead of unbinned
 ```
 
+## Minimizers
+
+These are alternative implementations of the minimizers in zfit meant to be used for special types of fits.
+
+### Anealing minimizer
+
+This minimizer is meant to be used for fits to models with many parameters, where multiple minima are expected in the
+likelihood. The minimizer use is illustrated in:
+
+```python
+from dmu.stats.minimizers  import AnealingMinimizer
+
+nll       = _get_nll()
+minimizer = AnealingMinimizer(ntries=10, pvalue=0.05)
+res       = minimizer.minimize(nll)
+```
+
+this will:
+
+- Take the `NLL` object.
+- Try fitting at most 10 times
+- After each fit, calculate the goodness of fit (in this case the p-value)
+- Stop when the number of tries has been exhausted or the p-value reached is higher than `0.05`
+- If the fit has not succeeded because of convergence, validity or goodness of fit issues, 
+randomize the parameters and try again.
+- If the desired goodness of fit has not been achieved, pick the best result.
+- Return the `FitResult` object and set the PDF to the final fit result.
+
+The $\chi^2/Ndof$ can also be used as in:
+
+```python
+from dmu.stats.minimizers  import AnealingMinimizer
+
+nll       = _get_nll()
+minimizer = AnealingMinimizer(ntries=10, chi2ndof=1.00)
+res       = minimizer.minimize(nll)
+```
+
 ## Fit plotting
 
 The class `ZFitPlotter` can be used to plot fits done with zfit. For a complete set of examples of how to use

{data_manipulation_utilities-0.1.7 → data_manipulation_utilities-0.1.9}/src/data_manipulation_utilities.egg-info/SOURCES.txt

@@ -22,6 +22,9 @@ src/dmu/rfile/rfprinter.py
 src/dmu/rfile/utilities.py
 src/dmu/stats/fitter.py
 src/dmu/stats/function.py
+src/dmu/stats/gof_calculator.py
+src/dmu/stats/minimizers.py
+src/dmu/stats/model_factory.py
 src/dmu/stats/utilities.py
 src/dmu/stats/zfit_plotter.py
 src/dmu/testing/utilities.py

{data_manipulation_utilities-0.1.7 → data_manipulation_utilities-0.1.9}/src/dmu/stats/fitter.py

@@ -4,6 +4,7 @@ Module holding zfitter class
 
 import pprint
 from typing                   import Union
+from functools                import lru_cache
 
 import numpy
 import zfit
@@ -100,8 +101,8 @@ class Fitter:
 
         return data
     #------------------------------
-    def _bin_pdf(self, nbins):
-        [[min_x]], [[max_x]] = self._pdf.space.limits
+    def _bin_pdf(self):
+        nbins, min_x, max_x = self._get_binning()
         _, arr_edg = numpy.histogram(self._data_np, bins = nbins, range=(min_x, max_x))
 
         size = arr_edg.size
@@ -117,23 +118,29 @@ class Fitter:
 
         return numpy.array(l_bc)
     #------------------------------
+    def _bin_data(self):
+        nbins, min_x, max_x = self._get_binning()
+        arr_data, _ = numpy.histogram(self._data_np, bins = nbins, range=(min_x, max_x))
+        arr_data    = arr_data.astype(float)
+
+        return arr_data
+    #------------------------------
+    @lru_cache(maxsize=10)
     def _get_binning(self):
         min_x = numpy.min(self._data_np)
         max_x = numpy.max(self._data_np)
         nbins = self._ndof + self._get_float_pars()
 
+        log.debug(f'Nbins: {nbins}')
+        log.debug(f'Range: [{min_x:.3f}, {max_x:.3f}]')
+
         return nbins, min_x, max_x
     #------------------------------
     def _calc_gof(self):
         log.debug('Calculating GOF')
-        nbins, min_x, max_x = self._get_binning()
 
-        log.debug(f'Nbins: {nbins}')
-        log.debug(f'Range: [{min_x:.3f}, {max_x:.3f}]')
-
-        arr_data, _ = numpy.histogram(self._data_np, bins = nbins, range=(min_x, max_x))
-        arr_data    = arr_data.astype(float)
-        arr_modl    = self._bin_pdf(nbins)
+        arr_data    = self._bin_data()
+        arr_modl    = self._bin_pdf()
         norm        = numpy.sum(arr_data) / numpy.sum(arr_modl)
         arr_modl    = norm * arr_modl
         arr_res     = arr_modl - arr_data

data_manipulation_utilities-0.1.9/src/dmu/stats/gof_calculator.py

@@ -0,0 +1,145 @@
+'''
+Module holding GofCalculator class
+'''
+from functools import lru_cache
+
+import zfit
+import numpy
+import pandas as pnd
+
+from scipy                  import stats
+from zfit.core.basepdf      import BasePDF   as zpdf
+from zfit.core.parameter    import Parameter as zpar
+from dmu.logging.log_store  import LogStore
+
+log = LogStore.add_logger('dmu:stats:gofcalculator')
+# ------------------------
+class GofCalculator:
+    '''
+    Class used to calculate goodness of fit from zfit NLL
+    '''
+    # ---------------------
+    def __init__(self, nll, ndof : int = 10):
+        self._nll     = nll
+        self._ndof    = ndof
+
+        self._pdf     = self._pdf_from_nll()
+        self._data_in = self._data_from_nll()
+        self._data_np = self._data_np_from_data(self._data_in)
+        self._data_zf = zfit.Data.from_numpy(obs=self._pdf.space, array=self._data_np)
+    # ---------------------
+    def _data_np_from_data(self, dat) -> numpy.ndarray:
+        if   isinstance(dat, numpy.ndarray):
+            return dat
+
+        if isinstance(dat, zfit.Data):
+            return zfit.run(zfit.z.unstack_x(dat))
+
+        if isinstance(dat, pnd.DataFrame):
+            return dat.to_numpy()
+
+        if isinstance(dat, pnd.Series):
+            dat    = pnd.DataFrame(dat)
+            return dat.to_numpy()
+
+        data_type = str(type(dat))
+        raise ValueError(f'Data is not a numpy array, zfit.Data or pandas.DataFrame, but {data_type}')
+    # ---------------------
+    def _pdf_from_nll(self) -> zpdf:
+        l_model = self._nll.model
+        if len(l_model) != 1:
+            raise ValueError('Not found one and only one model')
+
+        return l_model[0]
+    # ---------------------
+    def _data_from_nll(self) -> zpdf:
+        l_data = self._nll.data
+        if len(l_data) != 1:
+            raise ValueError('Not found one and only one dataset')
+
+        return l_data[0]
+    # ---------------------
+    def _get_float_pars(self) -> int:
+        npar     = 0
+        s_par    = self._pdf.get_params()
+        for par in s_par:
+            if par.floating:
+                npar+=1
+
+        return npar
+    # ---------------------
+    @lru_cache(maxsize=10)
+    def _get_binning(self) -> tuple[int, float, float]:
+        min_x = numpy.min(self._data_np)
+        max_x = numpy.max(self._data_np)
+        nbins = self._ndof + self._get_float_pars()
+
+        log.debug(f'Nbins: {nbins}')
+        log.debug(f'Range: [{min_x:.3f}, {max_x:.3f}]')
+
+        return nbins, min_x, max_x
+    # ---------------------
+    def _get_pdf_bin_contents(self) -> numpy.ndarray:
+        nbins, min_x, max_x  = self._get_binning()
+        _, arr_edg = numpy.histogram(self._data_np, bins = nbins, range=(min_x, max_x))
+
+        size = arr_edg.size
+
+        l_bc = []
+        for i_edg in range(size - 1):
+            low = arr_edg[i_edg + 0]
+            hig = arr_edg[i_edg + 1]
+
+            var : zpar = self._pdf.integrate(limits = [low, hig])
+            val = var.numpy()[0]
+            l_bc.append(val * self._data_np.size)
+
+        return numpy.array(l_bc)
+    #------------------------------
+    def _get_data_bin_contents(self) -> numpy.ndarray:
+        nbins, min_x, max_x = self._get_binning()
+        arr_data, _         = numpy.histogram(self._data_np, bins = nbins, range=(min_x, max_x))
+        arr_data            = arr_data.astype(float)
+
+        return arr_data
+    #------------------------------
+    @lru_cache(maxsize=30)
+    def _calculate_gof(self) -> tuple[float, int, float]:
+        log.debug('Calculating GOF')
+
+        arr_data    = self._get_data_bin_contents()
+        arr_modl    = self._get_pdf_bin_contents()
+
+        norm        = numpy.sum(arr_data) / numpy.sum(arr_modl)
+        arr_modl    = norm * arr_modl
+        arr_res     = arr_modl - arr_data
+
+        arr_chi2    = numpy.divide(arr_res ** 2, arr_data, out=numpy.zeros_like(arr_data), where=arr_data!=0)
+        sum_chi2    = numpy.sum(arr_chi2)
+
+        pvalue      = 1 - stats.chi2.cdf(sum_chi2, self._ndof)
+        pvalue      = float(pvalue)
+
+        log.debug(f'Chi2: {sum_chi2:.3f}')
+        log.debug(f'Ndof: {self._ndof}')
+        log.debug(f'pval: {pvalue:<.3e}')
+
+        return sum_chi2, self._ndof, pvalue
+    # ---------------------
+    def get_gof(self, kind : str) -> float:
+        '''
+        Returns good ness of fit of a given kind
+
+        kind: Type of goodness of fit, e.g. pvalue
+        '''
+
+        chi2, ndof, pval = self._calculate_gof()
+
+        if kind == 'pvalue':
+            return pval
+
+        if kind == 'chi2/ndof':
+            return chi2/ndof
+
+        raise NotImplementedError(f'Invalid goodness of fit: {kind}')
+# ------------------------

data_manipulation_utilities-0.1.9/src/dmu/stats/minimizers.py

@@ -0,0 +1,183 @@
+'''
+Module containing derived classes from ZFit minimizer
+'''
+import numpy
+
+import zfit
+from zfit.result                   import FitResult
+from zfit.core.basepdf             import BasePDF           as zpdf
+from zfit.minimizers.baseminimizer import FailMinimizeNaN
+from dmu.stats.gof_calculator      import GofCalculator
+from dmu.logging.log_store         import LogStore
+
+log = LogStore.add_logger('dmu:ml:minimizers')
+# ------------------------
+class AnealingMinimizer(zfit.minimize.Minuit):
+    '''
+    Class meant to minimizer zfit likelihoods by using multiple retries,
+    each retry is preceeded by the randomization of the fitting parameters
+    '''
+    # ------------------------
+    def __init__(self, ntries : int, pvalue : float = -1, chi2ndof : float = -1):
+        '''
+        ntries  : Try this number of times
+        pvalue  : Stop tries when this threshold is reached
+        chi2ndof: Use this value as a threshold to stop fits
+        '''
+        self._ntries   = ntries
+        self._pvalue   = pvalue
+        self._chi2ndof = chi2ndof
+
+        self._check_thresholds()
+
+        super().__init__()
+    # ------------------------
+    def _check_thresholds(self) -> None:
+        good_pvalue  = 0 <= self._pvalue < 1
+        good_chi2dof = self._chi2ndof > 0
+
+        if good_pvalue and good_chi2dof:
+            raise ValueError('Threshold for both chi2 and pvalue were specified')
+
+        if good_pvalue:
+            log.debug(f'Will use threshold on pvalue with value: {self._pvalue}')
+            return
+
+        if good_chi2dof:
+            log.debug(f'Will use threshold on chi2ndof with value: {self._chi2ndof}')
+            return
+
+        raise ValueError('Neither pvalue nor chi2 thresholds are valid')
+    # ------------------------
+    def _is_good_gof(self, ch2 : float, pvl : float) -> bool:
+        is_good_pval = pvl > self._pvalue   and self._pvalue   > 0
+        is_good_chi2 = ch2 < self._chi2ndof and self._chi2ndof > 0
+        is_good      = is_good_pval or is_good_chi2
+
+        if is_good_pval:
+            log.info(f'Stopping fit, found p-value: {pvl:.3f} > {self._pvalue:.3f}')
+
+        if is_good_chi2:
+            log.info(f'Stopping fit, found chi2/ndof: {ch2:.3f} > {self._chi2ndof:.3f}')
+
+        if not is_good:
+            log.debug(f'Could not read threshold, pvalue/chi2: {pvl:.3f}/{ch2:.3f}')
+
+        return is_good
+    # ------------------------
+    def _is_good_fit(self, res : FitResult) -> bool:
+        if not res.valid:
+            log.warning('Skipping invalid fit')
+            return False
+
+        if res.status != 0:
+            log.warning('Skipping fit with bad status')
+            return False
+
+        if not res.converged:
+            log.warning('Skipping non-converging fit')
+            return False
+
+        return True
+    # ------------------------
+    def _get_gof(self, nll) -> tuple[float, float]:
+        log.debug('Checking GOF')
+
+        gcl = GofCalculator(nll)
+        pvl = gcl.get_gof(kind='pvalue')
+        ch2 = gcl.get_gof(kind='chi2/ndof')
+
+        return ch2, pvl
+    # ------------------------
+    def _randomize_parameters(self, nll):
+        '''
+        Will move floating parameters of PDF according
+        to uniform PDF
+        '''
+
+        log.debug('Randomizing parameters')
+        l_model = nll.model
+        if len(l_model) != 1:
+            raise ValueError('Not found and and only one model')
+
+        model = l_model[0]
+        s_par = model.get_params(floating=True)
+        for par in s_par:
+            ival = par.value()
+            fval = numpy.random.uniform(par.lower, par.upper)
+            par.set_value(fval)
+            log.debug(f'{par.name:<20}{ival:<15.3f}{"->":<10}{fval:<15.3f}{"in":<5}{par.lower:<15.3e}{par.upper:<15.3e}')
+    # ------------------------
+    def _pick_best_fit(self, d_chi2_res : dict) -> FitResult:
+        nres = len(d_chi2_res)
+        if nres == 0:
+            raise ValueError('No fits found')
+
+        l_chi2_res= list(d_chi2_res.items())
+        l_chi2_res.sort()
+        chi2, res = l_chi2_res[0]
+
+        log.warning(f'Picking out best fit from {nres} fits with chi2: {chi2:.3f}')
+
+        return res
+    #------------------------------
+    def _set_pdf_pars(self, res : FitResult, pdf : zpdf) -> None:
+        '''
+        Will set the PDF floating parameter values as the result instance
+        '''
+        l_par_flt = list(pdf.get_params(floating= True))
+        l_par_fix = list(pdf.get_params(floating=False))
+        l_par     = l_par_flt + l_par_fix
+
+        d_val = { par.name : dc['value'] for par, dc in res.params.items()}
+
+        log.debug('Setting PDF parameters to best result')
+        for par in l_par:
+            if par.name not in d_val:
+                par_val = par.value().numpy()
+                log.debug(f'Skipping {par.name} = {par_val:.3e}')
+                continue
+
+            val = d_val[par.name]
+            log.debug(f'{"":<4}{par.name:<20}{"->":<10}{val:<20.3e}')
+            par.set_value(val)
+    # ------------------------
+    def _pdf_from_nll(self, nll) -> zpdf:
+        l_model = nll.model
+        if len(l_model) != 1:
+            raise ValueError('Cannot extract one and only one PDF from NLL')
+
+        return l_model[0]
+    # ------------------------
+    def minimize(self, nll, **kwargs) -> FitResult:
+        '''
+        Will run minimization and return FitResult object
+        '''
+
+        d_chi2_res : dict[float,FitResult] = {}
+        for i_try in range(self._ntries):
+            log.info(f'try {i_try:02}/{self._ntries:02}')
+            try:
+                res = super().minimize(nll, **kwargs)
+            except (FailMinimizeNaN, ValueError, RuntimeError) as exc:
+                log.warning(exc)
+                self._randomize_parameters(nll)
+                continue
+
+            if not self._is_good_fit(res):
+                continue
+
+            chi2, pvl = self._get_gof(nll)
+            d_chi2_res[chi2] = res
+
+            if self._is_good_gof(chi2, pvl):
+                return res
+
+            self._randomize_parameters(nll)
+
+        res = self._pick_best_fit(d_chi2_res)
+        pdf = self._pdf_from_nll(nll)
+        self._set_pdf_pars(res, pdf)
+
+        return res
+# ------------------------

data_manipulation_utilities-0.1.9/src/dmu/stats/model_factory.py

@@ -0,0 +1,207 @@
+'''
+Module storing ZModel class
+'''
+# pylint: disable=too-many-lines, import-error
+
+from typing import Callable, Union
+
+import zfit
+from zfit.core.interfaces   import ZfitSpace as zobs
+from zfit.core.basepdf      import BasePDF   as zpdf
+from zfit.core.parameter    import Parameter as zpar
+from dmu.logging.log_store  import LogStore
+
+log=LogStore.add_logger('dmu:stats:model_factory')
+#-----------------------------------------
+class MethodRegistry:
+    '''
+    Class intended to store protected methods belonging to ModelFactory class
+    which is defined in this same module
+    '''
+    # Registry dictionary to hold methods
+    _d_method = {}
+
+    @classmethod
+    def register(cls, nickname : str):
+        '''
+        Decorator in charge of registering method for given nickname
+        '''
+        def decorator(method):
+            cls._d_method[nickname] = method
+            return method
+
+        return decorator
+
+    @classmethod
+    def get_method(cls, nickname : str) -> Union[Callable,None]:
+        '''
+        Will return method in charge of building PDF, for an input nickname
+        '''
+        return cls._d_method.get(nickname, None)
+#-----------------------------------------
+class ModelFactory:
+    '''
+    Class used to create Zfit PDFs by passing only the nicknames, e.g.:
+
+    ```python
+    from dmu.stats.model_factory import ModelFactory
+
+    l_pdf = ['dscb', 'gauss']
+    l_shr = ['mu']
+    mod   = ModelFactory(obs = obs, l_pdf = l_pdf, l_shared=l_shr)
+    pdf   = mod.get_pdf()
+    ```
+
+    where one can specify which parameters can be shared among the PDFs
+    '''
+    #-----------------------------------------
+    def __init__(self, obs : zobs, l_pdf : list[str], l_shared : list[str]):
+        '''
+        obs:      zfit obserbable
+        l_pdf:    List of PDF nicknames which are registered below
+        l_shared: List of parameter names that are shared
+        '''
+
+        self._l_pdf           = l_pdf
+        self._l_shr           = l_shared
+        self._l_can_be_shared = ['mu', 'sg']
+        self._obs             = obs
+
+        self._d_par : dict[str,zpar] = {}
+    #-----------------------------------------
+    def _get_name(self, name : str, suffix : str) -> str:
+        for can_be_shared in self._l_can_be_shared:
+            if name.startswith(f'{can_be_shared}_') and can_be_shared in self._l_shr:
+                return can_be_shared
+
+        return f'{name}{suffix}'
+    #-----------------------------------------
+    def _get_parameter(self,
+                       name   : str,
+                       suffix : str,
+                       val    : float,
+                       low    : float,
+                       high   : float) -> zpar:
+        name = self._get_name(name, suffix)
+        if name in self._d_par:
+            return self._d_par[name]
+
+        par  = zfit.param.Parameter(name, val, low, high)
+
+        self._d_par[name] = par
+
+        return par
+    #-----------------------------------------
+    @MethodRegistry.register('exp')
+    def _get_exponential(self, suffix : str = '') -> zpdf:
+        c   = self._get_parameter('c_exp', suffix, -0.005, -0.05, 0.00)
+        pdf = zfit.pdf.Exponential(c, self._obs)
+
+        return pdf
+    #-----------------------------------------
+    @MethodRegistry.register('pol1')
+    def _get_pol1(self, suffix : str = '') -> zpdf:
+        a   = self._get_parameter('a_pol1', suffix, -0.005, -0.95, 0.00)
+        pdf = zfit.pdf.Chebyshev(obs=self._obs, coeffs=[a])
+
+        return pdf
+    #-----------------------------------------
+    @MethodRegistry.register('pol2')
+    def _get_pol2(self, suffix : str = '') -> zpdf:
+        a   = self._get_parameter('a_pol2', suffix, -0.005, -0.95, 0.00)
+        b   = self._get_parameter('b_pol2', suffix,  0.000, -0.95, 0.95)
+        pdf = zfit.pdf.Chebyshev(obs=self._obs, coeffs=[a, b])
+
+        return pdf
+    #-----------------------------------------
+    @MethodRegistry.register('cbr')
+    def _get_cbr(self, suffix : str = '') -> zpdf:
+        mu  = self._get_parameter('mu_cbr', suffix, 5300, 5250, 5350)
+        sg  = self._get_parameter('sg_cbr', suffix,   10,    2,  300)
+        ar  = self._get_parameter('ac_cbr', suffix,   -2,  -4.,  -1.)
+        nr  = self._get_parameter('nc_cbr', suffix,    1,  0.5,  5.0)
+
+        pdf = zfit.pdf.CrystalBall(mu, sg, ar, nr, self._obs)
+
+        return pdf
+    #-----------------------------------------
+    @MethodRegistry.register('cbl')
+    def _get_cbl(self, suffix : str = '') -> zpdf:
+        mu  = self._get_parameter('mu_cbl', suffix, 5300, 5250, 5350)
+        sg  = self._get_parameter('sg_cbl', suffix,   10,    2,  300)
+        al  = self._get_parameter('ac_cbl', suffix,    2,   1.,   4.)
+        nl  = self._get_parameter('nc_cbl', suffix,    1,  0.5,  5.0)
+
+        pdf = zfit.pdf.CrystalBall(mu, sg, al, nl, self._obs)
+
+        return pdf
+    #-----------------------------------------
+    @MethodRegistry.register('gauss')
+    def _get_gauss(self, suffix : str = '') -> zpdf:
+        mu  = self._get_parameter('mu_gauss', suffix, 5300, 5250, 5350)
+        sg  = self._get_parameter('sg_gauss', suffix,   10,    2,  300)
+
+        pdf = zfit.pdf.Gauss(mu, sg, self._obs)
+
+        return pdf
+    #-----------------------------------------
+    @MethodRegistry.register('dscb')
+    def _get_dscb(self, suffix : str = '') -> zpdf:
+        mu  = self._get_parameter('mu_dscb', suffix, 5300, 5250, 5400)
+        sg  = self._get_parameter('sg_dscb', suffix,   10,    2,   30)
+        ar  = self._get_parameter('ar_dscb', suffix,    1,    0,    5)
+        al  = self._get_parameter('al_dscb', suffix,    1,    0,    5)
+        nr  = self._get_parameter('nr_dscb', suffix,    2,    1,    5)
+        nl  = self._get_parameter('nl_dscb', suffix,    2,    0,    5)
+
+        pdf = zfit.pdf.DoubleCB(mu, sg, al, nl, ar, nr, self._obs)
+
+        return pdf
+    #-----------------------------------------
+    def _get_pdf_types(self) -> list[tuple[str,str]]:
+        d_name_freq = {}
+
+        l_type = []
+        for name in self._l_pdf:
+            if name not in d_name_freq:
+                d_name_freq[name] = 1
+            else:
+                d_name_freq[name]+= 1
+
+            frq = d_name_freq[name]
+            frq = f'_{frq}'
+
+            l_type.append((name, frq))
+
+        return l_type
+    #-----------------------------------------
+    def _get_pdf(self, kind : str, preffix : str) -> zpdf:
+        fun = MethodRegistry.get_method(kind)
+        if fun is None:
+            raise NotImplementedError(f'PDF of type {kind} is not implemented')
+
+        return fun(self, preffix)
+    #-----------------------------------------
+    def _add_pdf(self, l_pdf : list[zpdf]) -> zpdf:
+        nfrc = len(l_pdf)
+        if nfrc == 1:
+            log.debug('Requested only one PDF, skipping sum')
+            return l_pdf[0]
+
+        l_frc= [ zfit.param.Parameter(f'frc_{ifrc + 1}', 0.5, 0, 1) for ifrc in range(nfrc - 1) ]
+
+        pdf = zfit.pdf.SumPDF(l_pdf, fracs=l_frc)
+
+        return pdf
+    #-----------------------------------------
+    def get_pdf(self) -> zpdf:
+        '''
+        Given a list of strings representing PDFs returns the a zfit PDF which is
+        the sum of them
+        '''
+        l_type=   self._get_pdf_types()
+        l_pdf = [ self._get_pdf(kind, preffix) for kind, preffix in l_type ]
+        pdf   =   self._add_pdf(l_pdf)
+
+        return pdf
+#-----------------------------------------