粒子物理与核物理试验中的数据分析.PPT

Keep it. * * 查看Tree的信息 root -l tree1.root 打开root文件 root[1].ls 查看文件信息， 发现TTree t1 root[2]t1-Show(0); 显示第0个event的信息 root[3]t1-GetEntries() 总事例数 root[4]t1-Scan(); root[5]t1-Print(); root[6]t1-Draw(px); * 查看Tree的信息(续) 也可以 root -l 进入root root[0]TFile *f1=new TFile(tree1.root); root[1]t1-Draw( sqrt(px*px+py*py) ); root[2]TH1F *h1; root[3]t1-Draw(pxh1); root[4]t1-Draw(py,px0,sames); root[5]t1-Draw(py,,sames); root[5]t1-StartViewer() * 查看Tree的信息(2) TBrowser b TBrowser打开一个浏览器，从中可以选择root文件，并一层层进入其中的tree，branch以及leaf。类似于Windows下的Explorer。 这是当前目录，双击进入并选择要打开的root文件，以及文件中的tree，最后可以看到tree的各个leaf 双击leaf可以查看leaf的直方图 适合初步浏览，但不适合具体的数据分析处理。并不推荐使用。 * ... const Int_t kMaxTrack = 50; Int_t ntrack; Float_t px[kMaxTrack]; Float_t py[kMaxTrack]; Float_t zv[kMaxTrack]; Double_t pv[3]; TFile f(rootfile,recreate); TTree *t3 = new TTree(t3, Reconst events); t3-Branch(ntrack,ntrack,ntrack/I); t3-Branch(px,px,px[ntrack]/F); t3-Branch(py,py,py[ntrack]/F); t3-Branch(zv,zv,zv[ntrack]/F); t3-Branch(pv,pv,pv[3]/D); ... 运行:进入ROOT环境后 .L ex42.C ex42w() ex42r() void ex42r() {//读取数据，适用于简单分析 TFile *f = new TFile(rootfile); TTree *t3 = (TTree*)f-Get(t3); t3-Draw(sqrt(px*px+py*py)); htemp-SetLineColor(2); t3-Draw(sqrt(px*px+py*py), zv100,sames); } 如何读写含有不定长数组的tree(1)$EXAMPLEDIR/ex42.C 1)估计不定长数组的最大维数，以该维数定义数组；如float zv[kMaxTrack] 2)定义某变量，用于存放数组的实际维数。如int ntrack,表示一个事例中实际的径迹数。 3)定义tree，设置Branch。第三个参数给出数组的实际维数。如”zv[ntrack]/F” 很多时候不定长数组是必要的，比如正负电子对撞，记录末态粒子的信息，末态粒子数目是不固定的。 ！！如何获取root文件中的tree指针 直接画出Branch/Leaf， 可以加很多条件。 * void ex42r2() { TFile *f = new TFile(rootfile); TTree *t3 = (TTree*)f-Get(t3); //步骤1:定义好必要的变量 const Int_t kMaxTrack = 100; Int_t ntrack; Float_t px[kMaxTrack]; //[ntrack] Float_t py[kMaxTrack]; //[ntrack] Float_t zv[kMaxTrack]; //[ntrack] Double_t pv[3]; //步骤2: 用SetBranchAddress函数 //将tree的Branch与定义好的变量 /