TSFunc1DMarcoPdf.h

/*
 * TSFunc1DMarcoPdf.h
 *
 *  Created on: Aug 23, 2016
 *      Author: Maciej Lewicki
 *              malewick@cern.ch
 *
 */

#ifndef TSFUNC1DMARCOPDF_H_
#define TSFUNC1DMARCOPDF_H_


#include "TSFunc1DGaussPdf.h"

/* 
 * f(x) = 
 *  = \frac{1}{\sum_l n_l}\sum_l \frac{n_l}{\sqrt{2\pi}\sigma_{i,l}}\exp{-\frac{1}{2}\left( \frac{x-x_i}{(1 \pm \delta)\sigma_{i,l}} \right)^2}
 * , where:
 * \sigma_{i,l} = \frac{\sigma_0}{\sqrt{l}} \left( \frac{x_i}{x_1} \right)^{\alpha}
 */


class TSFunc1DMarcoPdf: public TSFunc1DGaussPdf {


public:

    TSFunc1DMarcoPdf() {};
    TSFunc1DMarcoPdf(TString name, TString title
                    ,const TSVariable& var_x, const TSParameter& mean, const TSParameter& sigma
                    ,const TSParameter& delta
                    ,const TSParameter& beta
                    ,std::map<int,int>& npoints );

    TSFunc1DMarcoPdf(TString name, TString title
                    ,const TSVariable& var_x, const TSParameter& mean, const TSParameter& sigma
                    ,const TSParameter& delta
                    ,std::map<int,int>& npoints );


    virtual ~TSFunc1DMarcoPdf();

    virtual double Evaluate() const;

    double GetDelta() const { return fDelta ? fDelta->GetValue(): 0; }
    double GetBeta() const { return fBeta ? fBeta->GetValue(): 0; }

    const TSParameter *GetParamDelta() const {return fDelta;}
    const TSParameter *GetParamBeta() const {return fBeta;}

private:

    const TSParameter *fDelta;
    const TSParameter *fBeta;
    std::map<int,int> fNPoints;

};

#endif /* TSFUNC1DMARCOPDF_H_ */