ROOT(a Data analysis Framework)-Note6: iSTEP day3-Random

最新推荐文章于 2023-01-04 18:57:41 发布
最新推荐文章于 2023-01-04 18:57:41 发布
阅读量486 收藏 5
点赞数 3
分类专栏： ROOT 文章标签： ROOT  CERN
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/weixin_43316938/article/details/96467969
版权
本文介绍了使用ROOT框架中的TRandom3生成随机数，并通过Monte-Carlo方法模拟特定分布。内容包括：1) 直接生成随机数并保留满足条件的数；2) 使用函数积分的反解方法；3) 应用Monte-Carlo方法模拟Gaussian分布。同时，讨论了直方图与PDF对比的问题，并展示了如何拟合数据。
摘要由CSDN通过智能技术生成

Note5: Random3

其实每天进展都还是听慢的，今天也只考虑了一个小问题，就是生成随机数。

Monte-Carlo方法产生特定分布的随机数

给定一个分布，或许不一定归一化：

P ( x = x ′ ) = f ( x ) P(x = x') = f(x) P(x=x′)=f(x)

方法一：

  • 生成一个平均分布的随机数 x ∈ [ a 1 , a 2 ] x \in [a_1,a_2] x∈[a1​,a2​]
  • 生成一个平均分布的随机数 y ∈ [ b 1 , b 2 ] y \in [b_1,b_2] y∈[b1​,b2​]
  • 如果 y < = P ( x = x
最低0.47元/天 解锁文章
X D
关注
  • 3
    点赞
  • 踩
  • 5
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
ROOT (a Data analysis Framework) TMatrix类的简单操作
Epsilom的博客
09-08 897
简单使用TMatrix类 ROOT中的TMatrix类可以实现矩阵运算。详细介绍见ROOTUsersGuide-14 Linear Algebra in ROOT 一个内含Double_t元素的矩阵可以这样定义： TMatrixD x = TMatrix(nRow, nCol,) ...
ROOT (a Data analysis Framework) - Note2：关于编译和运行
Epsilom的博客
05-25 1070
Note2：关于编译运行 写在前面 ROOT是基于C++的，即使有了cling的帮忙，也还是不能像python那样方便地交互运行。怎么说呢，就比如你写了一个脚本，在交互式界面中直接.x运行，ROOT还是把它当做一个C++脚本去看待的，诸如.L、.help之类的命令是不能出现在脚本中的。这样如果你自己写了一个类，在主函数中调用它之前，一般是include这个类的头文件，但是这在ROOT上并不好用，R...
ROOT(a Data analysis Framework)-Note4: iSTEP day1-TH1&TF1
Epsilom的博客
07-16 1123
ROOT(a Data analysis Framework)-Note4: iSTEP day1 最近参加了一个叫做iSTEP的暑期学校，这里教授一些基本的ROOT使用方法。曾经用root的过程忘的差不多了，也没有记录。现在正好记录一下。 TF1类 TF1类是一个数学函数的类，基本的用法就是画一些函数的图像。 root [0] TF1* f = new TF1("f","sin(x)",-10,...
SFrame - A ROOT data analysis framework-开源
04-28
SFrame是围绕ROOT库构建的C ++框架，用于分析粒子物理数据。 通过允许用户在分布式计算机场上运行他/她的代码，它为处理数据提供了非常高的性能。
ROOT (a Data analysis Framework) - Note1：先定一个小目标
Epsilom的博客
05-24 1881
ROOT-Note ROOT是一个基于C++的软件包，它被广泛地应用于高能物理实验的数据处理中。为了应付毕设， 近期开始学习ROOT的基本操作方法。 打开ROOT Root软件是基于C++的，但它嵌入了cling，因此可以通过命令行交互，还可以运行macro文件（实际上是一个C++脚本），可以解释运行也可以编译运行，这就很方便。 在设置好环境变量的前提下，打开ROOT交互界面可以输入命令： roo...
Pipr:.NET的简单责任链
04-13
用法简单的首先，创建一些实现IStep 的步骤。 public class ToStringStep : IStep<int>{ public string Execute(int input, PipelineContext context) { return input.ToString(); }}public class DoublerStep : ...
OLPC English Literacy Tool:OLPC软件提高英语素养-开源
07-19
该软件最初是作为乌拉圭 iSTEP 2011 实习计划的一部分创建的。 有关 iSTEP 2011 实习的更多信息，请访问 http://istep2011.techbridgeworld.org。 如果您对此版本或实习有任何疑问，请联系或访问 ...
Content Authoring Tool:OLPC 英语读写工具的内容创作工具-开源
06-08
该软件最初是作为乌拉圭 iSTEP 2011 实习计划的一部分创建的。 有关 iSTEP 2011 实习的更多信息，请访问 http://istep2011.techbridgeworld.org。 如果您对此版本或实习有任何疑问，请联系或访问 ...
iStep：用于移动传感器网络中扩散分析的可感知步骤的采样方法
02-25
iStep：用于移动传感器网络中扩散分析的可感知步骤的采样方法
matlab集成c代码-Streaming-DMD-integrated-within-NEK5000:流式DMD集成在NEK5000内
05-22
A，Qx，Qy和Gx矩阵 istep_n-以压缩数组的形式存储在时间步长n处每个点的uvw速度数据：[u1 u2 u3。 。 u_linterpol v1 v2 v3。 。 v_interpol w1 w2 w3。 。 w_interpol] istep文件在开始后以1000个时间步的时间间隔...
ROOT (a Data analysis Framework) - Note3：TCanvas类
Epsilom的博客
05-24 1298
TCanvas类 TCanvas是一个画布的类，任何一张图，无论是直方图、函数曲线还是其他图，最终显示，都需要在这样一个画布上完成。实际上，如上一节那样，当你没有建立TCanvas实例就使用h.Draw()这个方法时，系统总会默认给这个Draw()方法创建一个画布，并且默认给一个名字叫c1。 Info in <TCanvas::MakeDefCanvas>: created defa...
ROOT(a Data analysis Framework)-Note5: iSTEP day2-TH1::Fit()
Epsilom的博客
07-18 1250
Note5 ：数据拟合 今天还是关于TH1和TF1的内容，使用Fit方法，使用不同的函数对数据进行拟合。下面以问题的形式做个小的总结记录。 Section 1 目标： Read histogram “\h1” from file “\hists.root” that I give you. And Draw it. 读取一个.root文件： 基于TFile类,有这样几种方法： TFile f(...
【ROOT from CERN】——TH1绘制一维直方图
qq_20604319的博客
01-04 2441
TH1绘制一维直方图。
covariance matrix r语言_NONMEM中$COVARIANCE失败时的救济方法，谈一下参数uncertainty估算问题...
weixin_39811101的博客
11-30 358
最近我的好友朱校博士在他的wechat公主号（名称：PKPD，里面内容非常丰富）写了一篇关于NONMEM $COV中R矩阵的文章，我顺手在文章中给他加了一些其他的估算方法，把其中其他方法的部分搬过来与大家讨论一下。$COVARIANCE在NONMEM中意义很重要，当你完成建模后估算出参数值，这个参数值估的准不准肯定是modeler非常关心的部分，而这也是$COVARIANCE在NONMEM里的功能...
【ROOT from CERN】——传统的约定
qq_20604319的博客
08-21 563
介绍ROOT的一些约定的规则。
erf和erfc函数
热门推荐
格物致知的专栏 [音视频编解码 网络协议 计算机视觉 计算机图形学 图像理解 语音识别 机器学习 模拟电路 传感器]
12-07 3万+
erf即error function，erfc中的c是complementary的缩写，所以有erfc = 1 - erf. matlab里的定义描述：erf(x) = 2/sqrt(pi) * integral from 0 to x of exp(-t^2) dt. 其值域范围为(-1, 1)，区间长度为2.
math.erf- error function- 正太分布的CDF
weixin_45405128的博客
01-15 937
def cdf(std_dev, avg, x): inner = 1 + math.erf((x - avg) / (std_dev * math.sqrt(2))) return 0.5 * inner
ROOT学习——一维直方图的条形图示例(hbars.C)
M_M-铭的博客
04-17 1223
2021416
ROOT学习——高斯随机生成二维直方图的两种画法(candledecay.C)
M_M-铭的博客
04-27 1326
官方介绍：蜡烛衰变，说明了一定时间的发展。 下面给出具体代码步骤： 创建一个画板： auto c1 = new TCanvas("c1","Candle Decay",800,600); 将画板竖切为两半： c1->Divide(2,1); 新建一个随机变量： auto rng = new TRandom(); 新建一个二维直方图： auto h1 = new ...
%% 相关性分析 figure data2_limit_steady(:,end+1) = (data2_limit_steady(:,2) - data2_limit_steady(:,3))./data2_limit_steady(:,2)*100; n_data = length(data2_limit_steady(:,1)); n_cluster = 12; Q=zeros(n_cluster,1); step = floor(n_data/n_cluster); for i = 1:n_cluster subplot(n_cluster/2,2,i); index = find((data2_limit_steady((i-1)step+1:istep,10)==1) & (data2_limit_steady((i-1)step+1:istep,1)>480)); row = (i-1)step + index; plot(data2_limit_steady(row,6),data2_limit_steady(row,end),'s'); axis([0,500,60,100]); xlabel('喷氨量'); ylabel('效率'); title(num2str(i)); Q(i)=sum(data2_limit_steady(row,1).(data2_limit_steady(row,2)-data2_limit_steady(row,3)))./sum(data2_limit_steady(row,6));修改代码，使负荷大于480小于510
最新发布
05-29
以下是修改后的代码： ``` figure data2_limit_steady(:,end+1) = (data2_limit_steady(:,2) - data2_limit_steady(:,3))./data2_limit_steady(:,2)*100; n_data = length(data2_limit_steady(:,1)); n_cluster = 12; Q=zeros(n_cluster,1); step = floor(n_data/n_cluster); for i = 1:n_cluster subplot(n_cluster/2,2,i); index = find((data2_limit_steady((i-1)*step+1:i*step,10)==1) & (data2_limit_steady((i-1)*step+1:i*step,1)>480) & (data2_limit_steady((i-1)*step+1:i*step,1)<510)); row = (i-1)*step + index; plot(data2_limit_steady(row,6),data2_limit_steady(row,end),'s'); axis([0,500,60,100]); xlabel('喷氨量'); ylabel('效率'); title(num2str(i)); Q(i)=sum(data2_limit_steady(row,1).*(data2_limit_steady(row,2)-data2_limit_steady(row,3)))./sum(data2_limit_steady(row,6)); end ``` 主要的修改有： 1. 在第9行，将 `istep` 修改为 `i*step`。 2. 在第10行，增加一个条件 `& (data2_limit_steady((i-1)*step+1:i*step,1)<510)`，限制喷氨量小于510。 3. 在第11行，将 `(i-1)step` 修改为 `(i-1)*step`。 4. 在第13行，将 `data2_limit_steady(row,1).(data2_limit_steady(row,2)-data2_limit_steady(row,3))` 修改为 `data2_limit_steady(row,1).*(data2_limit_steady(row,2)-data2_limit_steady(row,3))`。
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值