2005 年度
P6 課題研究
マルチアノード光電子増倍管とヘッドアンプユニットを用いた動作実験
豊田武
安井健目次
1

目的と概要

22

実験で使用する装置について

32.1
2.2
2.3
3

ヘッドアンプユニット . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

34
5装置のセットアップについて
3.1 一連の装置の接続 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

66

PC でのデータ取得までの流れ . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

73.245シンチレータ . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
マルチアノード光電子増倍管 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

実験・解析

84.1
4.2前段階実験・自作プリアンプ . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

1ch Na22 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

894.34.44.516ch Na22 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
256ch x 30000 event の波高値 . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
256ch 分 Na22 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

1011124.64.7

256ch 分 Co57 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
256ch 分 Ba133 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

13144.84.9

Background について . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
キャリブレーション . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

1516

結果・まとめ5.1

これまでのまとめ . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

17
176

謝辞

187

参考文献

191

1

目的と概要
今回の実験では、マルチアノード光電子増倍管∗ (浜松フォトニクス製 H9500) にヘッドアンプユニット（読み出し回路）を直接とりつけデータを取得して処理する。
このヘッドアンプユニットは上記の MAPMT 用なので、接続が容易であり、コンパクトである。

MAPMT に接続したものを二次元平面に敷き詰めたりすることが可能である。†
また同時に、実験を通して回路を自作したり、解析手法を学ぶ。

∗ Multi Anode Photomultiplier 略：MAPMT† 「取扱説明書 80058 型 MAPMT ヘッドアンプユニット」(CLEAR-PULSE 社 2005) より

2

2

実験で使用する装置について

2.1

シンチレータ

荷電粒子は物質と相互作用すると、その物質の分子を電離・励起する。励起された分子は基底状態に戻るが、励起エネルギーが熱エネルギーに変換されてエネルギーを失い、基底状態に戻ることが、多いが、シンチレータは基底状態へ戻る際にエネルギーの一部を光として放出する。励起されて光を放出するまでの時間が短いものが放射線の検出に利用されている。
今回は CsI(Tl) シンチレータを用いた。CsI(Tl) の特徴としては原子番号、密度が NaI よりも大きく、γ線に対する吸収が大きく、検出効率は大きいが、シンチレーション効率は NaI の約半分である。潮解性はない。蛍光減衰時間は１μ sec である。
今回用いたのは縦 50mm × 横 50mm × 厚さ 2mm の大きさのものである。
結晶

蛍光減衰時間蛍光量最大波長

相対蛍光効率

密度NaI(Tl)

230ns

415nm

100%

3.67cm3

CsI(Tl)1000ns530nm

45%

4.51cm3

表 2.1: NaI(Tl) と CsI(Tl) の比較

図 2.1: CsI(Tl) シンチレータのイメージ (縦 × 横)

3

2.2

マルチアノード光電子増倍管

光電子増倍管はシンチレータからの光によって光電子を発生させる光電陰極、発生した光電子を効率よく第一ダイノードに集める収束電極、光電子の数を増倍させるダイノード電極、増倍された電子を集める陽極からなる。
今回使用するマルチアノード光電子増倍管は、256ch に対してそれぞれアノード・増幅部が対応し多数の信号の読み出しが可能になっている。

図 2.2: MAPMT の各 ch の配列

また、今回使用するマルチアノード光電子増倍管のダイノードは 12 段である。
ゲインは、各 ch ごとで異なっていて、中心付近の ch に大きいなゲインの ch が揃っている。
最小のゲイン：最大のゲイン=1:5.5 の比率の差がある。
ただ、今回は各ゲインの違いについて考慮する段階までは至らなかった。

電流増倍率 (ゲイン)

陽極暗電流1.05 × 102.85nA

6

表 2.2: 今回の MAPMT:H9500 のゲインと暗電流†増倍率 (ゲイン) につて
光電子増倍管の増倍率を M とすると、M = δ n で与えられる。δ は二次電子放出率、n はダイノード電極の段数。例えば、δ = 4, ダイノード 10 段のとき、M = δ n = 410 ≃ 106 。†
陽極暗電流について光電子増倍管を作動状態にしておくと、光が全然入らないくても陽極に流れるわずかな電流のこと。主な原因は光電陰極からの熱電子放出。ノイズの原因になる。

∗ 「H9500 試験成績表」(HAMAMATSU PHOTONICS 2004) より† log4≃0.6024

2.3ヘッドアンプユニット

今回使用するヘッドアンプユニットは CLEAR-PULSE 社製で、マルチアノード光電子増倍管に繋いで使用する。また、PC のボードとの接続には RJ-45 コネクタの LAN ケーブルを使用するなど、セットアップも容易。∗ また、外形は、縦 52mm× 横 52mm× 長さ 110mm。

ADC の分解能は 12bit、4096ch。主な役割としては Preamp（アナログパルス増幅）+ Shaper（波形整形）+ ADC（アナログパルスをデジタルデータに変換）+ トリガー（イベントの発生を知らせるのに使用）の機能を持つ。
今回は通常の全 256ch 読み出しモードでデータを取得した。1event で 512byte (2Byte x 256ch)（SDS モードというトリガーが鳴ったデータだけを取得するモードもあるが今回は使用しなかった。）

図 2.3: MAPMT とヘッドアンプの接続イメージ（矢印１本が 64 個の出力に対応）

64ch ごとにというのは、ヘッドアンプユニットと MAPMT の接続形態上、各アノードが 64 個ごとでまとまって接続されているだけで、256 個のアノード出力は全部情報としてきちんと出力されている。
トリガースレッショルドについては、ヘッドアンプ裏面のイモネジを外し、さらに内部にあるネジ† を回して、調整する。ネジは 4 つ存在し、それぞれ、64ch 分に相当しているが、この方法では64ch 分各 ch ごとの細やかな調整はできず、まとめて調整するということになる。¶∗ トリガー入出力用には LEMO ケーブルを使用† はじめ、イモネジの方を回していたのでスレッショルドの調節ができなかった。

¶ もっと細かく ch ごとにスレッショルドの調節をするには、データを読み出すプログラム上で、1336bit のレジスタマップという情報のスレッショルド情報に対応する部分に値を書きこんでやればよいということだがうまくいかなった。同じレジスタマップ上のマスクという部分の値を書き換え、その ch を反応させなくするといった単純な書き換えの動作確認はうまくいったが、結局レジスタマップ情報の扱いについては完全には理解できなかった。

53装置のセットアップについて

3.1

一連の装置の接続

まず、シンチレーターをマルチアノード光電子増倍管に接着した。接着面は全反射を避けるためにガラスに近い屈折率を持つグリスで密着させた。また、このとき光を遮蔽するために、シンチレータの周りを白のテープで覆い、さらに上から暗幕で覆った。黒のテープでもって、MAPMTとシンチレータの接合部分は黒のテープで覆った。最後に、全体を暗幕で覆って遮蔽した。図 3.2参照。

図 3.1: 検出器イメージ

図 3.2: 装置とデータの流れのイメージ

PC のボードは CLEAR-PULSE 社製の 80057A 型の PMT DATA Processor という VME バスに対応した物を使用した。RJ-45 コネクタに対応するポートが二つあり、そこのひとつにヘッドアンプを LAN ケーブルでもって繋いだ。†
おおまかに言って、図の で検出器にイベントが来た時にトリガーをヘッドアンプからボードに送る。そして、 でトリガーを返す。この間にデータを取り続ける。

† あとはトリガー入出力用の LEMO ケーブルを繋ぐところがある

6

3.2

PC でのデータ取得までの流れ

PC 上でプログラムを走らせた後の流れを以下に示す。

図 3.3: イメージ
はじめに、プログラム∗ を PC 上で走らせたとき、まだ、線源のイベントのデータは取リ始めない。まずは意図的に CPU トリガーをかけてペデスタルデータ† を取得する。そのペデスタルデー
タは 80057A 型ボード上のメモリに書き込まれる。‡

次に、イベントが来ればノーマルトリガーがかかり、データの取得が始まる。このときボード上でペデスタル分のデータが引き算される。また、1ch あたりのデータのフォーマットが 16 ビットになっていて実際の 12 ビットのデータにヘッダがくっついている¶ ので、ヘッダのデータを見てやることによって、256ch のデータのはじめのch のデータと終わりの ch のデータが分かる。それで 1 イベント・256ch ごとにデータをうまく区切って集めることが可能となる。

∗ 64chMAPMT 用のサンプルプログラムを 256ch 用に調整した† ペデスタルデータとは、検出器に何もイベントが来ていない状態のときにもある程度出てくる数値データのこと。イ
ベントが来た時のデータを補正するために使う。補正しない時との差は明らかに補正した方がよいと分かるものだった。ペ
デスタルデータは 1ch あたり 32 ビットのメモリを容量が確保されているが、下位 12 ビットのみを使用する。12 ビット
符号無し整数となる。
‡ 100 回分取得してそれを平均した。ペデスタルの補正はどれ位やればいいのか（回数や掛かる時間）についてはあいまいになってしまった。
¶ さらに符号データに 1 ビット分使われるので、13 ビットとなる。ペデスタルデータを引き算するためのもの。

74実験・解析

4.1

前段階実験・自作プリアンプ

まず、本実験の前段階として、MAPMT の生の信号を確認するために、プリアンプを自作し、256ch の内の一つの ch の生の信号をオシロスコープで見た。線源は Cs137。

図 4.1: 接続とイメージプリアンプの時定数 CR=51µsec。∗∗ 減衰時間が 1µsec 程度だから、51µsec で問題ない。

8

4.2

1ch Na22

ここからは、実際にヘッドアンプユニットを接続してデータを取得。
とりあえず、一つ一つの ch ごとのパルスハイトをグラフ化して見ていったが各 ch ごとのグラフで最大の波高値の位置が異なっていたので、ピークを決めることは無理と判断した。
以下は代表として ch1 のグラフ。

図 4.2: ch1 の Na22

標準偏差 [ch]最大 [ch]23.535

表 4.1: 1ch Na2294.3

16ch Na22

次に、ひとつの ch を見た場合と同じく、今度は、16 の ch ごとに区切ってグラフにしてみた。このとき、16ch 分のデータを足して、一つのデータとした。この場合も、さきほどと同じく、やはり、各 ch ごとのグラフで最大の波高値の位置が異なっていた∗ ので、ピークを決めることは無理と判断した。以下は代表としてそのひとつのグラフ。

図 4.3: MAPMT16ch 分の Na22

標準偏差 [ch]最大 [ch]99.4417

表 4.2: MAPMT16ch 分の Na22

∗ 標準偏差の値 σ =100∼300 位のばらつきがあった。

10

4.4

256ch x 30000 event の波高値

以下は、30000event のデータ† を、MAPMT 上の ch で区別せず、波高値を全て書き出して、カウント数を縦軸にとってグラフ化したものである。この場合は、256 × 30000 = 768000 個のデータになる。また、右側は縦軸を対数にしたもの。

図 4.4: 256ch x 30000 Na22
このグラフからも、結局、ある波高値のところにピークが見えるということが分からないかった。
以上より、結局 256ch 分を足して見ていくことにした。

† 1event は 256ch の各 ch それぞれの波高データでまとまっている。

11

4.5

256ch 分 Na22

線源 Na22。256ch 足して見てやった場合。

図 4.5: 256ch 分 Na22

Na22 のピークは 511keV(89.8%) で、これは約 15000 のピークに対応してると判断できる。10500程度にもピークが立っているがこれはエネルギーに換算すると、320keV 程度となって、この位置にピークが出てしまうのはよく分からなかった。∗

3400ch 程度にあるピークは Noise かどうかはっきりしなかった。†∗ コンプトン散乱の可能性を考えると、hν ′ の最大 =hν01+mc2 /2hν02hν=511keV, mc =511keV より、340keV 程度にはなるがはっきりしない。
† 4.8 バックグラウンドとの比較参照

12

4.6

256ch 分 Co57

線源 Co57。256ch 足して見てやった場合。

図 4.6: 全 chCo57

Co57 は一番大きなピークは 122keV(85.6%) なので図の 6212ch のピークに対応すると判断できた。他のエネルギーの値は 136keV(10.7%) があるが、これは 6400ch 程度のピークになるはずであるが、122keV のピークに近くて判別は無理だと思われる。また、14.4keV(9.2%) のときは 3700程度のところにピークが立つはずであるが、バックグラウンド† でも表れてしまうよく分からない

3000∼3400ch 付近のピークのせいで見られないと判断した。

† 4.8 のバックグラウンド比較を参照

13

4.7

256ch 分 Ba133

線源 Ba133。∗

256ch 足して見てやった場合。

図 4.8: Ba133 ガウシアンでフィッティング
図 4.7: Ba133

標準偏差 [ch]最大 [ch]130010692(356keV)
表 4.3:

Ba133 は 356keV のピークが 62%であるので、図でもっともピークが立っているのが 356keV のピークだと判断できるが、81keV(34.1%) のエネルギーに対応する ch は 3500 程度† になるが 3200ch位にピークがあって（この 3200 程度のピーク自体が他の線源の場合と同じく Ba133 によるものだとは言えないと判断したので）これによってはっきり見えているとは言えなかった。他にも、

303keV(18.3%) の場合は 10000 程度になるが、これは図の 356keV のピークの中に埋もれていてはっきり区別できなかったと思われる。

∗ この他に Cs137 も見てやったが分かりにくかったので、Na22 Ba133 Co57 の 3 つの線源に絞った。
† 次々ページのキャリブレーションの式より

14

4.8

Background について

線源 Na22 でバックグラウンドとの比較をしてみた。

図 4.9: Na22 vs Background

図 4.10: Co57(赤)&Ba133(緑) vs Background(青)Na22 の計数率を求める。測定時間は Na22 の 30000event が 15 分だったので、それを基準にBackgraound を測定した。
Na22Background

30000events

3847events

表 4.4: Na22 とバックグラウンドの比較
計数率=22 ± 0.12/sec となった。よって、Background は無視しても問題ないと思われる。やはり、2900ch 付近にピークが立ってしまっているが、この原因がよく分からなかった。†† 5.1 まとめ参照

15

4.9

キャリブレーション

以上三つの線源での値を使いキャリブレーション† して、エネルギーと波高値の関係を求めた。

図 4.11: NaBaCo
キャリブレーションによってキャリブレーションの直線の式を求めた。

ch = 22(±0.8)E(kev) + 3396(±300)† 実験の最後の方において、試しに SDS モードを使ってみた。すると、通常モードでのキャリブレーションよりもうまく直線がグラフの 3 点を通った。しかし、なぜそうなるのかなどよく分からなかった。あるいは、もっと何度もデータを取得して通常モードと比較する必要があると思われる。

165

結果・まとめ5.1

これまでのまとめ

各線源に対してキャリブレーションの式でエネルギー分解能を計算して表にしたものを以下に示す。

図 5.1: 分解能:FWHM で評価した各 ch ごとで評価せず全 ch 足し上げることで評価した。そのために、分解能が悪くなってしまった。また、Ba133，Co57 にピーク位置細かなピークがあるので Na22 よりも分解能が悪くなったと思われる。
課題としては、やはり、バックグランドで 2900ch 位に存在するピークの説明がうまくできなかった。しかも、微妙に各線源の同じものかと思われたピークとは一致しない。Na22≃ 3400, Co57≃

3000or3400† , Ba133≃ 3200 となり、完全に同じものかどうか判断しかねた。256ch 足してやっているためかもしれない。16ch ごとに足して見ていった場合にも同じようなピークが見えたりして、原因がはっきり確かめられなかった。また、ゲインの違いも今回は考慮していないので、そのこともあるかもしれない。さらに、本来は 256ch についての位置分解能についても考慮すべきと思われるが、その実験をするまでに至らなかった。

† 細かいピーク 3000 と 3400 に立っていた

17

6

謝辞
今回の実験に関しまして、御指導頂きました助手の松本さん、TA の乾さんに深く感謝を申し上げます。また、線源などを快く貸して頂いた宇宙線研究室の方々に感謝を申し上げます。P6 のメンバーにも同じく感謝申し上げます。最後に、レポートが遅れてしまったことをお詫びします。

18

7

参考文献
「放射線計測ハンドブック」 グレン F. ノル (日刊工業新聞社 2001)「アイソトープ手帳 10 版」 日本アイソトープ協会 (丸善 2002)「放射線計測学」 長哲二(南山堂 2002)「放射線計測」 加藤貞幸

(培風館 1994)「計測における誤差解析入門」John R. Taylor (東京化学同人 2000)「取扱説明書 80057A 型 PMT Data Processor ファームウェア Ver.2」 (CLEAR-PULSE 2005)

19