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

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

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

↑

技 †

↑

ROOTでSharedObject †

まずライブラリを作る。たとえば適当に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で乱数を使う †

ここを参考にさせていただいた。

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"を読み込む必要あり。

↑

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