猫でもわかるROOT
1.ヒストグラムを作成する
1−1:一次元のヒストグラムの作成
データファイルのある作業ディレクトリに移ってROOTを実行
今回はサンプルデータrun147.spcを利用
データ形式は
チャンネル カウント数
の順番に並んでいる
- //ROOT
- gROOT->Reset();
- gStyle->SetOptStat(0);
- gStyle->SetPalette(1);
- gStyle->SetCanvasColor(33);
- gStyle->SetFrameFillColor(10);
-
- TH1F hs1("hist","Scat+Abs",128,0,1.0);
- hs1->SetFillColor(kBlue);
-
- float energyS1, energyA1, energyA2, energyA3;
- TH1F h1("hist1","Nai-Cs",512,0,4095);
- h1->SetFillColor(kRed);
-
- ifstream data("run147.txt");
- while(data >>energy ) {
- h1.Fill(energyS1:energyA1);
- }
- data.close();
-
- TCanvas C1("C1","Nai Cs hist",10,10,800,600);
- c2h.Divide(2,2);
-
- h1->GetXaxis()->SetTitle("ch");
- h1->GetYaxis()->SetTitle("counts");
-
- C1.cd(1);
- hs1.Draw();
- C1.Update();