BCMath.h

#ifndef __BCMATH__H
#define __BCMATH__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 <stdexcept>
#include <vector>

class TRandom;

namespace BCMath
{

double LogGaus(double x, double mean = 0, double sigma = 1, bool norm = false);

double LogSplitGaus(double x, double mode, double sigma_below, double sigma_above, bool norm = false);

double LogPoisson(double x, double lambda);

double LogApproxBinomial(unsigned n, unsigned k, double p);

double LogBreitWignerNonRel(double x, double mean, double Gamma, bool norm = false);

double LogBreitWignerRel(double x, double mean, double Gamma);

double LogChi2(double x, int n);

double LogVoigtian(double x, double sigma, double gamma);

double LogGammaPDF(double x, double alpha, double beta);

double LogLogNormal(double x, double mean = 0, double sigma = 1);

double ApproxBinomial(unsigned n, unsigned k, double p);

double LogBinomFactor(unsigned n, unsigned k);

double ApproxLogFact(double x);

double LogBinomFactorExact(unsigned n, unsigned k);

double LogFact(unsigned n);

unsigned CacheFactorials(unsigned n);

int Nint(double x);

double CorrectPValue(const double& pvalue, const unsigned& npar, const unsigned& nobservations);

double FastPValue(const std::vector<unsigned>& observed, const std::vector<double>& expected,
                  unsigned nIterations = 1e5, unsigned seed = 0);

namespace Random
{
double Chi2(TRandom* rng, double dof);

double Gamma(TRandom* rng, double a, double b);
}
}

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

#endif