TSFunc1DMaciekPdf.h

/*
 * TSFunc1DMaciekPdf.h
 *
 *  Created on: Dec 4, 2018
 *      Author: Maciej Lewicki
 *      malewick@cern.ch
 *
 */

#ifndef TSFUNC1DMACIEKPDF_H_
#define TSFUNC1DMACIEKPDF_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 TSFunc1DMaciekPdf: public TSFunc1DGaussPdf {


public:

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


    virtual ~TSFunc1DMaciekPdf();

    virtual double Evaluate() const;

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

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

private:

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

};

#endif /* TSFUNC1DMACIEKPDF_H_ */