BCMath.h

Go to the documentation of this file.

#ifndef __BCMATH__H
#define __BCMATH__H

/*
 * Copyright (C) 2008, Daniel Kollar and Kevin Kroeninger.
 * All rights reserved.
 *
 * For the licensing terms see doc/COPYING.
 */

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

#define BCMATH_NFACT_ALIMIT 20

// ---------------------------------------------------------
//#include <cstring>
#include <vector>

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

   inline int Max(int x, int y)
      { return x >= y ? x : y; }

   inline double Max(double x, double y)
      { return x >= y ? x : y; }

   inline int Min(int x, int y)
      { return x <= y ? x : y; }

   inline double Min(double x, double y)
   { return x <= y ? x : y; }

   int Nint(double x);

   double rms(int n, const double * a);

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

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

// /**
//  * Fill vector with N random numbers distributed according to
//  * exponential distribution with rate lambda
//  * P(x|lambda)=lambda exp(-lambda*x)
//  * @param randoms will be filled with random exponentials
//  * @param lambda the rate of the exp. distribution
//  */
// void RandomExponential(std::vector<double> &randoms, double lambda, unsigned int seed=0);

}

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

#endif

---