ROOT のバックアップ(No.7)

ROOTの覚書き(ROOT様と仲良くなりたいんだよ[現在はとこくらいの関係]) †

UsersGuideはここにあり！
便利なクイックリファレンスはここ!
技
...NTupleに詰め込む
...線を引く
...Gauss分布でfitしたときの面積を算出
 ROOT環境
...root解析をSGEに投げる

技 †

NTupleに詰め込む †

たとえば、こんな感じ。

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");
TNtuple *hnew = (TNtuple*)f->FindObject("hnew");
hnew->Write();
tmproot->Close();

新しくつくったhnewのポインタをfというFileから探してきて、新しく作ったtmprootというところにwriteする。

注意
TNtuopleには一気に14個までしかFillできない。なので、それを超えるときには配列に詰め込んでやって、その配列を渡すしかない。配列の型はfloatにしないといけない。

↑

線を引く †

始点終点を指定すればOK。

TLine *line = new TLine(88,0,88,-300);

でいける。この場合は(88,0)から(88,-300)への直線を引いている。

↑

gauss分布でfitしたときの面積を算出 †

fitparameterの設定のときにgaussでなくgaussnで定義すれば、fitparameter0(P0)が面積になる。これはただ単に、fitする際の変数にしておくパラメータを高さから面積に変えただけ。

↑

ROOT環境 †

↑

root解析をSGEに投げる †

あんまり時間がかかる作業をさせるときに、自分の端末でやっていたら動作が遅くて大変。なので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