ROOT/対話型マクロ作成のための小技

対話型マクロ作成のための小技†

bin関係†

• bin数を調べる

  hist->GetXaxis()->GetNbins()

• ヒストグラムの定義範囲の最大値・最小値

  hist->GetXaxis()->GetXmax()
  hist->GetXaxis()->GetXmin()

  • TH2のy軸を調べる場合は hist->GetYaxis()->GetXmax() のようにする

現在表示しているヒストグラムのポインタを取得する†

TH1* h1 = (TH1*)gPad->GetListOfPrimitives()->At(1);

• これで、登録されているヒストグラム名にかかわらずh1として扱える

現在開いているファイルの一覧を取得する†

gROOT->GetListOfFiles()->Capacity() // これで現在開いているファイルの個数がわかる
gROOT->GetListOfFiles()->At(0)->GetName() // 0番目のファイル名
gFile->GetName() // 現在選択されているファイル

現在開いているキャンバスの一覧を取得する†

gROOT->GetListOfCanvases()->Capacity() // これで現在開いているキャンバスの個数がわかる
gROOT->GetListOfCanvases()->At(0)->GetTitle() // 0番目のキャンバス名 (ウィンドウに表示されるタイトル)
gPad->GetCanvas()->GetTitle() // 現在選択されているキャンバス
((TCanvas*)gROOT->GetListOfCanvases()->At(1))->cd() // 1番目のキャンバスに移動

ヒストグラム名を使わずにDrawする方法†

• .lsでn番目に表示されるヒストグラムのDraw

  n=0;
  gDirectory->Get(gDirectory->GetListOfKeys()->At(n)->GetName())->Draw();
  または
  _file0->Get(_file0->GetListOfKeys()->At(n)->GetName())->Draw();

  • もっとスマートな指定方法を知ってる人がいたら教えてください。
    • gDirectoryではなくgROOTを使えばいい -> 一部環境でしか使えなかった

      TH1* h1 = (TH1*)gROOT->GetList()->At(n);
      h1->Draw();
      TH1* h1 = (TH1*)gROOT->GetList()->After(gPad->GetListOfPrimitives()->At(1));
      h1->Draw();

    • (追記) これで行けるかも?

      ((TKey*)gDirectory->GetListOfKeys()->At(n))->ReadObj()->Draw()

• ANAPAWのhn(次のヒストグラムを表示),hb(前のヒストグラムを表示)的な使い方

  n=0;
  n++;_file0->Get(_file0->GetListOfKeys()->At(n)->GetName())->Draw();
  n--;_file0->Get(_file0->GetListOfKeys()->At(n)->GetName())->Draw();
  
  // nを使わないバージョン 
  gDirectory->Get(gDirectory->GetListOfKeys()->After(gDirectory->GetListOfKeys()->FindObject(gPad->GetListOfPrimitives()->At(1)->GetName()))->GetName())->Draw();
  // 長い。代わりに現在表示されているhistを取得してその次を調べているので
  // hist1->Draw(); と通常のやり方でDrawした後でも使える。

• 現在のヒストグラム番号nを取得する

  gDirectory->GetListOfKeys()->IndexOf(gDirectory->GetListOfKeys()->FindObject(gPad->GetListOfPrimitives()->At(1)->GetName()))

ヒストグラム上でマウスの位置の座標を取得する†

gPad->GetEventX() // TCanvas上での現在マウスポインタがあるx座標(多分ピクセル単位)
gPad->GetEventY() // TCanvas上での現在マウスポインタがあるy座標(多分ピクセル単位)

gPad->AbsPixeltoX(gPad->GetEventX()) // ヒストグラム上の座標

• これらの単独利用は現実的ではない。で、「何かしらのイベントがあったときに実行する関数」を設定する

gPad->AddExec( "myexec" , "myexec()" );
void myexec(){
  printf("%d %.1f %.1f \r",gPad->GetEvent(),  gPad->AbsPixeltoX(gPad->GetEventX()), gPad->AbsPixeltoY(gPad->GetEventY()));
}
gPad->DeleteExec( "myexec");

• gPad->GetEvent()の返り値は https://root.cern.ch/doc/master/Buttons_8h.html を参照。クリックしたときは1。
• 参考
  • https://root-forum.cern.ch/t/get-histogram-coordinates-with-mouse/21307
  • https://root.cern.ch/doc/master/mandelbrot_8C.html
  • https://root.cern.ch/doc/master/classTPad.html

現在設定中のx軸(SetRangeUserなどで設定した)の範囲を取得する†

hist->GetXaxis()->GetBinCenter(hist->GetXaxis()->GetFirst())
hist->GetXaxis()->GetBinCenter(hist->GetXaxis()->GetLast())

TPadに表示されているObjectの一覧を表示†

gPad->GetListOfPrimitives()->ls()

TCanvasの分割関係†

• 今のキャンバス

  gPad->GetCanvas()

• 今のキャンバスを何もなくす

  gPad->GetCanvas()->Clear()

• 今のキャンバスを横3 x 縦2に分割する

  gPad->GetCanvas()->Divide(3,2)

• その中のn番目のパッドに移動する

  gPad->GetCanvas()->cd(n)

• 分割数を取得するのは容易ではない
  • https://root-forum.cern.ch/t/get-number-of-divisions-in-a-tcanvas/41495
  • 多分一番確実なのは gPad->GetCanvas()->cd(n) して、する前gPadのポインタと変わってなければそのn番目のpadは存在しない、と判定できる。
  • gPad->GetCanvas()->GetListOfPrimitives()->GetEntries() はそれっぽい値を返すが、状況によっては違う値を返すらしい。

gDirectory, gPad, gROOT†

• これらについての説明はそのうち書く

ANAPAW風のコマンドを使う†

• ROOTのコマンドは長ったらしくて打ち込むのがめんどくさい。 それに対し、ANAPAWは"ht"でヒストグラム表示、"hn"で次のヒストグラムを表示...などコマンドがシンプルで覚えやすい。 これをROOTで使えるように実装するマクロを開発中。-> こちら ANAPAW Commands for ROOT
  • 参考 : PAW/ANAPAW コマンド - Takayuki Sako Webpage

なにかに使えるかもしれない†

• プロンプトを変更する

  ((TRint*)gApplication)->SetPrompt("hoge [%d] ");

  • root [0] ってところを任意の文字列に変えれる

うるせえ! そんな情報いちいち表示すんな! というとき†

Info in <TCanvas::Print>: pdf file c1.pdf has been created などを表示しない方法†

• ROOTでは情報表示の重要度レベルに応じてInfo, Warning, Errorなどが設定されてる。
• gErrorIgnoreLevel というグローバル変数に表示レベルを設定することで情報表示の冗長性を調整できる

  #include <TError.h>
  gErrorIgnoreLevel = 2000;

• Infoを表示しない場合は2000 (Warning以上のみを表示)、Warningを表示しない場合は3000 (Error以上を表示)

• 参考
  • How to set output verbosity level in ROOT? - ROOT - ROOT Forum
  • gErrorIgnoreLevel problem - ROOT - ROOT Forum
  • ROOT; core/base/inc/TError.h File Reference

Fit関係†

• Fitオプションに Q を指定することで Quiet mode になる。(端末上にフィット結果を表示しなくなる)
• N を指定するとキャンバスにフィット結果の関数を描画しなくなる。

Topへ戻る