杨振伟 清华大学 第三讲：ROOT在数据分析中的应用(1) 粒子物理与核物理实验中的数据分析 杨振伟 清华大学 第三讲：ROOT在数据分析中的应用(1)

上讲摘要 C++基本概念 int,float,double型变量,循环，类... Linux下用g++编译C++程序 g++ -o hello.exe -I<include> ./src/*.cc 当前目录下输出 指定include目录 源文件 可执行文件hello.exe 如-I./include Linux下编写makefile进行C++编译 gmake 进行编译 gmake clean 清除编译结果

本讲要点 ROOT 基本概念 登录ROOT ROOT的语法简介 ROOT的函数,直方图,随机数,文件,散点图 TF1,TH1I,TH1F,TH1D,gRandom TF2,TF3,TH2F,... TFile

什么是 ROOT ？ 面对对象的程序设计C++分析软件包 root 文件包括直方图和tree，用户可以存储目标模块，数值等等，具有很好的绘图与分析功能。 可进行C++解释和编译,并含拓展的语言CINT 是粒子与核物理实验的基本架构之一

安装ROOT 到ROOT主页下载相应版本到指定目录。比如要在SLC3系统的/projects/yangzw目录下安装5.16.00版本 cd /projects/$USER (对用户yangzw, $USER=yangzw;对用户libo, $USER=libo) wget ftp://root.cern.ch/root/root_v5.16.00.Linux.slc3.gcc3.2.3.tar.gz tar –zxvf roo_v5.16.00.Linux.slc3.gcc3.2.3.tar.gz 设置ROOT的环境变量 export ROOTSYS=/projects/$USER/root export PATH=$ROOTSYS/bin:$PATH export LD_LIBRARY_PATH=$ROOTSYS/lib:$LD_LIBRARY_PATH 可以把上面这3行放到$HOME/.bashrc或者.login文件中 这样每次登录到Linux系统，系统就自动设置ROOT的环境变量 如果是其它发行版的Linux，需要到http://root.cern.ch 下载ROOT的源代码，按照安装指南通过gmake编译安装。

安装root，在 自己home目录中的 .login或 .bashrc 或 .cshrc中设置 root 的环境 键入 help 指令，如 root[0]? 退出 root[1].q ROOT环境其它常用指令： .L macro.C Load文件macro.C .x macro.C 执行文件macro.C .ls 显示ROOT当前环境的所有信息 .! ls 显示Linux系统当前目录的所有信息 注：ROOT环境中，ROOT指令都以“.”开头 系统指令都以“.!”开头

ROOT语法 ROOT使用C++语法 一段C++程序可以直接在ROOT环境运行 类型重定义 int  Int_t float  Float_t double  Double_t ...... ROOT的类都以T开头 如TFile, TH1F, TTree, ... 详细规定参阅ROOT手册(5.16版)第18-20页，关于Convention和Global Variables部分。 可以直接在ROOT环境中运行macro文件(自动调用cint编译器)，也可以在makefile中设置好相关参数用g++编译运行。

ROOT脚本文件示例(1)：Macro文件 /home/yangzw/workdir/examples/Lec3/ex31.C { cout << "Hello ROOT" << endl; int Num=5; for (int i=0;i<Num;i++) { cout << "i=" << i << endl; } 纯粹C++语法，执行的时候只需要在命令提示行： cd /projects/$USER cp -r /home/yangzw/workdir/examples/Lec3 . (注意最后有个“.”) cd Lec3 root -l ex31.C

在root中制作函数曲线图 画图时采用 root[1]fun_name.Draw(); 制作一维函数曲线图 制作二维函数曲线图 root[0]TF1 *f1 = new TF1(“f1”,”x*sin(x)”,-5,5); 制作一维函数曲线图 TF1 *fun_name = new TF1(“fun_name”,”function”, x_low,x_high); root[0]TF2 *f2 = new TF2(“f2”,”x*sin(x)+y*cos(y)”, -5,5,-10,10); 制作二维函数曲线图 TF2 *fun_name = new TH2(“fun_name”,”function”, x_low,x_high, y_low,y_high); 画图时采用 root[1]fun_name.Draw(); root[0]TF3 *f3 = new TF2(“f3”,”x*sin(x)+y*cos(y) +z*exp(z)”,-5,5,-10,10,-20,20); 制作三维函数曲线图 TF3 *fun_name = new TF3(“fun_name”,”function”, x_low,x_high,y_low,y_high,z_low,z_high);

在root中制作统计图 定制一维直方图 定制二维图 定制三维图 填充统计图 绘图：root[0]hist_name.Draw(); TH1F *hist_name = new TH1F(“hist_name”,”hist_title”, num_bins,x_low,x_high); 定制二维图 TH2F *hist_name = new TH2F(“hist_name”,”hist_title”, num_bins_x,x_low,x_high,num_bins_y,y_low,y_high); 定制三维图 TH3F *hist_name = new TH3F(“hist_name”,”hist_title”, num_bins_x,x_low,x_high,num_bins_y,y_low,y_high, num_bins_z,z_low,z_high); 填充统计图 hist_name.Fill(x); hist_name.Fill(x,y); Hist_name.Fill(x,y,z); 绘图：root[0]hist_name.Draw();

ROOT脚本文件示例(2):数学函数定义 /home/yangzw/workdir/examples/Lec3/ex32.C //a simple ROOT macro, ex32.C //说明ROOT中数学函数的使用，如TF1 void ex32() { gROOT->Reset(); //定义函数 TF1 *f1 = new TF1("func1","sin(x)/x",0,10); f1->Draw();//画出函数图像 TF1 *f2 = new TF1("func1",“Gaus(x,0,1)",0,10); f2->SetLineColor(2);//设置颜色为红色 f2->Draw(“same”);//用参数”same”，把f1,f2画在同一个画布上 } 函数名称 函数表达式 函数区间 执行的时候只需要在命令提示行 root -l ex32.C 或者进入ROOT环境之后，运行 .x ex32.C 提示：1)脚本中void函数的名字必须与文件名相同(如ex32) 2)ROOT环境中定义类指针之后，如TF1 *f1，之后 输入“f1->”，然后按一下Tab键，可以自动列出 该类对象的成员函数和成员变量

ROOT脚本文件示例(3): 画布，保存文件 /home/yangzw/workdir/examples/Lec3/ex33.C //说明ROOT画布的使用，TCanvas，保存图形 void ex33() { gROOT->Reset(); //define a function sin(x)/x TF1 *f1 = new TF1("func1","sin(x)/x",0,10); //define a Gaussian function, mean=0, sigma=1 TF1 *f2 = new TF1("func2","Gaus(x,0,1)",-3,3); //定义一个画布, TCanvas TCanvas *myC1 = new TCanvas("myC1","A Canvas",10,10,800,600); //将画布分成两部分 myC1->Divide(2,1); myC1->cd(1); //进入第一部分 f1->Draw(); myC1->cd(2); //进入第二部分 f2->Draw(); myC1->SaveAs(“myex33.gif”); myC1->SaveAs(“myex33.eps”); } 像素坐标 (10,10):左上角 (800,600):右下角 名称 描述 执行的时候只需要在命令提示行 root -l ex33.C 或者进入ROOT环境之后，运行 .x ex33.C

ROOT脚本文件示例(4):直方图，随机数 /home/yangzw/workdir/examples/Lec3/ex34.C //说明ROOT直方图、随机数的使用，如TH1F, gRandom void ex34() { gROOT->Reset(); const Int_t NEntry = 10000 ; //创建一个root文件 TFile *file = new TFile(“hist1.root”,”RECREATE”); TH1F *h1 = new TH1F("h1","A simple histo",100,0,1); //填充直方图10000次，用(0,1)均匀分布 for (int i=0;i<NEntry;i++) h1->Fill( gRandom->Rndm() ); h1->Draw(); h1->GetYaxis()->SetRangeUser(0,150); h1->GetXaxis()->SetTitle("x"); h1->GetXaxis()->CenterTitle(); file->cd(); //进入文件file h1->Write();//将h1写入文件 } 名称 描述 No. of Bin 区间 调用均匀分布Rndm()，其它： Landau(mean,sigma); Binomial(ntot,prob); Poisson(mean); Exp(tau); BreitWigner (mean,sigma); 执行的时候只需要在命令提示行 root -l ex34.C 或者进入ROOT环境之后，运行 .x ex34.C

直方图、打开root文件 打开已有的root文件，如hist1.root: 终端提示行下: root –l hist1.root 直方图的名称 直方图统计信息 事例数：Entries 均 值：Mean 方 差：RMS 参见ROOT手册37页 “Statistics Display” X轴的名称 打开已有的root文件，如hist1.root: 终端提示行下: root –l hist1.root ROOT环境下： TFile f1(“hist1.root”); .ls h1->Draw();

ROOT脚本文件示例(5): 散点图 /home/yangzw/workdir/examples/Lec3/ex35.C //2维直方图TH2F，散点图，散点图的协方差 void ex35() { gROOT->Reset(); const Int_t NEntry = 10000 ; TH2F *hXY = new TH2F("hXY","2d histo",100,0,1,100,-3,3); for (int i=0;i<NEntry;i++) { float x = gRandom->Rndm() ; float y = gRandom->Gaus(0,1) ; hXY->Fill(x,y) ; //填充2维直方图 } hXY->Draw(); //2维直方图的散点图 hXY->GetXaxis()->SetTitle("X: Uniform" ); hXY->GetYaxis()->SetTitle("Y: Gaussian"); Float_t covar = hXY->GetCovariance(); //协方差 cout << "Covariance = " << covar << endl;

小结 ROOT简介 C++，面向对象，实验数据处理的强大工具 安装与登录 运行ROOT脚本 数学函数，画布，直方图，随机数，散点图 TF1，TCanvas，TH1F，gRandom， TH2F 新建root文件，查看root文件 TFile

练习 1. 写一个ROOT脚本，ex3_gaus.C, 产生高斯分布，区间(-10,10)，分100个bin，比较不同的参数的分布。参数组合为： (mean,sigma)=(0,1), (0,2), (0,3), (1,3),(2,3) 把这5个分布画在同一个图中进行比较。 2. 写一个ROOT脚本，ex3_pdf.C，作4个直方图，分别产生10000事例的Gauss,Poisson,Binomial,Landau分布。创建画布，分成3*2块，将4个直方图画在画布的1-4部分。调整参数，考察这几个分布在什么情况下分布大致相同，再将这四个直方图在画布的第5部分画到同一个图中。 3. 将练习2中产生的直方图储存到mypdf.root文件中。 将所画直方图的x/y轴添加上名称，不同分布用不同颜色。 将画布存成eps文件和gif文件 4. 将例题ex35.C中的事例数改为1000，屏幕打印出关联系数。 5. cp –r $ROOTSYS/tutorials /projects/$USER 运行以下几个文件，查看ROOT直方图的常用功能如何实现 twoscales.C, transpad.C, multicolor.C, logscales.C, hstack.C 6. 阅读ROOT手册第二章以及第三章(直方图) 熟悉ROOT语法惯例，直方图制作的各种参数，随机数的使用

参考资料 ROOT手册第2章，第3章 http://root.cern.ch http://root.cern.ch/root/Reference.html http://root.cern.ch/root/Tutorials.html http://root.cern.ch/root/HowTo.html $ROOTSYS/tutorials中的各个例子