BCTH1Prior.h

#ifndef __BCTH1PRIOR__H
#define __BCTH1PRIOR__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 <TH1.h>
#include <limits>

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

class BCTH1Prior : public BCPrior
{
public:
    BCTH1Prior(TH1& h, bool interpolate = false);

    BCTH1Prior(TH1* h, bool interpolate = false);

    BCTH1Prior(const BCTH1Prior& other);

    virtual ~BCTH1Prior();

    BCTH1Prior& operator=(BCTH1Prior rhs);

    friend void swap(BCTH1Prior& A, BCTH1Prior& B);

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

    virtual bool IsValid() const;

    virtual double GetPrior(double x, bool normalize = false)
    { return ((fInterpolate) ? fPriorHistogram->Interpolate(x) : fPriorHistogram->GetBinContent(fPriorHistogram->FindFixBin(x))) * ((normalize) ? exp(-fLogIntegral) : 1); }

    virtual double GetLogPrior(double x);

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

    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());

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

    virtual double GetVariance(double xmin = -std::numeric_limits<double>::infinity(), double xmax = std::numeric_limits<double>::infinity())
    { double s = GetStandardDeviation(xmin, xmax); return s * s; }

    virtual double GetStandardDeviation(double xmin = -std::numeric_limits<double>::infinity(), double xmax = std::numeric_limits<double>::infinity())
    { (void)xmin; (void)xmax; return fPriorHistogram->GetRMS(); }

    virtual BCH1D GetBCH1D(TH1* bins, const std::string& name = "prior");

    virtual void SetInterpolate(bool interpolate)
    { fInterpolate = interpolate; }

    virtual double GetRandomValue(double xmin, double xmax, TRandom* const R = NULL)
    {
        (void) xmin;
        (void) xmax;
        (void) R;
        return fPriorHistogram->GetRandom();
    }

    virtual TH1& GetHistogram()
    { return *fPriorHistogram; }

    virtual const TH1& GetHistogram() const
    { return *fPriorHistogram; }

    virtual bool GetInterpolate()
    { return fInterpolate; }

    void NormalizeHistogram();

protected:

    // We don't accept nullptr and used a reference up to bat 1.0-rc1
    // but unfortunately, TH1& operator=(const TH1&) is declared private
    // at least up root 5.34/30 so we cannot change it in the swap function, hence we need to use a pointer
    TH1* fPriorHistogram;  

    bool fInterpolate; 
};

#endif