進捗報告 M2 下村 140614 地上系会議.

Presentation on theme: "進捗報告 M2 下村 140614 地上系会議."— Presentation transcript:

1 進捗報告 M2 下村 地上系会議

2 内容 MDCSimulatorのMerge（Auxiliary） QLOrbISSの作成 ToDo

3 MDCSimulatorのMerge（Auxiliary）
動作テスト sock_ifでのデータ受け渡し：ok 力石ルーチンから取り出したLat,Lonと入力した値を比較 　→ 小数第6位の精度で一致 (altは第3位) 入力　Lat　　　　Lon　　　　Alt 出力　Lat　　　　Lon　　　　Alt < < < < < < < < < < < < < < < < < < < < > > > > > > > > > > > > > > > > > > > >

4 QLOrbISSの作成 QLTrendを流用して、高度400kmでのCALETの位置のモニターを作成した
機能：CALETの位置 　　　通過した軌跡（default:1000sec） 　　　通過したGridのRateをプロット（画像はSingle）

5 Todo QLOrbISSの改善 QLの全天モニターを作る 軌道予測の追加（SGP4）→探したが上手く行かず カウントレートの選択
カラースケールの追加 QLの全天モニターを作る

6 END

7 ToDo QLの作成案を練る Level0データの検証

8 Validate Level0 TcpServerEXMainを流用してLevel0データの検証用プログラムを作成した
UKey2/EventIDの連続性、チェックサムの確認等 実行例 Level0データ エラーがあるとここ にコメントが出る Telemetry Data File:HighLat/CR_MJD55501_00000_00050.bin.ics TcpServerEXMain ... secDiff= [CheckSum OK 0000] [CheckSum OK 0000] [CheckSum OK 0000] [CheckSum OK 0000] [CheckSum OK 0000] GSEヘッダ データ部 GSEヘッダ データ部 GSEヘッダ データ部 ・

9 QLSummary修正 ボタンを追加＋見栄えの修正

10 CALET ISS軌道データ作成 (EPICS)
/nfs/RAID/home/asaoka/calet/TelemetrySim にて環境変数をセット $source setenv.sh [Data set] /nfs/RAID/home/asaoka/calet/TelemetrySim/CALET_ISS にてISS軌道CALET 入力データを生成 $./telesim [-d MJD Date] [-s Start Sec] [-p Proc Sec] [-n Sec Div] データ生成ディレクトリ ~/work/TelemetrySim/data/CRSample/[Data Set]/M30/ 　→ CR_MJD[Date]_[Sec]_[Div].dat 入力データとして、データ・セットに対応するATMNC3データが必要。 　→ ~work/TelemetrySim/data/CRMap参照 /nfs/RAID/home/asaoka/calet/TelemetrySim/CALET_Ana/trunk にてEPICS データを生成 $./epicsSGE.py [-d Data set] [-j MJD Date] [-s Start Sec] [-p Proc Sec] [-f Sec Div] [-r Retry] データ生成ディレクトリ ~/work/TelemetrySim/data/EPICS/[Data Set]/M30/ → CR_MJD[Date]_[Sec]_[Div].root

11 ソケット/QL プログラム修正 ソケットプログラムのオプション追加 QLEvDisp、QLSummary、QLSumTextのルーチンの改善
“-c”:パケットに記述された時間を付け替えない QLEvDisp、QLSummary、QLSumTextのルーチンの改善 問題点：実行時に激しいメモリリーク 原因　：UniqueKey3が1,6以外のデータで、MDCSimulator のインスタンスをdeleteしてなかった 修正後：激しいメモリリークは解消 　　　　ただし、微妙に残っているので今後の課題。

12 MDC Simulator 修正 Writeルーチンの変更の必要性 Readルーチンの修正 (140517)
MSSのサンプルデータの解析によりシミュレーションデータに以下の相違点が判明 イベントデータ2ワード目以降の13ワード目：ユーザシーケンスカウンタの削除 CRCを削除 フラグメントの最大、最小がNASA/ICS経由で違う => 最小以下の場合、0埋めする必要あり NASA : Max Packet Size = 8(CCSDS Header) (内Event Data 730 words), Min Packet Size = 8(CCSDS Header) (内Event Data 30 words) ICS : Max Packet Size = (内Event Data 736 words), Min Packet Size = (内Event Data 12 words) Readルーチンの修正 (140517) Eventデータ ラストフレームのWord数 先頭フレームのTotal Word Lengthを使い導出

13 MDC Simulator 修正 テスト NASA/ICS経由のそれぞれについて10秒分データを作成
QLで確認 => データは問題なく読めた 　 過去のデータとの整合性は要検証 TcpClientQLMainのReadオプション “r”= 0 : 修正済みシミュレーションデータ “r”= 1 : MSSのサンプルデータ “r”= 2 : 過去のシミュレーションデータ

14 QL追加 QLSummary, QLSumText 定時データ、イベントデータのサマリを表示 定時データのサマリはMPCでのパルステスト用ソフトウェアのルーチンを流用 使用方法 TcpClientQLMainのオプションで、 “-S” : QLSummaryを起動 QTのWindowにサマリを表示 “-Q” : QLSumTextを起動 コマンドラインにサマリを表示 　　　　 キーボードによる入力受付

15 QLSummary 使用例 Very preliminary!

16 EPICS CADモデルの使用 現状 現在使用中のEPICSのGeometryもCADモデルである
-Bug -"EDepAll(_prim, Eff, Eff_prim)" were NOT initialized in each event. This is fixed. -Add -calculate the component number of first interaction point, and ouput it. -Update-1 -use "epq1stIntTrack" function(Epicsv9.161~) to find 1st interaction point -Update-2 -unify the numbering method of component number. the start number is "0", not "1" (ex.) SciFi number: (0~447) PWO-log number: (0~15)

17 課題リスト（140411） QLのHit情報バグ（途中まで調査） 100TeV問題 CRCの取り扱い
MDC Simulatorのメモリ管理（OK） データ生成　=>6周分で停止 広平面入射のHEトリガーレートが合わない トリガー効率の確認

18 QL Hit情報調査 Localで調査 git clone IEEE, GSE(ICS,NASA)ヘッダについてそれぞれ調査
最新の状態でテスト環境構築 test.sh実行でエラー発生(新環境) →バージョンを戻して確認 $ git checkout ‘ハッシュ’ (Date: Tue Jan 7 21:55: ) DeadTimeの取り扱い変更が原因、ソースに追記 →もう一度test.shを実行し、現行verの結果とdiffにて確認 最新のテストデータにアップデートする必要あり IEEE, GSE(ICS,NASA)ヘッダについてそれぞれ調査 HEでトリガーされた全イベントに対し、LD TASC HE Hit情報が0になっていないか確認 50秒分のデータをそれぞれ作成 →どちらもなっていなかった WCOCのマシンでも確認、ok →QLではおかしいので、QLのアルゴリズムの確認

19 QL Hit情報調査 QLのアルゴリズムの理解 内容
MDC Simulator 内部でパケットを読むルーチンが発見できなかった →解決済み（140419） 内容 ICS,NASA経由のデータを50秒分それぞれ生成 EPICS　all　55501　psingle QL実行 CHDのSingle,IMCY4のSingleがおかしい QLのインスタンス内でHit情報を確認 →上記バグが発生した MDCのインスタンス内でHit情報を確認 →上記バグが発生しない やった位置 ThreadClientRelay.cpp ReadGSE()の直下◯ CloneInstance()後のmdcQueue☓ MDCSimulator.cc CloneInstance()中☓ CopyValuables()中◯ QLEvDisp.cpp TQLDispDet.cc InitDraw()中☓ MDCSimulator::CloneInstance()で上手くインスタンスがコピーできてない 以上を踏まえ、原因としてmdcのインスタンスを渡す際にメモリを誤開放していることが考えられる。 使うネットワークによって、GSEヘッダが上手く読み込まれないときがある →帯域が問題 自前のMobile Wifi経由で発生 研究室内ネットワークでは動く

20 PWO 平均Energy Deposit PWO一本一本のΔEのレートを求める（軌道1周） 方法
使用データ：MJD55501、EPICS all particles teleout.ccを改造し、各PWO毎に1秒毎のΔEのsumを出力 getPwoEnergyDeposit.cc 出力を確認.　（ :06） 100TeVの粒子の影響により、甚大なΔEがちらほら確認できる →1TeV以上のイベントを一旦cut 5406秒分のデータを作成 レートの絵を描く GetPwoEDepRate.C MIPの値は2GeVのミューオンのΔEより（赤池さん） GeV / 1PWO 192本まとめてDraw

21 QL Hit情報調査 今日の目標：QLのアルゴリズムの理解 内容 まずICS経由のデータを50秒生成
EPICS　all　55501　psingle

22 Hit情報バグ QLでイベントビューアを呼び出す際に、トリガー情報と、LDHit情報が一致しない 原因調査
使用データ：EPICS all particle, 周分のデータ, Header:IEEE telereadによりイベント情報を目で確認 　　＝＞特に異常なイベントはなし HEでトリガーされた全イベントに対し、LD TASC HE Hit情報が0になっていないか確認 　　＝＞全て1になっていた。 　　　　　ここで確認した変数をQLで読み込んでい るので、バグの原因ではなさそうです。