|
采用CERN ROOT工具,可用以下命令:
const double PI=3.141592654;
theta=acos(2.*gRandom->Uniform(0,1)-1);//散射角各项同性
phi=gRandom->Uniform(0,2*PI);//方位角各项同性
注意不是:
theta=gRandom->Uniform(0,PI);//散射角各项同性
phi=gRandom->Uniform(0,2*PI);//方位角各项同性
因为:
dΩ=sinθdθdφ里还有个sinθ
参考资料
https://mathworld.wolfram.com/SpherePointPicking.html
Archiver|手机版|科学网 ( 京ICP备07017567号-12 )
GMT+8, 2024-11-16 23:31
Powered by ScienceNet.cn
Copyright © 2007- 中国科学报社