ポアソン分布で乱数生成したい（TRandomMixMax::Poisson）

1#include <TRandomMixMax.h>
2// 乱数生成器を初期化
3TRandomMixMax rng;
4// 平均5のポアソン分布に従う乱数を生成
5double randomValue = rng.Poisson(5);

TRandomMixMax::Poissonでポアソン分布にしたがう乱数を生成できます。 引数に平均値を指定することで、任意のポアソン分布の乱数を取得できます。

ポアソン分布のヒストグラムを作りたい

 1#include <TCanvas.h>
 2#include <TH1D.h>
 3#include <TRandomMixMax.h>
 4
 5TCanvas *c = new TCanvas("c", "Poisson Distribution", 800, 600);
 6TH1D *h = new TH1D("h", "Poisson Distribution", 100, 0, 20);
 7
 8TRandomMixMax rng;
 9for (int i = 0; i < 10000; i++) {
10    h->Fill(rng.Poisson(5));
11}
12
13h->Draw();

このサンプルでは、平均5のポアソン分布にしたがう乱数を10000個生成し、その分布をヒストグラムで可視化しています。