ANAROOT/要望

ANAROOT/要望について、ここに記述してください。

Contents

ANAROOT/要望

warning messesage

エラー等はく際は
- Info
- Error
- Warn (現段階ではDebugと同等)
- Debug
を用意しているが、Infoが多くでて来ると、start()した後にうざい。
その為にWarnを準備しているが、どこかにメッセージをつめておいて、
dmesg()
みたいな事をしたらメッセージをdumpする仕様にしようかと思っている。

getxml.pyの置き場

解析PC(c02)でgetxml.pyとかの置き場ですが、

$HOME/exp/scripts

でいいですか？HOMEが汚くなるのが嫌ですし。面倒なのでそこにPATHを通しました。

TArtDataObjectのSetDetectorNameについて

nai->SetDetectorName(para->GetDetectorName());

みたいに"*"演算子なしで受け渡ししたいので、

void SetDetectorName(TString *val){name = *val;}

をつけ加えませんか？ (2012 2/28 田中)

TClonesArrayのClear

TCLonesArrayでメモリリークが起こる理由がわかった気がする。公式をみると、

TClonesArray::Delete
 Clear the clones array. Use this routine when your objects allocate
 memory (e.g. objects inheriting from TNamed or containing TStrings
 allocate memory). If not you better use Clear() since if is faster.

と書いてあり、TStringsをクラス内で使っている場合はClearではなくDeleteを使えとのこと。基本的にClearではなくDeleteにしたほうがいいかも。 (2012 2/27 田中)

bug情報

これから先バグ報告スレッドがあると便利かも？あと、プログラムが落ちた場合は、

The lines below might hint at the cause of the crash.
...
that might help us fixing this issue.
===========================================================
#12 0x064ccc59 in TCanvas::RaiseWindow() () from /usr/local/root.5.32/lib/root/libGpad.so
...
#16 0x07a3eae6 in hn () at src/TArtAlias.cc:19
===========================================================

の部分を貼っつけてください。あとはどうやったら落ちたとか再現性とかまで言っていただけるとdebugが楽です。 (2012 2/25 田中)

Nadekoについて

rootを立ち上げていきなりhtp()すると

 *** Break *** segmentation violation



===========================================================
There was a crash (#7 0x0030e6cd in SigHandler ()).
This is the entire stack trace of all threads:
===========================================================
#0  0x00110416 in __kernel_vsyscall ()
#1  0x03d59953 in __waitpid_nocancel () from /lib/libc.so.6

で鬱陶しい。まあこんなことしないのだけど何にもヒストグラムねぇぞコラ的なメッセージを出してほしい。（2012 2/24　近藤）

--> バグ報告ありがとうございます。調べたところ、ポインタを初期化していないことにより起こるバグでした(なのでhnとか他のコマンドでも同じバグが潜んでいる)。近いうちに修正します。 Nadekoはいつどんな時にどのコマンドを使っても問題が起きないという設計をしているつもりですので、Nadekoのコマンドで落ちたらご連絡ください。何もhistが無いときレベルのことに対してメッセージを出すとなると結構な量のメッセージを書く必要があるので多分やりません。 (ちなみに、hn()でも一番最初は1個目のhistが表示されます。) (2012 2/25 田中)

branches/2012SAMURAI.tanaka では更新しました。 (2012 2/25 田中)

UIコマンドのページ

Startup Guideみたいな簡単な例があるといいね。（2012 02/13近藤）

UIコマンド以下にUI周りのドキュメント作りました。 (2012 2/19 田中)

スタートアップガイド、Nadeko、Analoopごとのページにするのが良いかも？少なくともページのコンテンツがどのライブラリ（Nadeko、Analoop）に関するもなのかは明記してほしい。（1012 2/25　近藤）

最適化

(2012 2/9 田中)FindなんちゃらでPlaとかParaを探すときに現在はシーケンシャルに検索しているけれど、これを、あらかじめmappingしておいて一発で検索出きるようにしたところNEBULAのencで大体二倍速になった。具体的にはTClonesArrayのnewをするところで同時にIDをkeyとし、何番目に入れたかやparaを値としてstd::mapにつめます。

 fIDNPlaMap.insert(std::pair<int, int>(para->GetID(), npla));

で、Findなんちゃらではこのmapから検索します。

 std::map<int,int>::const_iterator it = fIDNPlaParaMap.find(id);
 if(it != fIDNPlaParaMap.end()) return GetNEBULAPlaPara(it->second);
 else return 0;
(ClearDataでstd::map.clear()を呼のを忘れずに。)

Findとつくものに対して全てmapを作って上げれば完成です(NEBULAの場合は3箇所)。気が向いたら実装してみてください。

環境設定について

(2012 2/8 田中) thisroot.shはスクリプトがあるディレクトリで実行しないとおかしくなるので、僕は.bashrcに、

# root
export temp_dir=$PWD
cd /usr/local/root
source bin/thisroot.sh
cd $temp_dir
unset temp_dir

# anaroot
export TARTSYS=$HOME/anaroot
export LD_LIBRARY_PATH=$LD_LIBRARY_PATH:$TARTSYS/lib

こんなのを書いています。いずれディレクトリ構造や昔のanaloginみたいなのも作る必要があると思うので、そのときにこの問題も合わせて検討ください。

-> cdをcdlsという関数に張り替えているのが原因でした。普通の場合は

. /usr/local/root/bin/thisroot.sh

とかで問題ないのですね。お騒がせしました。

analoginみたいなもの

僕はこんなスクリプト使っています。（近藤） makeanadir.sh

(2012/2/14 田中) analoginをそのままにシンプルにしてみました(誤爆対策でy/n付き) analogin.sh (変更したいけど権限がないといって怒られる...)。

これを回す前に、

exp/commisioning_samurai
exp/commisioning_samurai/src(anaroot_vX.X.Xへのリンクを貼ることになる気がする)
exp/commisioning_samurai/users
exp/commisioning_samurai/users/default(新しいユーザーを作るときはこれをまるごとコピることになる)

を置いてもらい、あとは.bashrcに

alias analogin='source hogehoge/analogin.sh'

と書いてもらえば完成です。

TArtNEBULAについて

名前の付け方について

(近藤）

TArtNEBULAのfEnergy（ADCを較正したもの）はfEdepositとかにしたい。
TArtNEBULAParametersのqcalというパラメーターはqslopeとかにしたい。
TArtNEBULA（名前変えたい）のメンバ変数は何が必要？
- id,layer,ADCu/d,TDCu/d,Tcal(u/d),Tave,Aave,dT,x,y,z,...
- TOFを作るには別のplasticの時間が必要だけど...
今の感じだと例えば1つの検出器がどの層かを知ろうと思うと

  TArtNEBULA* pla = mynebula->GetNEBULA(ihit);
  Int_t layer = pla->GetPara()->GetLayer();

と面倒かつ直観的でない。layer=pla->GetLayer()で行けるようにしたい。ただしlayerはxmlで（エクセルで）書いておいたほうが編集しやすいのでTArtNEBULAで同じlayerというメンバを持っておいてparaからコピーするのか？2重で無駄？

->(2011 2/9 田中) 僕のページに新しく上げたソースでなんとなく(positionのoffsetは直書き)実装しました。名前はとりあえず自分規約で付けてます。

TArtCalibNEBULA

NEBULAの解析では、最終的にはソフトでthresholdを設定して、それを越えたものだけを解析したい。その場合にTArtNEBULAを継承したものを使うのか？それともTArtNEBULAPlaにそういうフラグみたいのを加えるのか？Online解析では単にゲートを掛ければよいけど。（2012/02/13　近藤）

TArtDALIについて

anaroot_v4.3.1についてです。(近藤)

DopplerCorrectはtofが引数になっているが、βとする方がわかりやすいと思います。しかもconstで良いか。
addback解析とかとつなげられるようになっていると良いか？TArtDALIArrayをとってきてマクロを書けばできるのでいらないかも。
PPACのTXsumの様な、ある時間ゲートに入ったものだけを真のイベントとみなして解析するようなものも必要（とりあえずは優先度低め）

このページへの要望

更新されてもしばらく気が付かないので、RSSとかで更新が出ると嬉しい。

-->　RSSはじめました。ANAROOTトップページのRSSアイコンをクリックしてください。出来はとってもいまいちです。更新履歴をRSSに変換しただけで、どのページが更新されたかがわかる程度で、ページのどこが更新されたかはRSSからは全然わからない。。。（2012 2/21近藤）

(2012 2/8 田中) 言うだけで何もしないですが、左側に小林さんのページみたいにリンク一覧があるとうれしいです。あと、何かしらの規則見たいな物があると書きやすいかも。

(2012/02/12 近藤) とりあえずSideBarつけました。

ERROR Messageについて

mapファイル（nebula.map）がなくてもデータを読み込むが、detidとかはもちろん正しくない。マップファイルがない場合にはエラーを表示してほしい。　（近藤）

回答？：mapファイル経由だとdevice,fpl,detector,geo,chから一度category detidをマップし、その後さらにADCへとの代入になる。ANAROOTではdevice,fpl,detector,geo,chから直接ADCへマッピングすることで計算速度を向上させようとしている。（このマッピング情報はキャリブレーションパラメータと同じxmlファイルに格納されている。）なので、mapファイルの使用はオプションと想定している。つまりユーザー側でエラー処理をして貰いたい。if(detid==-1)false;みたいな感じで。

（近藤）TArtCalibDALIのLoadRawDataでは書いてありますね。

        TArtRIDFMap mm(fpl,detector,geo,ch);
        TArtDALIPara *para = (TArtDALIPara*)setup->GetDALIPara(&mm);
        if(NULL == para){
          TArtCore::Info(__FILE__,"Could not find TArtDALIPara...: Dev:%d, %s", device, mm.GetMapInfo());
          continue;
        }

(2012 2/9 田中) TArtCalibNEBULAでは一時的にerror処理を切っていました。 loop内でエラーが発生した場合に無限に出続けて困るので、例えば、、、(以下いろいろ問題ありだったので消しました。)

TArtNEBULAについて

名前の付け方について

(近藤）

TArtNEBULAのfEnergy（ADCを較正したもの）はfEdepositとかにしたい。
TArtNEBULAParametersのqcalというパラメーターはqslopeとかにしたい。
TArtNEBULA（名前変えたい）のメンバ変数は何が必要？
- id,layer,ADCu/d,TDCu/d,Tcal(u/d),Tave,Aave,dT,x,y,z,...
- TOFを作るには別のplasticの時間が必要だけど...
今の感じだと例えば1つの検出器がどの層かを知ろうと思うと

  TArtNEBULA* pla = mynebula->GetNEBULA(ihit);
  Int_t layer = pla->GetPara()->GetLayer();

->(2011 2/9 田中) 僕のページに新しく上げたソースでなんとなく(positionのoffsetは直書き)実装しました。名前はとりあえず自分規約で付けてます。

TArtCalibNEBULA

TArtDALIについて

anaroot_v4.3.1についてです。(近藤)

DopplerCorrectはtofが引数になっているが、βとする方がわかりやすいと思います。しかもconstで良いか。
addback解析とかとつなげられるようになっていると良いか？TArtDALIArrayをとってきてマクロを書けばできるのでいらないかも。
PPACのTXsumの様な、ある時間ゲートに入ったものだけを真のイベントとみなして解析するようなものも必要（とりあえずは優先度低め）

このページへの要望

更新されてもしばらく気が付かないので、RSSとかで更新が出ると嬉しい。

(2012/02/12 近藤) とりあえずSideBarつけました。

ERROR Messageについて

（近藤）TArtCalibDALIのLoadRawDataでは書いてありますね。

        TArtRIDFMap mm(fpl,detector,geo,ch);
        TArtDALIPara *para = (TArtDALIPara*)setup->GetDALIPara(&mm);
        if(NULL == para){
          TArtCore::Info(__FILE__,"Could not find TArtDALIPara...: Dev:%d, %s", device, mm.GetMapInfo());
          continue;
        }

ANAROOT/要望

Menu

ANAROOT/要望

warning messesage

getxml.pyの置き場

TArtDataObjectのSetDetectorNameについて

TClonesArrayのClear

bug情報

Nadekoについて

UIコマンドのページ

最適化

環境設定について

analoginみたいなもの

TArtNEBULAについて

名前の付け方について

TArtCalibNEBULA

TArtDALIについて

このページへの要望

ERROR Messageについて

TArtNEBULAについて

名前の付け方について

TArtCalibNEBULA

TArtDALIについて

このページへの要望

ERROR Messageについて