如何根据概率密度函数生成随机分布

问题描述

根据 y=cos(theta)概率函数生成随机抽样

解决办法

对概率密度函数积分归一化得到概率函数，然后采用概率函数反函数生成随机分布；

代码实现

void GenerateDirection()
{
  TRandom3 rndm;
  
  TH1D *h1 = new TH1D("","",400,0.,3.1415/2);
  for(int i=0;i<100000;i++)
    {
      double y=rndm.Uniform();
      double Theta=asin(y);
      h1->Fill(Theta);
    }
  h1->Draw();
}

结果