GiriWiki ja:root

ja:root:axis

2018-10-20T18:38:56+00:00

軸の設定・装飾 ---------- はじめに ROOTでは軸の目盛文字のことを”Label”といい、[MeV]とか[cm]とかの軸に添わせた文字を”Title”という。また、目盛線のことを”Tick”という。 ----------

ja:root:color

2018-08-30T05:24:58+00:00

カラー・スタイル・マーカーカラーテーブル (TColor) TCanvas *c = new TCanvas("c","Fill Area colors",0,0,500,200); c.DrawColorTable(); * 使い方は単に SetFillColor(2) のように、番号で色を指定すればよい。

ja:root:fit

2017-08-30T07:42:55+00:00

Fitting ---------- 簡単な例 TH1F *hist = new TH1F("hist","hist",100,-10,10); hist->FillRandom("gaus",10000); hist->Fit("gaus","","",-10,10); としてやればガウシアンでのfittingの例が得られる。 “gaus”の他にも以下のような関数が使える。

ja:root:font

2016-06-20T05:09:26+00:00

フォント ---------- 基本的にはTAttTextを参照。 Text font code = 10*fontnumber + precision The precision can be: precision = 0 fast hardware fonts (steps in the size) precision = 1 scalable and rotatable hardware fonts (see below) precision = 2 scalable and rotatable hardware fonts precision = 3 scalable and rotatable hardware fonts. Text size is given in pixels.

ja:root:hist

2018-08-03T02:36:16+00:00

ヒストグラムの装飾 ---------- スケールする hist->Scale(2.5); //2.5倍に（縦に）スケールする * hist1 , hist2 , hist3 , ... とあって、hist1 のエントリー数でスケールしたい（規格化したい）ときは以下のようにすればよい（複数のヒストグラムの形を比べるときなど）

ja:root:install

2015-10-04T12:41:24+00:00

source からコンパイルする以下のようにしておけば、バージョン更新・戻しもシンボリックリンクの張り替えだけで済み便利 ---------- ダウンロード * 以下、ほとんどにsudoが要るかも（mac）

ja:root:stat

2019-04-18T16:23:40+00:00

統計情報（統計box） ---------- 統計情報の表示 * 表示される情報は、SetOptStat で制御される gStyle->SetOptStat(mode); * modeで指定するのは、尖度(k)、歪度(s)、ビン数(i)、オーバーフロー数(o)、アンダーフロー数(u)、RMS(

ja:root:tcanvas

2014-05-23T06:18:12+00:00

キャンバスの設定・装飾 ---------- Canvas の作成 * Canvas を作成（TCanvas クラスのインスタンスを作成） TCanvas *c = new TCanvas("c","c",600,500);//TCanvas(const char* name, const char* title, Int_t ww, Int_t wh)

ja:root:tcut

2014-04-14T04:34:28+00:00

TCut ---------- TCut * cutの条件を書いておける * 同じcutを使い回したりするときにかなり便利 * 使い方は、大体下のような感じ TCut c1 = "x<1"; //文字列を " " で囲んで書く TCut c2 = "y<0"; // TCut c3 = c1&&c2; //合成することもできる tree->Draw("x", c1);　　　　　//一番基本的な使い方 tree->Draw("x", c1||"x>0"); //他のcut条件(not TCut)と組み合わせるときは " " でくくって並べる tree->Draw("x", c1&&c2); //TCutを合わせるときは " " でくくらずそのまま tree->Draw("x", "(x+y)"*(c1&&c2)); // " * " て何？…

ja:root:tdatime

2017-04-06T00:19:03+00:00

TDatime ---------- 正直あんまり便利じゃないかも。時間を横軸にしたりするときに使えるかも。 #include TDatime::TDatime()// 引数なしだと現在時刻をとる。 TDatime::TDatime(Int_t year, Int_t month, Int_t day, Int_t hour, Int_t min, Int_t sec) TDatime::Set(Int_t year, Int_t month, Int_t day, Int_t hour, Int_t min, Int_t sec) 　UInt_t TDatime::Convert(Bool_t toGMT=kFALSE)// UNIX time で出力　char* TDatime::AsString()//出力　char* TDatime::AsSQLString()//SQLのフォーマットで出力…

ja:root:tf1

2017-08-30T14:38:22+00:00

TF1 ---------- 基本 TF1 *func1 = nes TF1("func1" , "x*sin(x)" , x_min , x_max);//x_min,x_maxは定義域 TF1 *func2 = nes TF1("func2" , "gaus" , x_min , x_max);// システム定義関数を使うことができる TF1 *func3 = nes TF1("func3" , "gaus(0)+pol1(3)" , x_min , x_max);// 複数使う場合は最初のパラメータの番号を指定する

ja:root:tgrapherrors

2017-02-02T14:25:24+00:00

TGraphErrors TGraph の上位互換（TGraphを継承している）なので、基本的にこちらを使っておけばよい。 ---------- コンストラクタ TGraphErrors(Int_t n, const Float_t* x, const Float_t* y, const Float_t* ex = 0, const Float_t* ey = 0) TGraphErrors(Int_t n, const Double_t* x, const Double_t* y, const Double_t* ex = 0, const Double_t* ey = 0) TGraphErrors(const TVectorF& vx, const TVectorF& vy, const TVectorF& vex, const TVectorF& vey) TGraphErrors(const TVectorD& vx, const TVectorD& vy, const TVectorD& vex, const TVectorD& vey) TG…

ja:root:tlegend

2018-10-16T12:09:36+00:00

凡例（TLegend） ---------- * ヒストグラム等の色分け等を表示 TLegend *legend = new TLegend( 0.8, 0.68, 0.99, 0.78) ; //（）の中は位置の指定（左下の x , y 、右上の x , y ） legend->AddEntry( h0, "histogram 0" , "l") ; // AddEntry( pointer , "interpretation" , "option" ) legend->AddEntry( h1, "histogram 1" , "l") ; // option は　"f"=box, "l"="L"=line, "p"=marker legend->AddEntry( h2, "histogram 2" , "l") ; legend->SetFillColor(0); legend->Draw() ;…

ja:root:tmultigraph

2014-08-19T14:05:32+00:00

TMultiGraph 複数のグラフを１つにして扱ったりできる。複数のグラフをマージして、Fitしたりという使い方ができる。 ---------- 使い方 TGraph *gr1 = new TGraph(... ; TGraphErrors *gr2 = new TGraphErrors(... ; TMultiGraph *mg = new TMultiGraph(); mg->Add(gr1,"lp"); mg->Add(gr2,"cp"); mg->Draw("a"); mg->Fit("pol1");

ja:root:tproof

2018-09-29T01:08:03+00:00

TProof ---------- マルチコアリソースを有効活用するのに便利なクラスらしい。 TChainでのDrawにしか使ったことがないのでその話。 TChain では複数のFileからtreeをつなぎ合わせてくるが、Drawの時にその処理を並列化してくれるというもの。

ja:root:tstopwatch

2014-06-11T08:18:15+00:00

TStopwatch ---------- * その名の通りストップウォッチの様に使える。 * #include 基本的な機能 TStopwatch sw; sw.Start(); //スタート　前に使っていたら、その値をresetして、０から計測する。 sw.Start(1);//上と同じ sw.Start(0);//前に使っていた値をresetせず、その続きから計測。 sw.Stop(); //ストップ sw.Continue(); //ストップしていたところから計測再開 sw.RealTime(); // （ストップをかけて、）running time を秒の単位で返す。返り値はDouble_t。 sw.CPUtime(); // （ストップをかけて、）CPU time を秒の単位で返す。返り値はDouble_t。 sw.Print(); //Real time 0:00:04, CP time 0.000 のように返す(running time 4[sec] ,CPU time 0 [sec]) sw.Prin…

ja:root:tstring

2018-02-06T07:21:43+00:00

TString ---------- TString型の変数に文字を入れる * 初期化 * やり方色々 TString str("hoge"); TString str = "hoge"; TString str; //まず宣言してから str.Form("hoge"); //Form関数を使って値を入れる

ja:root:ttree

2017-08-19T09:27:25+00:00

TTree ---------- 読み込むBranchを限定する方法１ --- SetBranchStatus デフォルトではGetEntryすると、SetBranchAddressで指定されていないBranchまでも読みにいく。というかすべてのbranchのデータをメモリに乗せるが、SetBranchAddressするとそのbranchに対応するメモリのポインタが指定した変数に割り当てられる。

ja:root:wantto1

2014-08-01T02:31:00+00:00

ヒストグラムが最大値を取るX座標を取ってきたい histo->GetXaxis()->GetBinCenter(histo->GetMaximumBin()); とすればよい。 * histo->GetMaximumBin() でヒストグラムが最大値を取るBin番号を取ってくる *

ja:root:wantto2

2018-09-05T17:42:34+00:00

ヒストグラムを規格化(Scale)して、さらにエラーバーをつけたい ---------- ヒストグラムを規格化するエントリー数の違うヒストグラムの形を比べる時などに活躍。 TH1::Scale を使って h->Scale(ScaleFactor);

ja:root:wantto3

2015-04-01T07:48:04+00:00

Shell command を実行させたい gSystem->Exec("COMMAND") を使えば良い。

ja:root:wantto4

2018-07-27T00:57:09+00:00

イベントループの進行状況を表示させたい Tree から1イベントずつ取ってきて処理するような解析をよくやるが、その時にイベントループの進行状況を表示させたい。 const Long64_t nentry = tree->GetEntries(); // Loopさせるイベント数を取得 const Long64_t n_100 = nentry/100; // イベント数÷100 を定義しておく cout << "# events : " << nentry << endl; for( Long64_t ientry = 0 ; ientry < nentry ; ++ientry ){ if( nentry>100 ){ if( (ientry%n_100==0) && (ientry!=0) ){ if( ientry % (n_100 * 10) == 0 ) std::cerr << (ientry / n_100) << "% processed. " << std::en…

ja:root:wantto5

2018-08-22T16:07:51+00:00

root fileをマージしたい ---------- root fileをマージするのに、haddというコマンドがある。(${ROOTSYS}/bin/ 以下) hadd output.root input1.root input2.root input3.root .... のように第一引数に新規の出力ファイル名、第二引数以降にマージしたいroot file名をとればよい。

ja:root:wantto6

2018-08-22T16:10:47+00:00

treeの中のbranchの有無を判定したい ---------- TTree::GetBranch(const char *name)を使う。返り値はそのbranchへのポインタ(TBranch*)だが、指定した名前のbranchがなかった場合、0が返る(sourceを見るとnullではなく0)ので、

ja:root:wantto7

2018-08-22T17:42:07+00:00

軸のラベルをtextにする ---------- 元のラベルを飛ばして、手動でtextを配置する方法 hist->GetXaxis()->SetLabelOffset(100); // 範囲外へ元の数字のラベルを飛ばしてしまう hist->GetXaxis()->SetNdivisions(3,1,1); TLatex *label = new TLatex(); label->SetTextSize(0.07); label->DrawLatexNDC(0.18,0.04,"K_{L}#rightarrow3#pi^{0}");// TLatexでtextを書いていく label->DrawLatexNDC(0.45,0.04,"K_{L}#rightarrow2#pi^{0}"); label->DrawLatexNDC(0.72,0.04,"K_{L}#rightarrow2#gamma");…

ja:root:wantto8

2017-01-23T02:27:35+00:00

ヒストグラムのbinの中のentry数を表示させたい ---------- TH1::Draw の“TEXT” optionを使う。ただし、“TEXT” だけだとhistogram自身が描画されず、数字だけ、しかも横向きに表示されてしまう。そのため、

ja:root:wantto9

2018-11-07T06:18:26+00:00

treeからグラフを作る ---------- TTree::GetSelectedRows()とかTTree::GetV1()とかの存在を初めて知ったのでメモ。 Rootのclass referenceの”How to obtain more info from TTree::Draw”のところに色々書いてある。 tree->Draw("X:Y",Cut);//一旦Drawしておく (普段の二次元ヒストグラムを描く時とは順番が違うのに注意) const Int_t N = tree->GetSelectedRows();//Drawでselectionを通ったイベントの数を取ってくる TGraph *graph = new TGraph(N, tree->GetV1(), tree->GetV2());…

ja:root:wantto10

2017-04-06T10:19:10+00:00

TGraphの統計Boxをいじりたい ---------- コード内で出てくる全てについて同じ処理をしたい時はTStyleの方のメソッドを使えばよい。（gStyle->SetOptFit(0)とか）以下は個別でいろいろ変えたい時用。

ja:root:wantto11

2018-08-22T17:27:00+00:00

TChainでループを回している時、参照しているイベントの属するTTreeが記録されたファイル名を知りたい ---------- TChainを使ってイベント毎に処理をさせている時、当該イベントがどのファイルに属しているtreeに含まれているイベントなのかを知るのには、

ja:root:wantto12

2017-07-04T05:45:31+00:00

TTreeのScanでカラムの幅、数値の精度などを指定したい ---------- TTree::Scan()の第3引数がオプション指定になっている。 tree->Scan("a:b:c","","colsize=30 precision=3 col=::20.10:#x:5ld"); のように指定する。

ja:root:wantto13

2017-09-23T11:40:01+00:00

Fit関数としてクラスのメンバ関数を使いたい ---------- static メンバ関数としておく必要があるみたい。

ja:root:wantto14

2018-08-03T02:24:39+00:00

ヒストグラムの最大値をとってきたい ---------- hist->GetMaximum() でOK。

ja:root:wantto15

2018-08-07T06:14:25+00:00

１次元ヒストグラムを描いた後、フレームの最大値・最小値を知りたい ---------- 例えばヒストグラムを描いた後に、1σのところにフレームいっぱいに縦に線をいれたい、とする。そうした時にヒストグラムのフレームの上限、下限の値が知りたいが、それはヒストグラムのエントリー数によっても変わるし、 Y軸がリニアかログかでも変わる。その辺をうまくやる方法が、…

ja:root:wantto16

2018-09-05T12:47:24+00:00

Eventごとにweightが付いている場合の統計誤差を得たい ---------- * weightなしの場合、Nイベントあった場合の統計誤差は√Nで表すことができる。(N ± √N) * weightが全イベントで同じ場合(W)は,中心値 W*N, 誤差W*√N (W*N ± W*√N)

ja:root:wantto17

2018-09-20T09:23:00+00:00

macroを実行したら自動的にROOTを終了するようにしたい ---------- マクロを実行後、そのままROOTを終了するようにするためには、 root -q macro.C のようにすればよい。毎度 -q を入れるのが面倒ならば、マクロの最後に以下のどれかを加えてやれば良い。

ja:root:wantto18

2018-09-29T08:42:38+00:00

macroに引数を渡したい ---------- たとえば、Int_t と TString の2変数をとるマクロがあるとする。 test.C void test( const Int_t N, const TString str ) { std::cout << "N = " << N << std::endl; std::cout << "str = " << str.Data() << std::endl; }

ja:root:wantto19

2018-11-03T07:40:35+00:00

ヒストグラムのbin、境界上はどっちに入る？ヒストグラムのbin IDは, ヒストグラムのbin数をnbinsとして bin = 0; underflow bin bin = 1; first bin with low-edge xlow INCLUDED bin = nbins; last bin with upper-edge xup EXCLUDED bin = nbins+1; overflow bin