BCMath.h

#ifndef __BCMATH__H
#define __BCMATH__H

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

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

#define BCMATH_NFACT_ALIMIT 20

// ---------------------------------------------------------
#include <stdexcept>
#include <vector>

class TH1D;

namespace BCMath
{

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

   double LogPoisson(double x, double par);

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

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

   double LogBinomFactor(int n, int k);

   double ApproxLogFact(double x);

   double LogNoverK(int n, int k);

   double LogFact(int n);

   int Nint(double x);

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

   void RandomChi2(std::vector<double> &randoms, int K);

   TH1D* ECDF(const std::vector<double>& data);

   std::vector<int> longestRuns(const std::vector<bool>& bitStream);

   std::vector<double> longestRunsChi2(const std::vector<double>& yMeasured,
      const std::vector<double>& yExpected, const std::vector<double>& sigma);

   double longestRunFrequency(unsigned int longestObserved, unsigned int nTrials);

   double SplitGaussian(double* x, double* par);

   void CacheFactorial(unsigned int n);

    double Rvalue(const std::vector<double> & chain_means, const std::vector<double> & chain_variances,
                  const unsigned & chain_length, const bool & strict = true) throw (std::invalid_argument, std::domain_error);

    double CorrectPValue(const double & pvalue, const unsigned & npar, const unsigned & nobservations) throw (std::domain_error);

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

}

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

#endif