*ROOTの覚書き(ROOT様と仲良くなりたいんだよ、現在多分はとこくらいの関係) [#p2eada52]
UsersGuideは[[ここ:http://root.cern.ch/root/]]にあり！~
便利なクイックリファレンスは[[ここ:http://root.cern.ch/root/html/ClassIndex.html]]!~
#contents
~
*Minuit [#x5d8a802]
Minuitを使うための備忘録。
**使用手順 [#y5bc88d5]
まずminuitを定義。
TMinuit minui(9);
この9は変数の数。~
次に関数を定義してやる。FCNという。~
それから、
minuit.SetFCN(function);
で関数をセット。関数はパラメータに対してフィッティングの正しさの基準となる値(χ^2,Likelihoodなど)を返す。これが最小化(時に最大化)するようなパラメータを探し、良いところで収束する。収束した時の値などは、minuit.GetParameter()でとってこれる。~
それから、各変数のパラメータをセットする。たとえば
minuit.DefineParameter(0, "x", 0.0, 1.0, -200.0, 200.0);
なら、0番目のxという変数を0.0を初期値で、ステップは1.0、範囲は-200.0〜200.まで、という意味になる。~
もし変数を固定したいときには、
minuit.FixParameter(0);
とすれば良い。0は変数の番号。~
実際に実行するときには
minuit.Migrad();
でOK。この場合はMigradというアルゴリズムを使う。他にもMinosとかいろいろあるらしい。
**収束を見分ける [#p9f8a522]
フィッティングの結果は、本当の値に近いところに収束するはず。それを判断するためのひとつの基準がEDM(Estimated Distance to the Minimum)というもの。これがある基準よりも小さくなったときに収束したとみなす。~
逆に、これが基準をみたしていないときには、収束しなかったといえる。~
*技 [#y1665aba]
**テキストを挿入する [#x8d9cb2d]
テキストを入れたいときには、TTextかTLatexを使う。たとえば、
TText *t = new TText(0.1,0.1,"ahoaho");
で定義するが、 このまま何もせずにDrawすればもともとCanvas内にかかれていたもの(たとえばTH1F)とかの座標の場所として、(0.1,0.1)の場所に表示される。~
これをCanvasの絶対座標にしたければ
t->SetNDC(1);
とすればよい。(ちなみに元に戻すのはSetNDC(0))~
こうすると、(0.1,0.1)はCanvasのx,y10%の場所に表示される。~
**ヒストグラムをコピー [#lb6f56d5]
ヒストグラムをコピーしたいときには、
TH1F *hnew = (TH1F*)h->Close();
hnew->SetName();
でOK。~
SetNameは名前がかぶるのを防ぐために必要。~
**1イベントずつデータをみる [#k22513e9]
roottupleはあるけど取り出したい情報がカットとかで出せなくて、結局1イベントずつみていきたいときがある。~
たとえばntuple.rootに
nt--adc
というツリー、ブランチがあるときに、adc(>200)のヒット数をとるマクロ。
TFile *f = new TFile("ntuple.root");
Double_t d_adc;
Int_t hit_adc = 0;
TTree *tr = =(TTree*)f->Get("nt");
tr->SetBranchAddress("adc" , &d_adc );
for(Int_t i = 0;i < tr->GetEntries();i++){
tr->GetEntry( i );
if( d_adc > 200){
nhit_adc++;
}
}
最初にd_adcのような適当な変数を定義しておいて、
TTree *tr = =(TTree*)f->Get("nt");
でツリーのアドレスをとってきて、
tr->SetBranchAddress("adc" , &d_adc );
というので、ntツリーの中のadcというのをd_adcという変数に関連付ける。他のもみたいときはこの行を増やせばいい。~
その後
for(Int_t i = 0;i < tr->GetEntries();i++){
tr->GetEntry( i );
------
}
しているが、GetEntryすることでd_adcに入る値がiイベント目のadcの値になるので、これをループすることで、1イベントずつみていける。~
**ROOTをC++とかで使う [#t10ab309]
この方が早くて便利だったりする。いろいろ楽だったりもする。
でも身に余る技術だもんで、使いこなせない。そんな技術に振り回された過程をメモります。~
~
GUIを呼び出したかったのでいろいろ調べた。~
[[山下研:http://www.icepp.s.u-tokyo.ac.jp/~kubota/pukiwiki/index.php?ROOT]]~
[[つくばの方:http://utkhii.px.tsukuba.ac.jp/~nagata/root/native.html]]~
~
結局おきまりとして、
#include "TApplication.h"
--main--
TApplication theApp("App", &argc, argv);
という感じでTApplication.hを読み込めば良いみたい。でも&argcとかは何でいるんだろう。なんとかしたい。
~
**TGraphErrorsを使いこなす [#k4b36382]
TGraphErrorsがむかつくので、好きになれるように覚書き。
-TGraphErrorsでレンジを設定~
Drawオプションで"A"を使えば自動的にレンジを決めてくれるが、複数重ねた途端に、最初のグラフに引っ張られてむかつく存在になる。~
そういうときのためにframeを使う。
TH1F *frame = gPad->DrawFrame(Xstrat,Ystart,Xend,Yend);
でレンジを決められる。このときは"A"オプションは不要。~
canvasを分割している場合でも、フォーカスしているcanvasに適応される。
-TGraphErrorsで軸を書く~
frameを使っていれば
frame->GetXaxis()->SetTitle("xaxis");
でOK。
-レジェンドを作る~
てか、なんでレジェンドっていうのかな？調べても出てこないです。わからん。->「凡例」のことみたいです。~
使い方は
TLegend *led = new TLegend(size,"");
leg->AddEntry(graph,name,optoin(ex."apl"));
で追加していける。~
~
(2011.5.22)~
TGraphにいれたデータを直接触りたいときには
for(int index = 0;index < g->GetN();index++){
x = g->GetX()[index];
xerr = g->GetErrorX(index);
y = g->GetY()[index];
yerr = g->GetErrorY(index);
}
という感じでいける。~
**グラフの情報ボックスをいじくる [#lccbb881]
たとえばfit情報を出したいときは
gStyle->SetOptFit();
で出てくる・・・が、そのままだと往々にして文字が小さすぎるので、大きさを変える方法のメモ。~
高さを変える
gStyle->SetStatH(3);
幅を変える
gStyle->SetStatW(0.15);
フォントのサイズを変える
gStyle->SetStatFontSize(30);
他にもいろいろできるみたい。
~
sameオプションで追加して書き込むときには"same"を"sames"にすることで、情報ボックスは再描画される。このときに、表示したものをオブジェクトに固めてやらないと統計ボックスは描画されない。
~
2011.05.06追記:~
たとえばTH1Fとかの情報ボックスを消したいときには
TH1F *h;
h->Draw();
h->SetStats( 0 );
で消せる。~
**ROOTでSharedObject [#je1b9c37]
まずライブラリを作る。たとえば適当にexample.ccを次のように書く。
void aho(void){
std::cout << "aho!" << std::endl;
}
そのあとrootを立ち上げて、
>.x example.cc+
とすれば、勝手にexample_C.soファイルを生成してくれる。これで簡単にダイナミックに関数なんかを読み込める。すげー！~
gccでやろうとしたが、普通のやり方では駄目みたい。できるらしいけど、僕は挫折。~
~
その後soファイルは残るので、マクロで
gSystem->Load("example_C.so");
としてやれば、オブジェクトを読み込める。これで楽々プログラミング。~
**ROOTで乱数を使う [#ac0d6d44]
[[ここ:http://blogs.yahoo.co.jp/eguchium/31456122.html]]を参考にさせていただいた。~
x = gRandom->Gaus(0,3);
でガウス分布に従う乱数をとってこれる。引数はmeanとsigma。この種類はいくつかあって、それぞれ引数が違う。種類は次のものがある。
Rndm() or Uniform(min,max), Gaus(mean,sigma), Exp(tau), BreitWigner(mean,sigma),
Landau(mean,sigma), Poisson(mean), Binomial(ntot,prob)
がしがし使っていこう。ちなみに、gRandomを使うには"TRandom.h"を読み込む必要あり。~
あと、種作りは忘れずに。
gRandom->SetSeed( time(NULL) )
**NTupleに詰め込む [#md817886]
たとえば、こんな感じ。
TNtuple *ntuple = new TNtuple("h3","castum1","iev:z:Lmu:nphot");
for(...){
.....
ntuple->Fill(i,Z[i],peak[i],NOP[i]);
}
TNtuple構造体を定義してやるときに詰め込む要素の名前を指定してやって、あとは詰め込むだけ。
この場合はh3というブランチに"iev","z","Lmu","nphot"というリーフをつけてやっている。こいつらは1イベントi毎に関連づけられた状態で詰め込まれていく。
サイズを設定したりしたほうがよりよい。あとはこの応用です。
~
新しく作ったntupleを保存するには次の手順を踏めばOK。
TFile *tmproot = new TFile("tmp.root","recreate");
(たまにここで f->cd(); しないとだめなときがある)
TNtuple *hnew = (TNtuple*)f->FindObject("hnew");
hnew->Write();
tmproot->Close();
新しくつくったhnewのポインタをfというFileから探してきて、新しく作ったtmprootというところにwriteする。
~
~
注意~
TNtuopleには一気に14個までしかFillできない。なので、それを超えるときには配列に詰め込んでやって、その配列を渡すしかない。配列の型はfloatにしないといけない。~
**線を引く [#ea501c56]
始点終点を指定すればOK。
TLine *line = new TLine(88,0,88,-300);
でいける。
この場合は(88,0)から(88,-300)への直線を引いている。
**gauss分布でfitしたときの面積を算出 [#r0d1ed87]
fitparameterの設定のときにgaussでなくgaussnで定義すれば、fitparameter0(P0)が面積になる。これはただ単に、fitする際の変数にしておくパラメータを高さから面積に変えただけ。~
*ROOT環境 [#u89e308a]
**root解析をSGEに投げる [#y72dc546]
あんまり時間がかかる作業をさせるときに、自分の端末でやっていたら動作が遅くて大変。なのでbsubで他のマシンに仕事をさせる。~
rootで気をつけることはDrawしたときに勝手にcanvaが表示されるので、
root -l -b -q aho.cc
のようにすると-b(バッチジョブフラグ)の効果で表示されなくなる。-qはマクロを実行したら勝手に終了させる引数で、これがないとjobが終わらない。~
実行シェルの例。~
#!/bin/csh
cd WORK_DIRECTORY
root -l -b -q ana4_jobs.cc
~
*エラー覚書 [#n10ce1f1]
でた！エラーがでたよ！わーい！
**Warning in <TROOT::Append>: Replacing existing TH1: h_aho (Potential memory leak). [#jf620b90]
TH1Fを触ったときに出るエラー。どうも何回も重ねてdeleteもせずにnewしていたのが問題みたい。~
ループの場所を変えて、newしないようにしたら、表示されなくなった。
