实际作图中,可能想要观察Entry前后分布的数值情况,可是使用动态绘制直方图的情况来了解,通过动态绘制直方图,能够大概知道前后Entry是随机分布,还是遵循何种规律。
下面介绍具体步骤:
- 创建一块画布:
TCanvas* c1 = new TCanvas("c1","Dynamic Filling Example",200,10,600,900);
- 创建一个TH1F直方图,和两个TH2K(高级直方图)直方图,TH1K类支持在聚类分析中广泛使用的最近K邻域方法,此方法对于小型统计尤其有用:
TH1 *hpx[3];
hpx[0] = new TH1F("hp0","Normal histogram",1000,-4,4);
hpx[1] = new TH1K("hk1","Nearest Neighbour of order 3",1000,-4,4);
hpx[2] = new TH1K("hk2","Nearest Neighbour of order 16",1000,-4,4,16);
c1->Divide(1,3);
Int_t j;
for (j=0;j<3;j++) {
c1->cd(j+1);
hpx[j]->SetFillColor(48);
hpx[j]->Draw();
}
- 随机填充直方图:
gRandom->SetSeed(12345); //设置随机种子
Float_t px, py, pz;
const Int_t kUPDATE = 10;
for (Int_t i = 0; i <= 300; i++) {
gRandom->Rannor(px,py);
for (j=0;j<3;j++) {hpx[j]->Fill(px);}
if (i && (i%kUPDATE) == 0) {
padRefresh(c1);
}
}
其中,padRefresh(c1)
函数定义如下:
void padRefresh(TPad *pad,int flag=0)
{
if (!pad) return;
pad->Modified();
pad->Update();
TList *tl = pad->GetListOfPrimitives();
if (!tl) return;
TListIter next(tl);
TObject *to;
while ((to=next())) {
if (to->InheritsFrom(TPad::Class())) padRefresh((TPad*)to,1);}
if (flag) return;
gSystem->ProcessEvents();
}
- 对一维直方图进行高斯拟合:
for (j=0;j<3;j++) hpx[j]->Fit("gaus");
padRefresh(c1);
代码地址:https://github.com/root-project/root/blob/master/tutorials/hist/hksimple.C