// A program to test generating multinomial random numbers

// Glen Cowan
// RHUL Physics
// January 2010

#include <iostream>
#include <fstream>
#include <cstdlib>
#include <cmath>
#include <string>
#include <sstream>
#include <vector>

#include <Math/Random.h>
#include <Math/GSLRndmEngines.h>

using namespace std;

int main() {

  unsigned int seed = 12345;
  ROOT::Math::Random<ROOT::Math::GSLRngGFSR4>*
    ran = new ROOT::Math::Random<ROOT::Math::GSLRngGFSR4>(seed);

  int nBins = 10;
  vector<double> p(nBins);
  for (int i=0; i<nBins; i++){
    p[i] = 1. / static_cast<double>(nBins);
  }

  unsigned int nTot = 100;
  vector<unsigned int> multDist(nBins);
  multDist = ran->Multinomial(nTot, p);

  for (int i=0; i<nBins; i++){
    cout << i << "  " << multDist[i] << endl;
  }

  return 0;

}