memo385 Forward <= Today => Previous memo382

run77でも残念ながらゲインシフト。

*timelikeスペクトルで16Mil以降が細い問題について => memo384 (解決)
*QDC, TDCのデータからイベントIDをぶっこ抜く      => memo384

*timelikeスペクトルの謎

timelikeに用いている変数と、実際の解析済みイベント数との関係：

16000000以降で伸びる。

=> timelikeとしてreal(4byte)を用いていることが原因。
=> real*8に変えた図：


DAQがバグっているということはない。安心。
この修正をanapawに適用してスペクトルを採る：

きれいである。

*rawデータからイベント数ぶっこ抜き

・rdf[LINK]そのものにはイベントカウントは無い。

->では個々のモジュールからイベントIDを得る?
->それではモジュールごとの足並みが揃わない。

memo370より
GEO=11のV775(UD101-115)のデータ:
5a00 4549 580d 45fb 581d 5708 5c96 5707

d11 = b01011
-> 01011 + 010 = b0101 b1010 = x5 xa
-> 01011 + 000 = b0101 b1000 = x5 x8
-> 01011 + 100 = b0101 b1100 = x5 xc

=>
5a00 4549 :head
580d 45fb :data  x0d = b00001101 -> 13ch => U114
581d 5708 :data  x1d = b00011101 -> 29ch => D114
5c96 5707 :end

V775では、enderの右24ビットがイベントカウンタになる。
-> x965707 = d9852679.

GEO=12(UD201-215)ではどうか?
->
6200 46d9 6000 478e 6010 46ec 6002 477e 6012 5709 6496 5708

d12 = b01100
b01100010 = x62
b01100000 = x60
b01100100 = x64
->
6200 46d9 :head
6000 478e :data
6010 46ec :data
6002 477e :data
6012 5709 :data
6496 5708 :end

x965708 = d9852680

*イベント数ぶっこ抜き 実践編

dec_vme.f にoutput用パラメータ iev を付加
--------------------------------------------------------------
c     end of block ! bit24-25=0 -> end of block
         ElseIf(ia .EQ. x'04000000' .AND. evtflag .EQ. 1)Then
            evtflag = 0
            iev = iand(idata(i),x'00ffffff')
         EndIf
      EndDo
--------------------------------------------------------------

iev / 1e6 in real*8 を61, 62番(TDC, QDC用)として使用。

使用例：


eventcountが24bitしかないので、2**24 = 16777216~約17Milイベントでループする。
2Mil付近が強いのはそのため。

*ゲインシフト

D213をLEMO-Tで分配して波形を見る。