BCSplitGaussianPrior.h

#ifndef __BCSPLITGAUSSIANPRIOR__H
#define __BCSPLITGAUSSIANPRIOR__H

/*
 * Copyright (C) 2007-2018, the BAT core developer team
 * All rights reserved.
 *
 * For the licensing terms see doc/COPYING.
 * For documentation see http://mpp.mpg.de/bat
 */

// ---------------------------------------------------------

#include "BCPrior.h"

#include <limits>

class TF1;

// ---------------------------------------------------------

class BCSplitGaussianPrior : public BCPrior
{
public:
    BCSplitGaussianPrior(double mode, double sigma_below, double sigma_above);

    BCSplitGaussianPrior(const BCSplitGaussianPrior& other);

    virtual ~BCSplitGaussianPrior() {};

    virtual BCPrior* Clone() const
    { return new BCSplitGaussianPrior(*this); }

    virtual bool IsValid() const
    {
        return std::isfinite(fMode)
               and std::isfinite(fSigmaBelow) and fSigmaBelow > 0
               and std::isfinite(fSigmaAbove) and fSigmaAbove > 0;
    }

    virtual double GetLogPrior(double x);

    virtual double GetMode(double xmin = -std::numeric_limits<double>::infinity(), double xmax = std::numeric_limits<double>::infinity())
    { if (fMode < xmin) return xmin; if (fMode > xmax) return xmax; return fMode; }

    virtual double GetRawMoment(unsigned n, double xmin = -std::numeric_limits<double>::infinity(), double xmax = std::numeric_limits<double>::infinity());

    virtual double GetIntegral(double xmin = -std::numeric_limits<double>::infinity(), double xmax = std::numeric_limits<double>::infinity());

    void SetMean(double mean)
    { SetMode(mean); }

    void SetMode(double mode)
    { fMode = mode; }

    void SetSigmaBelow(double sigma)
    { fSigmaBelow = sigma; }

    void SetSigmaAbove(double sigma)
    { fSigmaAbove = sigma; }

    void SetParameters(double mode, double sigma_below, double sigma_above)
    { SetMode(mode); SetSigmaBelow(sigma_below); SetSigmaAbove(sigma_above);}

protected:
    double fMode;                            
    double fSigmaBelow;              
    double fSigmaAbove;              
};

#endif