2013 May 28 10:00 Detchar Meeting 端山、譲原、山本、宗宮、関口、伊藤 アジェンダ |Announce |• 6月10日、11日 第4回Japan-Korea workshopでのKorea detcharグループとの議論について。 [端山説明] 現在進行中の共同研究トピック(添付資料参照)大きくわけて4つある。 (1) 韓国側がやりたいと言っているevent detection pipeline(韓国側がメインのプロジェクト) ・多チャンネル解析によるイベント解析とそれによるCBCイベントのVetoを考えている。 ・補助チャンネル解析のためのglitch探索パイプライン( = バースト探索の単一検出器版) を新規開発したい意向。 (2) future selection statistics(韓国) ・データから情報を要約する統計量の研究 ・彼らはt or z-statisticsと読んでいるが、何を使っているのか不明。 (3) 雑音ともっとも相関しているチャンネルを探す。(日本) ・noise classificationと呼ばれる分野。 ・多チャンネル解析によってノイズ源を特定する。 (4) 以上の手法でCLIO Hardware Injectionを、データのみから検出できるかどうかを調べる。 (5) 他に、LVCデータを用いて、LVCでノイズと分かっているイベントを検出できるかどうかを テストする。 日韓ワークショップ中のdetchar ミニワークショップ 今のところ10日13:00−15:00開催の予定だが、変更の可能性あり。 韓国側がHHTを用いて、山本君が作ったCLIO h(t)データを解析した結果を韓国側が示す。 韓国側が、以前LVCで使われていたKWトリガーと、今回のHHTトリガーのパフォーマンス比較を話すとのこと。 ただし、LVCではKWはすでに使われていない。Omicron( = Omegaの高速バージョン)が使われている。 |Special Topic |• LVCのメジャーライブラリへのKAGRA detcharび主体的貢献の可能性の検討(Haskellスタディ)(山本、端山) [端山説明] ・現在検討段階の事項。KとLVCでデータ共有された段階のことを考える。 ・メインで利用されるであろう解析ソフトはLALだろう。 ・日本からの貢献がまったく無い状況で共同研究をすすめるのは避けたい。 ・detcharとして貢献できるのは何か? ・スクリプト(Matlab, Python-PyLAL, ...)に対するLALラッパーの開発で貢献するのはどうか検討する。 ・LALが利用しているソフトウェアについては、端山資料参照。 ・Detcharに関係するのは、DMT( 資料2ページ目だとエンドユーザーソフトウェア ) ・ラッパーを通してC, C++ライブラリをMatlab/Pythonなどから利用するケースが非常に多い。 ・SWIG:ラッパーを構築するソフトウェア http://www.swig.org/ ・Haskellへのラッパーはどうか検討中。 ・京大の村主さんがHaskellのエキスパートで、興味を持ってくれている。 ・柴田大研究室とCBCテンプレートのHaskellによる実装について研究している。 ・大阪市立大でHaskell勉強会をしている。以下山本成果報告。 [山本報告] ・Haskellの特徴:純粋な関数型言語 [端山] 他人が書いたHaskellは読みやすいのか? [山本] CやC++よりも、簡潔に書けるので、読みやすい。 [端山] Haskellでないとできないこと、実装が大変なこと、ソフトウェアとして実装しにくいものが あると良い。 [端山] 遅延評価とは何? [山本] 計算が必要になったときに初めてその評価をおこなう。 すべてのデータを常に評価するわけではない。計算量の低減ができると期待している。 [伊藤] ・lispとの違い。 ・Haskellのユーザー数と将来性。 ・遅延評価は、Mathematicaでもできるが、何が違うのか? ・将来重力波データ解析も大量のデータを扱う。 大量のデータを効率的に処理する機構があるのか? [端山] ・lispよりも速い。(Cの1〜1/3倍程度の速さ) ・将来性は不明。 ・Haskellにどれくらいマンパワーを割くかも含めてまだ検討段階。 |• バイオリンモードスタディ(端山) [端山]資料提示、現状説明: ・それほど進展は無し。 ・関口が以前計算した感度曲線(緑)は、100Hz以下がviolin modeで支配されてしまう。 ・これは関口の以前のパラメータと宗宮のパラメータの相違による。 ・関口がQを3x10^5にしたときに計算してくれたのが赤い線。 ・鏡の大きさが宗宮オフィシャルよりも少し小さくなった影響が残る。 ・パワーを下げることはしていない。 [宗宮] 宗宮が やったのは、鏡の大きさが少し小さくなってパワーを下げること。 加えてコーティングロス、タングステンのロスなど何点か変更している。 それで230Mpcぐらいで、CQGに発表しているもの。 パワーを下げないと240Mpcぐらいになる。 宗宮・関口の計算している感度は、 今はパワー250で計算している。400までできるが、山元さんのところでどれくらい 冷やせるかによる。400なら240Mpcまでいける。 [端山] 宗宮からもらったMatlabスクリプトをそのまま使って 鏡を小さくするだけでよいのか? [宗宮] 今はそれで良いが、新しいものができたら送る。 [端山] 関口が当初計算したviolin modeが4 ~ 16に分裂する話。 [関口] 通常の半分のテンションの場合、通常の2倍テンションがかかった場合で計算した。 100数十Hz - 260Hzぐらいまで分裂する。 1つ1つ分裂を計算するのは大変なので、ピークをフィッティングしてもらって、ピークの取りうる範囲内で 適当に分裂させてもらった。 [宗宮] 赤い1本と、分裂した青いライン1本の振幅の関係は良いのか。 [端山] 4本に分裂した場合、4分の1にしている。 [端山] 16分裂すると、関口に許された範囲で適当にパラメータを振ってもS/Nは あまり変わらない。219Mpcぐらい。 [宗宮] 110Hzの縦振動変わらない? [関口] あまり変わらない。 ファイバーの太さにばらつきがある場合に変わるが、現状ファイバーの太さなどは分からない。 [端山] 100Hz - 300 Hzのviolin modeはバーストに影響を与える。 感度曲線と超新星爆発の特徴的周波数(数値計算結果)を示した資料提示。 [関口] 超新星爆発ということはバーストに近いもの。 ブロードバンドなのでは? [端山] 特徴周波数を示している。 検出のときに一番影響を受ける周波数を示している。 リコイルマスのQ値がどんな値になっているのかが心配。 [宗宮] 16本分裂は非常に最悪のケースだと思われる。 4本分裂ですでにあまりないのではないか? 分裂しても、4本の周波数は関口計算よりも近いところに 分布するのではないか? [関口] ワイヤーがものすごく硬いのでワイヤーの長さを数ミクロンであわせないと 周波数空間でかなり広がる可能性がある。現在の技術では、関口計算結果ぐらいは広がりうる。 [端山] 16本はありうるのか? [関口] ありうる。 [宗宮] 問題提起すべきである。 縦バネなど、解決策は考えられるはず。 [関口] 縦バネを入れることで改善する可能性ある。 [宗宮] 縦バネを入れることができない場合のことを考えて、 2本吊りも検討されているが、どうか? [関口] 2本吊りだと、ピッチのコントロールをすべてテストマスでやらなくては ならなくて、テストマス上のアクチュエータですべてやる必要がある。 アクチュエータのパワーを考えると難しいのでは。コントロールがやっかい。 [宗宮] 分裂だけで見ると、2本吊りの方が良いのではないか? (以下議論)